Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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,243 records · Page 69

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗

Adaptive, problem-tailored variational quantum eigensolver mitigates rough parameter landscapes and barren plateaus

Abstract Variational quantum eigensolvers (VQEs) represent a powerful class of hybrid quantum-classical algorithms for computing molecular energies. Various numerical issues exist for these methods, however, including barren plateaus and large numbers of local minima. In this work, we consider the Adaptive, Problem-Tailored Variational Quantum Eiegensolver (ADAPT-VQE) ansätze, and examine how they are impacted by these local minima. We find that while ADAPT-VQE does not remove local minima, the gradient-informed, one-operator-at-a-time circuit construction accomplishes two things: First, it provides an initialization strategy that can yield solutions with over an order of magnitude smaller error compared to random initialization, and which is applicable in situations where chemical intuition cannot help with initialization, i.e., when Hartree-Fock is a poor approximation to the ground state. Second, even if an ADAPT-VQE iteration converges to a local trap at one step, it can still “burrow” toward the exact solution by adding more operators, which preferentially deepens the occupied trap. This same mechanism helps highlight a surprising feature of ADAPT-VQE: It should not suffer optimization problems due to barren plateaus and random initialization. Even if such barren plateaus appear in the parameter landscape, our analysis suggests that ADAPT-VQE avoids such regions by design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THERMAL MODELING OF HANFORD CESIUM AND STRONTIUM CANISTERS DURING SIMULATED LOADING

A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies within a canister and overpack for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to make pretest numerical predictions for the behavior of the heated canister during the simulated radiolytic decay heat testing, which simulates the dry storage system during loading operations. The simulated radiolytic decay heat test is planned for mid-2024 in a configuration that includes the heater assembly, overpack, and canister, but with the lids removed to allow loading cesium and strontium capsules into the canister. One of the goals of the test is to evaluate the thermal behavior of the canister and overpack assembly in the ambient air of the test facility, which will provide data critical to validating the thermal models and understanding how the HLC will perform as a system once deployed. To best approximate real-world conditions, the CFD model includes the full air volume of the mock-up truck bay the heated canister test will be performed in, enabling detailed investigation of how the heated canister affects airflow around it. Rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2028 timeframe. This study presents the pre-test temperature predictions of the simulated radiolytic decay heat test. A description of the heater assembly, canister, and overpack system is presented. The model was developed with the commercial CFD software STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was preformed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.

Carpenter-Graffy, Dina E.↗

Preliminary Study on TRISO Fuel Cross Section Generation

Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A novel defoamer for processing nuclear waste: Testing and performance

Abstract Legacy radioactive waste from nuclear weapons material processing is the biggest environmental challenge in the state of South Carolina, and one of the biggest in the United States. Although substantial progress has been made in processing and vitrifying high‐level radioactive waste at the Savannah River site, approximately 35 million gallons remains to be treated and dispositioned. In this article, we show the development of a new defoamer for use in the processing of high‐level radioactive waste. The application of the new defoamer is not only more effective at controlling foam but will shorten the batch processing time leading to safer and more efficient processing. This has the potential to help shorten the overall site mission lifetime, saving the federal government hundreds of millions of dollars. The novel aspect of the new defoamer is its effectiveness and stability in the harsh conditions needed for processing high‐level radioactive waste: high temperatures, pH 4–13, strong oxidizing and reducing agents, and numerous metal catalysts. The defoamer will replace an existing antifoam agent that decomposes to form multiple flammable gases. Because the existing antifoam agent decomposes quickly, processing changes increased processing time. Statement of novelty We demonstrate that a new defoamer has been developed for use in processing high‐level radioactive waste. This novel defoamer is effective and stable in the harsh processing conditions with no detectable flammable by‐products. Because the existing antifoam agent decomposes quickly, processing changes, including lowering acid addition rates and evaporation rates, increased processing time. The new defoamer is not only more effective at controlling foam but will shorten the batch processing time and lead to safer and more efficient processing by eliminating the production of flammable gases. Statement of industrial relevance Foaming is of great concern in many industrial processes involving three‐phase gas/liquid/fine‐solids systems, like in water evaporation, mineral floatation, air sparging in situ remediation techniques, and in those found in the paper industry.

Lambert, Dan P.↗

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Energy performance augmentation of domestic refrigerators with microchannel condensers

Refrigerators are one of the most extensively used household appliances, accounting for approximately 4% of total household electricity consumption. Enhancing their energy performance can significantly reduce residential energy demand. As a result, this study focuses on reducing the energy consumption of a French-door, bottom-mount refrigerator using isobutane (R600a) as the working fluid. Here, the performance improvement of the refrigerators has been obtained mainly through improving the vapor compression refrigeration cycle, with a particular emphasis on condenser design. This study explores the impact of three-pass serpentine microchannel condenser on energy performance, along with the effects of two distinct fan combinations. Charge optimization experiments were conducted to identify the optimal refrigerant mass for the refrigerator unit with microchannel condenser. Additionally, a comprehensive numerical model was developed to analyze the behavior of vapor compression cycle under various refrigerant charges. The experimental findings revealed that the refrigerator unit with the three-pass microchannel condenser reduces the energy consumption of isobutane-based domestic refrigerators up to 16% while reducing the refrigerant charge by 11% in comparison to the refrigerator unit with conventional wire-and-tube condenser.

Charge optimization↗

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lattice quantum chromodynamics at large isospin density

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of n-pion correlation functions for n ϵ {1, … , N} using a single N × N matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 charged pions using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass m π = 170 MeV and spacetime volumes of (4.4 3 × 8.8) fm 4 and (5.8 3 × 11.6) fm 4 . We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well-approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential μ I ~ 1.5m π , signaling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound, c s , in the medium is seen to exceed the ideal-gas (conformal) limit ($c^{2}_{s} ≤ 1/3)$ over a wide range of chemical potential before falling towards the asymptotic expectation at μ I ~ 15m π . These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integration Study of Converter-Interfaced Combined Heat and Power Plants

Combined Heat and Power (CHP) systems, as a proven technology, provide numerous benefits to the plant owner by saving electricity bill, reducing carbon dioxide emissions, and improving electricity service reliability. However, their integration into the distribution grid faces many barriers including the complexity of grid code requirements and the lack of technical expertise in commercial and small industrial plants to interface with the utility and accelerate the interconnection process. To solve these issues, this study proposes a converter-interfaced solution for CHP integration. Simulations including various fault scenarios and load flow under different grid conditions are carried out to compare technical performances of the proposed solution with directly-coupled CHP. Results show that the generator size of converter-interfaced CHP can be reduced by approximately 25% and it also provides advantageous features such as reduced stress level under fault conditions and a better dynamic stability margin. In addition, a power factor controller is designed and reactive power from CHP is dispatched in such a way that the CHP plant is compliant with power factor requirements from utilities under all operating conditions.

Converter, Combined heat and power, Distributed en↗

Learning Stochastic Parametric Differentiable Predictive Control Policies

We present a scalable unsupervised learning-based method for obtaining explicit control policies for model predictive control problems for stochastic linear systems with additive uncertainties subject to nonlinear chance constraints. We call the proposed method stochastic parametric differentiable predictive control (SP-DPC), which extends the recently proposed deterministic DPC policy optimization algorithm. We formulate the SP-DPC as a deterministic approximation to the stochastic parametric constrained optimal control problem via independent sampling of the problem's parameters and uncertainties. This formulation allows us to directly compute the policy gradients via automatic differentiation of the problem's value function, evaluated over sampled parameters and uncertainties. In particular, the computed expectation of the problem's value function is backpropagated through the finite-time closed-loop system rollouts parametrized by a known nominal system dynamics model and neural control policy. We also provide theoretical probabilistic guarantees on closed-loop stability and chance constraints satisfaction for systems controlled by learned neural policies. We demonstrate the computational efficiency and scalability of the proposed policy optimization algorithm in three numerical examples, including systems with a large number of states or subject to nonlinear constraints.

Drgona, Jan↗

Experimental Evaluation of Interfacial Bonding Strength Between 304 Stainless Steel Substrate and Electrodeposited Nickel Coating Using Mesoscale Mechanical Testing Methods

Electrodeposition is a commonly used method for depositing a metal layer on a metallic substrate, to provide a decorative finish or for functional purposes including wear and corrosion resistance. The interfacial bonding strength is a critical factor in determining the quality of electrodeposition, as it ensures the adhesion of the deposited layer to the substrate. Despite its importance, there has not previously been a suitable mechanical testing method to quantitatively characterize the bonding strength between the electrodeposited layers and the substrate, due to the high strength of the materials and to restrictions imposed by the geometry (due to the small thickness of the plating layer) to manufacture tensile bars. In this study, we introduce a novel mesoscale mechanical testing method to overcome these limitations. This technique was applied to assess the bonding strength between a 304 stainless steel (SS) substrate and electrodeposited nickel with three different plating formulae (nickel sulfamate, Watts bath, and hard nickel). The thickness of the pure nickel coatings achieved for each of the three electrolyte solutions was approximately 1 mm on each side of the substrate. Mesoscale tensile bars of 1 mm length were manufactured by a femtosecond laser, with the material interface at the center of the gauge section, and then tensile-tested with digital image correlation. Further, the results proved that the electrodeposition technique is able to produce a very high bonding strength that is close to the yield strength of both the substrate material and the electrodeposited nickel layer. Additionally, in the case of the Watts bath formula, the interfacial bonding strength between the SS 304 and electrodeposited nickel can exceed that of the nickel layer. This method is expected to be useful for quantifying the interfacial bonding strength in numerous applications.

36 MATERIALS SCIENCE↗

Turbulent transport and reactions of plant-emitted hydrocarbons in an Amazonian rain forest

The processes governing the temporal and spatial patterns of isoprene and monoterpenes emitted by a rainforest in the central Amazon region of Brazil is investigated using a combination of field experiments and numerical simulations. Specifically, Large Eddy Simulations (LES) were used to resolve emissions of isoprene and monoterpenes, turbulent transport, and air chemistry. The coupled chemistry-transport LES included the effects of isoprene and monoterpenes reactivity due to reactions with hydroxyl radical and ozone. The LES results are used to compute vertically resolved budgets of isoprene and monoterpenes in the rainforest canopy in response to emissions, turbulent transport, surface deposition, and air chemistry. Our results indicated that emission and dispersion dominated the isoprene budget as the gases were transported out of the canopy space. In a region limited by nitrogen oxides (with prevailing nitric oxide levels of < 0.5 parts per billion), the in-canopy chemical destruction removed approximately 10% of locally emitted monoterpenes. Hydroxyl radical production rates from the ozonolysis of monoterpenes amounted to ≈ 2 x 10 6 radicals cm -3 s -1 and had similar magnitude to the light-dependent hydroxyl radical formation. One key conclusion was that the Amazonia rainforest abundantly emitted monoterpenes whose in-canopy ozonolysis yielded hydroxyl radicals in amounts similar to the magnitude of light-dependent formation. Reactions of monoterpenes and isoprene with hydroxyl radical and ozone were necessary for the maintenance of the Amazon rainforest canopy as a photochemically active environment suitable to generate oxidants and secondary organic aerosols.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep learning insights into cosmological structure formation

The evolution of linear initial conditions present in the early Universe into extended halos of dark matter at late times can be computed using cosmological simulations. However, a theoretical understanding of this complex process remains elusive; in particular, the role of anisotropic information in the initial conditions in establishing the final mass of dark matter halos remains a long-standing puzzle. Here, we build a deep learning framework to investigate this question. We train a three-dimensional convolutional neural network to predict the mass of dark matter halos from the initial conditions, and quantify in full generality the amounts of information in the isotropic and anisotropic aspects of the initial density field about final halo masses. We find that anisotropies add a small, albeit statistically significant amount of information over that contained within spherical averages of the density field about final halo mass. However, the overall scatter in the final mass predictions does not change qualitatively with this additional information, only decreasing from 0.9 dex to 0.7 dex. Given such a small improvement, our results demonstrate that isotropic aspects of the initial density field essentially saturate the relevant information about final halo mass. Therefore, instead of searching for information directly encoded in initial conditions anisotropies, a more promising route to accurate, fast halo mass predictions is to add approximate dynamical information based e.g. on perturbation theory. More broadly, our results indicate that deep learning frameworks can provide a powerful tool for extracting physical insight into cosmological structure formation. Published by the American Physical Society 2024

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗