The Lagrangian spectral relaxation model for differential diffusion in homogeneous turbulence

The Lagrangian spectral relaxation ~LSR! model is extended to treat turbulent mixing of two passive scalars (fa and fb) with different molecular diffusivity coefficients ~i.e., differential-diffusion effects!. Because of the multiscale description employed in the LSR model, the scale dependence of differential-diffusion effects is described explicitly, including the generation of scalar decorrelation at small scales and its backscatter to large scales. The model is validated against DNS data for differential diffusion of Gaussian scalars in forced, isotropic turbulence at four values of the turbulence Reynolds number (Rl538, 90, 160, and 230! with and without uniform mean scalar gradients. The explicit Reynolds and Schmidt number dependencies of the model parameters allows for the determination of the Re ~integral-scale Reynolds number! and Sc ~Schmidt number! scaling of the scalar difference z5fa2fb . For example, its variance is shown to scale like ^z2& ;Re20.3. The rate of backscatter (bD) from the diffusive scales towards the large scales is found to be the key parameter in the model. In particular, it is shown that bD must be an increasing function of the Schmidt number for Sc\u3c1 in order to predict the correct scalar-to-mechanical time-scale ratios, and the correct long-time scalar decorrelation rate in the absence of uniform mean scalar gradients

A study of optimum cowl shapes and flow port locations for minimum drag with effective engine cooling, volume 2

The listings, user's instructions, sample inputs, and sample outputs of two computer programs which are especially useful in obtaining an approximate solution of the viscous flow over an arbitrary nonlifting three dimensional body are provided. The first program performs a potential flow solution by a well known panel method and readjusts this initial solution to account for the effects of the boundary layer displacement thickness, a nonuniform but unidirectional onset flow field, and the presence of air intakes and exhausts. The second program is effectually a geometry package which allows the user to change or refine the shape of a body to satisfy particular needs without a significant amount of human intervention. An effort to reduce the cruise drag of light aircraft through an analytical study of the contributions to the drag arising from the engine cowl shape and the foward fuselage area and also that resulting from the cooling air mass flowing through intake and exhaust sites on the nacelle is presented. The programs may be effectively used to determine the appropriate body modifications or flow port locations to reduce the cruise drag as well as to provide sufficient air flow for cooling the engine

Automated design of minimum drag light aircraft fuselages and nacelles

The constrained minimization algorithm of Vanderplaats is applied to the problem of designing minimum drag faired bodies such as fuselages and nacelles. Body drag is computed by a variation of the Hess-Smith code. This variation includes a boundary layer computation. The encased payload provides arbitrary geometric constraints, specified a priori by the designer, below which the fairing cannot shrink. The optimization may include engine cooling air flows entering and exhausting through specific port locations on the body

Improved Lagrangian mixing models for passive scalars in isotropic turbulence

Lagrangian data for velocity, scalars, and energy and scalar dissipation from direct numerical simulations are used to validate Lagrangian mixing models for inert passive scalars in stationary isotropic turbulence. The scalar fluctuations are nearly Gaussian, and, as a result of production by uniform mean gradients, statistically stationary. Comparisons are made for Taylor-scale Reynolds numbers in the range 38 to about 240 and Schmidt numbers in the range 1/8 to 1. Model predictions for one-point, one-time Eulerian statistics ~Eulerian correspondence! and one-particle, two-time Lagrangian statistics ~Lagrangian correspondence! are examined. Two scalar mixing models, namely the Lagrangian Fokker–Planck model and the Lagrangian colored-noise ~LCN! model, are proposed and written in terms of stochastic differential equations ~SDE! with specified drift and diffusion terms. Both of these models rely on statistics of the scalar field conditioned upon the energy dissipation, as provided by the Lagrangian spectral relaxation ~LSR! model. With the exception of the scalar dissipation, the models are shown to capture the Reynolds and Schmidt-number dependence of the Lagrangian integral time scales. However, the LCN model provides a more realistic description of the Lagrangian scalar fluctuations as differentiable time series having the correct form of the scalar autocorrelation function. Further extensions of the new mixing models to non-Gaussian scalars are conceptually straightforward, but require a closure for the scalar-conditioned scalar dissipation rate matrix. Likewise, accurate prediction of joint statistics for differential diffusion between different scalars with unequal molecular diffusivities will require the formulation of a multiscale SDE similar to the LSR model

Flight test evaluation of predicted light aircraft drag, performance, and stability

A technique was developed which permits simultaneous extraction of complete lift, drag, and thrust power curves from time histories of a single aircraft maneuver such as a pull up (from V max to V stall) and pushover (to V max for level flight). The technique, which is an extension of nonlinear equations of motion of the parameter identification methods of Iliff and Taylor and includes provisions for internal data compatibility improvement as well, was shown to be capable of correcting random errors in the most sensitive data channel and yielding highly accurate results. Flow charts, listings, sample inputs and outputs for the relevant routines are provided as appendices. This technique was applied to flight data taken on the ATLIT aircraft. Lack of adequate knowledge of the correct full throttle thrust horsepower true airspeed variation and considerable internal data inconsistency made it impossible to apply the trajectory matching features of the technique

A multienvironment conditional probability density function model for turbulent reacting flows

The multienvironment conditional probability density function (MECPDF) model was first proposed by Fox [Computational Models for Turbulent Reacting Flows (Cambridge University Press, Cambridge, 2003)] as a simple extension of multienvironment probability density function models for turbulent reacting flows. Like the conditional moment closure (CMC) and the laminar flamelet model (LFM), the MECPDF model describes the reacting scalars conditioned on the value of the mixture fraction. However, unlike CMC and LFM, the new model provides a consistent description of conditional fluctuations in both the scalar dissipation rate and the reacting scalars, and hence can be used to model partial extinction and reignition in homogeneous turbulent reacting flows. In this work, a general derivation of the MECPDF model is presented for a single reaction-progress variable using the direct quadrature method of moments. Extensions of the model to multiple reaction-progress variables and conditioning on the mixture-fraction vector are also discussed. After deriving the model, the closure assumptions are validated using direct simulations for pure diffusion of two randomly distributed, initially correlated scalar fields. Two homogeneous applications are then considered: nonreactive mixing starting from nontrivial initial conditions, and reactive mixing with partial extinction and reignition

Application of the Fokker-Planck molecular mixing model to turbulent scalar mixing using moment methods

An extended quadrature method of moments using the beta kernel density function (beta-EQMOM) is used to approximate solutions to the evolution equation for univariate and bivariate composition probability distribution functions (PDFs) of a passive scalar for binary and ternary mixing. The key element of interest is the molecular mixing term, which is described using the Fokker-Planck (FP) molecular mixing model. The direct numerical simulations (DNSs) of Eswaran and Pope [ Direct numerical simulations of the turbulent mixing of a passive scalar, Phys. Fluids 31, 506 (1988)] and the amplitude mapping closure (AMC) of Pope [ Mapping closures for turbulent mixing and reaction, Theor. Comput. Fluid Dyn. 2, 255 (1991)] are taken as reference solutions to establish the accuracy of the FP model in the case of binary mixing. The DNSs of Juneja and Pope [ A DNS study of turbulent mixing of two passive scalars, Phys. Fluids 8, 2161 (1996)] are used to validate the results obtained for ternary mixing. Simulations are performed with both the conditional scalar dissipation rate (CSDR) proposed by Fox [Computational Methods for Turbulent Reacting Flows (Cambridge University Press, 2003)] and the CSDR from AMC, with the scalar dissipation rate provided as input and obtained from the DNS. Using scalar moments up to fourth order, the ability of the FP model to capture the evolution of the shape of the PDF, important in turbulent mixing problems, is demonstrated. Compared to the widely used assumed beta-PDF model [S. S. Girimaji, Assumed beta-pdf model for turbulent mixing: Validation and extension to multiple scalar mixing, Combust. Sci. Technol. 78, 177 (1991)], the beta-EQMOM solution to the FP model more accurately describes the initial mixing process with a relatively small increase in computational cost