Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “different dimensions”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cell dynamics simulations of coupled charge and magnetic phase transformation in correlated oxides

In this work, we present a comprehensive numerical study on the kinetics of phase transition that is characterized by two nonconserved scalar order parameters coupled by a special linear-quadratic interaction. This particular Ginzburg-Landau theory has been proposed to describe the coupled charge and magnetic transition in nickelates and the collinear stripe phase in cuprates. The inhomogeneous state of such systems at low temperatures consists of magnetic domains separated by quasimetallic domain walls where the charge order is reduced. By performing large-scale cell dynamics simulations, we find a two-stage phase-ordering process in which a short period of independent evolution of the two order parameters is followed by a correlated coarsening process. The long-time growth and coarsening of magnetic domains is shown to follow the Allen-Cahn power law. We further show that the nucleation-and-growth dynamics during phase transformation to the ordered states is well described by the Kolmogorov-Johnson-Mehl-Avrami theory in two dimensions. On the other hand, the presence of quasimetallic magnetic domain walls in the ordered states gives rise to a very different kinetics for transformation to the high-temperature paramagnetic phase. In this scenario, the phase transformation is initiated by the decay of magnetic domain walls into two insulator-metal boundaries, which subsequently move away from each other. Implications of our findings to recent nano-imaging experiments on nickelates are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative datasets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student-t process regression. We apply Student-t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via coloring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics, and classic machine learning clustering algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Stochastic parametric skeletal dosimetry model for humans: Anatomical-morphological basis and parameter evaluation

Radiation exposure of the hematopoietic system that results in a radiation dose to bone marrow of more than 100 mGy leads to an increase in the risk of leukemia in humans. Excess relative risk of leukemia was observed in cohorts whose members lived in the territories of the Southern Urals that were radioactively contaminated in the 1950s. As part of the dosimetric support of epidemiological studies of these cohorts, an original methodology for stochastic bone dosimetric modeling was developed, termed the Stochastic Parametric Skeletal Dosimetry (SPSD) model. The purpose of this work was to present the anatomical and morphological bases of the SPSD model, which includes an assessment of the parameters of the microstructure of the trabecular bone in the hematopoietic areas of the human skeleton, as well as a description of the macrostructural division (segmentation) of hematopoietic areas into simple bone segments. As a result, an anatomical-morphological basis of the SPSD model was created based on published data. Data collection work included the analysis of original articles, atlases, manuals, monographs and the formation of primary data files. Data on the duration of hematopoiesis in various parts of the skeleton; and data on age-related changes in the microstructure and linear dimensions of human bones and their segments were analyzed. The paper describes the full set of parameters to be used for the dosimetric model for newborns, children aged 1, 5 and 10 years, as well as for adolescents aged 15 years and adults; for the latter, sex differences in bone size were considered. In total, the SPSD model includes 289 unique basic (bone) phantom segments, each of which is described by 7 or more parameters describing the microstructure, thickness of the cortical layer and linear dimensions. Population variability was estimated for each parameter. The approach to SPSD modeling, i.e., the use of simple geometric shapes, was successfully verified using independent datasets on bone masses and volumes.

Science & Technology - Other Topics↗
DOE OSTI2025

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗
DOE OSTI2022

Automatic Segmentation of Building Envelope Point Cloud Data Using Machine Learning

About 50% of buildings in the US were constructed before energy codes were introduced. Modular overclad panel retrofits, in which a new envelope is constructed over the existing building, are a promising solution given that it minimizes occupant disruption and shortens construction time at the jobsite. Current state-of-the-art retrofit panel layout and dimensioning consists of three steps: 1) 3D point cloud data generation of the building envelope using commonly available surveying equipment, 2) manual segmentation of 3D point cloud data by a trained professional to identify and dimension window openings, door openings, and other architectural features, and 3) modular panel layout optimization and dimensioning by an architect or engineer. Among these steps, the second one remains the most difficult and costly because it is very labor-intensive. We propose a methodology to automatically label 3D point cloud data to reduce the time and expense spent in manual segmentation. Machine learning methods were employed to classify the point cloud data into distinct groups, each of which corresponds to different features of the building envelope. After classification, a segmentation algorithm was developed to perform boundary detection and separate the components of the façade. Finally, the algorithm returns the relative positions and dimensions of the features in the building envelope. The measurements obtained with the proposed automated method were compared against the actual dimensions to determine the overall algorithm accuracy. The proposed algorithm can then be used to reduce manual efforts for 3D point cloud labeling before modular panel layout optimization is performed.

Maldonado Puente, Bryan↗
DOE OSTI2024

Investigating Mechanisms Driving Differences in the Characteristics of Precipitation in the E3SM Multiscale Modeling Framework With 2D Versus 3D Cloud Resolving Model Configurations

Abstract In this study, we compare the Energy Exascale Earth Systems Model (E3SM) multiscale modeling framework (MMF) with the cloud resolving model (CRM) configured in two (2dMMF) and three (3dMMF) dimensions. We explore how CRM dimensionality impacts the representation of mean and extreme precipitation characteristics. Our results show that tropical mean precipitation patterns are better represented in 3dMMF compared to observations (Integrated Multi‐satellitE Retrivals for GPM and Global Precipitation Climatology Project One Degree Daily products), while 2dMMF better captures extreme precipitation intensity, with systematic land‐ocean differences in precipitation and cloud‐associated variables. These differences are attributed to the co‐occurrence of CRM throttling (i.e., suppressed convection in due to smaller numbers of CRM columns and domain size) and dilution (i.e., 3‐D cloud circulations with increased entrainment and lower precipitation efficiency) effects. Overall, throttling results in more low‐level humidity in 2dMMF and dilution contributes to more high clouds with less precipitation efficiency in 3dMMF. Since throttling occurs more strongly over the ocean than land, the 3dMMF tends to have less cloud liquid and precipitation over the ocean and more cloud ice and precipitation over land. These results may serve as a guide for choosing the CRM structure to reduce precipitation and cloud‐related biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Microstructure, Texture, and Mechanical Properties of Continuously Extruded and Rolled AZ31 Magnesium Alloy Sheets

It is difficult to fabricate magnesium alloy sheets at room temperature due to their hexagonal close-packed structure. In this paper, we proposed a new process that can be used to fabricate AZ31 magnesium alloy sheets with fine-grained microstructure and good mechanical properties. AZ31 magnesium alloy sheets with dimensions of 160 mm × 8 mm were obtained successfully by a continuous extrusion technique. The extruded sheets were then rolled by two different rolling routes, Route 1: the extruded sheet was heated at 350 °C and then rolled without heating of rolls and Route 2: the sheet was rolled without any previous heating by heating rolls at 280 °C. The microstructure and texture evolution as well as mechanical properties of the continuously extruded and rolled sheets were examined experimentally. The magnesium alloy sheets with relatively fine grains of 4.5-8.3 µm were obtained by continuous extrusion followed by rolling processes. The tensile strength of continuously extruded AZ31 Mg alloy sheets can be effectively increased by rolling using heating rolls (Route 2), while the ductility was enhanced by rolling heated sheets (Route 1). The mixed microstructure, consisting of refined DRX grains and deformed grains, contributes to enhancement in tensile strength of the sheet samples rolled by Route 2 due to grain-boundary strengthening and work hardening. Lower temperature and larger rolling reduction are two main factors for improving the tensile strength of AZ31magnesium alloy sheets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Synthesis and morphological characterization of linear and miktoarm star poly(solketal methacrylate)- block -polystyrene copolymers

The synthesis, molecular characterization, and morphological evaluation of AB n, (n = 2,3) miktoarm star block copolymers consisting of poly(glycerol monomethacrylate) (PGMA) and polystyrene (PS) with varying molecular weights and compositions is described. The system is known to demonstrate a remarkably high Flory-Huggins interaction parameter. The corresponding linear diblock copolymer analogues, were synthesized as well, and comparisons with regard to feature dimensions and morphologies are provided. Well-ordered nanostructures of various morphologies were formed with domain spacing as low as 7.2 nm. Different morphologies were attained by some of the topological isomers indicating that in miktoarm star block copolymers the phase boundaries were strongly shifted. Additionally, a triblock ABA analogue was studied to investigate the effect of triblock copolymer conformations. Noteworthy is that for copolymers with different macromolecular architecture leading to similar morphology; different domain spacings were obtained. The synthesis of all the samples was carried out by high-vacuum anionic polymerization techniques. Molecular characterization with Size Exclusion Chromatography (SEC) and Proton Nuclear Magnetic Resonance Spectroscopy ( 1 H NMR) confirmed well-defined copolymers obtained. The morphological characterization was accomplished by Small-Angle X-ray Scattering (SAXS). In conclusion, the observations from this study highlight the potential of incorporating macromolecular architecture in the self-assembly of strongly immiscible block copolymers to attain ultra-small nanofeatures with desired morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrated and differential fiducial cross-section measurements for the vector boson fusion production of the Higgs boson in the $H → WW^* → evμv$ decay channel at 13 TeV with the ATLAS detector

The vector-boson production cross-section for the Higgs boson decay in the H→WW*→eνμν channel is measured as a function of kinematic observables sensitive to the Higgs boson production and decay properties as well as integrated in a fiducial phase space. The analysis is performed using the proton-proton collision data collected by the ATLAS detector in Run 2 of the LHC at $\sqrt{s}$=13 TeV center-of-mass energy corresponding to an integrated luminosity of 139 fb -1 . The different flavor final state is studied by selecting an electron and a muon originating from a pair of W bosons and compatible with the Higgs boson decay. The data are corrected for the effects of detector inefficiency and resolution, and the measurements are compared with different state-of-the-art theoretical predictions. The differential cross-sections are used to constrain anomalous interactions described by dimension-six operators in an Effective Field Theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Scale and conformal invariance on (A)dS spacetimes

We examine the question of scale versus conformal invariance on maximally symmetric curved backgrounds and study general two-derivative conformally invariant free theories of vectors and tensors. For spacetime dimension D > 4 , these conformal theories can be diagonalized into standard massive fields in which unbroken conformal symmetry nontrivially mixes components of different spins. In D = 4 , the tensor case becomes a conformal theory mixing a partially massless spin-2 field with a massless spin-1 field. For massless linearized gravity in D = 4 , we confirm through direct calculation that correlation functions of gauge-invariant operators take the conformally invariant form, despite the absence of standard conformal symmetry at the level of the action. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Rapid Initial-State Preparation for the Quantum Simulation of Strongly Correlated Molecules

Studies on quantum algorithms for ground-state energy estimation often assume perfect ground-state preparation; however, in reality the initial state will have imperfect overlap with the true ground state. Here, we address that problem in two ways: by faster preparation of matrix-product-state (MPS) approximations and by more efficient filtering of the prepared state to find the ground-state energy. We show how to achieve unitary synthesis with a Toffoli complexity about 7 × lower than that in prior work and use that to derive a more efficient MPS-preparation method. For filtering, we present two different approaches: sampling and binary search. For both, we use the theory of window functions to avoid large phase errors and minimize the complexity. We find that the binary-search approach provides better scaling with the overlap at the cost of a larger constant factor, such that it will be preferred for overlaps less than about 0.003 . Finally, we estimate the total resources to perform ground-state energy estimation of Fe - S cluster systems, including the Fe Mo cofactor by estimating the overlap of different MPS initial states with potential ground states of the Fe Mo cofactor using an extrapolation procedure. With a modest MPS bond dimension of 4000 , our procedure produces an estimate of approximately 0.9 overlap squared with a candidate ground state of the Fe Mo cofactor, producing a total resource estimate of 7.3 × 10 10 Toffoli gates; neglecting the search over candidates and assuming the accuracy of the extrapolation, this validates prior estimates that have used perfect ground-state overlap. This presents an example of a practical path to prepare states of high overlap in a challenging-to-compute chemical system. Published by the American Physical Society 2025

Berry, Dominic W. (ORCID:0000000334461449)↗
DOE OSTI2020

Modeling Oxygen Transport in High Surface Area Carbon Supports for Polymer-Electrolyte Fuel Cells

Here, an analytical model for resistance to oxygen transport in air electrodes containing carbon black supports with high surface area was developed by combining a Thiele modulus—effectiveness factor approach at the agglomerate scale (~150 nm) with nanoscale diffusional resistance in carbon micropores/pits (~5 nm). This paper extends an earlier model for transport resistance to platinum nanoparticles on low surface area carbon. Differences in transport resistances between catalyst layers with high and low surface area carbon blacks predicted by the model with reasonable geometric dimensions and physical properties are consistent with experimental observations.

08 HYDROGEN↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2022

Renewable Energy and Efficiency Technologies in Scenarios of U.S. Decarbonization in Two Types of Models: Comparison of GCAM Modeling and Sector-Specific Modeling

Energy system projections from quantitative models inform actions ranging from short-term and local decisions, such as those about technology and infrastructure deployment, and global and long-term negotiations and targets. Computational limits require model designers to balance coverage and resolution (i.e., breadth versus depth). Some models, such as the Global Change Analysis Model (GCAM), represent all energy sources and uses with less resolution than models that focus on a single sector's energy use. GCAM balances global supply and demand of all energy carriers projecting prices using internal calculations for energy sources and costs of greenhouse gas mitigation while capturing interlinkages between the energy system, water, agriculture and land use, the economy, and the climate. This globally comprehensive model was used to frame the Long-Term Strategy of the United States: Pathways to Net-Zero Greenhouse Gas Emissions by 2050, which the White House released in 2021 and has been used to inform national and global economy-wide climate change mitigation discussions and strategy development for decades. Unlike GCAM, sectoral models focus on a portion of the energy sector and with greater detail and resolution. The Regional Energy Deployment System (ReEDS) electricity-sector model, for example, projects electricity system capacity expansion and operation with high-fidelity representation of emerging technologies for deep decarbonization, such as variable renewable energy and energy storage, and integration of these technologies into the electric grid. The Transportation Energy and Mobility Pathway Options (TEMPO) transportation-sector model enables analysis of household choices, with a focus on adoption, charging, and use of electric vehicles. The Scout buildings-sector model supports detailed consideration of the policies and markets that can accelerate the adoption of electrification and energy conservation measures in buildings. Such sector-specific models are instrumental in informing technology research, sectoral planning strategies, and sector-specific aspects of greenhouse gas mitigation strategies in the United States. The integrated multisector and sector-specific modeling approaches represented by GCAM and these sectoral models are complementary. The integrated multisector approach calculates energy pricing and resource allocation within the model, which is important for consistency when future conditions substantially diverge from current conditions in transformative scenarios. The sector-specific approach facilitates representation of granular details across spatial, temporal, technological, and market dimensions that enable exploration of particular interactions and trade-offs. This report presents the results of recent work to explore the differences and trade-offs between these approaches by comparing GCAM with the sector-specific ReEDS, TEMPO, and Scout models. The report compares both model structures and results, and it addresses their potential relevance and applications.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Renewable Energy and Efficiency Technologies in Scenarios of U.S. Decarbonization in Two Types of Models: Comparison of GCAM Modeling and Sector-Specific Modeling

Energy system projections from analytic models inform actions ranging from short-term and local decisions, such as technology and infrastructure deployment, to global and long-term negotiations and targets. Computational limits require the designers of these models to trade off between coverage and resolution. Some models, such as the Global Change Analysis Model (GCAM), represent all energy sources and uses but at a relatively coarse level of resolution. GCAM balances global supply and demand of all energy carriers by endogenously projecting prices for energy sources and costs of greenhouse gas mitigation while capturing interlinkages between the energy system, water, agriculture and land use, the economy, and the climate. This global model was used to frame the Long-Term Strategy released by the White House in 2021 and has been used to inform national and global economy-wide decarbonization discussions and strategy development for decades. Other models instead focus on a portion of the energy sector with greater detail and resolution. The Regional Energy Deployment System (ReEDS) electricity-sector model, for example, projects capacity expansion with an emphasis on integration of variable renewable energy into the grid of the future. The Transportation Energy and Mobility Pathway Options (TEMPO) transportation-sector model enables analysis of household choices in adoption, charging, and use of electric vehicles. The Scout buildings-sector model supports detailed consideration of the policies and markets that can accelerate the adoption of energy conservation measures in buildings. Such sector-specific models are instrumental in informing technology research, sectoral planning strategies, and sector-specific aspects of greenhouse gas (GHG) mitigation strategies in the United States. These global and sector-specific modeling approaches can complement each other. The global approach ensures consistent, endogenous energy pricing and resource allocation, which can substantially diverge from current conditions in transformative scenarios, while the sector-specific approach facilitates representation of granular details across spatial, temporal, technological, and market dimensions that enable exploration of particular interactions and trade-offs. This report presents the results of recent work to explore the differences and tradeoffs between these approaches by comparing GCAM with the sector-specific ReEDS, TEMPO, and Scout models. The report compares both model structures and results, and discusses their potential relevance and applications.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

The GGCMI Phase 2 emulators: global gridded crop model responses to changes in CO 2 , temperature, water, and nitrogen (version 1.0)

Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: that growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Extending and Smoothing Two-dimension Equation of State Simulation Data

The Multiphase Equation of State (MEOS) project works to produce high quality equation of state tables which are used in computationally intensive simulations of materials in different conditions. An equation of state table describes the state of matter under certain physical conditions such as pressure, volume, temperature, or internal energy. MEOS uses many different models and combinations thereof to produce accurate tables consisting of continuous smooth data, derivatives, and higher-order derivatives. Accurate and smooth data are important factors in producing precise simulations. This report focuses on improving the electron tables produced from Purgatorio and Thomas-Fermi data. Purgatorio data is extremely accurate but often jagged and discontinuous in the lower temperature and density region. Thomas-Fermi is a model that provides smooth data throughout. Together, this produces an EOS data table that is both accurate and smooth. This report describes the implementation of a new feature that allows users to define more detailed regions in the Purgatorio table to be replaced with Thomas-Fermi. As a result, MEOS can generate a table that is both more smooth and still accurate.

97 MATHEMATICS AND COMPUTING↗