Search NASA⌕ Search

SEARCH · Search NASA

Results for “process 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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
DOE OSTI2010

Nucleotide Docking: Prediction of Reactant State Complexes for Ribonuclease Enzymes

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools - MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0 - for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗
DOE OSTI2019

Vehicle energy consumption estimation using large scale simulations and machine learning methods

Vehicle energy consumption and powertrain operations for future vehicle powertrain technologies are predicted through full-vehicle simulations using Autonomie, a simulation tool developed at Argonne National Laboratory. Over 1 million vehicles, with different powertrain types (conventional and numerous electrified powertrains) and component technologies (advanced engines, batteries, electrical machines, light weighting) are simulated for outputs including vehicle fuel economy and cost. Despite the ability of Autonomic to predict fuel economy given various combination of vehicle inputs, the full-vehicle simulation is very time-consuming. It takes approximately 84 h to simulate for 33,060 vehicles even using distributed computing on a cluster of 128 worker nodes. Moreover, Autonomic is unable to populate a continuous output space through simulation since it only yields outputs with discrete values instead of an input-output relation. This paper proposes a novel large-scale learning and prediction process (LSLPP) via machine learning approaches to answer the need for a continuous space of vehicle fuel economy and a more efficient simulation process. By learning in a supervised fashion to discover structures in data, the machine learning approaches seek to train a model that can explain the data with which we can perform prediction over unseen data. As such, we are able to accomplish analytic continuation on the current discrete space of simulation results. In addition, we present a random sampling -based numerosity reduction algorithm incorporated in LSLPP that reduces simulation runs and improves efficiency of data procurement. We analytically derive a theoretical bound in order to guarantee the effectiveness of our algorithm. We also compare random sampling to stratified sampling and suggest sampling strategies. A user-friendly tool is developed to support the process. Our experimental results confirm that the proposed LSLPP is able to greatly accelerate prediction and analysis of vehicle fuel economy. The output of this study is used by the US government to evaluate the impact of its R&D funding on energy security and CO 2 emission as well as for the setting of Corporate Average Fuel Economy (CAFE) standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

HPC Molecular Simulation Tries Out a New GPU: Experiences on Early AMD Test Systems for the Frontier Supercomputer

Molecular simulation is an important tool for nu- merous efforts in physics, chemistry, and the biological sciences. Simulating molecular dynamics requires extremely rapid cal- culations to enable sufficient sampling of simulated temporal molecular processes. The Hewlett Packard Enterprise (HPE) Cray EX Frontier supercomputer installed at the Oak Ridge Leadership Computing Facility (OLCF) will provide an exascale resource for open science, and will feature graphics processing units (GPUs) from Advanced Micro Devices (AMD). The future LUMI supercomputer in Finland will be based on an HPE Cray EX platform as well. Here we test the ports of several widely used molecular dynamics packages that have each made substantial use of acceleration with NVIDIA GPUs, on Spock, the early Cray pre-Frontier testbed system at the OLCF which employs AMD GPUs. These programs are used extensively in industry for pharmaceutical and materials research, as well as academia, and are also frequently deployed on high-performance computing (HPC) systems, including national leadership HPC resources. We find that in general, performance is competitive and installation is straightforward, even at these early stages in a new GPU ecosystem. Our experiences point to an expanding arena for GPU vendors in HPC for molecular simulation.

Sedova, Ada↗
DOE OSTI2013

Recovery of yttrium from cathode ray tubes and lamps’ fluorescent powders: experimental results and economic simulation

Highlights: • Fluorescent powder of lamps. • Fluorescent powder of cathode ray rubes. • Recovery of yttrium from fluorescent powders. • Economic simulation for the processes to recover yttrium from WEEE. - Abstract: In this paper, yttrium recovery from fluorescent powder of lamps and cathode ray tubes (CRTs) is described. The process for treating these materials includes the following: (a) acid leaching, (b) purification of the leach liquors using sodium hydroxide and sodium sulfide, (c) precipitation of yttrium using oxalic acid, and (d) calcinations of oxalates for production of yttrium oxides. Experimental results have shown that process conditions necessary to purify the solutions and recover yttrium strongly depend on composition of the leach liquor, in other words, whether the powder comes from treatment of CRTs or lamp. In the optimal experimental conditions, the recoveries of yttrium oxide are about 95%, 55%, and 65% for CRT, lamps, and CRT/lamp mixture (called MIX) powders, respectively. The lower yields obtained during treatments of MIX and lamp powders are probably due to the co-precipitation of yttrium together with other metals contained in the lamps powder only. Yttrium loss can be reduced to minimum changing the experimental conditions with respect to the case of the CRT process. In any case, the purity of final products from CRT, lamps, and MIX is greater than 95%. Moreover, the possibility to treat simultaneously both CRT and lamp powders is very important and interesting from an industrial point of view since it could be possible to run a single plant treating fluorescent powder coming from two different electronic wastes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2011

Passive scalar transport in peripheral regions of random flows

We investigate statistical properties of the passive scalar mixing in random (turbulent) flows assuming its diffusion to be weak. Then at advanced stages of the passive scalar decay, its unmixed residue is primarily concentrated in a narrow diffusive layer near the wall and its transport to the bulk goes through the peripheral region (laminar sublayer of the flow). We conducted Lagrangian numerical simulations of the process for different space dimensions d and revealed structures responsible for the transport, which are passive scalar tongues pulled from the diffusive boundary layer to the bulk. We investigated statistical properties of the passive scalar and of the passive scalar integrated along the wall. Moments of both objects demonstrate scaling behavior outside the diffusive boundary layer. We propose an analytic scheme for the passive scalar statistics, explaining the features observed numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Non-oscillatory flux correlation functions for efficient nonadiabatic rate theory

There is currently much interest in the development of improved trajectory-based methods for the simulation of nonadiabatic processes in complex systems. An important goal for such methods is the accurate calculation of the rate constant over a wide range of electronic coupling strengths and it is often the nonadiabatic, weak-coupling limit, which being far from the Born-Oppenheimer regime, provides the greatest challenge to current methods. We show that in this limit there is an inherent sign problem impeding further development which originates from the use of the usual quantum flux correlation functions, which can be very oscillatory at short times. From linear response theory, we derive a modified flux correlation function for the calculation of nonadiabatic reaction rates, which still rigorously gives the correct result in the long-time limit regardless of electronic coupling strength, but unlike the usual formalism is not oscillatory in the weak-coupling regime. In particular, a trajectory simulation of the modified correlation function is naturally initialized in a region localized about the crossing of the potential energy surfaces. In the weak-coupling limit, a simple link can be found between the dynamics initialized from this transition-state region and an generalized quantum golden-rule transition-state theory, which is equivalent to Marcus theory in the classical harmonic limit. This new correlation function formalism thus provides a platform on which a wide variety of dynamical simulation methods can be built aiding the development of accurate nonadiabatic rate theories applicable to complex systems.

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

Atmospheric river representation in the Energy Exascale Earth System Model (E3SM) version 1.0

Abstract. The Energy Exascale Earth System Model (E3SM) project is an ongoing, state-of-the-science Earth system modeling, simulation, and prediction project developed by the US Department of Energy (DOE). With an emphasis on supporting the DOE's energy mission, understanding and quantifying how well the model simulates water cycle processes is of particular importance. Here, we evaluate E3SM version 1.0 (v1.0) for its ability to represent atmospheric rivers (ARs), which play significant roles in water vapor transport and precipitation. The characteristics and precipitation associated with global ARs in E3SM at standard resolution (1∘ × 1∘) are compared to the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA2). Global patterns of AR frequencies in E3SM show high degrees of correlation (≥0.97) with MERRA2 and low mean absolute errors (MAEs; <1 %) annually, seasonally, and across different ensemble members. However, some large-scale condition biases exist, leading to AR biases – most significant of which are the double intertropical convergence zone (ITCZ), a stronger and/or equatorward-shifted subtropical jet during boreal and austral winters, and enhanced Northern Hemisphere westerlies during summer. By comparing atmosphere-only and fully coupled simulations, we attribute the sources of the biases to the atmospheric component or to a coupling response. Using relationships revealed in Dong et al. (2021), we provide evidence showing the stronger North Pacific jet in winter and the enhanced Northern Hemisphere westerlies during summer, associated with E3SM's double ITCZ and related weaker Atlantic meridional overturning circulation (AMOC), respectively, which are significant sources of the AR biases found in the coupled simulations.

58 GEOSCIENCES↗
DOE OSTI2014

Atomic-scale investigation of point defects and hydrogen-solute atmospheres on the edge dislocation mobility in alpha iron

In this study, we present atomistic mechanisms of 1/2 [111](11{sup ¯}0) edge dislocation interactions with point defects (hydrogen and vacancies) and hydrogen solute atmospheres in body centered cubic (bcc) iron. In metals such as iron, increases in hydrogen concentration can increase dislocation mobility and/or cleavage-type decohesion. Here, we first investigate the dislocation mobility in the presence of various point defects, i.e., change in the frictional stress as the edge dislocation interacts with (a) vacancy, (b) substitutional hydrogen, (c) one substitutional and one interstitial hydrogen, (d) interstitial hydrogen, (e) vacancy and interstitial hydrogen, and (f) two interstitial hydrogen. Second, we examine the role of a hydrogen-solute atmosphere on the rate of local dislocation velocity. The edge dislocation simulation with a vacancy in the compression side of the dislocation and an interstitial hydrogen atom at the tension side exhibit the strongest mechanical response, suggesting a higher potential barrier and hence, the higher frictional stress (i.e., ∼83% higher than the pure iron Peierls stress). In the case of a dislocation interacting with a vacancy on the compressive side, the vacancy binds with the edge dislocation, resulting in an increase in the friction stress of about 28% when compared with the Peierls stress of an edge dislocation in pure iron. Furthermore, as the applied strain increases, the vacancy migrates through a dislocation transportation mechanism by attaining a velocity of the same order as the dislocation velocity. For the case of the edge dislocation interacting with interstitial hydrogen on the tension side, the hydrogen atom jumps through one layer perpendicular to the glide plane during the pinning-unpinning process. Finally, our simulation of dislocation interactions with hydrogen show first an increase in the local dislocation velocity followed by a pinning of the dislocation core in the atmosphere, resulting in resistance to dislocation motion as the dislocation moves though the hydrogen-solute atmospheres. With this systematic, atomistic study of the edge dislocation with various point defects, we show significant increase in obstacle strengths in addition to an increase in the local dislocation velocity during interaction with solute atmospheres. The results have implications for constitutive development and modeling of the hydrogen effect on dislocation mobility and deformation in metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Southeast Pacific atmospheric composition and variability sampled along 20˚S during VOCALS-REx

The VAMOS Ocean-Climate-Atmosphere-Land Regional Experiment (VOCALS-REx) was conducted from 15 October to 15 November 2008 in the South East Pacific region to investigate interactions between land, sea and atmosphere in this unique tropical eastern ocean environment and to improve the skill of global and regional models in representing the region. This study synthesises selected aircraft, ship and surface site observations from VOCALS-REx to statistically summarise and characterise the atmospheric composition and variability of the Marine Boundary Layer (MBL) and Free Troposphere (FT) along the 20{sup o} S parallel between 70{sup o} W and 85{sup o} W. Significant zonal gradients in mean MBL sub-micron aerosol particle size and composition, carbon monoxide, ozone and sulphur dioxide were seen over the campaign, with a generally more variable and polluted coastal environment and a less variable, more pristine remote maritime regime. Gradients are observed to be associated with strong gradients in cloud droplet number. The FT is often more polluted in terms of trace gases than the MBL in the mean; however increased variability in the FT composition suggests an episodic nature to elevated concentrations. This is consistent with a complex vertical interleaving of airmasses with diverse sources and hence pollutant concentrations as seen by generalised back trajectory analysis, which suggests contributions from both local and long-range sources. Furthermore, back trajectory analysis demonstrates that the observed zonal gradients both in the boundary layer and the free troposphere are characteristic of marked changes in airmass history with distance offshore - coastal boundary layer airmasses having been in recent contact with the local land surface and remote maritime airmasses having resided over ocean for in excess of ten days. Boundary layer composition to the east of 75{sup o} W was observed to be dominated by coastal emissions from sources to the west of the Andes, with evidence for diurnal pumping of the Andean boundary layer above the height of the marine capping inversion. The climatology presented here aims to provide a valuable dataset to inform model simulation and future process studies, particularly in the context of aerosol-cloud interaction and further evaluation of dynamical processes in the SEP region for conditions analogous to those during VOCALS-REx.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Heterogeneous chain dynamics and aggregate lifetimes in precise acid-containing polyethylenes: Experiments and simulations

Melt state dynamics for a series of strictly linear polyethylenes with precisely spaced associating functional groups were investigated. The periodic pendant acrylic acid groups form hydrogen-bonded acid aggregates within the polyethylene (PE) matrix. The dynamics of these nanoscale heterogeneous morphologies were investigated from picosecond to nanosecond timescales by both quasi-elastic neutron scattering (QENS) measurements and fully atomistic molecular dynamics (MD) simulations. Two dynamic processes were observed. The faster dynamic processes which occur at the picosecond timescales are compositionally insensitive and indicative of spatially restricted local motions. The slower dynamic processes are highly composition dependent and indicate the structural relaxation of the polymer backbone. Higher acid contents, or shorter PE spacers between pendant acid groups, slow the structural relaxation timescale and increase the stretching parameter (β) of the structural relaxation. Additionally, the dynamics of specific hydrogen atom positions along the backbone correlate structural heterogeneity imposed by the associating acid groups with a mobility gradient along the polymer backbone. At time intervals (<2 ns), the mean-squared displacements for the four methylene groups closest to the acid groups are up to 10 times smaller than those of methylene groups further from the acid groups. At longer timescales acid aggregates rearrange and the chain dynamics of the slow, near-aggregate regions and the faster bridge regions converge, implying a characteristic timescale for the passage of chains between aggregates. As a result, the characterization of the nanoscale chain dynamics in these associating polymer systems both provides validation of simulation force fields and provides understanding of heterogeneous chain dynamics in associating polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Final report for the DOE Early Career Award #DE-SC0003912

This DoE supported early career project was aimed at developing computational models, theory and simulation methods that would be then be used to predict assembly and morphology in polymer nanocomposites. In particular, the focus was on composites in active layers of devices, containing conducting polymers that act as electron donors and nanoscale additives that act as electron acceptors. During the course this work, we developed the first of its kind molecular models to represent conducting polymers enabling simulations at the experimentally relevant length and time scales. By comparison with experimentally observed morphologies we validated these models. Furthermore, using these models and molecular dynamics simulations on graphical processing units (GPUs) we predicted the molecular level design features in polymers and additive that lead to morphologies with optimal features for charge carrier behavior in solar cells. Additionally, we also predicted computationally new design rules for better dispersion of additives in polymers that have been confirmed through experiments. Achieving dispersion in polymer nanocomposites is valuable to achieve controlled macroscopic properties of the composite. The results obtained during the course of this DOE funded project enables optimal design of higher efficiency organic electronic and photovoltaic devices and improve every day life with engineering of these higher efficiency devices.

14 SOLAR ENERGY↗
DOE OSTI2017

The Turbulent Chiral Magnetic Cascade in the Early Universe

The presence of asymmetry between fermions of opposite handedness in plasmas of relativistic particles can lead to exponential growth of a helical magnetic field via a small-scale chiral dynamo instability known as the chiral magnetic effect. Here, we show, using dimensional arguments and numerical simulations, that this process produces through the Lorentz force chiral magnetically driven turbulence. A k{sup −2} magnetic energy spectrum emerges via inverse transfer over a certain range of wavenumbers k. The total chirality (magnetic helicity plus normalized chiral chemical potential) is conserved in this system. Therefore, as the helical magnetic field grows, most of the total chirality gets transferred into magnetic helicity until the chiral magnetic effect terminates. Quantitative results for height, slope, and extent of the spectrum are obtained. Consequences of this effect for cosmic magnetic fields are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗