Differential equation analytical solution

Mean temperature differences in such flow patterns are obtained by solving the differential equation. Analytical solutions have been found for the simpler cases, and numerical ones for many important complex patterns, whose results sometimes are available in generalized graphical form. 

While the boundary layer equations that were derived in the preceding sections are much simpler than the general equations from which they were derived, they still form a complex set of simultaneous partial differential equations. Analytical solutions to this set of equations have been obtained in a few important cases. For the majority of flows, however, a numerical solution procedure must be adopted. Such solutions are readily obtained today using modest modem computing facilities. This was, however, not always so. For this reason, approximate solutions to the boundary layer equations have in the past been quite widely used. While such methods of solution are less important today, they are still used to some extent. One such approach will, therefore, be considered in the present text. 

Mathematical models of flow processes are non-linear, coupled partial differential equations. Analytical solutions are possible only for some simple cases. For most flow processes which are of interest to a reactor engineer, the governing equations need to be solved numerically. A brief overview of basic steps involved in the numerical solution of model equations is given in Section 1.2. In this chapter, details of the numerical solution of model equations are discussed. 

In the previous three chapters, we described various analytical techniques to produce practical solutions for linear partial differential equations. Analytical solutions are most attractive because they show explicit parameter dependences. In design and simulation, the system behavior as parameters change is quite critical. When the partial differential equations become nonlinear, numerical solution is the necessary last resort. Approximate methods are often applied, even when an analytical solution is at hand, owing to the complexity of the exact solution. For example, when an eigenvalue expression requires trial-error solutions in terms of a parameter (which also may vary), then the numerical work required to successfully use the analytical solution may become more intractable than a full numerical solution would have been. If this is the case, solving the problem directly by numerical techniques is attractive since it may be less prone to human error than the analytical counterpart. 

The root time method of data analysis for diffusion coefficient determination was developed by Mohamed and Yong [142] and Mohamed et al. [153]. The procedure used for computing the diffusion coefficient utilizes the analytical solution of the differential equation of solute transport in soil-solids (i.e., the diffusion-dispersion equation) ... 

Analytical solutions of Fick s laws are most easily derived using Laplace transforms, a subject described in every undergraduate book on differential equations. The solution of diffusion equations has fascinated academic elec-troanalytical chemists for years, and they naturally have a tendency to expound on them at the slightest provocation. Fortunately, the chemist using electrode reactions can accomplish a great deal without more than a cursory appreciation of the mathematics. Our intention here is to provide this qualitative appreciation on a level sufficient to understand laboratory techniques. 

This equation defines directly the change in concentration of the spedes AB with given concentrations of the reactants A and B, and the product AB. This is a differential equation whose solution is an expression of the form cAB=f(t, c% eg). The solution involves a process of integration, which is often difficult, and sometimes impossible, at least analytically. In such cases, numerical integration can be used to simulate the time-dependent variation of cAB in an experiment, enabling theoretical data to be obtained even for complex systems. 

The formulation of heat conduction problems for the determination of the one-dimj .nsional transient temperature distribution in a plane wall, a cylinder, or a sphefeTesults in a partial differential equation whose solution typically involves irtfinite series and transcendental equations, wliicli are inconvenient to use. Bijt the analytical soluliop provides valuable insight to the physical problem, hnd thus it is important to go through the steps involved. Below we demonstrate the solution procedure for the case of plane wall. 

Polymerization, polycondensation and polymer modification reactions are dynamic processes that are often represented by sets of differential equations. Analytical integration of such equations is often difficult or impossible and the solutions that are obtained often have complex forms that provide limited insight concerning the nature of the processes. This has been particular-... 

As we learned in this chapter, the formulation of unsteady distributed problems leads to partial differential equations. The solution of these equations is much more involved than that of ordinary differential equations. Among the techniques available, the analytical and computational methods are most frequently referred to. Exact analytical methods such as separation of variables and transform calculus are beyond the scope of the text. However, the method of complex temperature and the use of charts based on exact analytical solutions, being useful for some practical problems, are respectively discussed in Sections 3.4 and 3.6. Among approximate analytical methods, the integral method, already introduced in Sections 2.4 and 3.1, is further discussed in Section 3.5. The analog solution technique is also briefly treated in Section 3.7. 

The conventional data analysis involves the fitting of data to an equation describing the time dependence of the reaction, leading to the best estimates for the constants defining the equations. Analytical solutions to most simple reaction sequences can be obtained (7, 5, 63). Solutions of differential equations describing the series of first-order (or pseudo-first-order) reactions will always be a sum of exponential terms [Eq. (22)]. Thus for a single exponential, the fitting process provides the amplitude (A), the rate of reaction (X), and the end point (C)... 

In the previous section we solved linear ordinary differential equations analytically, obtaining general solutions in terms of the parameters in the equations. Numerical methods can also be used to obtain solutions, using a computer. In Chapter 1 we looked at the dynamic responses of several processes by using numerical integration methods (Euler integration-see Table 1.2). 

The reason why the random flight model has proved so popular theoretically stems from its simplicity, which offers hope for the development of analytic solutions. The problem can usually be cast in the form of a diffusionlike or a Schrodinger-wave-equation-like differential equation, the solutions of which are reasonably well explored. A tendency has developed in recent times to apply extremely sophisticated mathematical procedures to what are really very primitive models for polymer chains (see, e.g. Levine et al., 1978). Whether the ends merit the means in such instances cannot yet be assessed objectively. A strategy that might be more productive in terms of the development of a practical theory for steric stabilization is to aim for a simpler mathematical description of more complex models of polymer chains. It should also be borne in mind in developing ab initio theories that a simple model that may well suffice in polymer solution thermodynamics may be quite inadequate for the simulation of the conformational properties of polymers. Polymer solution thermodynamics seem to be relatively insensitive to molecular architecture per se whereas the conformation of a polymer chain is extremely sensitive to it. 

For this problem, we can solve the original differential equation analytically. The exact solution is... 

Once the differential equations and mass balance have been written down, three approaches can be followed in order to model complex reaction schemes. These are (1) numerical integration of differential equations, (2) steady-state approximations to solve differential equations analytically, and (3) exact analytical solutions of the differential equations without using approximations. 

The finite element solution of differential equations requires function integration over element domains. Evaluation of integrals over elemental domains by analytical methods can be tedious and impractical and is not attempted in... 

The comparison between the finite element and analytical solutions for a relatively small value of a - 1 is shown in Figure 2.25. As can be seen the standard Galerkin method has yielded an accurate and stable solution for the differential Equation (2.80). The accuracy of this solution is expected to improve even further with mesh refinement. As Figmre 2.26 shows using a = 10 a stable result can still be obtained, however using the present mesh of 10 elements, for larger values of this coefficient the numerical solution produced by the standard... 

After the aussembly of elemental equations into a global set and imposition of the boundary conditions the final solution of the original differential equation with respect to various values of upwinding parameter jS can be found. The analytical solution of Equation (2.80) with a = 50 is found as... 

In applying quantum mechanics to real chemical problems, one is usually faced with a Schrodinger differential equation for which, to date, no one has found an analytical solution. This is equally true for electronic and nuclear-motion problems. It has therefore proven essential to develop and efficiently implement mathematical methods which can provide approximate solutions to such eigenvalue equations. Two methods are widely used in this context- the variational method and perturbation theory. These tools, whose use permeates virtually all areas of theoretical chemistry, are briefly outlined here, and the details of perturbation theory are amplified in Appendix D. 

---