Numerical Tools

The book is intended to provide a foundation for the senior-level or first-year graduate student interested in shock-compression science, and to expose the student to the basic experimental, theoretical, and numerical tools prevailing in the field. Students should find the introductory material presented to be useful in preparation for more advanced studies of shock compression. Several problems are also presented with the different chapters to aid the student in understanding the material presented. Those interested in pursuing this field are encouraged to understand these problems in depth. 

For example, the rate constant of the collinear reaction H -f- H2 has been calculated in the temperature interval 200-1000 K. The quantum correction factor, i.e., the ratio of the actual rate constant to that given by CLTST, has been found to reach 50 at T = 200 K. However, in the reactions that we regard as low-temperature ones, this factor may be as large as ten orders of magnitude (see introduction). That is why the present state of affairs in QTST, which is well suited for flnding quantum contributions to gas-phase rate constants, does not presently allow one to use it as a numerical tool to study complex low-temperature conversions, at least without further approximations such as the WKB one. ... 

These equations are nonlinear and cannot be solved analytically. They are included in this section because they are autocatalytic and because this chapter discusses the numerical tools needed for their solution. Figure 2.6 illustrates one possible solution for the initial condition of 100 rabbits and 10 lynx. This model should not be taken too seriously since it represents no known chemistry or... 

The Stribeck curve gives a general description for the transition of lubrication regime, but the quantitative information, such as the variations of real contact areas, the percentage of the load carried by contact, and changes in friction behavior, are not available due to lack of numerical tools for prediction. The deterministic ML model provides an opportunity to explore the entire process of transition from full-film EHL to boundary lubrication, as demonstrated by the examples presented in this section. 

In Section 1.4, the MCT was introduced in general terms as an important numerical tool for studying the relationship of variables in complex systems of equations. Here, the algorithm will be presented in detail, and a more complex example will be worked. 

The fact that self-interaction errors are canceled exactly in HF calculations suggests that a judicious combination of an HF-like approach for localized states with DFT for everything else may be a viable approach for strongly correlated electron materials. This idea is the motivation for a group of methods known as DFT+U. The usual application of this method introduces a correction to the DFT energy that corrects for electron self-interaction by introducing a single numerical parameter, U — J, where U and J involve different aspects of self-interaction. The numerical tools needed to use DFT+U are now fairly widely implemented in plane-wave DFT codes. 

In this central result the choice of the point q(0) is arbitrary. This means that at time t = 0 one can initiate trajectories anywhere and after a short induction time the reactive flux will reach a plateau value, which relaxes exponentially, but at a very slow rate, It is this independence on the initial location which makes the reactive flux method an important numerical tool. 

The time-independent and time-dependent approaches merely provide different views of the dissociation process and different numerical tools for the calculation of photodissociation cross sections. The time-independent approach is a boundary value problem, i.e., the stationary wavefunction... 

With the advent of nanomaterials, different types of polymer-based composites developed as multiple scale analysis down to the nanoscale became a trend for development of new materials with new properties. Multiscale materials modeling continue to play a role in these endeavors as well. For example, Qian et al. [257] developed multiscale, multiphysics numerical tools to address simulations of carbon nanotubes and their associated effects in composites, including the mechanical properties of Young s modulus, bending stiffness, buckling, and strength. Maiti [258] also used multiscale modeling of carbon nanotubes for microelectronics applications. Friesecke and James [259] developed a concurrent numerical scheme to evaluate nanotubes and nanorods in a continuum. 

The derivations made so far were done in terms of the wave functions (vectors). In fact all the basic tools of the quantum theory used throughout this book are covered by this brief account. However, it may be practical to have a variety of representations for the same set of basic techniques in different incarnations. Not adding too much either to pragmatic numerical tools or to a deeper understanding of what is going on, these tools are useful for getting general relations which are an important part of the present book and for a more economical representation of the variables of the problems considered here. For that reason we review them below. 

M is a singular Matrix. Zero entries on the main diagonal of this matrix identify the algebraic equations, and all other entries which have the value 1 represent the differential equation. The vector x describes the state of the system. As numeric tools for the solution of the DAE system, MATLAB with the solver odel5s was used. In this solver, a Runge Kutta procedure is coupled with a BDF procedure (Backward Difference Formula). An implicit numeric scheme is used by the solver. 

Three-dimensional models can handle complex cases with fewer empirical features. They require the solution of multidimensional, partial differential equations. These solutions are numerical tools and require the use of a supercomputer, or very long computation times on less sophisticated hardware. 

At this level, most of what can be done with numerical tools and simple laboratory equipment has been done. We must now check this design on a laboratory pilot plant. Experience shows that the numerical prediction enables us to get a set of flow rates which are very close to the adequate ones. 

A complete optimization of the process would require an intensive use of numerical tools to investigate the precise influence of the number of columns, number of columns per zone, column length, etc. on an defined economical function. 

The previous chapters have discussed examples of instabilities in combustors where the fuel was gaseous. However, in most practical devices, fuel is liquid (gasoline or kerosene) and the problem becomes much more difficult. In this chapter, the specificities of two-phase flow combustion will be discussed and the construction of a numerical tool to perform LES of liquid fuel combustion will be discussed. This chapter will also present the equations solved for gaseous and liquid combustion in more details than the previous ones. 

The EE approach has the important advantage to be straightforward to implement in a numerical tool, and immediately efficient as it allows the use of the same parallel algorithm than for the gas phase [301]. However it requires an initial modeling effort much larger than in the EL method [280] and faces difficulties in handling droplet clouds with extended size distributions. Moreover the resulting set of equations is numerically difficult to handle and requires special care [346]. 

---