Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “closure model”

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 973 records · Page 54

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

Testing and Implementation of Advanced Reynolds Stress Models

A research program was proposed for the testing and implementation of advanced turbulence models for non-equilibrium turbulent flows of aerodynamic importance that are of interest to NASA. Turbulence models that are being developed in connection with the Office of Naval Research ARI in Non-equilibrium are provided for implementation and testing in aerodynamic flows at NASA Langley Research Center. Close interactions were established with researchers at Nasa Langley RC and refinements to the models were made based on the results of these tests. The models that have been considered include two-equation models with an anisotropic eddy viscosity as well as full second-order closures. Three types of non-equilibrium corrections to the models have been considered in connection with the ARI on Nonequilibrium Turbulence: conducted for ONR.

Speziale, Charles G.↗

Dynamic Stochastic Superresolution of sparsely observed turbulent systems

Real-time capture of the relevant features of the unresolved turbulent dynamics of complex natural systems from sparse noisy observations and imperfect models is a notoriously difficult problem. The resulting lack of observational resolution and statistical accuracy in estimating the important turbulent processes, which intermittently send significant energy to the large-scale fluctuations, hinders efficient parameterization and real-time prediction using discretized PDE models. This issue is particularly subtle and important when dealing with turbulent geophysical systems with an vast range of interacting spatio-temporal scales and rough energy spectra near the mesh scale of numerical models. Here, we introduce and study a suite of general Dynamic Stochastic Superresolution (DSS) algorithms and show that, by appropriately filtering sparse regular observations with the help of cheap stochastic exactly solvable models, one can derive stochastically ‘superresolved’ velocity fields and gain insight into the important characteristics of the unresolved dynamics, including the detection of the so-called black swans. The DSS algorithms operate in Fourier domain and exploit the fact that the coarse observation network aliases high-wavenumber information into the resolved waveband. It is shown that these cheap algorithms are robust and have significant skill on a test bed of turbulent solutions from realistic nonlinear turbulent spatially extended systems in the presence of a significant model error. In particular, the DSS algorithms are capable of successfully capturing time-localized extreme events in the unresolved modes, and they provide good and robust skill for recovery of the unresolved processes in terms of pattern correlation. Moreover, we show that DSS improves the skill for recovering the primary modes associated with the sparse observation mesh which is equally important in applications. The skill of the various DSS algorithms depends on the energy spectrum of the turbulent signal and the observation time relative to the decorrelation time of the turbulence at a given spatial scale in a fashion elucidated here. The DSS technique exploiting a simple Gaussian closure of the nonlinear stochastic forecast model emerges as the most suitable trade-off between the superresolution skill and computational complexity associated with estimating the cross-correlations between the aliasing modes of the sparsely observed turbulent signal. Such techniques offer a promising and efficient approach to constraining unresolved turbulent fluxes through stochastic superparameterization and a subsequent improvement in coarse-grained filtering and prediction of the next generation atmosphere–ocean system (AOS) models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗

SootLib: A soot model library for combustion simulation

Soot formation in combustion is an important process that affects radiative heat transfer, flame temperatures, and emissions with health and environmental impacts. Soot formation involves complex chemistry for nucleation, growth, oxidation, and coagulation processes. The soot particles vary widely in size and accurate modeling requires representation of the particle size distribution (PSD). Modeling soot is not trivial, and is only one of several physical processes active in combustion systems. This paper presents a software package called SootLib, which is an open-source library for modeling soot formation and other aerosol systems. SootLib is written in C++, is documented with Doxygen, and is available on GitHub. The library includes several models for soot chemistry and coagulation, and it represents the PSD using either a sectional model or the method of moments (MOM). Four closure approaches for the MOM are implemented allowing up to eight moments: monodispersed, an assumed-shape lognormal distribution, the quadrature method of moments, and the method of moments with interpolative closure. SootLib provides an interface for inclusion in other combustion packages including CFD or reacting flow solvers. The range of models allows comparisons and sensitivity studies, and the modularity facilitates extension to other soot models.

97 MATHEMATICS AND COMPUTING↗

White Paper on Consequence Assessment of Transport of Hazardous Materials for PHMSA

This document describes a proposed methodology for assessing the consequences of accidents during transport of hazardous materials in the form of compressed gases as liquids in the United States. A quantitative population effect risk metric for a transport plan in terms of population numbers and effects is described. Quantitative results provide the Pipeline and Hazardous Materials Safety Administration (PHMSA) a means for evaluating routes, times, and other plans for material movement. Hazard releases are modeled with a consequence assessment tool maintained by the Defense Threat Reduction Agency (DTRA) that uses the Second-Order Closure Integrated Puff (SCIPUFF) atmospheric transport and dispersion model. LandScan TM 2016 USA Day and Night distributions are used for population data.

54 ENVIRONMENTAL SCIENCES↗

The vertical-velocity skewness in the inertial sublayer of turbulent wall flows

Empirical evidence is provided that within the inertial sublayer (i.e. logarithmic region) of adiabatic turbulent flows over smooth walls, the skewness of the vertical-velocity component 𝑆 𝑤 displays universal behaviour, being a positive constant and constrained within the range 𝑆 𝑤 ≈ 0.1–0.16, regardless of flow configuration and Reynolds number. A theoretical model is then proposed to explain this behaviour, including the observed range of variations of 𝑆 𝑤 . The proposed model clarifies why 𝑆 𝑤 cannot be predicted from down-gradient closure approximations routinely employed in large-scale meteorological and climate models. The proposed model also offers an alternative and implementable approach for such large-scale models.

turbulent boundary layers↗

Site-Wide Integrated Water Monitoring -- Defining and Implementing Sampling Objectives to Support Site Closure

The Underground Test Area (UGTA) activity is responsible for assessing and evaluating the effects of the underground nuclear weapons tests on groundwater at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), and implementing a corrective action closure strategy. The UGTA strategy is based on a combination of characterization, modeling studies, monitoring, and institutional controls (i.e., monitored natural attenuation). The closure strategy verifies through appropriate monitoring activities that contaminants of concern do not exceed the SDWA at the regulatory boundary and that adequate institutional controls are established and administered to ensure protection of the public. Other programs conducted at the NNSS supporting the environmental mission include the Routine Radiological Environmental Monitoring Program (RREMP), Waste Management, and the Infrastructure Program. Given the current programmatic and operational demands for various water-monitoring activities at the same locations, and the ever-increasing resource challenges, cooperative and collaborative approaches to conducting the work are necessary. For this reason, an integrated sampling plan is being developed by the UGTA activity to define sampling and analysis objectives, reduce duplication, eliminate unnecessary activities, and minimize costs. The sampling plan will ensure the right data sets are developed to support closure and efficient transition to long-term monitoring. The plan will include an integrated reporting mechanism for communicating results and integrating process improvements within the UGTA activity as well as between other U.S. Department of Energy (DOE) Programs.

54 ENVIRONMENTAL SCIENCES↗

Site-Wide Integrated Water Monitoring - Defining and Implementing Sampling Objectives to Support Site Closure - 13060

The Underground Test Area (UGTA) activity is responsible for assessing and evaluating the effects of the underground nuclear weapons tests on groundwater at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), and implementing a corrective action closure strategy. The UGTA strategy is based on a combination of characterization, modeling studies, monitoring, and institutional controls (i.e., monitored natural attenuation). The closure strategy verifies through appropriate monitoring activities that contaminants of concern do not exceed the SDWA at the regulatory boundary and that adequate institutional controls are established and administered to ensure protection of the public. Other programs conducted at the NNSS supporting the environmental mission include the Routine Radiological Environmental Monitoring Program (RREMP), Waste Management, and the Infrastructure Program. Given the current programmatic and operational demands for various water-monitoring activities at the same locations, and the ever-increasing resource challenges, cooperative and collaborative approaches to conducting the work are necessary. For this reason, an integrated sampling plan is being developed by the UGTA activity to define sampling and analysis objectives, reduce duplication, eliminate unnecessary activities, and minimize costs. The sampling plan will ensure the right data sets are developed to support closure and efficient transition to long-term monitoring. The plan will include an integrated reporting mechanism for communicating results and integrating process improvements within the UGTA activity as well as between other U.S. Department of Energy (DOE) Programs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Automated Symbolic Upscaling: 2. Model Generation for Extended Applicability Regimes

Abstract In this second part of the two paper series, we detail an algorithmic procedure for systematically implementing the generalized closure form strategy presented in Part 1. This strategy extends the applicability of homogenized models with respect to classical homogenization theory, as demonstrated in Part 1 where upscaled models are rigorously derived in moderately reactive physical regimes. After encoding the algorithm into Symbolica, an automated upscaling framework, we upscale two reactive mass transport problems and numerically validate the resulting nonlinear homogenized models by showing the absolute error estimates predicted by homogenization theory are satisfied. In both problems, nontrivial closure forms and closure problems are automatically formulated using the encoded strategy with no human interaction, nor prior knowledge regarding the closure required for the systems. We hope these demonstrations spark further interest in automated analytical frameworks for multiscale modeling, as such capabilities are invaluable for generating rigorous multiscale models of complex phenomena in porous media.

Pietrzyk, Kyle↗

Molecular model of TFIIH recruitment to the transcription-coupled repair machinery

Transcription-coupled repair (TCR) is a vital nucleotide excision repair sub-pathway that removes DNA lesions from actively transcribed DNA strands. Binding of CSB to lesion-stalled RNA Polymerase II (Pol II) initiates TCR by triggering the recruitment of downstream repair factors. Yet it remains unknown how transcription factor IIH (TFIIH) is recruited to the intact TCR complex. Combining existing structural data with AlphaFold predictions, we build an integrative model of the initial TFIIH-bound TCR complex. We show how TFIIH can be first recruited in an open repair-inhibited conformation, which requires subsequent CAK module removal and conformational closure to process damaged DNA. In our model, CSB, CSA, UVSSA, elongation factor 1 (ELOF1), and specific Pol II and UVSSA-bound ubiquitin moieties come together to provide interaction interfaces needed for TFIIH recruitment. STK19 acts as a linchpin of the assembly, orienting the incoming TFIIH and bridging Pol II to core TCR factors and DNA. Molecular simulations of the TCR-associated CRL4CSA ubiquitin ligase complex unveil the interplay of segmental DDB1 flexibility, continuous Cullin4A flexibility, and the key role of ELOF1 for Pol II ubiquitination that enables TCR. Collectively, these findings elucidate the coordinated assembly of repair proteins in early TCR.

Paul, Tanmoy↗

Final Technical Report

We have constructed and analyzed a series of idealized models of tropical convection interacting with large-scale circulations, with 25-50km resolution and with 1-2km cloud resolving resolution to set the stage for rigorous tests of convection closure schemes in high resolution global climate models. Much of the focus has been on the climatology of tropical cyclogenesis in rotating systems and the related problem of the spontaneous aggregation of convection in non-rotating systems. The PI (Held) will be delivering the honorary Bjerknes lecture at the Fall 2016 AGU meeting in December on this work. We have also provided new analyses of long-standing issues related to the interaction between convection and the large-scale circulation: Kelvin waves in the upper troposphere and lower stratosphere, water vapor transport into the stratosphere, and upper tropospheric temperature trends. The results of these analyses help to improve our understanding of processes, and provide tests for future high resolution global modeling. Our final goal of testing new convections schemes in next-generation global atmospheric models at GFDL has been left for future work due to the complexity of the idealized model results meant as tests for these models uncovered in this work and to computational resource limitations. 11 papers have been published with support from this grant, 2 are in review, and another major summary paper is in preparation.

54 ENVIRONMENTAL SCIENCES↗

Short Wavelength Electromagnetic Perturbations Excited Near the Solar Probe Plus Spacecraft in the Inner Heliosphere: 2.5D Hybrid Modeling

A 2.5D numerical plasma model of the interaction of the solar wind (SW) with the Solar Probe Plus spacecraft (SPPSC) is presented. These results should be interpreted as a basic plasma model derived from the SW-interaction with the spacecraft (SC), which could have consequences for both plasma wave and electron plasma measurements on board the SC in the inner heliosphere. Compression waves and electric field jumps with amplitudes of about 1.5 V/m and (12-18) V/m were also observed. A strong polarization electric field was also observed in the wing of the plasma wake. However, 2.5D hybrid modeling did not show excitation of whistler/Alfven waves in the upstream connected with the bidirectional current closure that was observed in short-time 3D modeling SPPSC and near a tether in the ionosphere. The observed strong electromagnetic perturbations may be a crucial point in the electromagnetic measurements planned for the future Solar Probe Plus (SPP) mission. The results of modeling electromagnetic field perturbations in the SW due to shot noise in absence of SPPSC are also discussed.

Lipatov, Alexander S.↗

Process‐Oriented Calibration of a Turbulence Scheme in the DOE's Global Storm‐Resolving Model Using Machine Learning

A process‐oriented calibration framework is developed for the Simplified Higher‐Order Closure (SHOC) turbulence scheme in DOE's Simple Cloud Resolving E3SM Atmospheric Model (SCREAM). This framework leverages machine learning surrogates and observational constraints to efficiently calibrate SHOC adjustable parameters across two convective regimes: clear‐sky dry convective boundary layer and fair‐weather shallow cumulus clouds from ARM observations. We use perturbed‐parameter ensembles of a doubly periodic version of SCREAM to train surrogates and apply Markov Chain Monte Carlo sampling guided by cost functions based on benchmarking large‐eddy simulations and observations to identify optimized parameter sets that perform well in both regimes. The calibrated SHOC parameters substantially improve boundary‐layer turbulence and cloud boundaries, and modeled cloud fraction and radiative effects align better with observations than the default. These results demonstrate that combining multiple process‐specific convective regimes with machine‐learning surrogates can reduce parametric uncertainties and yield a model more faithful to cloud–turbulence interactions.

58 GEOSCIENCES↗

E-Area Low-Level Waste Facility Vadose Zone Model: Confirmation of Water Mass Balance for Subsidence Scenarios

In preparation for the next revision of the E-Area Low-Level Waste Facility (LLWF) Performance Assessment (PA), a mass balance model was developed in Microsoft Excel to confirm correct implementation of intact- and subsided-area infiltration profiles for the proposed closure cap in the PORFLOW vadose-zone model. The infiltration profiles are based on the results of Hydrologic Evaluation of Landfill Performance (HELP) model simulations for both intact and subsided cases.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗