Delta function computation

In order to improve upon the mean-field approximation given in equation 7.112, we must somehow account for possible site-site correlations. Let us go back to the deterministic version of the basic Life rule (equation 7.110). We could take a formal expectation of this equation but we first need a way to compute expectation values of Kronecker delta functions. Schulman and Seiden [schul78] provide a simple means to do precisely that. We state their result without proof... 

Of all of the methods reviewed thus far in this book, only DNS and the linear-eddy model require no closure for the molecular-diffusion term or the chemical source term in the scalar transport equation. However, we have seen that both methods are computationally expensive for three-dimensional inhomogeneous flows of practical interest. For all of the other methods, closures are needed for either scalar mixing or the chemical source term. For example, classical micromixing models treat chemical reactions exactly, but the fluid dynamics are overly simplified. The extension to multi-scalar presumed PDFs comes the closest to providing a flexible model for inhomogeneous turbulent reacting flows. Nevertheless, the presumed form of the joint scalar PDF in terms of a finite collection of delta functions may be inadequate for complex chemistry. The next step - computing the shape of the joint scalar PDF from its transport equation - comprises transported PDF methods and is discussed in detail in the next chapter. Some of the properties of transported PDF methods are listed here. 

Figure 28. Geometry for computing the delta function response of the resist. (Reproduced with permission from Ref. 43)...

So much for the background. The question now is whether we can visualize the solution before going to the computer. This is particularly necessary as we have a delta function as initial condition. However, the equation looks rather like the equation for diffusion and flow in a straight tube, which, if to were the concentration of a solute, D( ) its spatially varying diffusion coefficient, and V(x) its spatially varying convective velocity, would be... 

Fig. 8. Generation of the form of the helical diffraction pattern. (A) shows that a continuous helical wire can be considered as a convolution of one turn of the helix and a set of points (actually three-dimensional delta-functions) aligned along the helix axis and separated axially by the pitch P. (B) shows that a discontinuous helix (i.e., a helical array of subunits) can be thought of as a product of the continuous helix in (A) and a set of horizontal density planes spaced h apart, where h is the subunit axial translation as in Fig. 7. This discontinuous set of points can then be convoluted with an atom (or a more complicated motif) to give a helical polymer. (C)-(F) represent helical objects and their computed diffraction patterns. (C) is half a turn of a helical wire. Its transform is a cross of intensity (high intensity is shown as white). (D) A full turn gives a similar cross with some substructure. A continuous helical wire has the transform of a complete helical turn, multiplied by the transform of the array of points in the middle of (A), namely, a set of planes of intensity a distance n/P apart (see Fig. 7). This means that in the transform in (E) the helix cross in (D) is only seen on the intensity planes, which are n/P apart. (F) shows the effect of making the helix in (E) discontinuous. The broken helix cross in (E) is now convoluted with the transform of the set of planes in (B), which are h apart. This transform is a set of points along the meridian of the diffraction pattern and separated by m/h. The resulting transform in (F) is therefore a series of helix crosses as in (E) but placed with their centers at the positions m/h from the pattern center. (Transforms calculated using MusLabel or FIELIX.)...

The Feynman diagram for the simplest annihilation event shows that annihilation is possible when the two particles are Ax h/mc 10 12 5 m apart, and that the duration of the event is At h/mc2 10-21 s. The distance is the geometric mean of nuclear and atomic dimensions, which is probably not significant. The distance is so much smaller than electronic wave functions that it may be assumed to be zero in computations of annihilation rates. The time is so short that, during it, a valence electron in a typical atom or molecule moves a distance of only ao/104, so that a spectator electron can be assumed to be stationary and the annihilating electron can be assumed to disappear in zero time. Thus the calculation of annihilation rates requires the evaluation of expectation values of the Dirac delta function, and the relaxation of the daughter system (post-annihilation remnant) can be understood with the aid of the sudden approximation [4], These are both relatively simple computations, providing an accurate wave function is available. 

The extent of gas dispersion can usually be computed from experimentally measured gas residence time distribution. The dual probe detection method followed by least square regression of data in the time domain is effective in eliminating error introduced from the usual pulse technique which could not produce an ideal Delta function input (Wu, 1988). By this method, tracer is injected at a point in the fast bed, and tracer concentration is monitored downstream of the injection point by two sampling probes spaced a given distance apart, which are connected to two individual thermal conductivity cells. The response signal produced by the first probe is taken as the input to the second probe. The difference between the concentration-versus-time curves is used to describe gas mixing. 

The expression for the anisotropic part of hyperfine coupling involves an integral over the spatial distribution of the unpaired electron, which is relatively easy to compute accurately even at a relatively low level of theory. The contact term, however, includes a delta-function that chips out the wave function amplitude at the nucleus point. The latter is quite difficult to compute both because standard Gaussian basis sets do not reproduce the wavefunction cusp at the nucleus point and because additional flexibility has to be introduced into the core part of the basis to account for the now essential core valence interaction. " ... 

The formulation that has been given here is not the only approach to the description of two-phase flows with nonequilibrium processes. Many different viewpoints have been pursued textbooks are available on the subject [43], [44], and a reasonably thorough review recently has been published [45]. Combustion seldom has been considered in this extensive literature. Most of the work that has addressed combustion problems has not allowed for a continuous droplet distribution function but instead has employed a finite number of different, discrete droplet sizes in seeking computer solution sets of conservation equations [5]. The present formulation admits discrete sizes as special cases (through the introduction of delta functions in fj) but also enables influences of continuous distributions to be investigated. A formulation of the present type recently has been extended to encompass thick sprays [25]. Some other formulations of problems of multiphase reacting flows have been mentioned in Sections 7.6 and 7.7. 

Before discussing the use of molecular dynamics to compute this quantity, one notes that by replacing the delta function by its integral representation ... 

The probability density function, written as pif), describes the fraction of time that the fluctuating variable/ takes on a value between/ and/ + A/. The concept is illustrated in Fig. 5.7. The fluctuating values off are shown on the right side while p(f) is shown on the left side. The shape of the PDF depends on the nature of the turbulent fluctuations of/. Several different mathematical functions have been proposed to express the PDF. In presumed PDF methods, these different mathematical functions, such as clipped normal distribution, spiked distribution, double delta function and beta distribution, are assumed to represent the fluctuations in reactive mixing. The latter two are among the more popular distributions and are shown in Fig. 5.8. The double delta function is most readily computed, while the beta function is considered to be a better representation of experimentally observed PDF. The shape of these functions depends solely on the mean mixture fraction and its variance. The beta function is given as... 

In this problem, elaborate on the Einstein model for structural change given in the chapter by replacing the delta function densities of states by rectangular densities of states. Assume that the phase with the lower internal energy has the broader vibrational density of states. Compute the transformation temperature as a function of the difference in the widths of the two rectangular bands. 

The reconstructed distribution function may be continuous (EQMOM) or discrete (QMOM), but we will assume that it is always realizable (i.e. nonnegative). For the case in which / is a set of weighted delta functions, the computation of the moments and is trivial. With EQMOM the integrals are evaluated using... 

The residence-time distribution of the PFR was shown to be arbitrarily sharp because all molecules spend identical times in the PFR. We introduced the delta function to describe this arbitrarily narrow RTD. We added a dispersion term to the PFR equations to model the spread of the RTD observed in actual tubular reactors. We computed the full, transient behavior of the dispersed plug-flow model, and displayed the evolution of the concentration profile after a step change to the feed concentration. 

For comparison, we show in the following and compute the case of a very narrow distribution and the actual values for the base case of a single holding time, that is for Dirac-Delta Function of residence times ... 

---