Partial Differential Equation systems

Note that this equation still retains the radial coordinate r. Therefore, unlike wedge case, there is not a unique ordinary differential that applies at any radius. Rather, there is an ordinary differential for every r position. Such local similarity behavior certainly represents a simplification compared to the original partial-differential-equation system. Nevertheless, the differential equation is more complex than that for the wedge case. 

The laminar flow assumption eliminates the non-linear term in the partial differential equations system (3.3), thus significantly reducing the computational cost. In addition, the present formulation often admits an exact solution. For example, in the case of an incompressible 2D laminar flow between two motionless parallel plates (i.e. planar SOFC configuration of Figure 3.1), Equation (3.29) reduces to ... 

Discretization of the partial differential equation system in axial (z) and radial (r) direction by means of the orthogonal collocation method (7) leads to the following system of ordinary differential equations. 

The conditions specified by Eq. (6.206) provide the conditions required to design the model, also called similarity requirements or modeling laws. The same analysis could be carried out for the governing differential equations or the partial differential equation system that characterize the evolution of the phenomenon (the conservation and transfer equations for the momentum). In this case the basic theorem of the similitude can be stipulated as A phenomenon or a group of phenomena which characterizes one process evolution, presents the same time and spatial state for all different scales of the plant only if, in the case of identical dimensionless initial state and boundary conditions, the solution of the dimensionless characteristic equations shows the same values for the internal dimensionless parameters as well as for the dimensionless process exits . 

The partial differential equations system for steady flows of Maxwell type (i.e., with = 0) is of composite type, neither elliptic, nor hyperbolic. This is not surprising, the same being true for instance for the stationary system of ideal incompressible fluids. The new feature, discovered in [8], is that some change of type may occur. In fact an easy but tedious calculation shows that three types of characteristics axe present ... 

Here we choose to apply the GITT on equations (23) in the partial transformation sfa-ategy, resulting in tiie parabolic partial differential equations system below ... 

It should be noted that in the case of the high intensive fields the CNT conduction band becomes highly populated and the contribution from the intraband motion of r-electrons to the total SWNT absorption can be essential and can lead to the so-called nonsaturated component observed in [1]. To study this case, the numerical solution of the partial differential equations system is required. 

Finally it should be mentioned that a direct numerical solution of the transport equations will allow us to obtain more exact quantitative results. This way of solving the complex partial differential equation system is not trivial. First results are obtained by MacLeod Radke... 

Generally speaking the mathematical problems tackled in voltammetry involve the resolution of partial differential equation systems by means of analytical, semi-analytical or numerical methods. The solutions of the problem are the concentration profiles of the different species, and from them the current-potential-time response of the system to a given electrical perturbation can be calculated. 

Reacting flow and its scalar field evolve in time according to an JVu-dimensional partial differential equation system, where Nu is the number of unknown chemical and physical variables (See Equation (1) for comparison) ... 

The numerical solution is performed by the method of lines. Spatial discretization of the partial differential-equation system using finite differences on statically adapted grids leads to large systems of ordinary differential and algebraic equations. This system of coupled equations is solved by an implicit extrapolation method using the software package LIMEX [14]. The code computes species mass-fraction and temperature profiles in the gas phase, fluxes at the gas-surface interface, and surface temperature and coverage as function of time. 

For a static equilibrium problem the partial differential equation system in Eulerian form together with the boundary conditions is given by... 

The total Lagrangian form of the equation of equilibrium is obtained by differentiating (2.116) directly. Thus the partial differential equation system together with the... 

Ortner et al. (20) attempted to extend the Lorbach and Marr model to countercurrent column operation. This resulted in a nonlinear partial differential equation system, which was complicated and was solved numerically. 

If the chemical system has a single unstable steady state and hence shows oscillatory behavior in the absence of starch, complex formation can stabilize the homogeneous steady state and make possible the appearance of Turing structures at parameters which would yield oscillatory kinetics in the complex-free system. Observe that in the above partial differential equation system (12) the effective ratio of diffusion coefficients is (1 -f K )c, which can be much greater than unity even if c < 1. Consequently, the presence of a species that forms an appropriate complex with the activator can allow Turing structures to form for, in principle, any ratio of the activator and inhibitor diffusion coefficients. 

The solutions of such partial differential equations require infomiation on the spatial boundary conditions and initial conditions. Suppose we have an infinite system in which the concentration flucPiations vanish at the infinite boundary. If, at t = 0 we have a flucPiation at origin 5C(f,0) = AC (f), then the diflfiision equation... 

In this paper, we discuss semi-implicit/implicit integration methods for highly oscillatory Hamiltonian systems. Such systems arise, for example, in molecular dynamics [1] and in the finite dimensional truncation of Hamiltonian partial differential equations. Classical discretization methods, such as the Verlet method [19], require step-sizes k smaller than the period e of the fast oscillations. Then these methods find pointwise accurate approximate solutions. But the time-step restriction implies an enormous computational burden. Furthermore, in many cases the high-frequency responses are of little or no interest. Consequently, various researchers have considered the use of scini-implicit/implicit methods, e.g. [6, 11, 9, 16, 18, 12, 13, 8, 17, 3]. 

Mathematical and Computational Implementation. Solution of the complex systems of partial differential equations governing both the evolution of pollutant concentrations and meteorological variables, eg, winds, requires specialized mathematical techniques. Comparing the two sets of equations governing pollutant dynamics (eq. 5) and meteorology (eqs. 12—14) shows that in both cases they can be put in the form ... 

Smith, I. M., J. L. Siemienivich, and I. Gladweh. A Comparison of Old and New Methods for Large Systems of Ordinary Differential Equations Arising from Parabolic Partial Differential Equations, Num. Anal. Rep. Department of Engineering, no. 13, University of Manchester, England (1975). 

An alternative procedure is the dynamic programming method of Bellman (1957) which is based on the principle of optimality and the imbedding approach. The principle of optimality yields the Hamilton-Jacobi partial differential equation, whose solution results in an optimal control policy. Euler-Lagrange and Pontrya-gin s equations are applicable to systems with non-linear, time-varying state equations and non-quadratic, time varying performance criteria. The Hamilton-Jacobi equation is usually solved for the important and special case of the linear time-invariant plant with quadratic performance criterion (called the performance index), which takes the form of the matrix Riccati (1724) equation. This produces an optimal control law as a linear function of the state vector components which is always stable, providing the system is controllable. 

Partieulate produets, sueh as those from eomminution, erystallization, preeipi-tation ete., are distinguished by distributions of the state eharaeteristies of the system, whieh are not only funetion of time and spaee but also some properties of states themselves known as internal variables. Internal variables eould inelude size and shape if partieles are formed or diameter for liquid droplets. The mathematieal deseription eneompassing internal eo-ordinate inevitably results in an integro-partial differential equation ealled the population balanee whieh has to be solved along with mass and energy balanees to deseribe sueh proeesses. 

Coupled-map Lattices. Another obvious generalization is to lift the restriction that sites can take on only one of a few discrete values. Coupled-map lattices are CA models in which continuity is restored to the state space. That is to say, the cell values are no longer constrained to take on only the values 0 and 1 as in the examples discussed above, but can now take on arbitrary real values. First introduced by Kaneko [kaneko83]-[kaneko93], such systems are simpler than partial differential equations but more complex than generic CA. Coupled-map lattices are discussed in chapter 8. 

LGs can also serve as powerful alternatives to PDEs themselves in modeling physical systems. The distinction is an important one. It must be remembered, however, that not all PDEs (and perhaps not all physical systems see chapter 12) are amenable to a LG simulation. Moreover, even if a candidate PDE is selected for simulation by a LG. there is no currently known cookbook recipe allowing a researcher to go from the PDE to a LG description (or vice versa). Nonetheless, by their very nature, LGs lend themselves to modeling any partial differential equation (PDE) for which the underlying physical basis for its construction involves a large number of particles with local interactions [wolf86c]. 

The overall set of partial differential equations that can be considered as a mathematical characterization of the processing system of gas-liquid dispersions should include such environmental parameters as composition, temperature, and velocity, in addition to the equations of bubble-size and residence-time distributions that describe the dependence of bubble nucleation and growth on the bubble environmental factors. A simultaneous solution of this set of differential equations with the appropriate initial and boundary conditions is needed to evaluate the behavior of the system. Subject to the Curie principle, this set of equations should include the possibilities of coupling effects among the various fluxes involved. In dispersions, the possibilities of couplings between fluxes that differ from each other by an odd tensorial rank exist. (An example is the coupling effect between diffusion of surfactants and the hydrodynamics of bubble velocity as treated in Section III.) As yet no analytical solution of the complete set of equations has been found because of the mathematical difficulties involved. To simplify matters, the pertinent transfer equation is usually solved independently, with some simplifying assumptions. 

In Schrodinger s wave mechanics (which has been shown4 to be mathematically identical with Heisenberg s quantum mechanics), a conservative Newtonian dynamical system is represented by a wave function or amplitude function [/, which satisfies the partial differential equation... 

Unlike stirred tanks, piston flow reactors are distributed systems with one-dimensional gradients in composition and physical properties. Steady-state performance is governed by ordinary differential equations, and dynamic performance is governed by partial differential equations, albeit simple, first-order PDEs. Figure 14.6 illustrates a component balance for a differential volume element. 

---