Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A continuous calibration of the ATLAS flavour-tagging classifiers via optimal transportation maps

A calibration of the ATLAS flavour-tagging algorithms using a new calibration procedure based on optimal transportation maps is presented. Simultaneous, continuous corrections to the b-jet, c-jet, and light-flavour jet classification probabilities from jet-tagging algorithms in simulation are derived for b-jets using $t\bar{t} \rightarrow e\mu \nu \nu bb$ data. After application of the derived calibration maps, closure between simulation and observation is achieved for jet flavour observables used in ATLAS analyses of Large Hadron Collider (LHC) Run 2 proton-proton collision data. This continuous calibration opens up new possibilities for the future use of jet flavour information in LHC analyses and also serves as a guide for deriving high-dimensional corrections to simulation via transportation maps, an important development for a broad range of inference tasks.

Aad, G. [Aix-Marseille Université] (ORCID:00000002↗
DOE OSTI2023

Neural networks for large eddy simulations of wall-bounded turbulence: numerical experiments and challenges

Here, we examine the application of neural network-based methods to improve the accuracy of large eddy simulations of incompressible turbulent flows. The networks are trained to learn a mapping between flow features and the subgrid scales, and applied locally and instantaneously—in the same way as traditional physics-based subgrid closures. Models that use only the local resolved strain rate are poorly correlated with the actual subgrid forces obtained from filtering direct numerical simulation data. We see that highly accurate models in a priori testing are inaccurate in forward calculations, owing to the preponderance of numerical errors in implicitly filtered large eddy simulations. A network that accounts for the discretization errors is trained and found to be unstable in a posteriori testing. We identify a number of challenges that the approach faces, including a distribution shift that affects networks that fail to account for numerical errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Substituent effect on napthodithiophene-fused porphyrins: Understanding the unusual trend of fluorescence quantum yield quenching

A series of π-extended porphyrins containing thiophene units was prepared to investigate an unusual trend observed in the fluorescence quantum yield of naphtho[2,1-b:3,4-b']dithiophene-fused porphyrins. Monobenzoporphyrins carrying 2-thiophenyl (vinyl thiophene porphyrin) groups (2VTP, 2VTBr and Br2VTP) (numbering 2 and 3 refers to the position of sulfur on the thiophene ring) were synthesized through a Heck-based cascade reaction followed by a newly developed bromination method. These were converted to naphtho[2,1-b:3,4-b']dithiophene-fused porphyrin derivatives (FBr2VTP and F2VTBr) via ring closure with an intramolecular Scholl reaction. Increasingly larger Stokes shifts with a higher number of bromo groups were observed in the unfused π-extended molecular systems, reflecting the heavy atom effect. Fluorescence spectroscopy further confirmed the unusual trend seen in the previous work: structural rigidification in naphtho[2,1-b:3,4-b']dithiophene-fused porphyrins leads to longer fluorescence lifetimes but unexpectedly lowers the quantum yield. Adding one bromo group to the naphtho[2,1-b:3,4-b']dithiophene unit does not change this trend. Furthermore, the presence of two bromo groups prevents the quantum yield from dropping. DFT, TDDFT, and NICS analyses suggest a drastic change in aromaticity in the pyrrole ring where the naphtho[2,1-b:3,4-b']dithiophene is fused in FBr2VTP, which might contribute to the quantum yield reductions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

ANDES (ARGONNE NUCLEAR DATA EXPLORATION SOFTWARE) - ndviz

The Argonne Nuclear Data Exploration Software (ANDES) is a collection of open-source C++ codes and ROOT scripts that enable the visualization of patterns emerging from nuclear data. From this global view, well-known phenomena such as major nuclear shell closures and pairing effects near N=Z nuclei, as well as more subtle features, such as quantum phase transition around the Zr isotopes or the charge radius staggering in the Hg isotopes, can be observed. This visualization tool apart from having pedagogic application can serve as a tool for discovering previously unknown phenomena in the nuclear chart.

Santiago-Gonzalez, Daniel↗
DOE OSTI2020

ANDES (ARGONNE NUCLEAR DATA EXPLORATION SOFTWARE) - Physics Tools

A collection of physics tools used by the ANDES. The Argonne Nuclear Data Exploration Software (ANDES) is a collection of open-source C++ codes and ROOT scripts that enable the visualization of patterns emerging from nuclear data. From this global view, well-known phenomena such as major nuclear shell closures and pairing effects near N=Z nuclei, as well as more subtle features, such as quantum phase transition around the Zr isotopes or the charge radius staggering in the Hg isotopes, can be observed. This visualization tool apart from having pedagogic application can serve as a tool for discovering previously unknown phenomena in the nuclear chart.

Santiago-Gonzalez, Daniel↗
DOE OSTI2024

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗
DOE OSTI2024

FAR3d

The FAR3d model calculates the linear and nonlinear stability properties of energetic particle driven Alfven instabilities for both tokamak and stellarator plasma confinement devices using gyro-landau closure methods. This is an important fundamental physics problem for existing fusion energy experiments and for future fusion reactors.

Varela, Jacobo↗
DOE OSTI2026

Code for the manuscript "Lagrangian Attention Tensor Networks for Velocity Gradient Statistical Mode

We disclose a python/pytorch implementation of the physics-informed machine learning algorithm described in "Lagrangian Attention Tensor Networks for Velocity Gradient Statistical Modeling", LA-UR-24-30678. Direct numerical simulation (DNS) of ubiquitous turbulence phenomena is computationally infeasible for realistic flows. As a result, reduced modeling for turbulent flows aim to reduce the number of resolved scales while retaining accurate representations of the small-scale physics. The dynamics of the velocity gradient tensor (VGT) is a key ingredient in reduced or subgrid turbulence models. The evolution equation for the VGT involves nonlocal terms, requiring closure modeling. This implementation of the novel methodology of Lagrangian Attention Tensor Networks (LATN), utilizes a structured representation of the history of the VGT to inform a physics-informed machine learning algorithm. This addition of structured memory terms is shown to outperform previous models when trained and evaluated on DNS data.

Livescu, Daniel [LANL]↗
DOE OSTI2026

elm-diagnostics

elm-diagnostics is a Python package for computing diagnostic analyses and visualizations for the E3SM Land Model (ELM) component and is meant to support new feature development in ELM. The tool reads model history files and performs quantitative analyses including budget-closure checking, variable transformations, temporal aggregations, and statistical summaries to support model evaluation, validation, and scientific interpretation. The framework is designed for extensibility, with modular architecture enabling straightforward addition of new diagnostic methods, derived variables, analysis types, visualization approaches, and model-specific adaptations

Hoffman, Matt [Los Alamos National Laboratory]↗
DOE OSTI2024

Spatially Local Surrogate Modeling of Subgrid-Scale Effects in Idealized Atmospheric Flows: A Deep Learned Approach Using High-Resolution Simulation Data

Abstract We introduce a machine learned surrogate model from high-resolution simulation data to capture the subgrid-scale effects in dry, stratified atmospheric flows. We use deep neural networks (NNs) to model the spatially local state differences between a coarse-resolution simulation and a high-resolution simulation. The setup enables the capture of both dissipative and antidissipative effects in the state differences. The NN model is able to accurately capture the state differences in offline tests outside the training regime. In online tests intended for production use, the NN-coupled coarse simulation has higher accuracy over a significant period of time compared to the coarse-resolution simulation without any correction. We provide evidence of the capability of the NN model to accurately capture high-gradient regions in the flow field. With the accumulation of the errors, the NN-coupled simulation becomes computationally unstable after approximately 90 coarse simulation time steps. Insights gained from these surrogate models further pave the way for formulating stable, complex, physics-based spatially local NN models which are driven by traditional subgrid-scale turbulence closure models. Significance Statement Flows in the atmosphere are highly chaotic and turbulent, comprising flow structures of broad scales. For effective computational modeling of atmospheric flows, the effects of the small- and large-scale structures need to be captured by the simulations. Capturing the small-scale structures requires fine-resolution simulations. Even with the current state-of-the-art supercomputers, it can be prohibitively expensive to simulate these flows when computed for the entire earth over climate time scales. Thus, it is necessary to focus on the larger-scale structures using a coarse-resolution simulation while capturing the effects of the smaller-scale structures using some parameterization (approximation) scheme and incorporating it into the coarse-resolution simulation. We use machine learning to model the effects of the small-scale structures (subgrid-scale effects) in atmospheric flows. Data from a fine-resolution simulation is used to compute the missing subgrid-scale effects in coarse-resolution simulations. We then use machine learning models to approximate these differences between the coarse- and fine-resolution simulations. We see improved accuracy for the coarse-resolution simulations when corrected using these machine learned models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Conserving Land–Atmosphere Synthesis Suite (CLASS)

Accurate estimates of terrestrial water and energy cycle components are needed to better understand climate processes and improve models’ ability to simulate future change. Various observational estimates are available for the individual budget terms; however, these typically show inconsistencies when combined in a budget. In this work, a Conserving Land–Atmosphere Synthesis Suite (CLASS) of estimates of simultaneously balanced surface water and energy budget components is developed. Individual CLASS variable datasets, where possible, 1) combine a range of existing variable product estimates, and hence overcome the limitations of estimates from a single source; 2) are observationally constrained with in situ measurements; 3) have uncertainty estimates that are consistent with their agreement with in situ observations; and 4) are consistent with each other by being able to solve the water and energy budgets simultaneously. First, available datasets of a budget variable are merged by implementing a weighting method that accounts both for the ability of datasets to match in situ measurements and the error covariance between datasets. Then, the budget terms are adjusted by applying an objective variational data assimilation technique (DAT) that enforces the simultaneous closure of the surface water and energy budgets linked through the equivalence of evapotranspiration and latent heat. Comparing component estimates before and after applying the DAT against in situ measurements of energy fluxes and streamflow showed that modified estimates agree better with in situ observations across various metrics, but also revealed some inconsistencies between water budget terms in June over the higher latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Understanding the Spatiotemporal Variability of Tropical Orographic Rainfall Using Convective Plume Buoyancy

Abstract Mechanical forcing by orography affects precipitating convection across many tropical regions, but controls on the intensity and horizontal extent of the orographic precipitation peak and rain shadow remain poorly understood. A recent theory explains this control of precipitation as arising from modulation of lower-tropospheric temperature and moisture by orographic mechanical forcing, setting the distribution of convective rainfall by controlling parcel buoyancy. Using satellite and reanalysis data, we evaluate this theory by investigating spatiotemporal precipitation variations in six mountainous tropical regions spanning South and Southeast Asia, and the Maritime Continent. We show that a strong relationship holds in these regions between daily precipitation and a measure of convective plume buoyancy. This measure depends on boundary layer thermodynamic properties and lower-free-tropospheric moisture and temperature. Consistent with the theory, temporal variations in lower-free-tropospheric temperature are primarily modulated by orographic mechanical lifting through changes in cross-slope wind speed. However, winds directed along background horizontal moisture gradients also influence lower-tropospheric moisture variations in some regions. The buoyancy measure is also shown to explain many aspects of the spatial patterns of precipitation. Finally, we present a linear model with two horizontal dimensions that combines mountain wave dynamics with a linearized closure exploiting the relationship between precipitation and plume buoyancy. In some regions, this model skillfully captures the spatial structure and intensity of rainfall; it underestimates rainfall in regions where time-mean ascent in large-scale convergence zones shapes lower-tropospheric humidity. Overall, these results provide new understanding of fundamental processes controlling subseasonal and spatial variations in tropical orographic precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2021

Characterization of an electrically conductive proppant for fracture diagnostics

Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.

Geochemistry & Geophysics↗
DOE OSTI2024

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗
DOE OSTI2024

Soil Desiccation Treatability Testing at BC Waste Disposal Cribs

During Hanford’s production period, low-level waste products generated from chemical processing of uranium fuel rods were discharged directly to the ground through a system of cribs and trenches located in the 200-BC-1 Operable Unit (OU). The site consists of 6 cribs and 20 trenches that received more than 117,000 m3 of radioactive liquid waste discharged to the soil. These unlined infiltration galleries held volumes of liquid waste while it seeped into the ground, with the understanding that the 100 m (330 ft) thick vadose zone in the area would effectively capture the effluent and prevent groundwater impacts. A conceptual model showing the operation of cribs and trenches is shown in Figure 1. Data show effluent from the 26 cribs and trenches containing about 410 curies of Technetium-99 (Tc-99) is primarily located between 30 m and 70 m (98 ft and 230 ft) depth (Corbin et al., 2005; Ward et al., 2004). Despite no evidence indicating that the contamination has reached the groundwater at BC cribs and trenches, the mobility of Tc-99 had been demonstrated in laboratory tests and was recognized as a threat to groundwater at the site. Using data from numerical models, laboratory analyses, field investigations, and information on historical discharges, the EPA and Ecology identified Tc-99 and U contamination of the vadose zone as a remediation priority. The U.S. DOE was notified by EPA and Ecology regarding risks associated with Tc-99 contamination in a letter requesting development of a strategy for improved methods to understand the nature and extent of vadose zone contamination, specifically Tc-99, and to develop remedial options for addressing such contamination. To develop the appropriate technology for characterizing, remediating, and monitoring the deep vadose zone Tc-99 contamination, the U.S. DOE worked with the EPA and Ecology to create a Treatability Test Plan under a Remedial Investigation/Feasibility Study (RI/FS) for the Hanford 200 Areas. Under this RI/FS, it was determined that a treatability test for soil desiccation should be carried out as it was identified as a promising in-situ treatment technology for mitigating risks posed by Tc-99 contamination to the groundwater table. The BC Cribs and Trenches site was identified as a representative site for Tc-99 and U contamination and selected for the soil desiccation treatability test. In this chapter, we summarize the overlying regulatory framework of RI/FS and treatability tests and illustrate how development and experimentation supported the evaluation of selected remedies. We briefly discuss the RI/FS for the 200 Areas of the Hanford Site and focus on the soil desiccation treatability testing performed at the BC cribs and trenches site under the Deep Vadose Zone Treatability Test Plan for the Hanford Central Plateau (DVZ-TT). The DVZ-TT is one component of the remedial investigation/feasibility study for the Hanford 200 Areas and represents the underlying regulatory framework that drives site operations towards records of decision and site closure.

Mangel, Adam R.↗