Temperature field, solutal model

The basic scheme for the numerical solution is the same as that used for the 1 -D model, except that in this case the solid temperature field used to solve the DAE system for each monolith channel must be calculated from the three-dimensional solid-phase energy balance equation. The three-dimensional energy balance equation can be solved by a nonlinear finite element solver (such as ABAQUS) for the solid-phase temperature field while a nonlinear finite difference solver for the DAE system calculates the gas-phase temperature and... 

The solution of the gas flow and temperature fields in the nearnozzle region (as described in the previous subsection), along with process parameters, thermophysical properties, and atomizer geometry parameters, were used as inputs for this liquid metal breakup model to calculate the liquid film and sheet characteristics, primary and secondary breakup, as well as droplet dynamics and cooling. The trajectories and temperatures of droplets were calculated until the onset of secondary breakup, the onset of solidification, or the attainment of the computational domain boundary. This procedure was repeated for all droplet size classes. Finally, the droplets were numerically sieved and the droplet size distribution was determined. 

Thus in the ideal solution model of a two phase field, a knowledge of the enthalpies of fusion of the pure components at their respective melting points allows simultaneous solution of these two equations for the two unknowns, NA(l> and NA(aj at the temperature of interest. 

To finalize the development of the aqueous CO2 force field parameters, the C02 model was used in free energy perturbation Monte Carlo (FEP/MC) simulations to determine the solubility of C02 in water. The solubility of C02 in water is calculated as a function of temperature in the development process to maintain transferability of the C02 model to different simulation techniques and to quantify the robustness of the technique used in the solubility calculations. It is also noted that the calculated solubility is based upon the change in the Gibbs energy of the system and that parameter development must account for the entropy/enthalpy balance that contributes to the overall structure of the solute and solvent over the temperature range being modeled [17]. 

The difference between RANS and LES is depicted in Figure 20.1, which shows the temperature fields of a pool fire flame. While the RANS result shows smooth variations and looks like a laminar flame, the LES result clearly illustrates the large-scale eddies. Both results are the correct solutions of the corresponding equations. However, the time accuracy of LES is also essential for the quantitative accuracy of the buoyancy-driven flows. As Rehm and Baum have shown [10], the dynamic motions or eddies are responsible for most of the air entrainment into the fire plumes. Because these motions cannot be captured by RANS, LES is usually better suited for fire-driven flow. LES typically requires a finer spatial resolution than RANS. Examples of RANS-based fire CFD models are JASMINE, KAMELEON [11], SMARTFIRE [12], SOFIE [13], ISIS [14], and ISIS-3D [15]. Examples of LES models are the FDS [4,5] and SMAFS [16], developed at Lund University. Fire simulations using LES have also been performed by Cheung et al. [17] and Gao et al. [18],... 

The mathematical model comprises a set of partial differential equations of convective diffusion and heat conduction as well as the Navier-Stokes equations written for each phase separately. For the description of reactive separation processes (e.g. reactive absorption, reactive distillation), the reaction terms are introduced either as source terms in the convective diffusion and heat conduction equations or in the boundary condition at the channel wall, depending on whether the reaction is homogeneous or heterogeneous. The solution yields local concentration and temperature fields, which are used for calculation of the concentration and temperature profiles along the column. 

As mentioned above, the nonequilibrium radiation code NEQAIR is employed for prediction of ultraviolet emission from the DSMC flow field solutions. The modeling of ultraviolet emission with this code is discussed for nitric oxide in Ref. 84 and for atomic oxygen in Ref. 87. A common assumption made in using the NEQAIR code is that a quasisteady state (QSS) exists for the number densities of the electronically excited species. The assumption requires that the time scale of chemical processes is much smaller than the time scales for diffusion and for changes in overall properties. Under these conditions, the local values of temperatures and ground state species number densities obtained from the DSMC computation may be used to compute the populations of the electronically excited states. 

The temperature fields induced by microwave heating can be modeled via the simultaneous solution of Maxwell s equations (for the electromagnetic component of the problem) and the heat equation. The modeling is very challenging, in part because the dielectric constants sj. and s" are functions of temperature and the microwave frequency as well as of the microstructural and chemical details of the ceramic [Eq. (lb)]. Also, the thermal conductivity, k (which is needed for the heat transfer calculations), typically is a function of temperature as well as of microstructural variables such as porosity. ... 

The solution of Eq. 8.25 is wavelike, suggesting that the temperature field propagates as a wave. The speed of propagation of this wave is equal to /klCt which also happens to be the speed of the energy carrier, for example, the speed of sound for phonons. So, this model is nonlocal in time but local in space since the temperature represents a spatially localized thermodynamic equilibrium. 

In the case of concentrated reactant solutions, we can observe a sharp temperature increase during acidic media neutralisation. The temperature field in the reaction zone can be adjusted by using small tubular turbulent reactors of cylinder or diffuser-confusor design. There are several options for temperature adjustment [16], e.g., the change of the device radius and reactant flow rate, allows the application of the zone model for fast chemical processes and the use of shell-and-tube reactors with pipe columns of small radius comprising the same reactor cross section in total [17],... 

FD method is applied to solve the temperature field. In traditional ordinary FD method, the 3D model is discretized into a 3D matrix, whose size depends on the smallest circumscribed box of all the geometry. Fig. 2 schematically shows that the vacuum elements in white occupy memory as well, but they are never used in the solution. 

---