Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Forster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to the native protein. Furthermore, numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of “local” Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

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

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Exploring Parameter Redundancy in the Unitary Coupled-Cluster Ansätze for Hybrid Variational Quantum Computing

One of the commonly used chemical-inspired approaches in variational quantum computing is the unitary coupled-cluster (UCC) ansatze. Despite being a systematic way of approaching the exact limit, the number of parameters in the standard UCC ansatze exhibits unfavorable scaling with respect to the system size, hindering its practical use on near-term quantum devices. Efforts have been taken to propose some variants of UCC ansatze with better scaling. In this paper we explore the parameter redundancy in the preparation of unitary coupled-cluster singles and doubles (UCCSD) ansatze employing spin-adapted formulation, small amplitude filtration, and entropy-based orbital selection approaches. Numerical results of using our approach on some small molecules have exhibited a significant cost reduction in the number of parameters to be optimized and in the time to convergence compared with conventional UCCSD-VQE simulations. Further, we also discuss the potential application of some machine learning techniques in further exploring the parameter redundancy, providing a possible direction for future studies.

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

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2020

Contrasting aerosol effects on longwave cloud forcing in South East Asia and Amazon simulated with Community Atmosphere Model version 5.3

Aerosols modify cloud microphysical and radiative properties and thus impact the shortwave (SW) and longwave cloud forcing (LWCF). This study first reports the finding of contrasting aerosol effects on LWCF in South East Asia and Amazon in the Community Atmosphere Model version 5.3 (CAM5.3), which corresponds to the sum of LW indirect and semi-direct effects investigated in Ghan et al., (2012). A series of numerical experiments is conducted to decompose the complex aerosol effects on LWCF. Our analysis indicates that the cooling (negative aerosol effects on LWCF) in Amazon is due mainly to the aerosol effects on warm clouds and the inhibition of vertical motion by the aerosol-induced radiative cooling. In contrast, the warming (positive aerosol effects on LWCF) in South East Asia is due mainly to the aerosol effect on homogeneous freezing, thus reducing the ice particle size and prolonging the existence of ice cloud. Our results emphasize that a comprehensive analysis of integrated aerosol effects on both warm and ice clouds is necessary for better understanding the aerosol effects on LWCF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2021

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Revised Kimberlina Groundwater Simulations

Kimberlina Groundwater Qaulity Simulations (Version 3: 2015 October) by Lawrence Livermore National Labratory. Submitted to the EDX team at the National Energy Technology Laboratory for research use. This dataset contains 562 updated groundwater quality simulations based on the Kimberlina site. It is a six-layer model including the aquifers and aquitards inbetween the ground surface and the reservoir. The input files and associated leakage information for these simulations are included in the zip file "nrap_kim_well31_fluxes and inputs.tar.gz". These 562 simulations are the successful results out of 1000 monte carlo simulations, with time length greater than 5 years. Each simulation consists 6 output files, therefore, there are 3372 files (6*562) in total. One single example simulation of full time length (200 year) is also provided in the zip file "single_simulation_200yr_sim0012.tar.gz". These 6 groundwater quality parameters are: - pH (use absolute value) - TDS (mg/kg) - Pressure (Pa) - CO2 Liquid Saturation (fraction range 0.0-1.0) - CO2 Gas Concentration (mol/L) - CO2 Liquid Concentration (mol/kg) Their correponding file names, take the first simulation for example, are: ph0001.ntab (for pH) , tds0001.ntab (for TDS), usnt0001.p.ntab (for pressure), usnt0001.s_liquid.ntab (for liquid saturation), usnt0001.c_co2_gas.ntab (for CO2 gas concentration) and usnt0001.c_co2_liquid.ntab (for CO2 liquid concentration). Each output file contains at least 19 columns (variables). The total number of columns depends on the simulation time length. The column data include the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the jth grid of y-axis 4. k: the kth grid of z-axis 5. element_ref: element type reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The rest (14, 15, 16...) variables are the simulated parameter values over time, take pH as an example, are: 14. ph000y: initial pH value. 15. ph001y: simulated pH value at the end of the 1st year. 16. ph002y: simulated pH value at the end of the 2nd year. ... (and so on)... The zip input files have the input fluxes of brine and CO2 (in kg/s) allocated to model layers 1,2,3 and 5 (layers 4 and 6 are assumed to be no-flow aquitards) through the following nodes: (src1 WELL 7 7 25 25 10 32) #Etchegoin (src2 WELL 7 7 25 25 33 72) #Macoma Chanac (src3 WELL 7 7 25 25 73 83) #Santa Margarita Mclure (src4 WELL 7 7 25 25 84 89) #Fruitvale Round Mountain (aquitard) (src5 WELL 7 7 25 25 90 96) #Olcese Sand (thief zone) (src6 WELL 7 7 25 25 97 99) #tremblor Freeman (aquitard) ex: brn1 is assigned to src 1 from i=(7,7), j=(25,25), k=(10,32) Edited by Ya-Mei Yang at National Energy Technology Laboratory

brine leakage↗
DOE OSTI2020

Implementing a microphysics model in hydrodynamic simulations to study the initial plasma formation in dielectric ablator materials for direct-drive implosions

A microphysics model to describe the photoionization and impact ionization processes in dielectric ablator materials like plastic has been implemented into the 1-D hydrodynamic code LILAC for planar and spherical targets. Currently, the initial plasma formation during the early stages of a laser drive are modeled in an ad hoc manner, until the formation of a critical surface. Implementation of the physics-based models predict higher values of electron density, electron temperature, and pressure than the ad hoc model. Moreover, the numerical predictions are consistent with previous experimental observations of the shinethrough mechanism in plastic ablators. For planar targets, a decompression of the rear end of the target was observed, that is similar to recent experiments. An application of this model is to understand the laser-imprint mechanism that is caused by nonuniform laser irradiation due to single beam speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗
DOE OSTI2024

The Effect of Normal Stress Oscillations on Fault Slip Behavior Near the Stability Transition From Stable to Unstable Motion

Tectonic fault zones are subject to normal stress variations with a wide range of spatio-temporal scales. Stress perturbations cover a wide range of frequencies and amplitudes from high frequency seismic waves generated by earthquakes to low frequency transients associated with solid Earth tides. These perturbations can reactivate critically stressed faults and trigger earthquakes. Here, we describe lab experiments to illuminate the physics of such changes in friction and the mechanics of earthquake triggering and fault reactivation. Friction tests were done in a double direct shear configuration for conditions near the stability transition from stable to unstable motion. We studied simulated fault gouge composed of quartz powder and conducted experiments at reference normal stress from 10 to 13.5 MPa. After shearing to steady state sliding, we applied sinusoidal normal stress oscillations of amplitude 0.5–2 MPa, and period of 0.5–50 s. We performed numerical simulations using measured values of rate and state friction (RSF) parameters to assess our data. Our results show that low frequency stress oscillations cause a Coulomb-like response of shear strength that transitions from stable slip to slow lab earthquakes as frequency increases. At the critical frequency predicted by RSF we observe periodic stick-slip behavior. Perturbations of high amplitude and short period weaken the fault, while lower amplitudes strengthen the fault. We find that a modified RSF formulation is able to accurately match our laboratory data. Our findings highlight the complex effects of stress perturbations on fault strength and the mode of fault slip.

58 GEOSCIENCES↗
DOE OSTI2020

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Nuclear Astrophysics through simulations of neutron star mergers using Monte‐Carlo neutrino radiation transport (DE-SC0020435 Final Technical Report)

Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Current measurements in the intermittent-contact mode of atomic force microscopy using the Fourier method: a feasibility analysis

Atomic force microscopy (AFM) is an important tool for measuring a variety of nanoscale surface properties, such as topography, viscoelasticity, electrical potential and conductivity. Some of these properties are measured using contact methods (static contact or intermittent contact), while others are measured using noncontact methods. Some properties can be measured using different approaches. Conductivity, in particular, is mapped using the contact-mode method. However, this modality can be destructive to delicate samples, since it involves continuously dragging the cantilever tip on the surface during the raster scan, while a constant tip–sample force is applied. In this paper we discuss a possible approach to develop an intermittent-contact conductive AFM mode based on Fourier analysis, whereby the measured current response consists of higher harmonics of the cantilever oscillation frequency. Such an approach may enable the characterization of soft samples with less damage than contact-mode imaging. To explore its feasibility, we derive the analytical form of the tip–sample current that would be obtained for attractive (noncontact) and repulsive (intermittent-contact) dynamic AFM characterization, and compare it with results obtained from numerical simulations. Although significant instrumentation challenges are anticipated, the modelling results are promising and suggest that Fourier-based higher-harmonics current measurement may enable the development of a reliable intermittent-contact conductive AFM method.

14 SOLAR ENERGY↗
DOE OSTI2024

Relativistic Alfvén Turbulence at Kinetic Scales

Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close to k −3 or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics.

79 ASTRONOMY AND ASTROPHYSICS↗