Finite difference formulation differential equations

Next, we develop the finite difference formulation of heat conduction problems by replacing the derivatives in the differential equations by differences. In the following section we do it using the energy balance method, which does not require any knowledge of differential equations. 

Note that the boundary conditions have no effect on the finite difference formulation of interior nodes of the medium. This is not surprising since the control volume used in the development of the formulation does not involve any part of the boundary. You may recall that the boundary conditions had no effect on the differential equation of heat conduction in the medium either. 

The finite difference formulation is given above to demonstrate how difference equations ate obtained from differential equations. However, we use the energy balance approach in the following section.s to obtain the numerical formulation because it is more intuitive and can handle boundary conditions more easily. Besides, the energy balance approach does not require having the differential equation before the analysis. 

Numerical Solution. In the numerical formulation of THCC, Equations (2) and (3) are substituted into Equation (1). The resulting set of Nf, partial differential equations is transcribed into Nb finite-difference equations, using central differencing in space and the Crank-Nicolson method to obtain second-order accuracy in time. The set of unknowns consists of i = 1,..., ATft, and Pjt, k = 1,..., A/p, at each finite-difference node. Residue equations for the basis species are formed by algebraically summing all terms in the finite-difference forms of the transport equations. The finite-difference analogs of Equation (1) provide Ni, residue equations at each node the remaining Np residue equations are provided by the solubility products for the reactive solids. 

The numerical solutions of ordinary and partial differential equations are based on the finite difference formulation of these differential equations. Therefore, the stability and convergence considerations of finite difference solutions have important implications on the numerical solutions of differential equations. This topic will be discussed in more detail in Chaps. 5 and 6. 

With the above formulism a method is now defined for forming a finite difference set of equations for a partial differential equation of the initial value type in time and of the boundary value type in a spatial variable. The method can be applied to both linear and nonlinear partial differential equations. The result is an implicit equation which must be solved for the spatial variation of the solution... 

The Poisson equation has been used for both molecular mechanics and quantum mechanical descriptions of solvation. It can be solved directly using numerical differential equation methods, such as the finite element or finite difference methods, but these calculations can be CPU-intensive. A more efficient quantum mechanical formulation is referred to as a self-consistent reaction field calculation (SCRF) as described below. 

These coupled second-order partial differential equations do not have a closed-form solution. Accordingly, the approximate numerical technique of finite differences is employed. First, however, the boundary conditions must be prescribed in order to complete the formulation of the problem. Symmetry of the laminate about several planes permits reduction of the region of consideration to a quarter of the laminate cross section in the y-z plane at any value of x as shown in Figure 4-52. There, along the stress-free upper surface. 

The unsteady model, originally formulated in terms of a partial differential equation, is thus transformed into N difference differential equations. As a result of the finite-differencing, a solution can be obtained for the variation with respect to time of the water concentration, for every segment, throughout the bed. 

In the text, however, the numerical problem is formulated using momentum and total-energy balances on a finite control volume. The intent of this problem is to write a numerical simulation that is based on a finite-difference representation of the differential equations. 

Discuss the relationship between the continuity equation (Eq. 7.44) and Eq. 7.60 that represents the relationship between the physical radial coordinate and the stream function. Note that one is a partial differential equation and that the other is an ordinary differential equation. Formulate a finite-difference representation of the continuity equation in the primative form. Be sure to respect the order of the equation in the discrete representation. 

The most common methodology when solving transient problems using the finite element method, is to perform the usual Garlerkin weighted residual formulation on the spatial derivatives, body forces and time derivative terms, and then using a finite difference scheme to approximate the time derivative. The development, techniques and limitations that we introduced in Chapter 8 will apply here. The time discretization, explicit and implicit methods, stability, numerical diffusion etc., have all been discussed in detail in that chapter. For a general partial differential equation, we can write... 

Up to this point we have shown how conduction problems can be solved by finite-difference approximations to the differential equations. An equation is formulated for each node and the set of equations solved for the temperatures throughout the body. In formulating the equations we could just as well have used a resistance concept for writing the heat transfer between nodes. Designating our node of interest with the subscript i and the adjoining nodes with subscript j, we have the general-conduction-node situation shown in Fig. 3-10. At steady state the net heat input to node i must be zero or... 

Finite-Difference Methods. The numerical analysis literature abounds with finite difference methods for the numerical solution of partial differential equations. While these methods have been successfully applied in the solution of two-dimensional problems in fluid mechanics and diffusion (24, 25), there is little reported experience in the solution of three-dimensional, time-dependent, nonlinear problems. Application of these techniques, then, must proceed by extending methods successfully applied in two-dimensional formulations to the more complex problem of solving (7). The various types of finite-difference methods applicable in the solution of partial differential equations and their advantages and disadvantages are discussed by von Rosenberg (26), Forsythe and Wasow (27), and Ames (2S). 

The system of linear equations originating from the difference equation (2.308) has to be supplemented by the difference equations for the points around the boundaries where the decisive boundary conditions are taken into account. As a simplification we will assume that the boundaries run parallel to the x- and y-directions. Curved boundaries can be replaced by a series of straight lines parallel to the x- and y-axes. However a sufficient degree of accuracy can only be reached in this case by having a very small mesh size Ax. If the boundaries are coordinate lines of a polar coordinate system (r, differential equation and its boundary conditions are formulated in polar coordinates and then the corresponding finite difference equations are derived. 

In the second chapter we consider steady-state and transient heat conduction and mass diffusion in quiescent media. The fundamental differential equations for the calculation of temperature fields are derived here. We show how analytical and numerical methods are used in the solution of practical cases. Alongside the Laplace transformation and the classical method of separating the variables, we have also presented an extensive discussion of finite difference methods which are very important in practice. Many of the results found for heat conduction can be transferred to the analogous process of mass diffusion. The mathematical solution formulations are the same for both fields. 

Finite difference methods (FDM) are directly derived from the space time grid. Focusing on the space domain (horizontal lines in Fig. 6.6), the spatial differentials are replaced by discrete difference quotients based on interpolation polynomials. Using the dimensionless formulation of the balance equations (Eq. 6.107), the convection term at a grid point j (Fig. 6.6) can be approximated by assuming, for example, the linear polynomial. 

It is obvious that we obtain a stability condition that is not much different from the stability condition of the initial value equation. If At is larger than 2 y/MfK, (the cosine is smaller than —1), the solution grows exponentially and is numerically unstable. Hence, in the straightforward boundary value formulation of classical mechanics, we gain very little in terms of stability and step size compared to the solution of the initial value differential equation. The difficulty is not in the philosophical view (global or local) but in the estimate of the time derivative, which is approximated by a local finite difference expression. 

The material presented earlier was confined to steady-state flows over simply shaped bodies such as flat plates, with and without pressure gradients in the streamwise direction, or stagnation regions on blunt bodies. The simplicity of these flow configurations allows reduction of the problems to the solution of steady-state ordinary differential equations. The evaluation of convective heat transfer to more complex three-dimensional configurations, characteristic of real aerodynamic vehicles, involves the solution of partial differential equations. Even when the latter are confined to steady-state problems, they require extensive use of computers in the solution of finite difference or finite element formulations Nonsteady flows further complicate the problems by introducing another dimension, namely, time. 

The Galerkin weighted residual method is employed to formulate the finite element discretisation. An implicit mid-interval backward difference algorithm is implemented to achieve temporal discretisation. With appropriate initial and boundary conditions the set of non-linear coupled governing differential equations can be solved. 

Maxwell s famous four equations are found in many versions, for example, in differential or integral form and with different parameters involved. Equations 9.1—9.4 show one example set (the Minkowski formulation). The differential form relates the time and space derivatives at a point (in an infinitesimal small volume) to the current density at that point. The integral form relates to a defined finite volume. Ideal charge distributions are often discontinuous, and so not differentiable therefore, the integral form is a more generally applicable form. 

---