Solving the governing equations

In this section we consider how Newton-Raphson iteration can be applied to solve the governing equations listed in Section 4.1. There are three steps to setting up the iteration (1) reducing the complexity of the problem by reserving the equations that can be solved linearly, (2) computing the residuals, and (3) calculating the Jacobian matrix. Because reserving the equations with linear solutions reduces the number of basis entries carried in the iteration, the solution technique described here is known as the reduced basis method.  

The computing time required to evaluate Equation 4.19 in a Newton-Raphson iteration increases with the cube of the number of equations considered (Dongarra et al., 1979). The numerical solution to Equations 4.3 1.6, therefore, can be found most rapidly by reserving from the iteration any of these equations that can be solved linearly. There are four cases in which equations can be reserved  

The basis entries corresponding to these two cases are given by the reduced basis,  

We will carry the subscript r (for reduced ) to indicate that a vector or matrix includes only entries conforming to one of the two nonlinear cases. 

The residual functions measure how well a guess (nw,m,i)r satishes the governing Equations 4.3 f.4. The form of the residuals can be written, 

The nonlinear portion of the problem, then, consists of just two parts  

---

The word deterministic" means that the model employs a specific surface geometry or prescribed roughness data as an input of the numerical procedure for solving the governing equations. The method was originally adopted in micro-EHL to predict local film thickness and pressure distributions over individual asperities, and it can be used to solve the mixed lubrication problems when properly combined with the solutions of asperity contacts. 

The computational code used in solving the hydrodynamic equation is developed based on the CFDLIB, a finite-volume hydro-code using a common data structure and a common numerical method (Kashiwa et al., 1994). An explicit time-marching, cell-centered Implicit Continuous-fluid Eulerian (ICE) numerical technique is employed to solve the governing equations (Amsden and Harlow, 1968). The computation cycle is split to two distinct phases a Lagrangian phase and a remapping phase, in which the Arbitrary Lagrangian Eulerian (ALE) technique is applied to support the arbitrary mesh motion with fluid flow. 

Compute the velocity field k" by solving the governing equation, Eqs. (5-6), using the cell-centered ICE technique and ALE technique (Kashiwa et al., 1994). 

Providing an additional piece of information about the size of each phase predicts that a total of Ni + N, or Nc, values is needed to constrain the system s state and extent. This total matches the number of variables we must supply in order to solve the governing equations. Hence, although we can make no claim that we have cast the governing equations in simplest form, we can say that we have reduced the number of independent variables to the minimum allowed by thermodynamics. 

The procedure for solving the governing equations parallels the technique described in Chapter 4, with the added complication of accounting for electrostatic effects. We begin as before by identifying the nonlinear portion of the problem to form the reduced basis,... 

The equations expressed in this stepwise manner are somewhat easier to integrate into certain reaction configurations, such as the flush or flow-through model described later in this chapter. We could also update Mm in this manner, but there is no need to do so. A gas species Am appears in the basis only when its fugacity fm is known, so the value of each Mm results from solving the governing equations, as described in Chapter 4. 

In this chapter we consider how to construct reaction models that are somewhat more sophisticated than those discussed in the previous chapter, including reaction paths over which temperature varies and those in which species activities and gas fugacities are buffered. The latter cases involve the transfer of mass between the equilibrium system and an external buffer. Mass transfer in these cases occurs at rates implicit in solving the governing equations, rather than at rates set explicitly by the modeler. In Chapter 16 we consider the use of kinetic rate laws, a final method for defining mass transfer in reaction models. 

To date, spray modeling has been largely dependent on solving the governing equations of multiphase flows and specifying initial and... 

The marginal stability envelopes are shifted when a bias voltage is applied, and recently Iwamoto et al. (1991) prepared a number of stability maps showing the effects of bias voltage. They solved the governing equations numerically. 

Ray et al. (1991b) wrote conservation equations for the two species in the droplet and solved the governing equations to yield the evaporation rate in terms of the square of the droplet radius. If the outer material is relatively nonvolatile, the core material must diffuse through the coating of constant... 

We will introduce the product rule through demonstrating its use in an example problem. The product rule can be used to expand a solution without source and sink terms to the unsteady, one-dimensional diffusion equation to two and three dimensions. It does not work as well in developing solutions to all problems and therefore is more of a technique rather than a rule. Once again, the final test of any solution is (1) it must solve the governing equation(s) and (2) it must satisfy the boundary conditions. 

Full nodes (b). These nodes are the ones that have already filled during the mold filling process (/6 = 1). When solving the governing equations at any given time step, the pressure is unknown inside these control volumes. 

It is not, at the present time, possible on a routine basis to solve the governing equations to obtain the variation of the flow variables with time in turbulent flow. In most analyses of turbulent flow it is therefore usual to express the variables in terms of a time-averaged mean value plus a time-varying deviation from this mean value, i.e., as discussed in Chapter 2 [1],[2],[3],[4],[5], to express the variables in the following way ... 

The present section is not really concerned with flow over a flat plate. It is instead concerned with flows in which the frce-stream velocity is varying wit. . x. The solution in such a case usually has to be obtained by numerically solving the governing equations. For this purpose, it is convenient to introduce the following dimensionless variables ... 

Turbulent natural convective flows can also be analyzed by numerically solving the governing equations together with some form of turbulence model. This is... 

Finding the source location and the time history of the solute in ground-water can be categorized as a problem of time inversion. This means that we have to solve the governing equations backward in time. Modeling contaminant transport using reverse time is an ill-posed problem since the process, being dispersive is irreversible. Because of this ill-posedness, the problems have discontinuous dependence on data and are sensitive to the errors in data. 

To analyze the transport and retention of chemical contaminants in groundwater flowing through soils, experimental and theoretical studies generated several reliable models. Diverse numerical methods have been applied to solve the governing equations efficiently. Some computer models include the simulation of physical and chemical processes. 

The numerical solution to the advection-dispersion equation and associated adsorption equations can be performed using finite difference schemes, either in their implicit and/or explicit form. In the one-dimensional MRTM model (Selim et al., 1990), the Crank-Nicholson algorithm was applied to solve the governing equations of the chemical transport and retention in soils. The web-based simulation system for the one-dimensional MRTM model is detailed in Zeng et al. (2002). The alternating direction-implicit (ADI) method is used here to solve the three-dimensional models. 

A thermodynamic quantity of considerable importance in many combustion problems is the adiabatic flame temperature. If a given combustible mixture (a closed system) at a specified initial T and p is allowed to approach chemical equilibrium by means of an isobaric, adiabatic process, then the final temperature attained by the system is the adiabatic flame temperature T. Clearly depends on the pressure, the initial temperature and the initial composition of the system. The equations governing the process are p = constant (isobaric), H = constant (adiabatic, isobaric) and the atom-conservation equations combining these with the chemical-equilibrium equations (at p, T ) determines all final conditions (and therefore, in particular, Tj). Detailed procedures for solving the governing equations to obtain Tj> are described in [17], [19], [27], and [30], for example. Essentially, a value of Tf is assumed, the atom-conservation equations and equilibrium equations are solved as indicated at the end of Section A.3, the final enthalpy is computed and compared with the initial enthalpy, and the entire process is repeated for other values of until the initial and final enthalpies agree. 

---