Search NASA⌕ Search

SEARCH · Search NASA

Results for “model structural error”

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 235 records · Page 13

Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-β (Final Report)

This is the final report for the DOE award DE-SC0019239, titled “Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-b” for the period July 1, 2018 – July 2, 2023. This award supported research into computational models for plasma equilibrium and 3D conducting structures in the Lithium Tokamak eXperiment-b (LTX-b) using the PSI-Tri/TokaMaker and PSI-Tet codes. Additional detail on progress in support of these tasks is provided below.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies

KIC 08626021 is a pulsating DB white dwarf (WD) of considerable recent interest, and the first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based WD models choose a luminosity profile where the luminosity is proportional to the enclosed mass, ${L}_{r}\,\propto \,{M}_{r}$, independent of the effective temperature T eff . Evolutionary models of young WDs with T eff ≳ 25,000 K suggest that neutrino emission gives rise to luminosity profiles with L r ${/}\!\!\!\!\!\!{\propto }$ M r . We explore this contrast by comparing the oscillation frequencies between two nearly identical WD models: one with an enforced ${L}_{r}\propto {M}_{r}$ luminosity profile, and the other with a luminosity profile determined by the star's previous evolution history. We find that the low-order g-mode frequencies differ by up to ≃70 μHz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an ${L}_{r}\propto {M}_{r}$ profile may have uncorrected, effectively random errors at the level of tens of μHz. A mean frequency difference of 30 μHz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of ≃12% in WD mass, ≃9% in the radius, and ≃3% in the central oxygen mass fraction.

79 ASTRONOMY AND ASTROPHYSICS↗

Data and code from: Multivariate bayesian regression model for predicting disposed ash composition at U.S. coal fired power stations

This dataset contains the code and data files needed for implementation of a Multivariate Bayesian Regression model, described in Jin et al. (2025), for the historical prediction of the chemical composition of disposed coal ash at U.S. coal fired power plants as a function of annualized coal purchase data. The integrated coal supply data file (CoalSupplyDataset.csv) represents a compilation of monthly fuel purchase records for the period 1973-2022 at major U.S. power stations. These records were obtained from the U.S. Energy Information Administration. The CSV file also contains, for each coal purchase record, the coal region of the mine as defined by the U.S. Geological Survey. Data entry errors and data gaps in the EIA records were corrected as described in Jin et al. This CSV file represents the integrated coal supply data after corrections were made. The model structure and fitting parameters are encoded in pickle file format (Bayesian.pkl). The model was developed with the coal supply data and coal ash composition data, apportioned according to the Stratified Shuffle Split for training and testing subsets. The model was built using Python and the PyMC library. Reference Publication: Jin, Z.; Huang, J.; Hower, J.C.; Hsu-Kim, H.(2025). Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites. Environmental Science & Technology.

Coal ash composition↗

Two-level structural sparsity regularization for identifying lattices and defects in noisy images

Here, this paper presents a regularized regression model with a two-level structural sparsity penalty applied to locate individual atoms in a noisy scanning transmission electron microscopy image (STEM). In crystals, the locations of atoms is symmetric, condensed into a few lattice groups. Therefore, by identifying the underlying lattice in a given image, individual atoms can be accurately located. We propose to formulate the identification of the lattice groups as a sparse group selection problem. Furthermore, real atomic scale images contain defects and vacancies, so atomic identification based solely on a lattice group may result in false positives and false negatives. To minimize error, model includes an individual sparsity regularization in addition to the group sparsity for a within-group selection, which results in a regression model with a two-level sparsity regularization. We propose a modification of the group orthogonal matching pursuit (gOMP) algorithm with a thresholding step to solve the atom finding problem. The convergence and statistical analyses of the proposed algorithm are presented. The proposed algorithm is also evaluated through numerical experiments with simulated images. The applicability of the algorithm on determination of atom structures and identification of imaging distortions and atomic defects was demonstrated using three real STEM images. In conclusion, we believe this is an important step toward automatic phase identification and assignment with the advent of genomic databases for materials.

36 MATERIALS SCIENCE↗

ZeoNet: 3D convolutional neural networks for predicting adsorption in nanoporous zeolites

Zeolites are one of the most widely used materials in the chemical industry due to their nanometer-sized pores that can adsorb and react upon molecules selectively. With hundreds of known framework topologies and hundreds of thousands of computationally predicted structures, the ability to rapidly predict zeolite performance allows researchers to prioritize their efforts on the most promising structures for a given application. Although the accuracy of forcefield-based atomistic simulations has advanced significantly in the past two decades, these simulations can be computationally expensive, especially for long-chain, complex molecules. Here, we present ZeoNet, a representation learning framework using convolutional neural networks (ConvNets) and 3D volumetric representations for predicting adsorption in zeolites. ZeoNet was trained on the task of predicting Henry's constants for adsorption, k H , of n-octadecane in more than 330 000 known and predicted zeolite materials. Employing a 3D grid based on the distances to solvent-accessible surfaces, a volumetric representation that can be generated efficiently, the best-performing ZeoNet achieved a correlation coefficient r 2 = 0.977 and a mean-squared error MSE = 3.8 in ln k H , which corresponds to an error of 9.3 kJ mol -1 in adsorption free energy. In comparison, a model based on hand-designed geometric features has values of r 2 = 0.783 and MSE = 35.7. ZeoNet is also relatively efficient and can process ≈8 structures per second on an Nvidia RTX 2080TI GPU, orders of magnitude faster than forcefield-based simulations. A systematic analysis was conducted to investigate how the choice of ConvNet architectures, the linear dimension (L) and spatial resolution (Δd) of the distance grids, batch size, optimizer, and learning rate impact the model performance. We found that ConvNets based on the ResNet architecture offer the best tradeoff between expressiveness and efficiency. The performance for all models reaches a plateau at L = 30–45 Å and depends less sensitively on grid resolution, with a small benefit around Δd = 0.30–0.45 Å. Finally, saliency maps were visualized to identify which regions of the materials contributed the most to model predictions. It was found, interestingly, that the predictions are driven primarily by the accessible pore volume rather than the region occupied by the framework atoms.

36 MATERIALS SCIENCE↗

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗

Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.

The remoteness and extreme conditions of the Arctic make it a very difficult environment to investigate. In these polar regions covered by sea ice, the wind is relatively strong due to the absence of obstructions and redistributes a large part of the deposited snow mass, which complicates estimates for precipitation hardly distinguishable from blowing or drifting snow. Moreover, the snow mass balance in the sea ice system is still poorly understood, notably due to the complex structure of its surface. Quantitatively assessing the snow distribution on sea ice and its connection to the sea ice surface features is an important step to remove the snow mass balance uncertainties (i.e., snow transport contribution) in the Arctic environment. In this work we introduce snowBedFoam 1.0., a physics-based snow transport model implemented in the open-source fluid dynamics software OpenFOAM. We combine the numerical simulations with terrestrial laser scan observations of surface dynamics to simulate snow deposition in a MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate) sea ice domain with a complicated structure typical for pressure ridges. The results demonstrate that a large fraction of snow accumulates in their vicinity, which compares favorably against scanner measurements. However, the approximations imposed by the numerical framework, together with potential measurement errors (precipitation), give rise to quantitative inaccuracies, which should be addressed in future work. The modeling of snow distribution on sea ice should help to better constrain precipitation estimates and more generally assess and predict snow and ice dynamics in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

A machine learning approach to quantify degradation of nuclear fuels and the effects of fission products

Nuclear fuel performance is critically dependent on understanding the evolution of fuel properties under operational conditions, a complex challenge driven by chemical changes and substantial radiation damage during fission. Traditionally, property evolution has been determined via empirical data collected following irradiation. However, these empirical correlations are limited in their applicability beyond the specific conditions in which they were obtained. This study explores a novel approach to address this challenge by applying materials informatics to develop a machine learning random forest (ML-RF) model that captures the effects of fission products on fuel compounds. The model predicts formation enthalpy (ΔH f ) by leveraging extensive quantum materials property data and correlating it with material descriptors such as composition, atomic and site features, and crystal lattice properties. This ML-RF model enables rapid interpolation across the compositional and structural spaces covered by the training data, thus supporting high-throughput screening and energetic ranking of candidate phases. The model demonstrates the ability to predict ΔH f with a mean absolute error (MAE) of approximately 0.1 to 0.2 eV/atom across a wide range of compounds, including key nuclear fuel systems (U-O, U-N, U-C, U-Si, and U-Mo). For example, it was used to assess shifts in stoichiometry for UO 2 (O/M) and UN (N/M) fuels, revealing their distinct tendencies in chemical potential variation and enabling preliminary convex hull analyses. Furthermore, the model provides insights into how individual fission products affect fuel properties. Results indicate that larger fission products (e.g., Nd, Pu, Ce) have a more pronounced impact on UO 2 , while lighter ones (e.g., Zr) strongly influence UN. Here, the model developed in this work can be used to support the Accelerated Fuel Qualification approach by facilitating preliminary evaluations prior to extensive materials modeling and experimentation. To this end, the trained model has been made available to the fuel community to support ongoing fuel development efforts.

Accelerated fuel qualification↗

Data for: Climatic Imprint on Interfacially-Controlled Platinum-Palladium Resources

Data package for manuscript "Climatic Imprint on Interfacially-Controlled Platinum-Palladium Resources" by Emily G. Wright, Ivey Wang, Yihang Fang, Elaine D. Flynn, and Jeffrey G. Catalano. This dataset contains adsorption results from experiments designed to investigate the effect of chloride on Pd(II) adsorption to goethite and Pt(II) adsorption to hematite and goethite, including lab experiments, X-ray absorption fine structure spectroscopy, and models of retention within a laterite. See the associated manuscript for full methods information. The file "Wright2025_PtAds_data.csv" contains the target starting Pt concentration (uM), final aqueous Pt and associated error (in uM), calculated adsorbed Pt and associated error (in umol/m2), target and measured aqueous chloride (mM), target aqueous nitrate (mM), final pH, and mineral concentration/loading (g/L). Associated mineral-free controls (mineral loading = 0 g/L) are included; the aqueous Pd error was not calculated and chloride was not measured in every sample. These data appear in Figures 1, S3, S4, S5, S20, and S22 in the associated manuscript. The file "Wright2025_PdAds_data.csv" contains the target starting Pd concentration (uM), final aqueous Pd and associated error (in uM), calculated adsorbed Pd and associated error (in umol/m2), target and measured aqueous chloride (mM), and mineral concentration/loading (g/L). Associated mineral-free controls (mineral loading = 0 g/L) are included; the aqueous Pd error was not calculated and chloride was not measured in every sample. These data appear in Figures 1, S3, S4, S5, and S20 in the associated manuscript. The file "Wright2025_MineralBatches_data.csv" contains the mineral identity and BET specific surface area (m2/g) for every mineral batch synthesized and used in experiments. The annealing time used is listed for hydrothermally annealed goethite. These data appear in Table S2 in the associated manuscript. The file "Wright2025_XRD_data.csv" contains the XRD patterns for every mineral batch synthesized as the counts as a function of two theta (in degrees). See "Wright2025_MineralBatches_data.csv" for more details on specific mineral batches. These data appear in Figure S2 in the associated manuscript. The file "Wright 2025_ZetaPotential_data.csv" contains the measured zeta potentials for samples of goethite (batch G2) at pH 4 the presence of varying amounts of sodium chloride. These data appear in Table S3 in the associated manuscript. The file "Wright2025_XAFSSamples_data.csv" contains the specific mineral batch, measured final aqueous Pd or Pt (uM), measured final aqueous chloride (mM), and estimated adsorbed Pd or Pt (umol/m2) of all XAFS samples. These data appear in Tables S4, S7, S8, and S10 in the associated manuscript. The files "Wright2025_PdXAFS_data.csv" and "Wright2025_PtXAFS_data.csv" contain the normalized spectra of Pd and Pt, respectively, adsorbed to minerals at varying chloride concentrations. See "Wright2025_XAFSSamples_data.csv" for a guide to sample names. Note that "05" in a sample name is equivalent to "0.5". These data appear in Figures 2, S6, S7, S8, S12, S13, and S14 in the associated manuscript. The file "Wright2025_LateriteProfileProfileModelParameters_data.csv" include the ratio of hematite to hematite and goethite in two synthetic, modeled profiles, as well as the modeled surface areas of goethite and hematite as a function of relative depth within the modeled weathering zone. These data were used, in conjunction with equations presented in the paper, to calculate the theoretical concentrations of Pd and Pt (and the resulting Pt/Pd ratio) within the profiles. These data appear in Figure 3 in the associated manuscript. The file "Wright2025_Imagery_data.zip" is a zipped folder containing the TEM and STEM images appear in Figures S18 and S19. Individual files are labeled as either STEM (Fig. S18) or TEM (Fig. S19) with a letter representing the part of the multipart figure.

58 GEOSCIENCES↗

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transformer Leakage Inductance Design Methodology

The leakage inductance exhibited by a transformer depends on its winding geometry, which generally involves the selection of several key design parameters in addition to the winding structure and the interleaving configuration. With few resources explaining the effects of these design choices on the observed leakage inductance, numerous trial-and-error iterations become necessary to realize the desired leakage inductance. This paper explores more than a hundred winding geometries feasible in a 2-winding transformer comprising the same magnetic core, number of turns, and wire gauge, and finds the leakage inductance for each unique design using 2-D finite element method (FEM) simulations in association with the semi-analytical double- 2-D model. These leakage inductances are plotted and further analyzed to understand the effects of different design parameters on the effective leakage inductance. The results presented herein and the conclusions drawn from this research can serve as a valuable resource for future design practitioners from both industry and academia.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A discussion on various experimental methods of impact ionization coefficient measurement in GaN

Impact ionization coefficients play a critical role in semiconductors. In addition to silicon, silicon carbide and gallium nitride are important semiconductors that are being seen more as mainstream semiconductor technologies. As a reflection of the maturity of these semiconductors, predictive modeling has become essential to device and circuit designers, and impact ionization coefficients play a key role here. Recently, several studies have measured impact ionization coefficients. We dedicated the first part of our study to comparing three experimental methods to estimate impact ionization coefficients in GaN, which are all based on photomultiplication but feature characteristic differences. The first method inserts an InGaN hole-injection layer, the accuracy of which is challenged by the dominance of ionization in InGaN, leading to possible overestimation of the coefficients. The second method utilizes the Franz–Keldysh effect for hole injection but not for electrons, where the mixed injection of induced carriers would require a margin of error. The third method uses complementary p–n and n–p structures that have been at the basis of this estimation in Si and SiC and leans on the assumption of a constant electric field, and any deviation would require a margin of error. In the second part of our study, we evaluated the models using recent experimental data from diodes demonstrating avalanche breakdown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generalization error guaranteed auto-encoder-based nonlinear model reduction for operator learning

Many physical processes in science and engineering are naturally represented by operators between infinite-dimensional function spaces. The problem of operator learning, in this context, seeks to extract these physical processes from empirical data, which is challenging due to the infinite or high dimensionality of data. An integral component in addressing this challenge is model reduction, which reduces both the data dimensionality and problem size. In this paper, we utilize low-dimensional nonlinear structures in model reduction by investigating Auto-Encoder-based Neural Network (AENet). AENet first learns the latent variables of the input data and then learns the transformation from these latent variables to corresponding output data. Our numerical experiments validate the ability of AENet to accurately learn the solution operator of nonlinear partial differential equations. Furthermore, we establish a mathematical and statistical estimation theory that analyzes the generalization error of AENet. Finally, our theoretical framework shows that the sample complexity of training AENet is intricately tied to the intrinsic dimension of the modeled process, while also demonstrating the robustness of AENet to noise.

Auto-encoder↗

General Protocol for the Accurate Prediction of Molecular 13 C/ 1 H NMR Chemical Shifts via Machine Learning Augmented DFT

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. In this work, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aemulus ν: precision halo mass functions in wνCDM cosmologies

Precise and accurate predictions of the halo mass function for cluster mass scales in wνCDM cosmologies are crucial for extracting robust and unbiased cosmological information from upcoming galaxy cluster surveys. Here, we present a halo mass function emulator for cluster mass scales (≳ 1013 M ⊙/h) up to redshift z = 2 with comprehensive support for the parameter space of wνCDM cosmologies allowed by current data. Based on the Aemulus ν suite of simulations, the emulator marks a significant improvement in the precision of halo mass function predictions by incorporating both massive neutrinos and non-standard dark energy equation of state models. This allows for accurate modeling of the cosmology dependence in large-scale structure and galaxy cluster studies. We show that the emulator, designed using Gaussian Process Regression, has negligible theoretical uncertainties compared to dominant sources of error in future cluster abundance studies. Our emulator is publicly available (https://github.com/DelonShen/aemulusnu_hmf), providing the community with a crucial tool for upcoming cosmological surveys such as LSST and Euclid.

cluster counts↗

Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling

Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.

36 MATERIALS SCIENCE↗