Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface 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.

46 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

SECARB Regional Project Assessment

The SECARB Cranfield Early Test Project was designed to safely demonstrate large-scale, long-term CO 2 injection and storage in a CO 2 -enhanced oil recovery (EOR) operation and associated saline reservoir. Storage during EOR operations has the potential to offset commercial costs, holding a significant promise for future large-scale deployment in the southeastern United States. Denbury Onshore LLC began CO 2 -EOR operations in July 2008 at the Cranfield oilfield located east of Natchez, Mississippi (Figure 1). The location offered an ideal setting to investigate storage monitoring at a commercial CO 2 -EOR operation. The field team, led by the Gulf Coast Carbon Center of the Bureau of Economic Geology at the University of Texas at Austin, deployed the largest suite of monitoring, verification, accounting (MVA) technologies to collect data to assess long-term CO 2 utilization and storage. The CO 2 injection interval is within the 3,000-meter-deep fluvial lower Tuscaloosa Formation, a large and regionally extensive saline reservoir, at a broad four-way structural closure that defines the Cranfield oilfield. Tests were conducted in the oil-producing area as well as in the non-oil-producing brine aquifer. The CO 2 used during EOR operation was transported via pipeline from the Jackson Dome, a natural source of CO 2 . The study focused on four integrated research areas: (1) the High-Volume Injection Test area (HiVIT); (2) the Detailed Area of Study (DAS); (3) the Geomechanical area; and (4) the near-surface observatory, also called the “P-site.” CO 2 injection occurred at the HiVIT and the DAS (Figure 2 and Figure 3). The DAS was located down-dip of the oil-water interface of the Cranfield oilfield. The well layout included one injection well and two down-dip observation wells. The team was the first of the regional partnerships to monitor more than 1 million metric tons of CO 2 . At the end of the project, the team monitored more than 11 million metric tons of CO 2 injected and 5 million metric tons of CO 2 stored. The Cranfield Project was the fifth project worldwide, and the first such project in the United States, to reach this CO 2 injection volume while being monitored to demonstrate storage effectiveness.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Neural Network-Enhanced Reproducing Kernel Particle Method for Image-Based Multiphysics Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant stresses during charge/discharge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), image-based modeling techniques are used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electrochemical-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in capturing crack propagation due to mesh dependency. Additionally, commonly used damage models, such as the continuous damage model and the cohesive zone model, often have steep tradeoffs between discontinuous field accuracy and computational expense. In this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022.

damage modeling↗
DOE OSTI2020

ZPPR-15 Criticality Measurements

A reference critical configuration was established for each phase of the ZPPR-15 experimental program. The reference critical configuration for each phase corresponds to a configuration in which the measured excess reactivity was adjusted to a standard set of reporting conditions with an average reactor temperature of 293 K, a zero interface gap and all operational control rods parked 30 inches from the axial centerline of the core. The ZPPR staff published the measured excess reactivity and the statistical uncertainty in the measurement. In reality, the published uncertainty for each reference critical configuration is only a small component of the total uncertainty for that configuration. Additional uncertainties related measurement technique, configuration geometry and material compositions exist. The geometry and material composition uncertainties are the dominant uncertainties for ZPPR-15. A full uncertainty analysis was performed for each ZPPR-15 reference critical configuration. The significant uncertainties related to measurement technique, configuration geometry and material compositions were quantified, and a total uncertainty was derived for each reference critical configuration. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models for the reference critical configurations established in ZPPR-15A loading 15, ZPPR-15B loading 88, ZPPR-15C loading 166 and ZPPR-15D loading 185. The as-built models neglect the small, non-uniform gap between matrix halves at full closure, any structures outside the matrix tubes and minor components, e.g., impurities, in some materials. The biases associated with the neglected items were quantified. The experimental value of keff and its associated uncertainty for each reference critical configuration were adjusted to account for biases introduced by items that are neglected in the as-built model for that configuration.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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 OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2025

Crystal structures of Salmonella enterica FraB deglycase reveal a conformational heterodimer with remarkable structural plasticity at the active site

Abstract Thefralocus ofSalmonella entericaencodes five genes for metabolism of fructose‐asparagine, an Amadori product formed by condensation of asparagine with glucose. In the last step of this pathway, the FraB deglycase cleaves 6‐phospho‐fructose‐aspartate into glucose‐6‐phosphate and aspartate. In homology models, FraB forms a homodimer with two equivalent active sites located at the dimer interface. E214 and H230, two invariant residues essential for catalysis, project into each active site cleft from opposing subunits of the dimer. Here, we have determined six crystal structures of FraB, three of a variant containing an N‐terminal His 6 tag and two mutations needed for crystallization (hereafter referred to as WT′), two with additional mutations to active site residues (E214A and P232A), and one of a variant with C‐terminal residues 313–325 deleted. Surprisingly, in the WT′ FraB structure, the two catalytic residues, E214 (general base) and H230 (general acid), are positioned ~22 Å apart. In the E214A and C‐terminus‐truncated FraB variants, however, a conformational change in the E214‐residing helix brings E214 and H230* to ~7 Å (* indicates residue from the second protomer that creates the inter‐subunit catalytic center). The loop bearing H230 also exhibits significant variation, ranging from being completely disordered to adopting open or closed states, with the nearby P232* residue being eithercisortrans. The C‐terminal residues 313–325 form a flexible “C‐tail” that can be fully disordered, bind in the active site to block access of substrate, or angle across the active site to wrap across the other subunit of the dimer and potentially close over substrate. Collectively, these structures reveal that FraB is a conformational heterodimer with two chemically identical subunits that are constrained to adopt different structures as they come together for catalysis. This plasticity likely involves correlated opening and closure of the two active sites for their respective binding and release of substrates and ligands.

Biochemistry & Molecular Biology↗