Computer, solving differential equations with

There were not many computer-related developments before World War II. However, in the 1930s V. Bush, an electrical engineer, developed the first analog computer that would solve differential equations with up to 18 independent variables. This was a forerunner of the electronic computers developed in the late 1940s after World War II. 

G. Fasshauer. Solving differential equations with RBFs Multilevel methods and smoothing . Advances in Computational Mathematics 11 (1999), 139-159. 

It is impossible to solve these partial differential equations with the varied input conditions resulting from the California cycle. A computer was used in the study, where distance is divided into increments of 0.08 to 0.30 in. in thickness, and time is divided into increments of 0.4 sec. They obtained the temperature profile in the bed as a function of time from a cold start. This procedure could result in a very long computation, as the California cycle has a duration of 16 min, requiring 2400 time intervals and computations. 

Dynamic simulations are also possible, and these require solving differential equations, sometimes with algebraic constraints. If some parts of the process change extremely quickly when there is a disturbance, that part of the process may be modeled in the steady state for the disturbance at any instant. Such situations are called stiff, and the methods for them are discussed in Numerical Solution of Ordinary Differential Equations as Initial-Value Problems. It must be realized, though, that a dynamic calculation can also be time-consuming and sometimes the allowable units are lumped-parameter models that are simplifications of the equations used for the steady-state analysis. Thus, as always, the assumptions need to be examined critically before accepting the computer results. 

CFD may be loosely thought of as computational methods applied to the study of quantities that flow. This would include both methods that solve differential equations and finite automata methods that simulate the motion of fluid particles. We shall include both of these in our discussions of the applications of CFD to packed-tube simulation in Sections III and IV. For our purposes in the present section, we consider CFD to imply the numerical solution of the Navier-Stokes momentum equations and the energy and species balances. The differential forms of these balances are solved over a large number of control volumes. These small control volumes when properly combined form the entire flow geometry. The size and number of control volumes (mesh density) are user determined and together with the chosen discretization will influence the accuracy of the solutions. After boundary conditions have been implemented, the flow and energy balances are solved numerically an iteration process decreases the error in the solution until a satisfactory result has been reached. 

In this text all numerical problems involve integration of simultaneous ordinary differential equations or solution of simultaneous algebraic equations. You should have no trouble finding ways to solve algebraic equations with a calculator, a spreadsheet, a personal computer, etc. 

In this section we deal with estimating the parameters p in the dynamical model of the form (5.37). As we noticed, methods of Chapter 3 directly apply to this problem only if the solution of the differential equation is available in analytical form. Otherwise one can follow the same algorithms, but solving differential equations numerically whenever the computed responses are needed. The partial derivations required by the Gauss - Newton type algorithms can be obtained by solving the sensitivity equations. While this indirect method is... 

A reviewer has kindly provided references to other programs for solving differential equations by numerical methods (13-21). These will be useful to individuals that do not have IBM equipment but would like to do calculations of the sort outlined above. The general topic has been discussed (13). Sebastian, et al. (14) discuss DPS (Dynamic Process Simulator) and provide references to descriptions of MIMIC (15), ACSL (16), ISIS (17), BEDSOCS (18), DSL/77 (19), DARE (20) and PMSP (21). References to other programs can be found in papers that accompany reference 2 and in other Proceedings of Conferences on Applications of Continuous System Simulation. Because of the simplicity and power of system simulation programs such as are discussed herein, one could expect that every major computer producer would have available software with the capability of CSMP. 

Owing to the possibility to solve the governing differential equations with numerical methods, we can in principle compute future system reactions, provided that initial... 

Solutions of assemblies of ordinary differential equations with time as the independent variable are ideally suited for solution by analog computation. Hence complex kinetics equations of the type considered in this section may conveniently be solved with an analog computer. This is illustrated in the problems at the end of this chapter. 

Note that if the concentrations and reaction rates have already been computed and stored, then this is a set of linear differential equations (with non-constant coefficients) for the sensitivity coefficients. Note also that each sensitivity coefficient depends only on the sensitivity coefficients of the other species to the same parameter, but not on sensitivity coefficients with respect to other parameters. Thus, the sensitivity coefficients with respect to a given parameter (K of them) are coupled and must be computed simultaneously, but it is not necessary to solve for all the sensitivity coefficients (2 x / x of them) simultaneously. 

Wang, M.L. Liou, C.T., and Chang, R.Y., Numerical technique for solving partial differential equations with applications to adsorption process, Comput. Chem. Eng., 4(2), 85-92 (1980). 

There are many books on numerical methods available that contain exercises that allow you to practice writing yoin own algorithms to solve differential equations. It is equally important to learn how to use the existing software products for numerical analysis, and select appropriate numerical methods for stability and efficiency reasons. With this in mind, some of the problems in this section require that you work with software for munerical computing, e.g. MATLAB. 

Computer simulation of the reactor kinetic hydrodynamic and transport characteristics reduces dependence on phenomenological representations and idealized models and provides visual representations of reactor performance. Modem quantitative representations of laminar and turbulent flows are combined with finite difference algorithms and other advanced mathematical methods to solve coupled nonlinear differential equations. The speed and reduced cost of computation, and the increased cost of laboratory experimentation, make the former increasingly usehil. 

In general, comprehensive, multidimensional modeling of turbulent combustion is recognized as being difficult because of the problems associated with solving the differential equations and the complexities involved in describing the interactions between chemical reactions and turbulence. A number of computational models are available commercially that can do such work. These include FLUENT, FLOW-3D, and PCGC-2. 

With the introduction of Gear s algorithm (25) for integration of stiff differential equations, the complete set of continuity equations describing the evolution of radical and molecular species can be solved even with a personal computer. Many models incorporating radical reactions have been pubHshed. 

Rigorous error bounds are discussed for linear ordinary differential equations solved with the finite difference method by Isaacson and Keller (Ref. 107). Computer software exists to solve two-point boundary value problems. The IMSL routine DVCPR uses the finite difference method with a variable step size (Ref. 247). Finlayson (Ref. 106) gives FDRXN for reaction problems. 

The development of mathemafical models is described in several of the general references [Giiiochon et al., Rhee et al., Riithven, Riithven et al., Suzuki, Tien, Wankat, and Yang]. See also Finlayson [Numerical Methods for Problems with Moving Front.s, Ravenna Park, Washington, 1992 Holland and Liapis, Computer Methods for Solving Dynamic Separation Problems, McGraw-Hill, New York, 1982 Villadsen and Michelsen, Solution of Differential Equation Models by... 

---