Finite difference equation

To solve the PB equation for arbitrary geometries requires some type of discretization, to convert the partial differential equation into a set of difference equations. Finite difference methods divide space into a cubic lattice, with the potential, charge density, and ion accessibility defined at the lattice points (or grid points ) and the permittivity defined on the branches (or grid lines ). Equation [1] becomes a system of simultaneous equations referred to as the finite difference Poisson-Boltzmann (FDPB) equation ... 

Difference equations (finite changes, steady stale)... 

Since this equation introduces a new variable, Cq, another equation is needed and is obtained by writing the finite difference equation for = 1, too. 

If average diffusion coefficients are used, then the finite difference equation is as follows. 

We denote by C the value of c(x , t) at any time. Thus, C is a function of time, and differential equations in C are ordinary differential equations. By evaluating the diffusion equation at the ith node and replacing the derivative with a finite difference equation, the following working equation is derived for each node i, i = 2,. . . , n (see Fig. 3-52). 

Elliptic Equations Elhptic equations can be solved with both finite difference and finite element methods. One-dimensional elhptic problems are two-point boundary value problems. Two- and three-dimensional elliptic problems are often solved with iterative methods when the finite difference method is used and direct methods when the finite element method is used. So there are two aspects to consider howthe equations are discretized to form sets of algebraic equations and howthe algebraic equations are then solved. 

Later work by Hashem and Shepcevich [Chem. E/ig. Prog., 63, Symp. Sen 79, 35, 42 (1967)] offers more accurate second-order finite-difference equations. 

If the stress is at the primary time step loeation and the veloeities are at the middle of the time step, then the resulting finite-difference equation is second-order accurate in space and time for uniform time steps and elements. If all quantities are at the primary time step, then a more complicated predictor-corrector procedure must be used to achieve second-order accuracy. A typical predictor-corrector scheme predicts the stresses at the middle of the time step and uses them to calculate the divergence of the stress tensor. 

Smith, G.D., 1985. Numerical Solution of Partial Differential Equations Finite Difference Methods, 3rd edition. Clarendon Press. 

The resulting finite difference equations constitute a set of nonho-mogeneous linear algebraic equations. Because there are three dependent variables, the number of equations in the set is three times the number of material points. Obviously, if a large number of points is required to accurately represent the continuous elastic body, a computer is essential. 

Numerical methods have been developed by replacing the differential equation by a finite difference equation. Thus in a problem of unidirectional flow of heat ... 

The most general method of tackling the problem is the use of the finite-element technique 8 to determine the temperature distribution at any time by using the finite difference equation in the form of equation 9.40. 

Finite difference equation, unsteady hear transfer 398... 

From a physical point of view, the finite difference method is mostly based based on the further replacement of a continuous medium by its discrete model. Adopting those ideas, it is natural to require that the principal characteristics of a physical process should be in full force. Such characteristics are certainly conservation laws. Difference schemes, which express various conservation laws on grids, are said to be conservative or divergent. For conservative schemes the relevant conservative laws in the entire grid domain (integral conservative laws) do follow as an algebraic corollary to difference equations. 

Lax, P. and Richtmyer, R. (1956) A servey of stability of linear finite difference equations. Comm. Pure Appl. Mathem., 9, 267-293. 

Alternatively, the difference-equation model form can be derived directly by dividing the length of the heat exchanger into N finite-difference elements or segments, each of length AZ, as shown in Fig. 4.27. 

This then provides a physical derivation of the finite-difference technique and shows how the solution to the differential equations can be propagated forward in time from a knowledge of the concentration profile at a series of mesh points. Algebraic derivations of the finite-difference equations can be found in most textbooks on numerical analysis. There are a variety of finite-difference approximations ranging from the fully explicit method (illustrated above) via Crank-Nicolson and other weighted implicit forward. schemes to the fully implicit backward method, which can be u.sed to solve the equations. The methods tend to increase in stability and accuracy in the order given. The difference scheme for the cylindrical geometry appropriate for a root is... 

Smith and Brinkley developed a method for determining the distribution of components in multicomponent separation processes. Their method is based on the solution of the finite-difference equations that can be written for multistage separation processes, and can be used for extraction and absorption processes, as well as distillation. Only the equations for distillation will be given here. The derivation of the equations is given by Smith and Brinkley (1960) and Smith (1963). For any component i (suffix i omitted in the equation for clarity)... 

The computational procedure can now be explained with reference to Fig. 19. Starting from points Pt and P2, Eqs. (134) and (135) hold true along the c+ characteristic curve and Eqs. (136) and (137) hold true along the c characteristic curve. At the intersection P3 both sets of equations apply and hence they may be solved simultaneously to yield p and W for the new point. To determine the conditions at the boundary, Eq. (135) is applied with the downstream boundary condition, and Eq. (137) is applied with the upstream boundary condition. It goes without saying that in the numerical procedure Eqs. (135) and (137) will be replaced by finite difference equations. The Newton-Raphson method is recommended by Streeter and Wylie (S6) for solving the nonlinear simultaneous equations. In the specified-time-... 

Export processes are often more complicated than the expression given in Equation 7, for many chemicals can escape across the air/water interface (volatilize) or, in rapidly depositing environments, be buried for indeterminate periods in deep sediment beds. Still, the majority of environmental models are simply variations on the mass-balance theme expressed by Equation 7. Some codes solve Equation 7 directly for relatively large control volumes, that is, they operate on "compartment" or "box" models of the environment. Models of aquatic systems can also be phrased in terms of continuous space, as opposed to the "compartment" approach of discrete spatial zones. In this case, the partial differential equations (which arise, for example, by taking the limit of Equation 7 as the control volume goes to zero) can be solved by finite difference or finite element numerical integration techniques. 

Models of the above have been presented by various researchers of the U.S. Geological Survey (USGS) and the academia. The above equation has been solved principally (a) numerically over a temporal and spatial discretized domain, via finite difference or finite element mathematical techniques (e.g., 11) (b) analytically, by seeking exact solutions for simplified environmental conditions (e.g., 12) or (c) probabilistically (e.g., 13). 

Strictly speaking, finite difference or finite element solutions to differential equations are simply multiplying the number of comparments many times, but the mathematical rules for linking cells in difference calculations are rigorously set by the form of the equations. 

Rendering a transport equation in finite difference form is a straightforward and well known process (e.g., Peaceman, 1977 Smith, 1986). The derivatives of C in time and space are replaced with finite difference equivalents. The resulting difference equation written at a nodal block (7, J) now includes concentration... 

Finite-difference equations will be formulated with five radial increments m = 0,1, 2, 3, 4, and 5. These terms occur in the equations... 

For a finite, (rather than an infinite) chain under the action of an applied field, it turns out that the energy spectrum forms only an approximate WSL (Fleinrichs and Jones 1972). To see this result, we start with the difference equations for a finite TB chain of length m ... 

In summary, DQMOM is a numerical method for solving the Eulerian joint PDF transport equation using standard numerical algorithms (e.g., finite-difference or finite-volume codes). The method works by forcing the lower-order moments to agree with the corresponding transport equations. For unbounded joint PDFs, DQMOM can be applied... 

---