Continuity equation numerical solution

Computational fluid dynamics (CFD) emerged in the 1980s as a significant tool for fluid dynamics both in research and in practice, enabled by rapid development in computer hardware and software. Commercial CFD software is widely available. Computational fluid dynamics is the numerical solution of the equations or continuity and momentum (Navier-Stokes equations for incompressible Newtonian fluids) along with additional conseiwation equations for energy and material species in order to solve problems of nonisothermal flow, mixing, and chemical reaction. 

Practical problems involving variable-density PFRs require numerical solutions, and for these it is better to avoid expanding Equation (3.4) into separate derivatives for a and u. We could continue to use the molar flow rate, Na, as the dependent variable, but prefer to use the molar flux,... 

Numerical Solution. The momentum Equation 5 is solved simultaneously along with the energy Equation 6 to obtain axial velocity, v, and temperature fields. The continuity equation with the known axial velocity is used to... 

The objective of most of the theories of transport in porous media is to derive analytical or numerical functions for the effective diffusion coefficient to use in the preceed-ing averaged species continuity equations based on the structure of the media and, more recently, the structure of the solute. 

While offering a more inherently realistic method of solution, however, the technique may cause some additional problems in the numerical solution, since high values of Kl can lead to increased stiffness in the differential equations. Thus in using this technique, a compromise between the approach to equilibrium and the speed of numerical solution may have to be adopted. Continuous single-stage extraction is treated in the simulation example EQEX. Reaction with integrated extraction is demonstrated in simulation example REXT. 

Another kind of situation arises when it is necessary to take into account the long-range effects. Here, as a rule, attempts to obtain analytical results have not met with success. Unlike the case of the ideal model the equations for statistical moments of distribution of polymers for size and composition as well as for the fractions of the fragments of macromolecules turn out normally to be unclosed. Consequently, to determine the above statistical characteristics, the necessity arises for a numerical solution to the material balance equations for the concentration of molecules with a fixed number of monomeric units and reactive centers. The difficulties in solving the infinite set of ordinary differential equations emerging here can be obviated by switching from discrete variables, characterizing macromolecule size and composition, to continuous ones. In this case the mathematical problem may be reduced to the solution of one or several partial differential equations. 

In this case, an exact analytical solution of the continuity equations for A and B does not exist. An approximate solution has been developed by Van Krevelen and Hoftijzer (1948) in terms of E. Results of a numerical solution could be fitted approximately by the implicit relation... 

Figure 5. Exact (numerical solution, continuous line) and linearised (equation (24), dotted line) velocity profile (i.e. vy of the fluid at different distances x from the surface) at y = 10-5 m in the case of laminar flow parallel to an active plane (Section 4.1). Parameters Dt = 10 9m2 s-1, v = 10-3ms-1, and v = 10-6m2s-1. The hydrodynamic boundary layer thickness (<50 = 5 x 10 4 m), equation (26), where 99% of v is reached is shown with a horizontal double arrow line. For comparison, the normalised concentration profile of species i, ct/ithe linear profile of the diffusion layer approach (continuous line) and its thickness (<5, = 3 x 10 5m, equation (34)) have been added. Notice that the linearisation of the exact velocity profile requires that <5,

Equation (7.19) is a self-consistent equation for AEn, in the form of a sum of a pair of continued fractions (CFs). Although numerical solutions to (7.19) are feasible, we are only concerned with its qualitative features. In particular, we note that an exact WSL occurs when AEn = 0, which happens only if both CFs contain the same number of terms (apart from the trivial case 0 = 0). For the infinite chain, this situation is the case for every allowed energy, so an exact WSL is indeed found. But, for a finite chain, AEn = 0 only for the center state, which thus possesses the exact WSL energy. Therefore, the set of energies for a finite chain form only an approximate WSL. 

In the MO approach molecular orbitals are expressed as a linear combination of atomic orbitals (LCAO) atomic orbitals (AO), in return, are determined from the approximate numerical solution of the electronic Schrodinger equation for each of the parent atoms in the molecule. This is the reason why hydrogen-atom-like wavefunctions continue to be so important in quantum mechanics. Mathematically, MO-LCAO means that the wave-functions of the molecule containing N atoms can be expressed as... 

In this chapter we consider the performance of isothermal batch and continuous reactors with multiple reactions. Recall that for a single reaction the single differential equation describing the mass balance for batch or PETR was always separable and the algebraic equation for the CSTR was a simple polynomial. In contrast to single-reaction systems, the mathematics of solving for performance rapidly becomes so complex that analytical solutions are not possible. We will first consider simple multiple-reaction systems where analytical solutions are possible. Then we will discuss more complex systems where we can only obtain numerical solutions. 

The fundamental physical laws governing motion of and transfer to particles immersed in fluids are Newton s second law, the principle of conservation of mass, and the first law of thermodynamics. Application of these laws to an infinitesimal element of material or to an infinitesimal control volume leads to the Navier-Stokes, continuity, and energy equations. Exact analytical solutions to these equations have been derived only under restricted conditions. More usually, it is necessary to solve the equations numerically or to resort to approximate techniques where certain terms are omitted or modified in favor of those which are known to be more important. In other cases, the governing equations can do no more than suggest relevant dimensionless groups with which to correlate experimental data. Boundary conditions must also be specified carefully to solve the equations and these conditions are discussed below together with the equations themselves. 

The system of hyperbolic and parabolic partial differential equations representing the ID or 2D model of monolith channel is solved by the finite differences method with adaptive time-step control. An effective numerical solution is based on (i) discretization of continuous coordinates z, r and t, (ii) application of difference approximations of the derivatives, (iii) decomposition of the set of equations for Ts, T, c and cs, (iv) quasi-linearization of... 

The stationary system, e.g., (5.3.1) (5.3.5) is replaced by its time-dependent counterpart. In this counterpart, the Poisson equation is replaced by the total current continuity equation (1.5), obtained as a linear combination of the original equations. The resulting system is then solved by quasilinearization [9] with a simultaneous solution of quasilinearized equations and subsequent Newton s iterations at each time step. Integration is continued in time until the steady state is reached. This numerical procedure is a modification of that suggested by Mock in [10]. 

The first group of methods then manipulate a very small subset of vector elements Vi at a time, and a direct method continually updates the affected elements rj. Such methods are collectively known as relaxation methods, and they are primarily used in situations where, for each elements >< to be changed, the set of affected rjt and the matrix elements Ay, are immediately known. This applies in particular to difference approximations and also to the so-called Finite Element Method for obtaining tabular descriptions of the wave function, i.e. a list of values of the wave function at a set of electron positions. (Far some reason, such a description is commonly referred to as a numerical solution to the Schrodinger equation). Relaxation methods have also been applied to the Cl problem in the past, but due to their slow convergence they have been replaced by analytical methods. Even for numerical problems, the relaxation methods are slowly yielding to analytical methods. 

In all the above derivations in this section, the influence of viscosity is neglected so that analytical solutions for velocity and pressure profiles can be obtained. When the viscosity of fluid is taken into account, it is difficult to obtain any analytical solution. Kuts and Dolgushev [35] solved numerically the flow field in the impingement of two axial round jets of a viscous impressible liquid ejected at the same velocity from conduits with the same diameter and located very close to each other. The mathematical formulation incorporated the complete Navier-Stokes equations transformed into stream and velocity functions in cylindrical coordinates r and z, with the assumption that the velocity profiles at the entrance and the exit of the conduit were parabolic. The continuity equation is given by Eq. (1.22) and the equations for motion in dimensionless form are ... 

The integrated continuity equation is a weaker form of the full continuity equation. This is noticed in numerical solutions of mold filling problems, where continuity is never fully satisfied. However, this violation of continuity is insignificant and will not hinder the solution of practical mold filling problems. The integrated continuity equation reduces to... 

The LGA is a variant of a cellular automaton, introduced as an alternative numerical approximation to the partial differential equation of Navier-Stokes and the continuity equations, whose analytical solution leads to the macroscopic approach of fluid dynamics. The microscopic behavior of the LGA has been shown to be very close to the Navier-Stokes (N-S) equations for incompressible fluids at the macroscopic level. 

---