General Numerical Simulation

The second classification of simulation methods is the set of Monte Carlo techniques, which refer to very general numerical simulations that have been applied to a wide variety of problems. The Monte Carlo methods are not able... 

In general, numerical simulation using full Boltzmann equation suggests that the BGK model is accurate for isothermal flows. However, for nonisothermal flows, corrections for the Prandtl number (for collision frequency) need to be introduced. 

Although direct numerical simulations under limited circumstances have been carried out to determine (unaveraged) fluctuating velocity fields, in general the solution of the equations of motion for turbulent flow is based on the time-averaged equations. This requires semi-... 

Numerical simulations are designed to solve, for the material body in question, the system of equations expressing the fundamental laws of physics to which the dynamic response of the body must conform. The detail provided by such first-principles solutions can often be used to develop simplified methods for predicting the outcome of physical processes. These simplified analytic techniques have the virtue of calculational efficiency and are, therefore, preferable to numerical simulations for parameter sensitivity studies. Typically, rather restrictive assumptions are made on the bounds of material response in order to simplify the problem and make it tractable to analytic methods of solution. Thus, analytic methods lack the generality of numerical simulations and care must be taken to apply them only to problems where the assumptions on which they are based will be valid. 

In general, discontinuities constitute a problem for numerical methods. Numerical simulation of a blast flow field by conventional, finite-difference schemes results in a solution that becomes increasingly inaccurate. To overcome such problems and to achieve a proper description of gas dynamic discontinuities, extra computational effort is required. Two approaches to this problem are found in the literature on vapor cloud explosions. These approaches differ mainly in the way in which the extra computational effort is spent. 

Blast effects can be represented by a number of blast models. Generally, blast effects from vapor cloud explosions are directional. Such effects, however, cannot be modeled without conducting detailed numerical simulations of phenomena. If simplifying assumptions are made, that is, the idealized, symmetrical representation of blast effects, the computational burden is eased. An idealized gas-explosion blast model was generated by computation results are represented in Figure 4.24. Steady flame-speed gas explosions were numerically simulated with the BLAST-code (Van den Berg 1980), and their blast effects were calculated. 

In this section we briefly outline a general mean-field theory approach to arbitrary PCA and then apply the formalism to a particular class of one-parameter rules. We then compare the theoretical predictions to numerical simulations on lattices of dimension 1 < d < 4. 

The analytic theory outlined above provides valuable insight into the factors that determine the efficiency of OI.EDs. However, there is no completely analytical solution that includes diffusive transport of carriers, field-dependent mobilities, and specific injection mechanisms. Therefore, numerical simulations have been undertaken in order to provide quantitative solutions to the general case of the bipolar current problem for typical parameters of OLED materials [144—1481. Emphasis was given to the influence of charge injection and transport on OLED performance. 1. Campbell et at. [I47 found that, for Richardson-Dushman thermionic emission from a barrier height lower than 0.4 eV, the contact is able to supply... 

Proc. of the International summer school on Experimental physics of gravitational waves, (Barone, M. et al. Eds., World Scientihc, London 2000). Contains a valuable chapter on General relativity by P. Tourrenc (contains a precise description of the various coordinates systems and their use, OBLIGATORY), a chapter by S. Bonazzola and E. Gourgoulhon on compact sources, in particular neutron stars, and a chapter by Jean-Yves Vinet on numerical simulations of interferometric gw detectors. 

The modern discipline of Materials Science and Engineering can be described as a search for experimental and theoretical relations between a material s processing, its resulting microstructure, and the properties arising from that microstructure. These relations are often complicated, and it is usually difficult to obtain closed-form solutions for them. For that reason, it is often attractive to supplement experimental work in this area with numerical simulations. During the past several years, we have developed a general finite element computer model which is able to capture the essential aspects of a variety of nonisothermal and reactive polymer processing operations. This "flow code" has been Implemented on a number of computer systems of various sizes, and a PC-compatible version is available on request. This paper is intended to outline the fundamentals which underlie this code, and to present some simple but illustrative examples of its use. 

The partial differential equations describing the catalyst particle are discretized with central finite difference formulae with respect to the spatial coordinate [50]. Typically, around 10-20 discretization points are enough for the particle. The ordinary differential equations (ODEs) created are solved with respect to time together with the ODEs of the bulk phase. Since the system is stiff, the computer code of Hindmarsh [51] is used as the ODE solver. In general, the simulations progressed without numerical problems. The final values of the rate constants, along with their temperature dependencies, can be obtained with nonlinear regression analysis. The differential equations were solved in situ with the backward... 

Generally, the results obtained through the numerical simulation showed good agreement with the experimental data leading to the conclusion that CFD techniques can be effectively used in consequence assessment procedures concerning toxic/flammable dispersion scenarios in real terrains, where box models have limited capabilities. 

Chapter 8 presented the last of the computational approaches that I find widely useful in the numerical simulation of environmental properties. The routines of Chapter 8 can be applied to systems of several interacting species in a one-dimensional chain of identical reservoirs, whereas the routines of Chapter 7 are a somewhat more efficient approach to that chain of identical reservoirs that can be used when there is only one species to be considered. Chapter 7 also presented subroutines applicable to a generally useful but simple climate model, an energy balance climate model with seasonal change in temperature. Chapter 6 described the peculiar features of equations for changes in isotope ratios that arise because isotope ratios are ratios and not conserved quantities. Calculations of isotope ratios can be based directly on calculations of concentration, with essentially the same sources and sinks, provided that extra terms are included in the equations for rates of change of isotope ratios. These extra terms were derived in Chapter 6. 

In general, it has to be noticed that data from experiments have to be compared with data from numerical simulation with great care and prudence. As a matter of fact, most of the numerical codes describe propagation of the ultrashort laser pulse (1) without any pre-pulse and (2) in a medium already ionized. These two conditions are usually far from being fulfilled in a real experiment. In the following, we will show how precursors of the main laser pulse and ionization of the medium can deeply change the propagation of an ultrashort intense laser pulse. 

The control strategies for determining the feed policies were decided on the basis of a numerical solution of the terpolymerisations described by equations 1 - 3 using a microcomputer and a general purpose simulation package, BEEBSOC (10). Where necessary, these data were acquired in the course of this study, otherwise literature values were used. The apparent first order rate constants in terpolymerisations have been shown to be composition dependent. The variation in rate constants with... 

The Stratonovich SDEs for either generalized or Cartesian coordinates could be numerically simulated by implementing the midstep algorithm of Eq. (2.238). Evaluation of the required drift velocities would, however, require the evaluation of sums of derivatives of B or whose values will depend on the decomposition of the mobility used to dehne these quantities. This provides a worse starting point for numerical simulation than the forward Euler algorithm interpretation. 

---