Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Insight on electrolyte infiltration of lithium ion battery electrodes by means of a new three-dimensional-resolved lattice Boltzmann model

Electrolyte filling takes place between sealing and formation in Lithium Ion Battery (LIB) manufacturing process. This step is crucial as it is directly linked to LIB quality and affects the subsequent time consuming electrolyte wetting process. Although having fast, homogeneous and complete wetting is of paramount importance, this process has not been sufficiently examined and fully understood. For instance, experimentally available data is insufficient to fully capture the complex interplay upon filling between electrolyte and air inside the porous electrode. We report here for the first time a 3D-resolved Lattice Boltzmann Method (LBM) model able to simulate electrolyte filling upon applied pressure of LIB porous electrodes obtained both from experiments (micro X-ray tomography) and computations (stochastic generation, simulation of the manufacturing process using Coarse Grained Molecular Dynamics and Discrete Element Method). The model allows obtaining advanced insights about the impact of the electrode mesostructures on the speed of electrolyte impregnation and wetting, highlighting the importance of porosity, pore size distribution and pores interconnectivity on the filling dynamics. Furthermore, we identify scenarios where volumes with trapped air (dead zones) appear and evaluate the impact of those on the electrochemical behavior of the electrodes.

25 ENERGY STORAGE↗
DOE OSTI2020

Crossover from Rouse to Reptation Dynamics in Salt-Free Polyelectrolyte Complex Coacervates

Considerable interest in the dynamics and rheology of polyelectrolyte complex coacervates has been motivated by their industrial application as viscosity modifiers. A central question is the extent to which classical Rouse and reptation models can be applied to systems where electrostatic interactions play a critical role on the thermodynamics. By relying on molecular simulations, we present a direct analysis of the crossover from Rouse to reptation dynamics in salt-free complex coacervates as a function of chain length. This crossover shifts to shorter chain lengths as electrostatic interactions become stronger, which corresponds to the formation of denser coacervates. To distinguish the roles of Coulomb interactions and density, here we compare the dynamics of coacervates to those of neutral, semidilute solutions at the same density. Both systems exhibit a universal dynamical behavior in the connectivity-dominated (subdiffusion and normal diffusion) regimes, but the monomer relaxation time in coacervates is much longer and increases with increasing Bjerrum length. This is similar to the cage effect observed in glass-forming polymers, but the local dynamical slowdown is caused here by strong Coulomb attractions (ion pairing) between oppositely charged monomers. Our findings provide a microscopic framework for the quantitative understanding of coacervate dynamics and rheology.

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

The Present and Future of QCD: QCD Town Meeting White Paper – An Input to the 2023 NSAC Long Range Plan

It is currently understood that there are four fundamental forces in nature: gravitational, electromagnetic, weak and strong forces. The strong force governs the interactions between quarks and gluons, elementary particles whose interactions give rise to the vast majority of visible mass in the universe. The mathematical description of the strong force is provided by the non-Abelian gauge theory Quantum Chromodynamics (QCD). While QCD is an exquisite theory, constructing the nucleons and nuclei from quarks, and furthermore explaining the behavior of quarks and gluons at all energies, remain to be complex and challenging problems. Such challenges, along with the desire to understand all visible matter at the most fundamental level, position the study of QCD as a central thrust of research in nuclear science. Experimental insight into the strong force can be gained using large particle accelerator facilities, which are necessary to probe the very short distance scales over which quarks and gluons interact. The Long Range Plans (LRPs) exercise of 1989 and 1996 led directly to the construction of two world-class facilities: the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) that is focused on studying how the structure of hadrons emerges from QCD (cold QCD research), and the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL) that aims at the discovery and study of a new state of matter, the quark-gluon plasma (QGP), at extremely high temperatures (hot QCD research). These past investments have produced major advances. Nucleons and nuclei are being studied with increasing precision with a unified description of the partonic structure utilizing multi-dimensional imaging. Significant progress has been made, paving the way towards a complete picture of how quarks and gluons give rise to the mass, spin, and momentum of the nucleon. In hot QCD, the QGP is created in the collisions of nuclei at RHIC and the Large Hadron Collider (LHC) and is observed to behave like a fluid with very low specific shear viscosity; the current goals are to understand how the fluid behavior emerges from QCD and to characterize the temperature (and chemical potential) dependence of the properties of the QGP. As this White Paper is written, current experimental programs at CEBAF, RHIC and the LHC continue to provide exciting near term opportunities to capitalize on the investments in experimental equipment and accelerator operations. Most importantly, the QCD community looks forward to the construction of the Electron Ion Collider (EIC) as a major new facility to push forward QCD research in the next decades, with significant focus on exploring the properties of gluons, the mediators of the strong force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Developing Methods to Assess Changes in Mechanical Properties of Shale Modified by Engineered Mineral Precipitation

Fractures in subsurface shale formations serve multiple purposes, for example, in the recovery of resources in hydraulic fracturing or as potential harmful leakage passages through caprocks that may contribute undesired fluids to the atmosphere or functional groundwater aquifers. A proposed method to seal or influence fracture properties is Ureolysis-Induced Calcium Carbonate Precipitation (UICP), a bio-mineralization technology driven by the enzymatic hydrolysis of urea, resulting in the formation of calcium carbonate. The resulting calcium carbonate can bridge the gaps in fractured shale and reduce fluid flow through fractures. This study represents the first step toward determining the influence of UICP treatment on shale material and its subsequent mechanical strength properties. The goal of this preliminary work is twofold: first, we aim to identify a method to test tensile strength along a core axis and second, we seek to assess the effect of temperature on the tensile strength of intact, unfractured shale cores (2.54 cm (1 in) diameter, 5.08 cm (2 in) long for comparison with future fractured and UICP-treated cores. A modified Brazilian indirect tensile strength test successfully measured splitting tensile strength of shale cores from Eagle Ford and Wolfcamp formations at room temperature and 60°C.

clastic rock↗
DOE OSTI2023

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Hotel Room Computational Fluid Dynamics to Investigate Airborne Pathogen Dispersal Patterns

A hotel room unit consisting of a bedroom and bathroom was modelled using computational fluid dynamics (CFD) to investigate airborne pathogen dispersal patterns. The full-scale model includes a ‘typical’ hotel room configuration, furniture, and vents. The air sources and sinks include a bathroom vent, a heating, ventilation, and cooling (HVAC) unit located in the bedroom, and a ½” gap at the bottom of the entry door. In addition, the entry door and window can be opened or closed, as desired. Three key configuration simulations were conducted: 1) both the bathroom vent and HVAC were on, 2) only the HVAC was on, and 3) only the bathroom vent was on. If the HVAC air is from a fresh, clean source, or passes through a high-efficiency filter/UV device, then the first configuration is the safest, as contaminated air is highly reduced. The second configuration is also safe, but does not benefit from the outsourcing of potentially-infected air, such as contaminated air flowing through an ineffective filter. The third configuration should be avoided, as the bathroom vent causes air to flow from the hallway, which can be of dubious origin. The CFD simulations also showed that recirculation and swirling regions tend to accumulate the largest concentrations of heavier airborne particles, pathogens, dust, etc. These regions are associated with the largest turbulence kinetic energy (TKE) , and tend to occur in areas with flow recirculation and corners. Therefore, TKE presents a reasonable metric to guide the strategic location of pathogen mitigation devices. The simulations show complex flow patterns with distinct upper and lower flow regions, swirling flow, and significant levels of turbulent mixing. These simulations provide intriguing insights that can be applied to help mitigate pathogen aerosol dispersal, generate building design guidelines, as well as provide insights for the strategic placement of mitigation devices, such as ultraviolet (UV) light, supplemental fans, and filters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Single-Cell Genomics of Novel Actinobacteria With the Wood–Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (~50°C; pH~10.5–11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), “Parcubacteria” (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98–1.00; estimated completeness: 33–93%; estimated genome size: ~2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO 2 fixation via the Wood–Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H + /Na + pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group “OPB41.” We propose to name this bacterium ‘Candidatus Hakubanella thermoalkaliphilus.'

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Rheological Properties of Enzymatically Hydrolyzed Corn Stover Pretreated via Deacetylation and Mechanical Refining

Lignocellulosic biomass is a feedstock for fuels and chemicals that does not compete with food resources and has less contaminants than refuse-derived biomass feedstocks. To convert lignocellulosics to biofuels or value-added products, multiple processing steps are typically necessary. One method of producing biofuels from lignocellulosic biomass utilizes a deacetylation and mechanical refining pretreatment and an enzymatic hydrolysis reaction to produce fermentable sugars from cellulose and hemicellulose. The rheological properties of biomass, such as yield stress and plastic viscosity, change during enzymatic hydrolysis and alter the energy requirements for pumping and mixing, an important consideration for the design of processing equipment. The dynamic changes in rheological properties that occur in a corn stover feedstock undergoing enzymatic hydrolysis are characterized in this work, and the influence on pressure losses in piping systems is estimated. Two rheometer geometries were fabricated with stereolithography 3D printing to reduce wall slip and sample ejection. The slurries have complex rheological behaviors that include shear-thinning behavior. Shear stress ramps were performed on samples at 20 and 50 degrees C using the custom geometries, and the Herschel-Bulkley model was fit to the data. The dynamic nature of the rheological properties is correlated with changes in the average fiber length at various extents of reaction, and the influence of solids concentration on the observed rheology and piping pressure losses is discussed.

09 BIOMASS FUELS↗
DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
DOE OSTI2024

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2025

DeFault: DEep‐Learning‐Based FAULT Delineation Using the IBDP Passive Seismic Data at the Decatur CO2 Storage Site

Abstract The carbon capture, utilization, and storage (CCUS) framework is an essential component in reducing greenhouse gas emissions, with its success hinging on the comprehensive knowledge of subsurface geology and geomechanics. Passive seismic event relocation and fault detection offer vital insights into subsurface structures and the ability to monitor fluid migration pathways. Accurate identification and localization of seismic events, however, face significant challenges, including the necessity for high‐quality seismic data and advanced computational methods. To address these challenges, we introduce a novel deep learning method, , specifically designed for passive seismic source relocation and fault delineating for passive seismic monitoring projects. By leveraging data domain‐adaptation, allows us to train a neural network with labeled synthetic data and apply it directly to field data. Using , the passive seismic sources are automatically clustered based on their recording time and spatial locations, and subsequently, faults and fractures are delineated accordingly. We demonstrate the efficacy of on a field case study involving injection related microseismic data from Decatur, Illinois area. Our approach accurately and efficiently relocated passive seismic events, identified faults and could aid in potential damage induced by seismicity. Our results highlight the potential of as a valuable tool for passive seismic monitoring, emphasizing its role in ensuring CCUS project safety. This research bolsters the understanding of subsurface characterization in CCUS, illustrating machine learning’s capacity to refine these methods. Ultimately, our work has significant implications for CCUS technology deployment, an essential strategy in combating climate change. Plain Language Summary In our quest to tackle climate change, we use a strategy known as carbon capture, utilization, and storage (CCUS) to keep greenhouse gases out of the atmosphere. This strategy relies heavily on our ability to understand what's happening deep under the earth's surface. To make sure we store super critical safely, we need to accurately map out the geological structure, especially faults, but this is tough without high‐quality data and complex computer programs. We've developed a new tool called “DeFault,” which uses advanced machine learning to improve how we find and map these underground features. “DeFault” is smart enough to learn from numerically simulated data and then apply what it’s learned to real‐world situations. It groups together seismic activity—tiny tremors and shifts in the earth—based on when and where they happen, which helps us spot where there might be cracks or faults. We tested “DeFault” in Illinois, where CO 2 is injected underground, and it successfully pinpointed where these tremors occurred and mapped out the faults, helping to prevent accidents accurately in the future. Our study shows that “DeFault” will be a powerful ally in making CCUS safer and more effective, especially for the Illinois Basin Decatur Project. Key Points Faults and fractures introduced by carbon storage can be monitored by passive seismicity DeFault algorithm enables an automatic process for accurate and efficient passive seismic event locating and clustering

58 GEOSCIENCES↗
DOE OSTI2022

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗