Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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

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

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development of Efficient Process for Manufacturing of Thermoplastic Composites with Tailored Properties (CRADA 511)

This is a collaborative effort between Battelle Memorial Institute as manager and operator of Pacific Northwest National Laboratory (PNNL) and ESI North America Inc. (“ESI” or “Participant”) to apply computation and data analytics to the challenge of light weighting with a focus on the battery enclosures of electric vehicles (EVs). EVs use heavy batteries to increase range and power. A complex-shaped battery enclosure is required to meet a host of challenging performance requirements. The ability to virtually develop composite parts such as battery enclosure with tailored properties to meet required performance will be highly valuable to the automotive industry. However, efficient simulation of composite-manufacturing processes remains a challenging issue since simulation involves multiscale models in space and time, highly non-linear and anisotropic behavior, strongly coupled multi-physics, and complex geometries. This work will advance the state of the art by reducing the computational burden of composite optimization by using simulation data from a limited number of configurations off-line and then developing a reduced order model (ROM) using data analytics and machine learning (ML). Develop a data driven approach to link features of the material and manufacturing processes to the mechanical properties of thermoplastic composite parts.

42 ENGINEERING↗
DOE OSTI2022

Thiocyanate Ions Form Antiparallel Populations at the Concentrated Electrolyte/Charged Surfactant Interface

Anions play significant roles in the separation of lanthanides and actinides. The molecular-scale details of how these anions behave at aqueous interfaces are not well understood, especially at high ionic strengths. Here, we describe the interfacial structure of thiocyanate anions at a soft charged interface up to 5 M bulk concentration with combined classical and phase-sensitive and molecular dynamics (MD) simulations. At low concentrations thiocyanate ions are mostly oriented with their sulfur end pointing toward the charged surfactants. In this work, the VSFG signal reaches a plateau at around 100 mM bulk concentration, followed by significant changes above 1 M. At high concentrations a new thiocyanate population emerges with their sulfur end pointing toward the bulk liquid. The –CN stretch frequency is different for up and down oriented SCN – ions, indicating different coordination environments. These results provide key molecular-level insights for the interfacial behavior of complex anions in highly concentrated solutions.

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

Quantifying Temperature Dependence of Electronic Band Gaps and Optical Properties in SnO 2 and SnO via First-Principles Simulations

Tin metal oxides SnO x ( x = 1, 2) have gained interest as gas-sensing materials. For their applications as high-temperature sensors, however, a better understanding of their temperature dependence sensing responses is needed. Here, we comparatively quantify the temperature-dependent electronic band gaps and optical properties of SnO 2 and SnO using first-principles calculations. Without considering the temperature effect, SnO 2 and SnO are predicted to have direct and indirect band gaps of 2.18 and 1.66 eV, respectively, at the PBE + U -GGA level that we employed. The temperature effect on the electronic and optical properties is captured by taking account of the electron–phonon interaction. Band gap renormalization with temperature is calculated via the Allen–Heine–Cardona theory. For both oxides, we find a monotonic decrease in the electronic band gap such that renormalization at zero point (0 K) is ~–0.17 and ~–0.52 eV at 1000 K. These results are also analyzed by employing an analytical equation that helps characterize the band gap shift with temperature. In addition, the optical properties at finite temperatures are simulated using the frozen-phonon method that combines electron–phonon coupling with the momentum matrix. As temperature increases, the optical property spectra are smoothed because of the smearing effect, which diminishes optical constants at shorter wavelengths. Our results are of interest for high-temperature functional materials in applications of optical detection.

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

Biomembrane Structure and Material Properties Studied With Neutron Scattering

Cell membranes and their associated structures are dynamical supramolecular structures where different physiological processes take place. Detailed knowledge of their static and dynamic structures is therefore needed, to better understand membrane biology. The structure–function relationship is a basic tenet in biology and has been pursued using a range of different experimental approaches. In this review, we will discuss one approach, namely the use of neutron scattering techniques as applied, primarily, to model membrane systems composed of lipid bilayers. An advantage of neutron scattering, compared to other scattering techniques, is the differential sensitivity of neutrons to isotopes of hydrogen and, as a result, the relative ease of altering sample contrast by substituting protium for deuterium. This property makes neutrons an ideal probe for the study of hydrogen-rich materials, such as biomembranes. In this review article, we describe isotopic labeling studies of model and viable membranes, and discuss novel applications of neutron contrast variation in order to gain unique insights into the structure, dynamics, and molecular interactions of biological membranes. We specifically focus on how small-angle neutron scattering data is modeled using different contrast data and molecular dynamics simulations. We also briefly discuss neutron reflectometry and present a few recent advances that have taken place in neutron spin echo spectroscopy studies and the unique membrane mechanical data that can be derived from them, primarily due to new models used to fit the data.

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

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗
DOE OSTI2022

The basics of quantum computing for chemists

The rapid and successful strides in quantum chemistry in the past decades can be largely credited to a conspicuous synergy between theoretical and computational advancements. However, the architectural computer archetype that enabled such a progress is approaching a state of more stagnant development. One of the most promising technological avenues for the continuing progress of quantum chemistry is the emerging quantum computing paradigm. This revolutionary proposal comes with several challenges, which span a wide array of disciplines. In chemistry, it implies, among other things, a need to reformulate some of its long established cornerstones in order to adjust to the operational demands and constraints of quantum computers. Due to its relatively recent emergence, much of quantum computing may still seem fairly nebulous and largely unknown to most chemists. It is in this context that here we review and illustrate the basic aspects of quantum information and their relation to quantum computing insofar as enabling simulations of quantum chemistry. As a result, we consider some of the most relevant developments in light of these aspects and discuss the current landscape when of relevance to quantum chemical simulations in quantum computers.

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

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical “slaving onset” stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

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

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

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

Validation of Fast Reactor Depletion Tools Using EBR-II Measured Data

The validation of simulation tools for calculating fuel depletion and evolution in fast reactors is vital for design, licensing, deployment, operations, and material accountancy. The Physics Analysis Database (PADB) and Analytical Laboratory (AL) database contain measured data collected from Experimental Breeder Reactor II (EBR-II) and were used to validate the most recent versions of the Argonne Reactor Computation (ARC) tool suite and ORIGEN-S for calculating isotopic compositions in irradiated fast reactor fuel. The PADB contains important modeling and operational information about the EBR-II core design, fuel cycle, and analytical results from the legacy versions of the ARC tool suite. The AL database contains the measured isotopic compositions of irradiated samples taken from core subassemblies. A new procedure was developed for the ARC tool suite to perform the EBR-II depletion simulation, as well as to perform more detailed isotopic calculations using ORIGEN-S calculations by coupling it with the ARC suite. Both the ARC and ARC-ORIGEN results were compared with the AL measured data for all relevant samples and showed good agreement for the major actinides. Good agreement with measured data was also achieved using the ARC-ORIGEN approach for several fission products.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2020

Non-cancellation of the parity anomaly in the strong-field regime of QED 2+1

Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in 2 + 1 space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

New Simulations for Ion-Production and Back-Bombardment in GaAs Photo-guns

GaAs-based DC high voltage photo-guns used at accelerators with extensive user programs must exhibit long photocathode operating lifetime. Achieving this goal represents a significant challenge for proposed high average current facilities that must operate at tens of milliamperes or more. Specifically, the operating lifetime is dominated by ion back-bombardment of the photocathode from ionized residual gas. While numerous experiments have been performed to characterize the operating lifetime under various conditions, detailed simulations of the ion back-bombardment mechanism that explains these experiments are lacking. Recently, a new user routine was implemented using the code General Particle Tracer (GPT) to simulate electron impact ionization of residual beam line gas and simultaneously track the incident electron, the ejected electron, and the newly formed ion. This new routine was benchmarked against analytical calculations and then applied to experiments performed at the CEBAF injector at the Thomas Jefferson National Accelerator Facility to study the effectiveness of limiting ions from entering the cathode-anode gap using a positively biased anode. These simulations were performed using detailed 3D field maps produced with CST Microwave Studio describing the photo-gun electrostatics. Discussion of the experiment and the application of this new GPT routine to model the experiments will be presented at the workshop.

Yoskowitz, Joshua↗
DOE OSTI2022

Quantitative analysis of Cu XANES spectra using linear combination fitting of binary mixtures simulated by FEFF9

We report x-ray absorption near edge structure (XANES) is a powerful tool to probe the fingerprint of local structures, and when coupled with X-ray microscopy, the small spot size enables one to probe very specific regions of interest in a device or material, e.g. interfaces/bulk, different grains, good/bad electrical areas. In this work, we investigate the use of linear combination fitting (LCF) of XANES spectra for the particular case of Cu doping in CdTe. We show that the experimental data seem to be accurately represented by standards of Cu 2 Te and its substoichiometric counterpart, Cu 1.43 Te. We use Cu in CdTe as a case study to evaluate the accuracy of linear combination fitting using simulated standards, given that experimental standards for certain phases (e.g. Cu 1.43 Te) or defect structures (e.g. vacancies) cannot be readily obtained. We discuss how spectral features of the FEFF9-simulated standards, fitting ranges, and noise levels all dictate the accuracy of this type of analysis. We show that the greater the spectral difference between the two standards, the better the LCF is able to differentiate between the two structures and to tolerate experimental noise. Finally, we estimate the error of the fitted weights for different spectral features and noise levels and propose a framework to study local structures semi-quantitatively by using binary mixtures of FEFF9-simulated standards.

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

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

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

Role of Molecular Architecture on Ion Transport in Ethylene oxide-Based Polymer Electrolytes

This work aims to develop a detailed mechanistic understanding of the role of a graft polymer architecture on lithium ion (Li + ) transport in poly(ethylene oxide)-based polymer electrolytes. Specifically, we compare Li + transport in poly-(ethylene oxide) (PEO) versus poly(oligo oxyethylene methacrylate) (POEM) polymers doped with lithium bis( trifluoromethanesulfonyl) (LiTFSI) salts, using both experimental electrochemical characterization and molecular dynamics (MD) simulations. Furthermore, our results indicate that POEM exhibits a range of relaxation processes that cannot be interpreted solely in terms of glass-transition temperature (T g ) effects. Due to its side-chain architecture, the segmental relaxation of POEM is nonuniform across ether oxygens (EOs) and shows a more pronounced sensitivity to temperature above T g compared to PEO. Moreover, POEM also exhibits a nonuniform Li + coordination behavior, in which Li + is primarily solvated by two different chains in POEM, compared to a single chain in PEO. Li + transport in POEM occurs via two events with distinct characteristic times: a fast intrachain hopping along side chains and a slow interchain hopping between side chains. Taken together, the relaxation processes and ion transport mechanisms identified in POEM provide useful insights into design of more effective solid polymer electrolytes.

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

Tidal Stripping of a White Dwarf by an Intermediate-mass Black Hole

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10 2-5 M ⊙ ), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z 5/2 , where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M h /10 4 M ⊙ ) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗