Nonlinear Partial Differential Equations

Errors are proportional to At for small At. When the trapezoid rule is used with the finite difference method for solving partial differential equations, it is called the Crank-Nicolson method. The implicit methods are stable for any step size but do require the solution of a set of nonlinear equations, which must be solved iteratively. The set of equations can be solved using the successive substitution method or Newton-Raphson method. See Ref. 36 for an application to dynamic distillation problems. 

If the reaetion rate is a funetion of pressure, then the momentum balanee is eonsidered along with the mass and energy balanee equations. Both Equations 6-105 and 6-106 are eoupled and highly nonlinear beeause of the effeet of temperature on the reaetion rate. Numerieal methods of solution involving the use of finite differenee are generally adopted. A review of the partial differential equation employing the finite differenee method is illustrated in Appendix D. Eigures 6-16 and 6-17, respeetively, show typieal profiles of an exo-thermie eatalytie reaetion. 

In the previous sections, we briefly introduced a number of different specific models for crystal growth. In this section we will make some further simplifications to treat some generic behavior of growth problems in the simplest possible form. This usually leads to some nonlinear partial differential equations, known under names like Burgers, Kardar-Parisi-Zhang (KPZ), Kuramoto-Sivashinsky, Edwards-Anderson, complex Ginzburg-Landau equation and others. 

Mathematical physics deals with a variety of mathematical models arising in physics. Equations of mathematical physics are mainly partial differential equations, integral, and integro-differential equations. Usually these equations reflect the conservation laws of the basic physical quantities (energy, angular momentum, mass, etc.) and, as a rule, turn out to be nonlinear. 

The partial differential equations describing the catalyst particle are discretized with central finite difference formulae with respect to the spatial coordinate [50]. Typically, around 10-20 discretization points are enough for the particle. The ordinary differential equations (ODEs) created are solved with respect to time together with the ODEs of the bulk phase. Since the system is stiff, the computer code of Hindmarsh [51] is used as the ODE solver. In general, the simulations progressed without numerical problems. The final values of the rate constants, along with their temperature dependencies, can be obtained with nonlinear regression analysis. The differential equations were solved in situ with the backward... 

Figure 8 depicts our view of an ideal structure for an applications program. The boxes with the heavy borders represent those functions that are problem specific, while the light-border boxes represent those functions that can be relegated to problem-independent software. This structure is well-suited to problems that are mathematically either systems of nonlinear algebraic equations, ordinary differential equation initial or boundary value problems, or parabolic partial differential equations. In these cases the problem-independent mathematical software is usually written in the form of a subroutine that in turn calls a user-supplied subroutine to define the system of equations. Of course, the user must write the subroutine that defines his particular system of equations. However, that subroutine should be able to make calls to problem-independent software to return many of the components that are needed to assemble the governing equations. Specifically, such software could be called to return in-... 

The combined fiuld fiow, heat transfer, mass transfer and reaction problem, described by Equations 2-7, lead to three-dimensional, nonlinear, time dependent partial differential equations. The general numerical solution of these goes beyond the memory and speed capabilities of the current generation of supercomputers. Therefore, we consider appropriate physical assumptions to reduce the computations. 

The modeling of steady-state problems in combustion and heat and mass transfer can often be reduced to the solution of a system of ordinary or partial differential equations. In many of these systems the governing equations are highly nonlinear and one must employ numerical methods to obtain approximate solutions. The solutions of these problems can also depend upon one or more physical/chemical parameters. For example, the parameters may include the strain rate or the equivalence ratio in a counterflow premixed laminar flame (1-2). In some cases the combustion scientist is interested in knowing how the system mil behave if one or more of these parameters is varied. This information can be obtained by applying a first-order sensitivity analysis to the physical system (3). In other cases, the researcher may want to know how the system actually behaves as the parameters are adjusted. As an example, in the counterflow premixed laminar flame problem, a solution could be obtained for a specified value of the strain... 

Let us consider the genera class of systems described by a system of n nonlinear parabolic or hyperbolic partial differential equations. For simplicity vve assume that we have only one spatial independent variable, z. 

The material and energy balances of a tubular vessel are based on the conservation law, Eq 2.42, applied to a differential ring between r and r+dr and z and z+dz. A material balance is derived, for example, in problem P5.08.01, and is quoted in Table 2.6 along with the heat balance. The result is a pair of second order partial differential equations, usually nonlinear, that must be solved numerically. Table 2.6 indicates one possible procedure, but computer software is plentiful. 

The partial differential equations representing material and energy balances of a reaction in a packed bed are rarely solvable by analytical means, except perhaps when the reaction is of zero or first order. Two examples of derivation of the equations and their analytical solutions are P8.0.1.01 and P8.01.02. In more complex cases analytical, approximations can be made (by "Collocation" or "Perturbation", for instance), but these usually are quite sophisticated to apply. Numerical solutions, on the other hand, are simple in concept and are readily implemented on a computer. Two such methods that are suited to nonlinear kinetics problems will be described. 

The nonlinearity of the system of partial differential equations (51) and (52) poses a serious obstacle to finding an analytical solution. A reported analytical solution for the nonlinear problem of diffusion coupled with complexation kinetics was erroneous [12]. Thus, techniques such as the finite element method [53-55] or appropriate change of variables (applicable in some cases of planar diffusion) [56] should be used to find the numerical solution. One particular case of the nonlinear problem where an analytical solution can be given is the steady-state for fully labile complexes (see Section 3.3). However, there is a reasonable assumption for many relevant cases (e.g. for trace elements such as... 

While nonlinear in g2 (the coefficients 7,.r and b are lengthy integral expressions), the partial differential equation is linear in the derivatives. It can thus be solved by the method of characteristics, with the boundary conditions given by the coupling at /x = 0, as obtained from finite-T lattice QCD. 

The problem in obtaining a state space model for the dynamics of the CSD from this physical model is that the population balance is a (nonlinear) first-order partial differential equation. Consequently, to obtain a state space model the population balance must be transformed into a set of ordinary differential equations. After this transformation, the state space model is easily obtained by substitution of the algebraic relations and linearization of the ordinary differential equations. 

The process inputs are defined as the heat input, the product flow rate and the fines flow rate. The steady state operating point is Pj =120 kW, Q =.215 1/s and Q =.8 1/s. The process outputs are defined as the thlrd moment m (t), the (mass based) mean crystal size L Q(tK relative volume of crystals vr (t) in the size range (r.-lO m. In determining the responses of the nonlinear model the method of lines is chosen to transform the partial differential equation in a set of (nonlinear) ordinary differential equations. The time responses are then obtained by using a standard numerical integration technique for sets of coupled ordinary differential equations. It was found that discretization of the population balance with 1001 grid points in the size range 0. to 5 10 m results in very accurate solutions of the crystallizer model. 

The simulation of a continuous, evaporative, crystallizer is described. Four methods to solve the nonlinear partial differential equation which describes the population dynamics, are compared with respect to their applicability, accuracy, efficiency and robustness. The method of lines transforms the partial differential equation into a set of ordinary differential equations. The Lax-Wendroff technique uses a finite difference approximation, to estimate both the derivative with respect to time and size. The remaining two are based on the method of characteristics. It can be concluded that the method of characteristics with a fixed time grid, the Lax-Wendroff technique and the transformation method, give satisfactory results in most of the applications. However, each of the methods has its o%m particular draw-back. The relevance of the major problems encountered are dicussed and it is concluded that the best method to be used depends very much on the application. 

First-order error analysis is a method for propagating uncertainty in the random parameters of a model into the model predictions using a fixed-form equation. This method is not a simulation like Monte Carlo but uses statistical theory to develop an equation that can easily be solved on a calculator. The method works well for linear models, but the accuracy of the method decreases as the model becomes more nonlinear. As a general rule, linear models that can be written down on a piece of paper work well with Ist-order error analysis. Complicated models that consist of a large number of pieced equations (like large exposure models) cannot be evaluated using Ist-order analysis. To use the technique, each partial differential equation of each random parameter with respect to the model must be solvable. 

The reason for constructing this rather complex model was that even though the mathematical equations may be easily set up using the dispersion model, the numerical solutions are quite involved and time consuming. Deans and Lapidus were actually concerned with the more complicated case of mass and heat dispersion with chemical reactions. For this case, the dispersion model yields a set of coupled nonlinear partial differential equations whose solution is quite formidable. The finite-stage model yields a set of differential-double-difference equations. These are ordinary differential equations, which are easier to solve than the partial differential equations of the dispersion model. The stirred-tank equations are of an initial-value type rather than the boundary-value type given by the dispersion model, and this fact also simplifies the numerical work. 

The principal use of Eq. (173) is in conjunction with a similar heat dispersion equation. Unfortunately, a system of coupled nonlinear partial differential equations then has to be solved, which is very difficult even with the aid of computers. In the oxidation of sulfur dioxide. Hall and Smith (HI) found relatively good agreement between theory and experiment near the center of the reactor. Their calculations were based on the heat-dispersion equation, and they did not take detailed mass dispersion into account. Baron (B2) later solved the mass and heat dispersion equations simultaneously by a novel graphical method, and found better agreement between his calculations and the data of Hall and Smith. 

Equation (E4.4.1) is a nonlinear partial differential equation, because of the velocity u that appears in front of the velocity gradient du/dx. The boundary layer thickness is generally defined as the distance from the plate where the momentum reaches 99% of the free-stream momentum. We will assign (Blasius, 1908)... 

The mathematical model developed in the preceding section consists of six coupled, three-dimensional, nonlinear partial differential equations along with nonlinear algebraic boundary conditions, which must be solved to obtain the temperature profiles in the gas, catalyst, and thermal well the concentration profiles and the velocity profile. Numerical solution of these equations is required. 

---