Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “range prediction”

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,063 records · Page 59

High-pressure melt dynamics in shock-compressed titanium

In this work we study the high-pressure melting behavior of titanium using laser-driven shock compression with in situ femtosecond x-ray diffraction and molecular-dynamics simulations based on a machine-learned interatomic potential. The MD simulations predict the solid-liquid coexistence on the Hugoniot in the ∼111−124GPa range. Experimentally, we observe the first evidence of liquid at 86 GPa. We also observe pronounced microstructural changes with pressure, with strong grain refinement associated with the emergence of liquid, within the solid-liquid coexistence (∼110−126GPa). Above 126 GPa, we observe the persistence of residual levels of highly textured crystalline Ti to ∼180GPa, well above the expected melt completion pressure. We discuss the accuracy that current laser-shock experimental platforms have at determining the melt onset and completion pressures.

36 MATERIALS SCIENCE↗

Generalized parton distributions from symbolic regression

AI/ML informed Symbolic Regression is the next stage of scientific modeling. We utilize a highly customizable symbolic regression package “PySR” to model the x and t dependence of the flavor isovector combination Hu-d(x, t, ) at = 0. These PySR models were trained on GPD results provided by both Lattice QCD and phenomenological sources GGL, GK, and VGG. We demonstrate, for the first time, the consistency and systematic convergence of Symbolic Regression by quantifying the disparate models through their Taylor expansion coefficients. In addition to PySR penalizing models with higher complexity and mean-squared error, we implement schemes that test specific physics hypotheses, including force-factorized x and t dependence and Regge behavior in PySR GPDs. We show that PySR can identify factorizing GPD sources based on their response to the Force-Factorized model. Here, knowing the precise behavior of the GPDs, and their uncertainties in a wide range in x and t, crucially impacts our ability to concretely and quantitatively predict hadronic spatial distributions and their derived quantities.

Dotson, Andrew [New Mexico State University, Las C↗

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Community size-spectra applied to recreational freshwater fisheries in Puerto Rican reservoirs

We developed a management framework based on community size-spectra (CSS) to establish targets to maximize desirable gamefish. We evaluated the influence of relative abundance of fish functional groups and watershed characteristics on the slope of the CSS or ecological efficiency. Both ecological efficiency and capacity were strongly influenced by relative abundance of fish functional groups and the degree of water level variation. Steeper slopes indicated overfished gamefish populations, whereas flatter slopes indicated systems dominated by undesirable invasive large-bodied low trophic position fish. We hypothesize that a desirable fishery would have a target range of CSS slopes that provides a deterministic community-based framework of predictable responses to environmental and food web changes that can be used to evaluate outcomes of management actions and a means to assess annual progress toward desirable conditions.

60 APPLIED LIFE SCIENCES↗

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗

Improved Parametric Models for Explosion Pressure Signals Derived From Large Datasets

Accurate recording and characterization of explosion-induced pressure signals are key components of the forensic analysis of explosion events in the atmosphere. Parametric overpressure models based on several key waveform features (peak overpressure, positive pulse duration, and impulse) are widely used to estimate explosion energy in terms of trinitrotoluene equivalent yield. However, those models are often developed by a limited dataset, including only a few events or recordings at relatively short propagation distances. Here, we develop empirical waveform-parameter models based on a regression analysis of a large set of data curated from four chemical explosion experiments including 16 detonations. We measured peak overpressure and impulse for positive and negative phases from 1000 pressure signals recorded at local ranges (<20 km) with scaled distance up to 8000 m/kg 1/3 . Additionally, the measured waveform parameters showed large variation with respect to observing distances indicating the effects of atmospheric propagation. In this study, a second-order polynomial model was used in a least-squares regression to account for those propagation effects and to improve data fitting. In addition to model parameters for waveform features, we also determined range-dependent model uncertainties based on data variance. The model uncertainty represents the prediction error of our models and can be critical to evaluating the uncertainty of yield estimate.

58 GEOSCIENCES↗

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗

ASCR@40: Highlights and Impacts of ASCR’s Programs

A report compiled by the ASCAC Subcommittee on the 40-year history of ASCR for the U.S. Department of Energy’s Office of Advanced Scientific Computing Research. The Office of Advanced Scientific Computing Research (ASCR) sits within the Office of Science in the Department of Energy (DOE). Per their web pages, “the mission of the ASCR program is to discover, develop, and deploy computational and networking capabilities to analyze, model, simulate, and predict complex phenomena important to the DOE.” This succinct statement encompasses a wide range of responsibilities for computing and networking facilities; for procuring, deploying, and operating high performance computing, networking, and storage resources; for basic research in mathematics and computer science; for developing and sustaining a large body of software; and for partnering with organizations across the Office of Science and beyond. While its mission statement may seem very contemporary, the roots of ASCR are quite deep—long predating the creation of DOE. Applied mathematics and advanced computing were both elements of the Theoretical Division of the Manhattan Project. In the early 1950s, the Manhattan Project scientist and mathematician John von Neumann, then a commissioner for the AEC (Atomic Energy Commission), advocated for the creation of a Mathematics program to support the continued development and applications of digital computing. Los Alamos National Laboratory (LANL) scientist John Pasta created such a program to fund researchers at universities and AEC laboratories. Under several organizational name changes, this program has persisted ever since, and would eventually grow to become ASCR.

97 MATHEMATICS AND COMPUTING↗

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗

PRIMED for the Future: Purposing Raw Intake for Machine Learning-Enabled Detection (Final Report)

The COVID-19 pandemic demonstrated how a novel, elusive, and diffuse biological threat can engender uncertainty and misinformation, and it underscored the need for flexible analytical modalities agnostic to the identity of biological material. Yet even before the pandemic, recognition of the limitations of the current, list-based approach, which focuses on known pathogens and biotoxins, and of the importance of agent-agnostic biodetection was growing within the biosecurity community. In a 2018 report on “Biodefense in the Age of Synthetic Biology,” for example, the National Academy of Sciences stated that “an overreliance on the Select Agent List is a systematic weakness affecting many aspects of the United States’ current biodefense mitigation capability”. More recently, a group of biodefense researchers proposed the identification and adoption of “bioagent-agnostic signatures (BASs)” as a way of detecting and characterizing not only existing agents but also novel ones, an approach they believe will “enable a more flexible and resilient biodefense posture”. Indeed, the future of biodetection requires us to begin developing novel analytics that can identify anomalies and/or characteristics that indicate a potential threat, whether known or unknown, without looking for a specific signature that has been identified previously. To assess potential threats more rapidly, it is critical to develop agnostic artificial intelligence (AI)/machine learning (ML) systems that can be employed for real-time assessment of the nature and source of a perturbation. Such systems should be multi scale and multi-dimensional, integrating sensor data from a range of biological, chemical, and physical application spaces. Emerging deep learning (DL) models demonstrate exceptional promise for identification of discriminatory features within multi-dimensional datasets. DL models have the capacity to recognize and encode highly complex patterns in a wide range of input data modalities, including images, text, and biological/chemical/physical spectra. As such, they can execute a wide range of assessments and determinations that have traditionally required a human operator. The promise of advances in DL is apparent in the realm of human health and medicine. DL models have been validated for evaluating a variety of clinical threats to human health in a range of contexts, including infection and cancer, and they demonstrated improved performance in predicting stroke relative to human neurologists in some categories of data. Continuously evolving advances in AI/ML are expected to support more efficient evaluation of raw sequence, spectroscopy, and spectrometry data. For instance, recent advances and deployment of large language models (LLM) such as Generative Pre training Transformer (GPT) and Bidirectional Encoder Representations from Transformers (BERT) have already motivated application of these models for biological function prediction. As frameworks such as LLMs become larger and more complex in their representations, their capacity to serve as pre-trained models that can be fine-tuned for biological/biodetection purposes will similarly be amplified. While existing and emerging AI/ML have found broad applicability and use cases in the clinical sciences, development for environmental evaluation and biodetection has been limited. Functionalizing such capabilities for this purpose requires an understanding of the existing technical landscape and how the respective tools and algorithms are currently being employed. This landscape awareness then allows an assessment of the current practical capabilities of existing models and the anticipated requirements and development efforts that will be needed to adapt available algorithms for biodetection applications relevant to DHS. Leveraging expertise in biodetection, ML, and operational biodetection, the effort described in this report is comprised of a systematic landscape assessment (Subtask 2.1), comparative evaluation (Subtask 2.2), and formulation of a value proposition (Subtask 2.3) for the prospect of ML-enabled, agnostic biodetection from raw, or minimally-processed, datasets.

59 BASIC BIOLOGICAL SCIENCES↗

Spread Spectrum Time Domain Reflectometry (SSTDR) and Frequency Domain Reflectometry (FDR) for Detection of Cable Anomalies Using Machine Learning

Cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, continued use of these cables must shift to a performance-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. A variety of cable tests are available and are commonly applied during outages when the cables can be taken out of service. Frequency domain reflectometry (FDR) is one of these test methods that is being more broadly accepted and used because it not only detects anomalies along the cable with a low-voltage signal that does not stress the cable insulation, but the technique also locates the anomalies. This supports follow-up local inspection and local repair or partial replacement of a damaged cable segment. Currently, FDR testing is only applied to cables that are taken out of service since the test instrument would be damaged by operational voltages. A related technology that has found some acceptance in the aircraft and rail industry is spread spectrum time domain reflectometry (SSTDR). This technology has been implemented with a custom commercial instrument by LiveWire Innovation Inc. that is designed to operate on live cables up to 1000 volts. One of the main conclusions of a previous effort was that cable reflectometry plots can be difficult for humans to analyze due to baseline noise, low or noisy anomaly response peaks, or large responses from cable ends. Detection of cable anomalies for many of these frequencies and test conditions was challenging for manual analysis. This presented an ideal opportunity for ML analysis to distinguish undamaged cable indications from anomalous cable indications. This research discusses application of machine learning (ML) to reflectometry cable test methods. The goal was to assess feasibility to distinguish undamaged cable reflectometry responses from damaged or anomalous cable reflectometry responses. The assessment considered the 3 instruments, multiple frequency bandwidths from each instrument, multiple cable anomalies and test conditions, and both supervised and unsupervised ML approaches. Although approaches and analysis methods were not identical or directly comparable, both outputs were encouraging. The unsupervised prediction weighted accuracy was assessed by instrument and by frequency. It performed better at high frequencies with the highest prediction accuracy of 0.84 for the higher frequency FDR, 0.79 for the 48-MHz LiveWire SSTDR, and 0.77 for 300-MHz PNNL SSTDR. The initial weighted accuracy average across all frequencies for using supervised ML was 0.56 to 0.68. The supervised analysis was repeated with noisier training data removed resulting in weighted accuracies of 0.69 to 0.87. These weighted accuracies are not directly comparable due to differences in the supervised and unsupervised analysis details but do indicate an encouraging trend. Even with limited and unbalanced data, strong prediction accuracies seem encouraging for further work including more data under a wider range of conditions.

42 ENGINEERING↗

Phase-curve Pollution of Exoplanet Transmission Spectra

The occurrence of a planet transiting in front of its host star offers the opportunity to observe the planet’s atmosphere filtering starlight. The fraction of occulted stellar flux is roughly proportional to the optically thick area of the planet, the extent of which depends on the opacity of the planet’s gaseous envelope at the observed wavelengths. Chemical species, haze, and clouds are now routinely detected in exoplanet atmospheres through rather small features in transmission spectra, i.e., collections of planet-to-star area ratios across multiple spectral bins and/or photometric bands. Technological advances have led to a shrinking of the error bars down to a few tens of parts per million (ppm) per spectral point for the brightest targets. The upcoming James Webb Space Telescope (JWST) is anticipated to deliver transmission spectra with precision down to 10 ppm. The increasing precision of measurements requires a reassessment of the approximations hitherto adopted in astrophysical models, including transit light-curve models. Recently, it has been shown that neglecting the planet’s thermal emission can introduce significant biases in the transit depth measured with the JWST/Mid-InfraRed Instrument, integrated between 5 and 12 μm. In this paper, we take a step forward by analyzing the effects of the approximation on transmission spectra over the 0.6–12 μm wavelength range covered by various JWST instruments. We present open-source software to predict the spectral bias, showing that, if not corrected, it may affect the inferred molecular abundances and thermal structure of some exoplanet atmospheres.

47 OTHER INSTRUMENTATION↗

Potential bioavailability of representative pyrogenic organic matter compounds in comparison to natural dissolved organic matter pools

Abstract. Pyrogenic organic matter (PyOM) from wildfires impacts river corridors globally and is widely regarded as resistant to biological degradation. Though recent work suggests PyOM may be more bioavailable than historically perceived, estimating bioavailability across its chemical spectrum remains elusive. To address this knowledge gap, we assessed potential bioavailability of representative PyOM compounds relative to ubiquitous dissolved organic matter (DOM) with a substrate-explicit model. The range of potential bioavailability of PyOM was greater than natural DOM; however, the predicted thermodynamics, metabolic rates, and carbon use efficiencies (CUEs) overlapped significantly between all OM pools. Compound type (e.g., natural versus PyOM) had approximately 6-fold less impact on predicted respiration rates than simulated carbon and oxygen limitations. Within PyOM, the metabolism of specific chemistries differed strongly between unlimited and oxygen-limited conditions – degradations of anhydrosugars, phenols, and polycyclic aromatic hydrocarbons (PAHs) were more favorable under oxygen limitation than other molecules. Notably, amino sugar-like, protein-like, and lignin-like PyOM had lower carbon use efficiencies relative to natural DOM of the same classes, indicating potential impacts in process-based model representations. Overall, our work illustrates how similar PyOM bioavailability may be to that of natural DOM in the river corridor, furthering our understanding of how PyOM may influence riverine biogeochemical cycling.

59 BASIC BIOLOGICAL SCIENCES↗

Measurement of the Deeply Virtual Compton Scattering Cross Section from the Proton at 10.6 GeV using the CLAS12 Detector

Deeply Virtual Compton Scattering (DVCS) is an exclusive process that produces a real photon when a lepton scatters from a quark inside a nucleon or a nucleus. Measurement of the DVCS cross section enables the study of the Generalized Parton Distributions (GPD), which plays a central role in understanding the QCD dynamics inside a hadron. Thus, the quark and gluon origin of the nucleon spin and mass can be probed and three-dimensional images of the target nucleon or nucleus can be realized. This thesis presents a cross section analysis of DVCS from the proton in the presence of its background, Bethe-Heitler (BH) process. The CEBAF Large Acceptance Spectrometer for operation at 12 GeV beam energy (CLAS12) collaboration has taken electron-proton scattering data in fall 2018 using a liquid hydrogen target and the 10.6 GeV polarized electron beam from the Continuous Electron Beam Accelerator Facility (CEBAF). The CLAS12 detector is a nearly hermetic fixed-target detector, located in Hall B, Jefferson Lab at Newport News, Virginia. The experimentally determined BH-DVCS cross section is in good agreement with a phenomenological-model based theoretical prediction. The kinematic dependence of the cross section is reported over a wide range. The short-term plan to utilize the results presented here for a thorough tomography study and the long-term plan for GPD studies at future facilities such as the Electron-Ion Collider (EIC) are discussed.

Lee, Sangbaek↗

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

Here, we analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call “complexity reduction”. In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature T g of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2–6 nm at T g . Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

cancer↗

Two-Nucleon Short-Range Correlations in Light Nuclei

Understanding the nucleon-nucleon interaction is a fundamental task in nuclear physics, as NN interaction models are a crucial input to modern nuclear structure calculations. While great progress has been made toward understanding this interaction, the available state-of-the-art models predict significantly different behaviors at short distances and high momenta (scale-and-scheme dependence), where two-nucleon Short-Range Correlations (SRCs) dominate the nuclear wave function. Thus, SRCs are a unique tool to constrain the NN interaction and vice versa. SRCs are naturallyoccurring high-local-density NN pairs that, as a result of their short-distance (approx. 0.1 fm) repulsive interaction, fly apart with high momenta, hence populating momentum states above the Fermi level. The study of SRCs also has significant implications for other fields, such as the astrophysics of neutron stars and the behavior of cold atomic gasses. This thesis describes experimental and phenomenological studies of the short-distance / high-momentum structure of the NN interaction through the study of SRCs and vice versa. Experimentally, I report the first measurement of the 3He and 3H(e,e'p) reactions in Hall A of the Thomas Jefferson National Accelerator Facility in kinematics in which the measured cross sections should be sensitive to the underlying nucleon momentum distributions in the range 40 to 500 MeV/c. The resulting cross-section ratios and absolute cross sections were compared to momentum-distribution ratios and precise cross-section calculations respectively. Phenomenologically, I report the generalization of the Contact Formalism (GCF) to nuclear systems, which exploits scale separation and universality to describe nucleons at short distances and high momenta.

Cruz Torres, Reynier↗

Impedance-Based Prediction of SSR-Generated Harmonics in Doubly-Fed Induction Generators

Doubly-fed induction generators (DFIGs) are prone to subsynchronous resonance (SSR) with series-compensated transmission lines and experiencing SSR-generated harmonics of wide-ranging magnitudes. This paper presents a large-signal impedance-based approach for the prediction of the magnitude of SSR-generated harmonics in DFIGs. Note that the large-signal impedance of a network shows variation in its impedance response with the magnitude of perturbation at its terminals. It is discovered in this paper that the impedance of a series-compen- sated transmission line starts changing with the magnitude of SSR because of the metal-oxide varistor-based protection of series capacitors in the line. The changing line impedance halts the growth of the SSR beyond a certain point, where the DFIG forms a limit-cycle mode and emits harmonics at the SSR frequency. The discovered role of the protection of series capacitors in SSR not only enables the prediction of SSR-generated harmonics but also can be used for the timely detection of an SSR event to avoid system disruption and sudden loss of bulk generation. The large- signal impedance-based prediction of SSR-generated harmonics in DFIG-based wind and hydro generators is demonstrated using numerical simulations of a 2.5-MVA, 690-V DFIG supplying to a series-compensated transmission line.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Additionally, Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

79 ASTRONOMY AND ASTROPHYSICS↗