Solution Methods for Second Order Nonlinear Equations

we observe that the arbitrary constant of integration is implicit, which is quite usual for nonlinear systems. 

4 SOLUTION METHODS FOR SECOND ORDER NONUNEAR EQUATIONS 

As stated earlier, much is known about linear equations of higher order, but no general technique is available to solve the nonlinear equations that arise frequently in natural and man-made systems. When analysis fails to uncover the analytical solution, the last recourse is to undertake numerical solution methods, as introduced in Chapters 7 and 8. 

We begin this section by illustrating the types of nonlinear problems that can be resolved using standard methods. Some important nonlinear second order 

The general strategy for attacking nonlinear equations is to reduce them to linear form. Often, inspection of the equation suggests the proper approach. The two most widely used strategies are as follows. 

---

Solution Methods for Second Order Nonlinear Equations 51... 

Solution Methods for Second Order Nonlinear Equations 57 and the arbitrary second constant of integration is obtained from y = 1, jc = 0,... 

The spurious or satellite term in the solution is introduced by using a second-order difference equation to approximate a first-order differential equation. An extra condition is needed to fix the solution of the second-order equation, and this condition must be that the coefficient of the spurious part of the solution is zero. In the general case of a nonlinear difference equation, no method is available for meeting this condition exactly. 

For nonlinear reaction kinetics, a numerical solution of the balance Equation 4.121 is carried out. For example, for second-order kinetics, R = kcACB, with an arbitrary stoichiometry, the generation rate expressions, ta = —va CaCb and tb = —vb caCb, are inserted into the mass balance expression, which is solved numerically using, for example, a polynomial approximation (orthogonal collocation method). The performances of the normal dispersion model and its segregated or maximum-mixed variants are compared in Figure 4.34. The symbols are explained in the figure. The comparison reveals that the differences between the segregated, maximum-mixed, and normal axial dispersion models are notable at moderate Damkohler numbers R = Damkohler number). 

Steady state mass or heat transfer in solids and current distribution in electrochemical systems involve solving elliptic partial differential equations. The method of lines has not been used for elliptic partial differential equations to our knowledge. Schiesser and Silebi (1997)[1] added a time derivative to the steady state elliptic partial differential equation and applied finite differences in both x and y directions and then arrived at the steady state solution by waiting for the process to reach steady state. [2] When finite differences are applied only in the x direction, we arrive at a system of second order ordinary differential equations in y. Unfortunately, this is a coupled system of boundary value problems in y (boundary conditions defined at y = 0 and y = 1) and, hence, initial value problem solvers cannot be used to solve these boundary value problems directly. In this chapter, we introduce two methods to solve this system of boundary value problems. Both linear and nonlinear elliptic partial differential equations will be discussed in this chapter. We will present semianalytical solutions for linear elliptic partial differential equations and numerical solutions for nonlinear elliptic partial differential equations based on method of lines. 

To guide model development, the observed data were first examined graphically to determine general characteristics and to look for trends with respect to dose, time, and the impact of anti-mAb antibodies. Models were developed using NONMEM (Version 5). Two different model types were developed the first model (MODEL 1, see Appendix 45.1) used an analytical solution (closed-form) where the nonlinearity was accounted for by allowing the model parameters to be a function of mAb dose and the titer of anti-mAb antibody, while the second model (MODEL 2, see Appendix 45.2) used differential equations to allow a more mechanistic approach to characterize the nonlinearity. For each model, three estimation methods were evaluated first-order (FO), first-order conditional estimation (FOCE), and FOCE with interaction. Various forms of between-subject variability models were evalu-... 

Methods for solving mass and heat transfer problems. The convective diffusion equation (3.1.1) is a second-order linear partial differential equation with variable coefficients (in the general case, the fluid velocity depends on the coordinates and time). Exact closed-form solutions of the corresponding problems can be found only in exceptional cases with simple geometry [79,197, 270, 370, 516]. This is especially true of the nonlinear equation (3.1.17). Exact solutions are important for adequate understanding of the physical background of various phenomena and processes. They can serve as test solutions to verify whether the problem is well-posed or to estimate the accuracy of the corresponding numerical, asymptotic, and approximate methods. 

Solution of Mathematical Model for Case 1. For the Case 1 solution iterative techniques were ruled unacceptable owing to the excessive time requirements of such methods. Several investigators (27, 28, 29, 30) working with similar noncoupled systems found that the Crank-Nicholson 6-point implicit differencing method (31) provided an excellent solution. For the solution of Equation (8) we decided to apply the Crank-Nicholson method to the second-order partials and corresponding explicit methods to the first-order partials. Nonlinear coefficients were treated in a special manner outlined by Reneau et al (5). 

An analysis of chemical desorption has recently been published (Chem.Eng.Sci., 21 0980)), which is based on a number of simplifying assumptions the film theory model is assumed, the diffusivities of all species are taken to be equal to each other, and in the solution of the differential equations an approximation which is second order with respect to distance from the gas-liquid interface is used this approximation was introduced as early as 1948 by Van Krevelen and Hoftizer. However, the assumptions listed above are not at all drastic, and two crucial elements are kept in the analysis reversibility of the chemical reactions and arbitrary chemical mechanisms and stoichiometry.The result is a methodology for developing, for any given chemical mechanism, a highly nonlinear, implicit, but algebraic equation for the calculation of the rate enhancement factor as a function of temperature, bulk-liquid composition, interface gas partial pressure and physical mass transfer coefficient The method of solution is easily gene ralized to the case of unequal diffusivities and corrections for differences between the film theory and the penetration theory models can be calculated. 

B. If the reaction is second order and the numerical value of Da (Damkohler number) is the same as in part 1 (although its definition is slightly different), find the exit conversion using the Fox s iterative method (explain your formulation of adjoint equations for the iterative solution of the nonlinear two-point boundary-value differential equation). 

In [165] a study of a new methodology for development of efficient methods for the numerical solution of second-order periodic initial value problems (IVPs) of ordinary differential equations is presented. The methodology is based on the development of nonlinear numeircal methods. In this paper the authors study the following nonlinear scheme ... 

Equation (1), with the associated boundary conditions, is a nonlinear second-order boundary-value ODE. This was solved by the method of collocation with piecewise cubic Hermite polynomial basis functions for spatial discretization, while simple successive substitution was adequate for the solution of the resulting nonlinear algebraic equations. The method has been extensively described before [9], and no problems were found in this application. 

PDEs can be classified in different ways. The classification is important because the solution methods often apply only to a specified class of PDEs. To start with, PDEs can be classified by the number of variables, e.g. Uf = which contain two independent variables, t and x. The order of a PDE is the order of the highest order derivative that appears in the PDE. For example, Ut = is a first-order PDE, whereas Ut = Uxx is a second-order PDE. In addition, it is important to make a distinction between nonlinear and linear PDEs. An example of a well-known non-linear PDE is the Navier-Stokes equation, which describes the motion of fluids. In a linear PDE, the dependent variable and its derivatives appear in a linear fashion. The linear second-order PDE in Equation (6.69) can be classified as... 

---