Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Small angle scattering of diblock copolymers profiled by machine learning

We outline a machine learning strategy for quantitively determining the conformation of AB-type diblock copolymers with excluded volume effects using small angle scattering. Complemented by computer simulations, a correlation matrix connecting conformations of different copolymers according to their scattering features is established on the mathematical framework of a Gaussian process, a multivariate extension of the familiar univariate Gaussian distribution. We show that the relevant conformational characteristics of copolymers can be probabilistically inferred from their coherent scattering cross sections without any restriction imposed by model assumptions. This work not only facilitates the quantitative structural analysis of copolymer solutions but also provides the reliable benchmarking for the related theoretical development of scattering functions.

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

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Machine Learning of Performance Space Mapping for the DPD Simulation of Drug Delivery to Endothelial Cells

Despite huge effort over the years, the design of functionalised nanocarriers (NCs) for targeted drug delivery to endothelial cells is still to be completely unveiled. Dissipative particle dynamics (DPD) simulation is used to study the adhesion of NCs to endothelial cells under the influence of series of parameters such as shape, size, and ligand density of the NCs. However, preparing a performance space mapping that illustrates the penetration depths of NCs as a function of variations in their properties requires simulations of all possible NCs with the above-mentioned properties, which is not feasible. This challenge was addressed by leveraging a Gaussian process regression (GPR)-informed active learning strategy and an extensive exploration of numerous samples, each representing different properties of NCs. The performance space mapping reveals that NCs with rod and disc shapes exhibit superior penetration capabilities compared to those with a spherical shape. Furthermore, it demonstrates that smaller-sized rod-shaped NCs and larger-sized disc-shaped NCs tend to achieve better penetration. When considering smaller NCs, the influence of ligand density appears to be limited. On the contrary, for larger NCs, an increase in ligand density correlates with greater penetration depth, underscoring its substantial role in shaping their penetration capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of phonon softening induced by anisotropic fluctuations in the enhanced mobility at free glassy surfaces

The surface of a glassy material exhibits enhanced mobility compared to the bulk counterpart, however the underlying mechanism for this remains elusive. Herein, we present studies of the dynamical properties of a prototypical glass-forming metallic liquid Zr 50 Cu 50 as a function of the distance from both the free surface and pinned surface using molecular dynamics simulations. We found that the surface mobility increases gradually on approaching the free surface, with a concomitant increase of the non-Gaussianity. The phonon density of states at the free surface exhibits lower characteristic frequencies than in the bulk and pinned surface. These findings suggest phonon softening caused by anisotropic fluctuations at free surfaces as an alternative physical mechanism leading to the enhanced dynamics at free glassy surfaces from the perspective of collective excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2022

Fast increased fidelity samplers for approximate Bayesian Gaussian process regression

Gaussian processes (GPs) are common components in Bayesian non-parametric models having a rich methodological literature and strong theoretical grounding. The use of exact GPs in Bayesian models is limited to problems containing several thousand observations due to their prohibitive computational demands. We develop a posterior sampling algorithm using H-matrix approximations that scales at O(n log 2 n). We show that this approximation’s Kullback-Leibler divergence to the true posterior can be made arbitrarily small. Though multidimensional GPs could be used with our algorithm, d-dimensional surfaces are modeled as tensor products of univariate GPs to minimize the cost of matrix construction and maximize computational efficiency. We illustrate the performance of this fast increased fidelity approximate GP, FIFA-GP, using both simulated and non-synthetic data sets

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Curvature Energetics Determined by Alchemical Simulation on Four Topologically Distinct Lipid Phases

The relative curvature energetics of two lipids are tested using thermodynamic integration (TI) on four topologically distinct lipid phases. Simulations use TI to switch between choline headgroup lipids (POPC; that prefers to be flat) and ethanolamine headgroup lipids (POPE; that prefer, for example, the inner monolayer of vesicles). Here, the thermodynamical moving of the lipids between planar, inverse hexagonal (H II ), cubic (Q II ; Pn3m space group), and vesicle topologies reveals differences in material parameters that were previously challenging to access. The methodology allows for predictions of two important lipid material properties: the difference in POPC/POPE monolayer intrinsic curvature (ΔJ 0 ) and the difference in POPC/POPE monolayer Gaussian curvature modulus (Δ$\bar{κ_m}$), both of which are connected to the energetics of topological variation. Analysis of the TI data indicates that, consistent with previous experiment and simulation, the J 0 of POPE is more negative than POPC (ΔJ 0 = –0.018 ± 0.001 Å –1 ). The theoretical framework extracts significant differences in $\bar{κ_m}$ of which POPE is less negative than POPC by 2.0 to 4.0 kcal/mol. The range of these values is determined by considering subsets of the simulations, and disagreement between these subsets suggests separate mechanical parameters at very high curvature. Finally, the fit of the TI data to the model indicates that the position of the pivotal plane of curvature is not constant across topologies at high curvature. Overall, the results offer insights into lipid material properties, the limits of a single HC model, and how to test them using simulation.

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

Optimization and Evaluation of Stochastic Unified Convection Using Single-Column Model Simulations at Multiple Observation Sites

We extend the previously developed stochastic unified convection scheme (UNICON) for shallow convection to deep convection by parameterizing the impact of mesoscale organized flow on updraft properties. The extended stochastic UNICON parameterizes thermodynamic properties of updrafts at the near-surface as a multivariate Gaussian distribution, where the variances of the distribution are the summation of variances from non-organized turbulence and mesoscale organized flow. The distribution of updraft radius is parameterized as a power-law distribution with a scale break which is parameterized as a linear function of the strength of mesoscale organized flow. The proposed parameterization is validated using a series of large-eddy simulations of deep convection. The free parameters introduced in the formulation of stochastic UNICON are optimized using 10 cases of single-column model simulations over the ocean. Stochastic UNICON with the optimized parameters significantly reduces the biases of thermodynamic profiles and surface precipitation rates simulated in the original UNICON for tropical convection cases. The simulation of the variation in anomalies of temperature and moisture associated with the Madden-Julian oscillation is also improved. The overall improvements in simulated thermodynamic profiles are found to be due to the increased heating and drying tendencies by convective processes in stochastic UNICON. An additional simulation of an idealized deep convection case shows that stochastic UNICON produces enhanced cloud variabilities with dependency on updraft radius, indicating its ability to represent the coexistence of shallow and deep convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Stochastic nonlinear analysis of unidirectional fiber composites using image-based microstructural uncertainty quantification

Here, we present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

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

Fresnel zone plate point spread function approximation for zeroth order mitigation in millimetric field of view x-ray imaging

High spatial and temporal resolution x-ray radiography images are required at the National Ignition Facility (NIF) for high-energy density experiments. One technique that is in development to achieve the required resolution uses Fresnel zone plate (FZP) optics to image an object that is backlit by an x-ray source. The multiple FZP diffraction orders do not focus on the same plane, which increases the background and reduces the contrast. Understanding the point spread function of the different diffraction orders will allow the prediction of the expected background using simulations. We find that the two-dimensional point spread function of the FZP can be approximated by the addition of a sharp Gaussian with a disk. This allowed for the estimation of the background in NIF experimental images of Rayleigh–Taylor spikes and their interpretation. An alternative design of FZP is discussed to allow the inclusion of a zeroth order blocker to reduce the background.

Do, A. (ORCID:0000000174372147)↗
DOE OSTI2025

High-Temperature Gas Sensor Materials with Properties Predicted via First-Principles Calculations with Machine Learning Modeling and Experimental Corroboration

Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.

bandgap renormalization↗