Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixed Compression”

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.

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

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Computer program for solving laminar, transitional, or turbulent compressible boundary-layer equations for two-dimensional and axisymmetric flow

A numerical algorithm and computer program are presented for solving the laminar, transitional, or turbulent two dimensional or axisymmetric compressible boundary-layer equations for perfect-gas flows. The governing equations are solved by an iterative three-point implicit finite-difference procedure. The software, program VGBLP, is a modification of the approach presented in NASA TR R-368 and NASA TM X-2458, respectively. The major modifications are: (1) replacement of the fourth-order Runge-Kutta integration technique with a finite-difference procedure for numerically solving the equations required to initiate the parabolic marching procedure; (2) introduction of the Blottner variable-grid scheme; (3) implementation of an iteration scheme allowing the coupled system of equations to be converged to a specified accuracy level; and (4) inclusion of an iteration scheme for variable-entropy calculations. These modifications to the approach presented in NASA TR R-368 and NASA TM X-2458 yield a software package with high computational efficiency and flexibility. Turbulence-closure options include either two-layer eddy-viscosity or mixing-length models. Eddy conductivity is modeled as a function of eddy viscosity through a static turbulent Prandtl number formulation. Several options are provided for specifying the static turbulent Prandtl number. The transitional boundary layer is treated through a streamwise intermittency function which modifies the turbulence-closure model. This model is based on the probability distribution of turbulent spots and ranges from zero to unity for laminar and turbulent flow, respectively. Several test cases are presented as guides for potential users of the software.

Harris, J. E.↗
DOE OSTI2017

Combustion energy frontier research center (CEFRC) final report (August 1, 2009 – July 31, 2016)

The Combustion Energy Frontier Research Center (CEFRC) was established to tackle the single overarching grand challenge of energy sustainability, energy security and global warming: to develop a “validated, predictive, multi-scale, combustion modeling capability to optimize the design and operation of evolving fuels in advanced engines for transportation applications,” as identified in the DOE report on “Basic Energy Needs for Clean and Efficient Combustion of 21st Century Transportation Fuels”. The challenge is particularly daunting since energy conversion efficiencies and exhaust emissions are governed by coupled chemical and transport processes at multiple length scales ranging from electron excitation to molecular rearrangements to nanoscale particulate formation to turbulent fuel/air mixing. To tackle this challenge, the CEFRC assembled a world-class team of 15 principal investigators, with the objectives to: 1) develop and test theoretical models to predict elementary reaction rates, molecule thermalization rates, chemical bond dissociation energies, and nonequilibrium transport properties using quantum chemistry calculations that account for strong electron correlation and multiple electronic potential energy surfaces; 2) develop automated kinetic mechanism generation, reduction, and error control methods for predicting alternative fuel including biofuel oxidation, fuel droplet decomposition, and NOx and particulate formation; 3) validate and improve the predictions of these models by measuring ignition delay times, species profiles, flame structures, burning limits, turbulence-kinetic coupling, and NOx and soot emissions at high-pressures and near-limit conditions, by using advanced experimental diagnostic techniques including multiple laser techniques, molecular beam sampling and synchrotron photoionization, and by conducting the measurements in high-pressure shock tubes, jet-stirred and flow reactors, flame bombs, counterflow flames, and advanced-design rapid compression ignition instruments; and 4) develop a suite of validated petascale high-fidelity simulation and modeling capabilities to understand and predict chemistry-turbulence-radiation coupling for new fuels in new regimes, including the high pressure, low-temperature combustion in advanced engine and turbine designs, and 5) establish a knowledge highway between researchers and engineers in academia, national laboratories, and industry to facilitate the dissemination and exchange of knowledge on national and international levels, and enrich the talent pool and capabilities of the next generation of combustion scientists and engineers. The technical activities of the CEFRC were conducted through three Disciplinary Working Groups – Chemistry Theory, Experiment and Mechanism, and Reacting Flows, which coordinated the Center’s research on the development of combustion chemistry of Foundation Fuels (C0–C4 hydrocarbons), Alcohols, and Biodiesel through three corresponding Mechanism Thrust Groups. Such a two-dimensional coordinated and tightly interwoven research structure has been proven to be highly effective in assuring the interplay between the developments of the fundamentals of combustion science and the utilization of the various categories of fuels. The Center has accomplished the above goals over the five year period (August 1, 2009 – July 31, 2014) with appropriated funding, followed by two additional no-cost-extension (NCE) years. The research results are documented in 230 journal publications, with six legacy review papers on the study of combustion chemistry using shock tubes, flow reactors, rapid compression machines, and flames, on uncertainty quantification, and on theoretical reaction dynamics and chemical modeling of combustion. A robust outreach program complemented these PI-led research activities, consisting of: 1) a roving post-doc program comprised of a corps of Center-appointed, co- or multi-sponsored post-doctoral fellows with rotating assignments to conduct seed projects initiated by at least two PIs, in residence with these sponsoring PIs, to rapidly pursue new and high-risk, high-payoff interdisciplinary ideas; 2) an annual summer school on combustion heavily attended (~200) by senior graduate students and practicing researchers covering advanced topics on chemical kinetics, fluid mechanics, turbulent combustion, engine combustion, new technologies, etc.; 3) a robust open web-site providing Center and community information as well as the lecture videos and notes of the summer school; and 4) widely distributed biannual newsletters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRSDate not supplied

Mitigation of Silver Ion Loss from Solution by Polymer Coating of Metal Surfaces, Part IV

Ionic silver (Ag+) biocide is a leading candidate to provide residual microbial control in spacecraft potable water systems, but suffers from rapid concentration loss due to interactions with the metallic storage containers and tubing. One approach to resolve this problem is the coating of metal surfaces with an inert barrier film. In previous reports, we have described our investigations addressing Ag+ loss mitigation and adhesion performance of parylene barrier coatings on coupons of several metal alloys and 316Ltubing under static immersion. In such conditions, parylene-C and -AF4 coatings have shown excellent Ag+loss mitigation and mixed long-term adhesion performance, depending on parylene species and substrate surface chemistry/structure. In Part IV of this series, we report on our work to investigate the performance of parylene-C barrier coatings, under more challenging and realistic conditions, in order to evaluate potential suitability for use. The resilience and associated Ag+ loss with Parylene-C coating on½” (12.7 mm)316L tubing and fittings under medium-term immersion and repeated fitting dis/reassembly were investigated. Potential mechanical challenges to barrier coatings in spacecraft potable water systems were investigated, with two main focuses being identified: liquid flow/pressure cycling in tubing and cyclic operation of bellows in positive-expulsion storage tanks. Two corresponding testbeds for experimental characterization of coating performance were developed. In the Flow/Pressure Testbed System, Ag+ solution will be pumped through internally coated 316L tubing, with optional pressure cycling. In the Bellows Tank Testbed System, an internally coated edge-welded 316L bellows containing Ag+ solution will be cyclically extended and compressed, analogous to the operation of bellows tanks used for potable water storage and delivery on the International Space Station. The design and operation of these testbeds are described. Finally, the Ag+ adsorption and saturation behaviors of Kalrez and Viton, for potential use as seal materials in fittings and other connections, were characterized in limited experiments.

mitigation↗
NASA NTRS1990

Three-dimensional compressible convection at low Prandtl numbers

Numerical simulations are used to study fully compressible thermal convection at large Rayleigh numbers. Results are presented from a sequence of three-dimensional simulations that reveal a transition from gradually-evolving laminar convection to nearly turbulent convection as the Prandtl number is reduced from a value of unity to one-tenth. The convective flows form irregular cellular patterns near the upper surface, possesing a network of fast downflow at cell peripheries and gentler upflow at cell centers. At greater depths the curving sheets of downflow collapse into plumes which may twist and possess substantial vertical vorticity. For the lowest Prandtl number, the convection near the bottom of the layer appears to be turbulent, yet the rapidly varying small-scale flow structure there is accompanied by more ordered sites of wavering upflow, with the latter able to penetrate all the way to the upper boundary. Thus a significant component of the flow is able to extend over multiple density scale heights, in contrast to what is argued in formulating mixing-length models for stellar convection. Results are also shown from two-dimensional simulations carried out with very high spatial resolution, which reveal that supersonic convection with fluttering shock systems can be realized.

Toomre, Juri↗
NASA NTRS1994

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence technique in the complex three-dimensional compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4% for pressure, temperature, and velocity and 3% for injectant mole fraction. Comparisons with numerical simulations using the SPARK three-dimensional Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector. Calculations and measurements are in good agreement throughout the flowfield, with deviations on the order of 5%; however, in specific regions, such as in the base of the ramp, deviations are larger. A weak asymmetry in the incoming flowfield appears to be amplified by boundary-layer separation occurring when the ramp-generated shock reflects off the tunnel walls. Ramp-generated vortices are weaker in the calculated results due to the effects of numerical viscosity in the vortex cores. This leads to less turning and mixing of the jet plume than observed in the experiments. The rate of decay of the maximum injectant mole fraction with streamwise distance is greater for the present ramp injection scheme than for previously measured transverse injection schemes. In this recirculation region at the base of the injector, laminar calculations show better agreement with the measurements than turbulent calculations.

Donohue, James M.↗
DOE OSTI2017

Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method

Here we investigate an efficient preconditioning of iterative methods (such as GMRES) for solving dense linear systems Ax = b that follow from a boundary element method (BEM) for the 3D Helmholtz equation, focusing on the low-frequency regime. While matrix–vector products in GMRES can be accelerated through the low-frequency fast multipole method (LFFMM), the BEM often remains computationally expensive due to the large number of GMRES iterations. We propose the application of the inverse fast multipole method (IFMM) as a preconditioner to accelerate the convergence of GMRES. The IFMM is in essence an approximate direct solver that uses a multilevel hierarchical decomposition and low-rank approximations. The proposed IFMM-based preconditioning has a tunable parameter ε that balances the cost to construct a preconditioner M, which is an approximation of A -1 , and the cost to perform the iterative process by means of M. Namely, using a small (respectively, large) value of ε takes a long (respectively, short) time to construct M, while the number of iterations can be small (respectively, large). A comprehensive set of numerical examples involving various boundary value problems with complicated geometries and mixed boundary conditions is presented to validate the efficiency of the proposed method. We show that the IFMM preconditioner (with a nearly optimal ε of 10 -2 ) clearly outperforms some common preconditioners for the BEM, achieving 1.2–10.8 times speed-up of the computations, in particular when the scale of the underlying scatterer is about five wavelengths or more. In addition, the IFMM preconditioner is capable of solving complicated problems (in a reasonable amount of time) that BD preconditioner can not.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗
NASA NTRS1999

Towards Integrated Pulse Detonation Propulsion and MHD Power

The interest in pulse detonation engines (PDE) arises primarily from the advantages that accrue from the significant combustion pressure rise that is developed in the detonation process. Conventional rocket engines, for example, must obtain all of their compression from the turbopumps, while the PDE provides additional compression in the combustor. Thus PDE's are expected to achieve higher I(sub sp) than conventional rocket engines and to require smaller turbopumps. The increase in I(sub sp) and the decrease in turbopump capacity must be traded off against each other. Additional advantages include the ability to vary thrust level by adjusting the firing rate rather than throttling the flow through injector elements. The common conclusion derived from these aggregated performance attributes is that PDEs should result in engines which are smaller, lower in cost, and lighter in weight than conventional engines. Unfortunately, the analysis of PDEs is highly complex due to their unsteady operation and non-ideal processes. Although the feasibility of the basic PDE concept has been proven in several experimental and theoretical efforts, the implied performance improvements have yet to be convincingly demonstrated. Also, there are certain developmental issues affecting the practical application of pulse detonation propulsion systems which are yet to be fully resolved. Practical detonation combustion engines, for example, require a repetitive cycle of charge induction, mixing, initiation/propagation of the detonation wave, and expulsion/scavenging of the combustion product gases. Clearly, the performance and power density of such a device depends upon the maximum rate at which this cycle can be successfully implemented. In addition, the electrical energy required for direct detonation initiation can be significant, and a means for direct electrical power production is needed to achieve self-sustained engine operation. This work addresses the technological issues associated with PDEs for integrated aerospace propulsion and MHD power. An effort is made to estimate the energy requirements for direct detonation initiation of potential fuel/oxidizer mixtures and to determine the electrical power requirements. This requirement is evaluated in terms of the possibility for MHD power generation using the combustion detonation wave. Small scale laboratory experiments were conducted using stoichiometric mixtures of acetylene and oxygen with an atomized spray of cesium hydroxide dissolved in alcohol as an ionization seed in the active MHD region. Time resolved thrust and MHD power generation measurements were performed. These results show that PDEs yield higher I(sub sp) levels than a comparable rocket engine and that MHD power generation is viable candidate for achieving self-excited engine operation.

Litchford, Ron J.↗
NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗
DOE OSTI2020

The effect of loading direction on the fracture behaviors of cortical bone at a dynamic loading rate

Here, the dynamic cracking processes in porcine cortical bone were visualized in real-time using the high-speed synchrotron X-ray phase-contrast imaging (PCI) technique in three osteon orientations: in-plane transverse, out-of-plane transverse and in-plane longitudinal. The dynamic flexural loading applied on the pre-notched bone specimens was introduced by a modified Kolsky compression bar. High-speed X-ray images of the entire loading events were documented with a high-speed camera. Three-dimensional X-ray micro-computed tomography was conducted to examine the intact microstructures and obtain the basic material properties of the bone material used for mechanical characterizations. The onset location, where crack initiated, and the subsequent direction, along which the incipient crack propagated, were measured quantitatively using the high-speed X-ray images and the latter was found dependent on the osteon direction significantly. The crack propagation velocities were dependent on crack extension over the entire crack path significantly for all the three directions while the initial velocity for in-plane longitudinal direction was lower than the other two directions. Straight-through crack paths were observed for in-plane longitudinal specimens while the cracks were deflected and twisted in the in-plane transverse direction. For out-of-plane transverse direction, the cracks follow paths with tortuosity fall in between the other two directions, showing a mixed mode of fractures of the former two extreme cases. The toughening mechanisms, visualized by the high-speed X-ray images, and the corresponding fracture toughness, evaluated in terms of fracture initiation toughness and crack growth resistance curve (R-curve), were also found significantly different among the three osteon directions, suggesting an overall transition from brittle to ductile-like fracture behaviors at the dynamic displacement rate (5.4 m/s) as the osteon orientation varies from in-plane longitudinal to out-of-plane transverse, and to in-plane transverse eventually.

42 ENGINEERING↗
DOE OSTI2025

Ultrahigh-pressure crystallographic passage towards metallic hydrogen

The structural evolution of molecular hydrogen H 2 under multi-megabar compression and its relation to atomic metallic hydrogen is a key unsolved problem in condensed-matter physics. Although dozens of crystal structures have been proposed by theory, only one, the simple hexagonal-close-packed (hcp) structure of only spherical disordered H 2 , has been previously confirmed in experiments. Through advancing nano-focused synchrotron X-ray probes, here we report the observation of the transition from hcp H 2 to a post-hcp structure with a six-fold larger supercell at pressures above 212 GPa, indicating the change of spherical H 2 to various ordered configurations. Theoretical calculations based on our XRD results found a time-averaged structure model in the space group $P\bar{6}2c$ with alternating layers of spherically disordered H 2 and new graphene-like layers consisting of H 2 trimers (H 6 ) formed by the association of three H 2 molecules. Here, this supercell has not been reported by any previous theoretical study for the post-hcp phase, but is close to a number of theoretical models with mixed-layer structures. The evidence of a structural transition beyond hcp establishes the trend of H 2 molecular association towards polymerization at extreme pressures, giving clues about the nature of the molecular-to-atomic transition of metallic hydrogen. Considering the spectroscopic behaviours that show strong vibrational and bending peaks of H 2 up to 400 GPa, it would be prudent to speculate the continuation of hydrogen molecular polymerization up to its metallization.

Ji, Cheng [Center for High Pressure Science and Te↗
NASA NTRS2004

Urethane/Silicone Adhesives for Bonding Flexing Metal Parts

Adhesives that are blends of commercially available urethane and silicone adhesives have been found to be useful for bonding metal parts that flex somewhat during use. These urethane/silicone adhesives are formulated for the specific metal parts to be bonded. The bonds formed by these adhesives have peel and shear strengths greater than those of bonds formed by double-sided tapes and by other adhesives, including epoxies and neat silicones. In addition, unlike the bonds formed by epoxies, the bonds formed by these adhesives retain flexibility. In the initial application for which the urethane/silicone adhesives were devised, there was a need to bond spring rings, which provide longitudinal rigidity for inflatable satellite booms, with the blades that provide the booms axial strength. The problem was to make the bonds withstand the stresses, associated with differences in curvature between the bonded parts, that arose when the booms were deflated and the springs were compressed. In experiments using single adhesives (that is, not the urethane/ silicone blends), the bonds were broken and, in each experiment, it was found that the adhesive bonded well with either the ring or with the blade, but not both. After numerous experiments, the adhesive that bonded best with the rings and the adhesive that bonded best with the blades were identified. These adhesives were then blended and, as expected, the blend bonded well with both the rings and the blades. The two adhesives are Kalex (or equivalent) high-shear-strength urethane and Dow Corning 732 (or equivalent) silicone. The nominal mixture ratio is 5 volume parts of the urethane per 1 volume part of the silicone. Increasing the proportion of silicone makes the bond weaker but more flexible, and decreasing the proportion of silicone makes the bond stronger but more brittle. The urethane/silicone blend must be prepared and used quickly because of the limited working time of the urethane: The precursor of the urethane adhesive is supplied in a two-part form, comprising a resin and a hardener that must be mixed. The resulting urethane adhesive has a working time of 3 to 5 minutes. To prepare the urethane/silicone blend, one must quickly add the silicone to the urethane adhesive and mix it in thoroughly within the working time of the urethane. Once the urethane/silicone blend has been mixed and applied to the bond surfaces, it takes about 2 hours for the adhesive to cure under pressure. However, it takes about 24 hours for the adhesive to reach full strength.

Edwards, Paul D.↗
NASA NTRS2011

Coordinated Analyses of Diverse Components in Whole Stardust Cometary Tracks

Analyses of samples returned from Comet 81P/Wild-2 by the Stardust spacecraft have resulted in a number of surprising findings that show the origins of comets are more complex than previously suspected. However, these samples pose new experimental challenges because they are diverse and suffered fragmentation, thermal alteration, and fine scale mixing with aerogel. Questions remain about the nature of Wild-2 materials, such as the abundances of organic matter, crystalline materials, and presolar grains. To overcome these challenges, we have developed new sample preparation and analytical techniques tailored for entire aerogel tracks. We have successfully ultramicrotomed entire "carrot" and "bulbous" type tracks along their axis while preserving their original shapes. This innovation allowed us to examine the distribution of fragments along the track from the entrance hole all the way to the terminal particle (TP). We will present results of our coordinated analysis of the "carrot" type aerogel tracks #112 and #148, and the "bulbous" type aerogel tracks #113, #147 and #168 from the nanometer to the millimeter scale. Scanning TEM (STEM) was used for elemental and detailed mineralogy characterization, NanoSIMS was used for isotopic analyses, and ultrafast two-step laser mass spectrometry (ultra L2MS) was used to investigate the nature and distribution of organic phases. The isotopic measurements were performed following detailed TEM characterization for coordinated mineralogy. This approach also enabled spatially resolving the target sample from fine-scale mixtures of compressed aerogel and melt. Eight of the TPs of track #113 are dominated by coarse-grained enstatite (En90) that is largely orthoenstatite with minor, isolated clinoenstatite lamellae. One TP contains minor forsterite (Fo88) and small inclusions of diopside with % levels of Al, Cr and Fe. Two of the TPs contain angular regions of fine-grained nepheline surrounded by enstatite. Their O isotopic compositions are in the range of meteoritic materials, implying that they originated in the inner Solar System. Complex aromatic hydrocarbons are distributed along aerogel tracks and in TPs. These organics are likely cometary but were affected by shock heating. Three TPs of track #147 and two of track 168 have completely different mineralogy. TP2 of track #147 entirely consists of Fe-Ni alloy (5 at% Ni) and TP3 contains Fa28 with partial olivine-pyroxene intergrowth and minor albite. TP4 contains pentlandite, Fe-olivine, albite and high Ca pyroxene with Na and Cr (kosmochlor component). TP1 of #168 contains Fe-olivine, albite and pentlandite, and the concentric TP2 has a core of olivine grains with co-existing indigenous amorphous SiO2 surrounded by a carbon mantle, which in turn is surrounded by a layer of compressed aerogel. The TP of the carrot track #112 is a (16)O-rich forsteritic olivine grain that likely formed in the inner Solar System. The track also contains submicron-sized diamond grains of likely Solar System origin.

Nakamura-Messenger, Keiko↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
NASA NTRS1998

Application of a Planar Doppler Velocimetry System to a High Reynolds Number Compressible Jet

A Planar Doppler Velocimetry (PDV) system has been constructed and used to investigate the instantaneous turbulent velocity structure of a round high-speed compressible air jet with a low-speed co-flow. The exit condition was Mach=0.85 at ambient pressure, yielding a Reynolds number of about 650,000 on diameter. The PDV system was installed at NASA Langley Research Center in the Small Anechoic Jet Facility (SAJF), a chamber in which both the acoustic and aerodynamic properties of jets can be studied. For this test, the goal was to gather data which can be used to relate the turbulence structure of the jet to the levels and character of the acoustic noise produced by the jet. The current PDV system can acquire single-velocity-component, single-shot, planar images (15ns exposures) at 30 Hz. For this paper, the primary data set consists of 240 frames of velocity data acquired with both the jet and the low-speed co-flow seeded with light-scattering articles. Thus, velocities could be measured everywhere in the jet shear layer, both in the jet fluid and in the entrained co-flow. Some data were also taken with only the jet flow seeded. These provided mixing concentration images along with the reduced velocity fields. Other images were taken with only the co-flow seeded. These produced unique quantitative images of high speed entrainment. Optical "laser speckle" noise is the largest source of random noise in pulsed PDV systems. Components for the PDV imaging system were specifically selected to minimize speckle noise. To reduce systematic velocity errors due to laser drift, a frequency monitoring reference leg with a temperature-tuned reference iodine cell, was employed. In the course of this study, a novel flow seeder was developed. It enabled continuously variable seeding of the flow with particles of Sheared Pyrogenic Amorphous Hydrophobic Silica (SPAHS). The seeder comprised a dry fluidized bed hopper and a supersonic nozzle "pickup." Shearing action in the pickup dispersed the seed material in an exceptionally fine cloud (approximately 0.3 micron). These particles followed the flow well, did not clump or cake on screens or model surfaces, and were not susceptible to evaporation. Because of the refractory nature of the particles, SPAHS seeding should also be applicable to anticipated future testing at high temperatures.

Smith, Michael W.↗
DOE OSTI2018

Prediction of Cyclic Variability and Knock-Limited Spark Advance (KLSA) in Spark-Ignition (SI) Engine

Engine knock remains one of the major barriers to further improve thermal efficiency of Spark Ignition (SI) engines. Knock can be suppressed by lowering the compression ratio, or retarding the spark ignition timing, however, at an expense of efficiency penalty. SI engine is usually operated at knock-limited spark advance (KLSA) to achieve possibly maximum efficiency with given engine hardware and fuel properties, such as Research Octane Number (RON), Motor Octane Number (MON), and heat of vaporization, etc. Co-optimization of engine design and fuel properties is promising to improve the engine efficiency and predictive CFD models can be used to facilitate this optimization process. However, difficulties exist in predicting KLSA in CFD simulations. First, cyclic variability of SI engine demands that multi-cycle results are required to capture the extreme conditions. Secondly, Mach Courant-Friedrichs-Lewy (CFL) number of 1 is desired to accurately predict the knock intensity (KI), resulting in unaffordable computational cost, especially for multi-cycle simulations. In this study, a new approach to numerically predict KLSA using large Mach CFL number of 50 is proposed. This approach is validated against experimental data for a boosted Direct Injection Spark Ignition (DISI) engine at multiple loads and spark timings. G-equation combustion model coupled with well-mixed chemical kinetic model are used to predict the turbulent flame propagation and end-gas auto-ignition, respectively. Simulations run for 10 consecutive engine cycles at each condition. The results show good agreement between model predictions and experiments in terms of cylinder pressure, combustion phasing and cyclic variation. Engine knock is predicted with early spark ignition timing, indicated by significant pressure wave oscillation and end-gas heat release. Maximum Amplitude of Pressure Oscillation (MAPO) analysis is performed to quantify the KI, and the slope change point in KI extrema is used to indicate the KLSA accurately. Using a smaller Mach CFL number of 5 also results in the same conclusions thus demonstrating that this approach is insensitive to the Mach CFL number. The use of large Mach CFL number allows us to achieve fast turn-around time for multi-cycle engine CFD simulations.

Yue, Zongyu↗
NASA NTRS1994

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗
DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗