Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure correlations”

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 271 records · Page 15

The coupling between flame surface dynamics and species mass conservation in premixed turbulent combustion

Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Trouve, A.↗

Sensible-Enthalpy-Based Conditional Moment Closure Model for Homogeneous Charge Compression Ignition with Temperature Inhomogeneity

Abstract This paper presents a sensible-enthalpy-based conditional moment closure (CMC) model for homogeneous charge compression ignition (HCCI) with temperature inhomogeneity. The focus of study is placed on the effects of a conditioning variable on the prediction of heat release rates and developing the mixing models for sensible enthalpy, a reacting scalar. The sensible-enthalpy-based CMC model is validated a priori with 2-D direct numerical simulation (DNS) data for HCCI, and compared with the total-enthalpy-based formulation. Sensible enthalpy is shown to be strongly correlated with temperature, which makes sensible enthalpy well suited for a conditioning variable in CMC for HCCI. Results show that sensible-enthalpy-based CMC performs well for the prediction of heat release rates for all temperature stratification levels investigated in this study, while total-enthalpy-based CMC leads to the overprediction of heat release rates when the temperature stratification level is high. While conditional fluctuations of species mass fractions are substantial, sensible-enthalpy-based CMC is found to accurately estimate heat release rates during ignition, as well as reaction rates of key species, due to substantially suppressed conditional temperature fluctuations. The mapping closure models for the probability density function (PDF) and conditional dissipation rates of sensible enthalpy are proposed. The proposed models are found to capture the characteristics of reactive scalar mixing.

42 ENGINEERING↗

Effect of chemistry and turbulence on NO formation in oxygen-natural gas flames

The effects of chemistry and turbulence on NO formation in oxygen-natural turbulent diffusion flames gas flames have been investigated. The chemistry of nitric oxides has been studied numerically in the counterflow configuration. Systematic calculations with the GRI 2.11 mechanism for combustion of methane and NO chemistry were conducted to provide a base case. It was shown that the 'simple' Zeldovich mechanism accounts for more than 75% of N2 consumption in the flame in a range of strain-rates varying between 10 and 1000 s-l. The main shortcomings of this mechanism are: 1) overestimation (15%) of the NO production rate at low strain-rates because it does not capture the reburn due to the hydrocarbon chemistry, and 2) underestimation (25%) of the NO production rate at high strainrates because it ignores NO production through the prompt mechanism. Reburn through the Zeldovich mechanism alone proves to be significant at low strain-rates. A one-step model based on the Zeldovich mechanism and including reburn has been developed. It shows good agreement with the GRI mechanism at low strain-rates but underestimates significantly N2 consumption (about 50%) at high strain-rates. The role of turbulence has been assessed by using an existing 3-D DNS data base of a diffusion flame in decaying turbulence. Two PDF closure models used in practical industrial codes for turbulent NO formation have been tested. A simpler version of the global one-step chemical scheme for NO compared to that developed in this study was used to test the closure assumptions of the PDF models, because the data base could not provide all the necessary ingredients. Despite this simplification, it was possible to demonstrate that the current PDF models for NO overestimate significantly the NO production rate due to the fact that they neglect the correlations between the fluctuations in oxygen concentration and temperature. A single scalar PDF model for temperature that accounts for such correlations based on laminar flame considerations has been developed and showed excellent agreement with the values given by the DNS.

Samaniego, J. -M.↗

Spatio‐temporally Resolved Methane Fluxes From the Los Angeles Megacity

Abstract We combine sustained observations from a network of atmospheric monitoring stations with inverse modeling to uniquely obtain spatiotemporal (3‐km, 4‐day) estimates of methane emissions from the Los Angeles megacity and the broader South Coast Air Basin for 2015–2016. Our inversions use customized and validated high‐fidelity meteorological output from Weather Research Forecasting and Stochastic Time‐Inverted Lagrangian model for South Coast Air Basin and innovatively employ a model resolution matrix‐based metric to disentangle the spatiotemporal information content of observations as manifested through estimated fluxes. We partially track and constrain fluxes from the Aliso Canyon natural gas leak and detect closure of the Puente Hills landfill, with no prior information. Our annually aggregated fluxes and their uncertainty excluding the Aliso Canyon leak period lie within the uncertainty bounds of the fluxes reported by the previous studies. Spatially, major sources of CH 4 emissions in the basin were correlated with CH 4 ‐emitting infrastructure. Temporally, our findings show large seasonal variations in CH 4 fluxes with significantly higher fluxes in winter in comparison to summer months, which is consistent with natural gas demand and anticorrelated with air temperature. Overall, this is the first study that utilizes inversions to detect both enhancement (Aliso Canyon leak) and reduction (Puente Hills) in CH 4 fluxes due to the unintended events and policy decisions and thereby demonstrates the utility of inverse modeling for identifying variations in fluxes at fine spatiotemporal resolution.

Yadav, Vineet↗

Development of a second order closure model for computation of turbulent diffusion flames

A typical eddy box model for the second-order closure of turbulent, multispecies, reacting flows developed. The model structure was quite general and was valid for an arbitrary number of species. For the case of a reaction involving three species, the nine model parameters were determined from equations for nine independent first- and second-order correlations. The model enabled calculation of any higher-order correlation involving mass fractions, temperatures, and reaction rates in terms of first- and second-order correlations. Model predictions for the reaction rate were in very good agreement with exact solutions of the reaction rate equations for a number of assumed flow distributions.

Varma, A. K.↗

Measurement and analysis of critical crack tip processes associated with variable amplitude fatigue crack growth

Crack growth retardation following overloads can result in overly conservative life predictions in structures subjected to variable amplitude fatigue loading when linear damage accumulation procedures are employed. Crack closure is believed to control the crack growth retardation, although the specific closure mechanism is debatable. Information on the relative contributions to crack closure from: (1) plasticity left in the wake of the advancing crack and (2) crack tip residual stresses is provided. The delay period and corresponding crack growth rate transients following overloads are systematically measured as a function of load ratio (R) and overload magnitude. These responses are correlated in terms of the local 'driving force' for crack growth as measured by crack tip opening loads and delta K sub eff. The latter measurements are obtained using a scanning electron microscope equipped with a cyclic loading stage; measurements are quantified using a relatively new stereoimaging technique. Combining experimental results with analytical predictions suggests that both plastic wake and residual stress mechanism are operative, the latter becoming predominate as R increases.

Hudak, S. J., Jr.↗

An Enriched Shell Element for Delamination Simulation in Composite Laminates

A formulation is presented for an enriched shell finite element capable of delamination simulation in composite laminates. The element uses an adaptive splitting approach for damage characterization that allows for straightforward low-fidelity model creation and a numerically efficient solution. The Floating Node Method is used in conjunction with the Virtual Crack Closure Technique to predict delamination growth and represent it discretely at an arbitrary ply interface. The enriched element is verified for Mode I delamination simulation using numerical benchmark data. After determining important mesh configuration guidelines for the vicinity of the delamination front in the model, a good correlation was found between the enriched shell element model results and the benchmark data set.

McElroy, Mark↗

Statistical model for premixed turbulent flames

Statistical concepts and physical arguments are used to develop closure models for the mean energy equation governing one-dimensional, premixed, intense turbulent flames. It is proposed that the turbulent transport terms can be modeled by a turbulent conductivity based on turbulence kinetic energy (k) and energy dissipation rate (epsilon). The models obtained are used to derive correlations for the turbulent flame speed. Results show that there are three different flame regimes: a low-intensity turbulence, wrinkled flame; a high-intensity, fast-chemistry flame; and a high-intensity, finite-rate-chemistry flame.

Gouldin, F. C.↗

Production and Study of Neutron-rich Nuclei Using the LICORNE Directional Neutron Source

We have recently successfully demonstrated a new technique for production and study of many of the most exotic neutron-rich nuclei at moderate spins. LICORNE, a newly developed directional inverse-kinematic fast neutron source at the IPN Orsay, was coupled to the MINIBALL high resolution γ-ray spectrometer to study nuclei the furthest from stability using the 238 U(n,f)reaction. This reaction and 232 Th(n,f)are the most neutron-rich fission production mechanisms achievable and can be used to simultaneously populate hundreds of neutron-rich nuclei up to spins of ≈ 16h. High selectivity in the experiment was achieved via tripleγ-ray coincidences and the use of a 400 ns period pulsed neutron beam,a technique which is unavailable to other population mechanisms such as 235 U(nth,f) and 252 Cf(SF). The pulsing allows time correlations to be exploited to separate delayedγrays from isomeric states in the hundreds of nuclei produced, which are then used to cleanly select a particular nucleus and its exotic binary partners. In the recent experiment, several physics cases are simultaneously addressed such as shape coexistence, the evolution of shell closures far from stability, and the spectroscopy of nuclei in the r-process path near N= 82. Overall, preliminary physics results on anomalies in the 238 U(n,f) fission yields and the structure of the 138 Te and100Sr nuclei will soon be published. A future project,ν-ball, to couple LICORNE with a hybrid escape-suppressed spectrometer to refine further the technique and achieve a large increase in the observational limit is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stimulation of isoprene emissions and electron transport rates are a key mechanism of thermal tolerance in the tropical species Vismia guianensis

Tropical forests absorb large amounts of atmospheric CO 2 through photosynthesis, but high surface temperatures suppress this absorption while promoting isoprene emissions. While mechanistic isoprene emission models predict a tight coupling to photosynthetic electron transport (ETR) as a function of temperature, direct field observations of these phenomenon are lacking in the tropics and are necessary to assess the impact of a warming climate on global isoprene emissions. In this paper, we demonstrate that in the early successional species Vismia guianensis in the central Amazon, ETR rates increased with temperature in concert with isoprene emissions, even as stomatal conductance ($g _s$) and net photosynthetic carbon fixation ($P_n$) declined. We observed the highest temperatures of continually increasing isoprene emissions yet reported (50°C). While $P_n$ showed an optimum value of 32.6 ± 0.4°C, isoprene emissions, ETR, and the oxidation state of PSII reaction centers (q L ) increased with leaf temperature with strong linear correlations for ETR (ρ = 0.98) and q L (ρ= 0.99) with leaf isoprene emissions. In contrast, other photoprotective mechanisms, such as non-photochemical quenching (NPQ), were not activated at elevated temperatures. Inhibition of isoprenoid biosynthesis repressed $P_n$ at high temperatures through a mechanism that was independent of stomatal closure. While extreme warming will decrease $g_s$ and $P_n$ in tropical species, our observations support a thermal tolerance mechanism where the maintenance of high photosynthetic capacity under extreme warming is assisted by the simultaneous stimulation of ETR and metabolic pathways that consume the direct products of ETR including photorespiration and the biosynthesis of thermoprotective isoprenoids. Our results confirm that models which link isoprene emissions to the rate of ETR hold true in tropical species and provide necessary "ground-truthing" for simulations of the large predicted increases in tropical isoprene emissions with climate warming.

54 ENVIRONMENTAL SCIENCES↗

Second order closure for variable density free shear layer

A full second order closure turbulence model for the prediction of free shear flows with variable density is developed. The importance of density fluctuations due to inhomogeneities in fluid composition is emphasized. Some simplicity is achieved by using density weighted averaging to account implicitly for density fluctuation effects. Transport equations for various second order correlations - such as Reynolds stress, turbulent mass fluxes or density velocity correlations - are derived, modelled and solved, together with the mean equations for momentum and species conservation. Solutions of these equations are compared with the two different sets of experimental data for low speed free mixing layers (helium into nitrogen and freon 12 into air).

Vandromme, D.↗

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

Machine learning method for enforcing variable independence in background estimation with LHC data: ABCDisCoTEC

A novel solution is presented for the problem of estimating the backgrounds of a signal search using observed data while simultaneously maximizing the sensitivity of the search to the signal. The 'ABCD method' provides a reliable framework for background estimation by partitioning events into one signal-enhanced region (A) and three background-enhanced control regions (B, C, and D) via two smoothly varying, statistically independent variables. In practice, even slight correlations between the two variables can significantly undermine the method's performance. Thus, choosing appropriate variables by hand can present a formidable challenge, especially when background and signal differ only subtly. To address this issue, the ABCD with distance correlation (ABCDisCo) method was developed to construct two learned variables via a neural network trained to provide strong signal-background discrimination with small values of the distance correlation (DisCo) measure between the two learned variables. However, relying solely on minimizing the DisCo can result in learned variables that may not have distributions of background events that are smoothly varying and localized at extreme values, as necessary for the validity of the background estimation. The ABCDisCo training enhanced with closure (ABCDisCoTEC) method is introduced to solve this issue by directly minimizing the nonclosure, expressed as a dedicated differentiable loss term. This extended method is applied to a data set of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS detector at the CERN Large Hadron Collider. Additionally, given the complexity of the minimization problem with constraints on multiple loss terms, the modified differential method of multipliers is applied and shown to greatly improve the stability and robustness of the ABCDisCoTEC method, compared to grid search hyperparameter optimization procedures.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

Temporal entanglement transitions in the periodically driven Ising chain

Periodically driven quantum systems can host non-equilibrium phenomena without static analogs, including in their entanglement dynamics. Here, we discover (TET) in a Floquet spin chain, which correspond to a quantum phase transition in the spectrum of the entanglement Hamiltonian and are signaled by dynamical spontaneous symmetry breaking. We identify the symmetry principles underlying these transitions: they appear when the driven Hamiltonian preserves global symmetry (here, ℤ2), the initial state respects this symmetry, and the reduced density matrix carries weight in both subsystem-parity sectors, with TET occurring precisely when the sector weights become equal (given the previous two conditions are also satisfied). Intriguingly, we find these transitions across a broad range of driving frequencies (from adiabatic to high-frequency regime) and independently of drive details, where they manifest as periodic, sharp entanglement spectrum reorganizations marked by the Schmidt-gap closure, a vanishing entanglement echo, and symmetry-quantum-number flips, while remaining invisible to conventional local observables. At high frequencies, the entanglement Hamiltonian acquires an intrinsic timescale decoupled from the drive period, rendering the transitions genuine steady-state features. Finite-size scaling reveals universal critical behavior with correlation-length exponent 𝜈 = 1, matching equilibrium Ising universality despite its emergence from purely dynamical mechanisms decoupled from static criticality. Furthermore, our work establishes TET as novel features in Floquet quantum matter.

36 MATERIALS SCIENCE↗

Simulating the Clamped Tapered Beam Specimen Under Quasi-Static and Fatigue Loading Using Floating Node Method

As part of the NASA Advanced Composites Project (ACP), a sub-element has been designed to provide validation data for progressive damage analysis models. The clamped tapered beam is a cross-ply laminated composite specimen designed to validate the simulation of the onset of matrix cracks and their interaction with delaminations, including delamination migration. A tapered geometry was used to localize the first damage occurrence in the tapered region, without prescribing an initial crack. The boundary and loading conditions were chosen to favor delamination growth and subsequent migration after the first damage occurrence. The typical sequence of events consists of a matrix crack located at the tapered region, leading to delamination onset, followed by delamination growth and subsequent delamination migration to a different interface via a dominant matrix crack. The Clamped Tapered Beam (CTB) was tested in both quasi-static and fatigue regimes. The results obtained are used in this study to assess and validate a methodology based on the Floating Node Method (FNM) implemented as an Extended Interface Element. In this methodology, quasi-static and fatigue damage formation and development are modeled by combining FNM to represent crack networks, with Directional Cohesive Zone Elements (DCZE) and Virtual Crack Closure Technique (VCCT), respectively. Qualitatively, the methodology is capable of predicting the sequence of events and overall failure morphology. Quantitatively, the simulation results generally bound the experimental data, based on the range of the characterization data used. In this paper, the results from quasi static and fatigue simulations are compared and correlated with experimental data.

Seshadri, B. R.↗

On the correlation between ‘non-local’ effects and intrinsic rotation reversals in Alcator C-Mod

Contemporary predictive models for heat and particle transport in tokamak plasmas are based on the assumption that local fluxes can be described in terms of local plasma parameters, where electromagnetic drift-wave-type turbulence is driven by local gradients and results in cross-field transport. The question of whether or not transport could be dominated by non-local terms in certain circumstances is essential for our understanding of transport in magnetically confined plasmas, and critical for developing predictive models for future tokamaks, such as ITER. Perturbative transport experiments using cold-pulse injections at low density seem to challenge the local closure of anomalous transport: a rapid temperature increase in the core of the plasma following a sharp edge cooling is widely observed in tokamaks and helical devices. Past work in Ohmic plasmas in Alcator C-Mod and in ECH plasmas in KSTAR found that the temperature inversions disappear at higher densities, above the intrinsic toroidal rotation reversal density. These observations suggested that the so-called 'non-local' heat transport effects were related to the intrinsic rotation reversal, and therefore to changes in momentum transport. In this work, new experiments and analysis at Alcator C-Mod show that intrinsic rotation reversals and disappearance of temperature inversions are not concomitant in Ohmic plasmas at high plasma current and in ICRH L-modes. This new data set shows that the correlation between transient temperature inversions and intrinsic rotation reversals is not universal, suggesting that 'non-local' heat transport and momentum transport effects may be affected by different physical mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the correlation between ‘non-local’ effects and intrinsic rotation reversals in Alcator C-Mod

Contemporary predictive models for heat and particle transport in tokamak plasmas are based on the assumption that local fluxes can be described in terms of local plasma parameters, where electromagnetic drift-wave-type turbulence is driven by local gradients and results in cross-field transport. The question of whether or not transport could be dominated by non-local terms in certain circumstances is essential for our understanding of transport in magnetically confined plasmas, and critical for developing predictive models for future tokamaks, such as ITER. Perturbative transport experiments using cold-pulse injections at low density seem to challenge the local closure of anomalous transport: a rapid temperature increase in the core of the plasma following a sharp edge cooling is widely observed in tokamaks and helical devices. Past work in Ohmic plasmas in Alcator C-Mod and in ECH plasmas in KSTAR found that the temperature inversions disappear at higher densities, above the intrinsic toroidal rotation reversal density. These observations suggested that the so-called 'non-local' heat transport effects were related to the intrinsic rotation reversal, and therefore to changes in momentum transport. In this work, new experiments and analysis at Alcator C-Mod show that intrinsic rotation reversals and disappearance of temperature inversions are not concomitant in Ohmic plasmas at high plasma current and in ICRH L-modes. This new data set shows that the correlation between transient temperature inversions and intrinsic rotation reversals is not universal, suggesting that 'non-local' heat transport and momentum transport effects may be affected by different physical mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aeolian Shear Stress Ratio Measurements within Mesquite-Dominated Landscapes of the Chihuahuan Desert, New Mexico, USA

A field study was conducted to ascertain the amount of protection that mesquite-dominated communities provide to the surface from wind erosion. The dynamics of the locally accelerated evolution of a mesquite/coppice dune landscape and the undetermined spatial dependence of potential erosion by wind from a shear stress partition model were investigated. Sediment transport and dust emission processes are governed by the amount of protection that can be provided by roughness elements. Although shear stress partition models exist that can describe this, their accuracy has only been tested against a limited dataset because instrumentation has previously been unable to provide the necessary measurements. This study combines the use of meteorological towers and surface shear stress measurements with Irwin sensors to measure the partition of shear stress in situ. The surface shear stress within preferentially aligned vegetation (within coppice dune development) exhibited highly skewed distributions, while a more homogenous surface stress was recorded at a site with less developed coppice dunes. Above the vegetation, the logarithmic velocity profile deduced roughness length (based on 10-min averages) exhibited a distinct correlation with compass direction for the site with vegetation preferentially aligned, while the site with more homogenously distributed vegetation showed very little variation in the roughness length. This distribution in roughness length within an area, defines a distribution of a resolved shear stress partitioning model based on these measurements, ultimately providing potential closure to a previously uncorrelated model parameter.

King, James↗