SEARCH · Engineering Papers

Results for “Point Kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗

DOE OSTI2020

Pt-Assisted Carbon Remediation of Mo2C Materials for CO Disproportionation

Using the CO disproportionation (Boudouard) reaction as a probe reaction, an in-depth analysis of temperature-programmed pulse response data shows that the addition of Pt to Mo2C mitigates deactivation of Mo active sites by acting as a carbon collector. CO2 production on Mo2C and Pt/Mo2C materials is dependent on both the activation energy and the CO surface concentration. Detailed plane-wave density functional theory calculations of the CO adsorption and disproportion reactions on Mo2C-supported Pt nanoparticles (NPs) are reported. The Mo2C was modeled by the ß-Mo2C (100) surface, and the Pt/Mo2C interface was modeled by the addition of 12 Pt atoms to the Mo2C (100) surface (12Pt@Mo2C). The potential energy surfaces of the Boudouard reaction were calculated on pure Mo2C, 12Pt@Mo2C, and Pt (111) surfaces. CO dissociation readily occurs on the Mo2C (100) surface, but not on the Pt (111) surface, with the former being exothermic and the latter being endothermic. At the Pt/Mo2C interface, CO dissociation is still exothermic, but with a larger energy barrier. The Boudouard reaction takes place on the Mo2C region, where CO2 is formed from a surface O atom dissociated from one CO molecule in reaction with another CO molecule, leaving one C atom on the surface. C adsorption is preferential on the Pt site in comparison to the Mo site. The supported Pt domains can collect the remaining C atoms, facilitating further CO2 formation on the active Mo sites. A Bader charge analysis shows that the surface metal-carbon bond is a mixture of covalent and ionic bonds, whereas the surface metal-oxygen bond is ionic. Electron localization function (ELF) and partial charge density calculations agree well with the Bader charge analysis. These computational results are consistent with experimental observations of the interaction of CO with Mo2C nanotube supported Pt domains in the transient regime under far from equilibrium conditions. The Boudouard reaction is an important side reaction, and the unexpected role found for Pt as a carbon collector, with Mo serving as a disproportionation site, provides a unique vantage point for understanding carbon and coke formation on catalytic materials.

Pt-assisted carbon remediation, Pt/Mo2C interface,↗

DOE OSTI2024

Rapid-charging aluminium-sulfur batteries operated at 85 °C with a quaternary molten salt electrolyte

Molten salt aluminum-sulfur batteries are based exclusively on resourcefully sustainable materials, and are promising for large-scale energy storage owed to their high-rate capability and moderate energy density; but the operating temperature is still high, prohibiting their applications. Here we report a rapid-charging aluminium-sulfur battery operated at a sub-water-boiling temperature of 85°C with a tamed quaternary molten salt electrolyte. The quaternary alkali chloroaluminate melt – possessing abundant electrochemically active high-order Al-Cl clusters and yet exhibiting a low melting point – facilitates fast Al 3+ desolvation. A nitrogen-functionalized porous carbon further mediates the sulfur reaction, enabling the battery with rapid-charging capability and excellent cycling stability with 85.4% capacity retention over 1400 cycles at a charging rate of 1 C. Importantly, we demonstrate that the asymmetric sulfur reaction mechanism that involves formation of polysulfide intermediates, as revealed by operando X-ray absorption spectroscopy, accounts for the high reaction kinetics at such temperature wherein the thermal management can be greatly simplified by using water as the heating media.

25 ENERGY STORAGE↗

DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1976

A comparison of measured and inferred temperatures from AEROS-B

The neutral composition and the temperature data obtained from the AEROS-B Neutral Atmosphere Temperature Experiment (NATE) and the Neutral and Ion Mass Spectrometer (NIMS) are compared, and the general validity of inferring gas temperatures from N2 and Ar density profiles is examined by comparing them with the in situ measured values of the neutral kinetic temperature (NATE). At times serious discrepancies are noted between the inferred and the measured temperature. This is particularly evident during periods of increased magnetic activity when the normally observed latitudinal variations are apparently modulated by waves propagating from the polar region to low latitudes. Under these conditions the N2 and Ar densities and temperature oscillations are usually out of phase, and the temperatures inferred from N2 and Ar at a given point become meaningless.

Chandra, S.↗

NASA NTRS2016

On the Connection Between Microbursts and Nonlinear Electronic Structures in Planetary Radiation Belts

Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave-particle interactions. We show that wave parameters, consistent with large amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (greater than 1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles Theta (sub k)B greater than 50 degrees and phase-speeds v(sub phi) greater than or equal to c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E greater than or equal to 100 keV) on kinetic timescales, which is much faster than previously inferred.

plasmas↗

NASA NTRS1991

Chemical-kinetic problems of future NASA missions

Thermochemical nonequilibrium in the shock layer surrounding vehicles entering the atmospheres of earth and Mars at superescape velocities is studied, deriving reaction rate coefficients that reproduce experimental data obtained in shock tubes. Thermodynamic properties and emitted radiation intensities are obtained for shock tube flow and flow in a shock layer over a blunt body. The results indicate that the viscous layer of the ablation product over an ablating heat shield is likely to be in chemical nonequilbrium. For earth entry flight, the thickness of the nonequilbrium region is between and 2 cm at the expected peak radiation point in the aerobraking trajectory, For Martian entry flight it is between 8 and 23 cm. For the earth entry case, nonequilibrium phenomena reduce radiative heating rate, while the opposite occurs for the Martian case. The radiative heat transfer rates are significant for the Mars entry conditions at entry velocities equal to or greater than 7 km/s.

Park, Chul↗

DOE OSTI2026

Effects of Polymer Morphology on Solvent and Catalyst Accessibility during Polyethylene and Polystyrene Autoxidation

Efficient catalytic deconstruction of plastics requires facile solvent and catalyst access to polymer substrates to minimize mass transfer effects. Autoxidation using Co(II) acetate, Mn(II) acetate, and a radical carrier in acetic acid is a promising strategy to deconstruct mixed plastic waste, yet the role of polymer morphology in governing solvent and catalyst accessibility remains poorly understood. Here, in situ simultaneous small- and wide-angle X-ray scattering (SAXS/WAXS), complemented by X-ray fluorescence (XRF) imaging and high-pressure differential scanning calorimetry (DSC), were used to elucidate interactions between acetic acid, a Co/Mn catalyst solution, and semicrystalline polyethylene (PE) and amorphous polystyrene (PS) from room temperature to 160 °C. In PE, acetic acid and catalyst access were confined to amorphous regions and cryomilled particle interfaces at room temperature, while crystalline lamellae remained intact after soaking for up to 34 h. Increasing temperature enabled solvent uptake into PE, followed by solvent-assisted softening above 100 °C, and a modest melting-point depression that removed lamellar transport barriers upon melting. Conversely for PS, acetic acid penetrated the glassy polymer without inducing chain mobility until the glass transition was reached, above which the observed structural changes were consistent with enhanced segmental mobility which enabled bulk penetration. These results suggest that polymer morphology and thermally activated physical transitions arising from diffusion and polymer–solvent interactions can influence whether autoxidation of plastics is transport-limited or kinetically controlled, providing a framework for aligning reaction conditions with reaction outcomes.

09 BIOMASS FUELS↗

DOE OSTI2024

Direct observation of ultrafast cluster dynamics in supercritical carbon dioxide using X-ray Photon Correlation Spectroscopy

Supercritical fluids exhibit distinct thermodynamic and transport properties, making them of particular interest for a wide range of scientific and engineering applications. These anomalous properties emerge from structural heterogeneities due to the formation of molecular clusters at conditions above the critical point. While the static behavior of these clusters and their effects on the thermodynamic response functions have been recognized, the relation between the ultrafast cluster dynamics and transport properties remains elusive. By measuring the intermediate scattering function in carbon dioxide at conditions near the critical point with X-ray photon correlation spectroscopy, we directly capture the cross-over dynamics between 4 and 13 picoseconds, revealing the transition between ballistic and diffusive motion. Complementary analysis using large-scale molecular dynamics simulations reveals that this behavior arises from collisions between unbound molecules and clusters. This study provides direct evidence of the ultrafast momentum exchange between clusters, which has significant impact on transport properties, solvation processes, and reaction kinetics in supercritical fluids.

carbon capture and storage↗

DOE OSTI2025

Cryosectioning-enhanced super-resolution microscopy for single-protein imaging across cells and tissues

DNA-points accumulation for imaging in nanoscale topography (DNA-PAINT) enables nanoscale imaging with virtually unlimited multiplexing and molecular counting. Here, we address challenges, such as variable imaging performance and target accessibility, that can limit its broader applicability. Specifically, we enhance its capacity for robust single-protein imaging and molecular counting by optimizing the integration of total internal reflection fluorescence microscopy with physical sectioning, in particular, Tokuyasu cryosectioning. Our method, tomographic and kinetically enhanced DNA-PAINT (tkPAINT), achieves 3 nm localization precision across diverse samples, enhanced imager binding, and improved cellular integrity. tkPAINT can facilitate molecular counting with DNA-PAINT inside the nucleus, as demonstrated through its quantification of the in situ abundance of RNA Polymerase II in both HeLa cells as well as mouse tissues. Anticipating that tkPAINT could become a versatile tool for the exploration of biomolecular organization and interactions across cells and tissues, we also demonstrate its capacity to support multiplexing, multimodal targeting of proteins and nucleic acids, and three-dimensional (3D) imaging.

Science & Technology - Other Topics↗

DOE OSTI2024

Single-molecule tracking reveals dual front door/back door inhibition of Cel7A cellulase by its product cellobiose

Degrading cellulose is a key step in the processing of lignocellulosic biomass into bioethanol. Cellobiose, the disaccharide product of cellulose degradation, has been shown to inhibit cellulase activity, but the mechanisms underlying product inhibition are not clear. We combined single-molecule imaging and biochemical investigations with the goal of revealing the mechanism by which cellobiose inhibits the activity of Trichoderma reesei Cel7A, a well-characterized exo-cellulase. We find that cellobiose slows the processive velocity of Cel7A and shortens the distance moved per encounter; effects that can be explained by cellobiose binding to the product release site of the enzyme. Cellobiose also strongly inhibits the binding of Cel7A to immobilized cellulose, with a K i of 2.1 mM. The isolated catalytic domain (CD) of Cel7A was also inhibited to a similar degree by cellobiose, and binding of an isolated carbohydrate-binding module to cellulose was not inhibited by cellobiose, suggesting that cellobiose acts on the CD alone. Finally, cellopentaose inhibited Cel7A binding at micromolar concentrations without affecting the enzyme’s velocity of movement along cellulose. Together, these results suggest that cellobiose inhibits Cel7A activity both by binding to the “back door” product release site to slow activity and to the “front door” substrate-binding tunnel to inhibit interaction with cellulose. These findings point to strategies for engineering cellulases to reduce product inhibition and enhance cellulose degradation, supporting the growth of a sustainable bioeconomy.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1988

Kinetic isolation tether experiment

Progress was made on the analysis of tether damping and on experimentation of the control system on the laboratory simulator. The damping analysis considers the dynamics of a long tether connecting two spacecraft in Earth orbit, one of the spacecraft having dominant mass. In particular, it considers the material damping of the tether. The results show that, with properly chosen tether material and braiding structure, longitudinal vibration of the tethered system is well damped. A particularly effective method of implementing attitude control for tethered satellites is to use the tether tension force to generate control torques by moving the tether attach point relative to the satellite center of mass. A scaled, one dimensional laboratory simulation of the KITE mission was built and preliminary experiments of the proposed attitude control system were performed. The simulator was built to verify theoretical predictions of attitude controllability, and to investigate the technological requirements in order to implement this concept. A detailed description of the laboratory apparatus is provided, and in addition, the results of the preliminary experiments are presented and discussed.

Powell, J. David↗

DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗

DOE OSTI2023

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1977

Surface charge kinetics near metal-dielectric interfaces exposed to kilovolt electron flux

Interfaces between dielectric films and grounded metallic boundaries were exposed, in vacuum, to monoenergic electron fluxes having energies up to 22 keV. Two principal concerns were the measuring of the charge distributions on dielectrics and the determining of causes of flashovers, events where dielectric surface charges abruptly transfer to the metallic structures. Surface charges are perturbed within 10 mm of interfaces. Perturbations are relatively small except within about 3 mm of the interface. The probability of flashover was found to be related to microscopic imperfections in the interfaces. As flashovers occur in an exposed metal substrate, points become burned into the dielectric along the slit. As these points develop, the probability of flashover increases greatly. An interface which is highly immune to flashover was formed by covering a dielectric film with a 1.5-mm-thick aperture plate which exposes the film through a machined opening.

Robinson, J. W.↗

NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗