Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel 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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Compact Ferroelectric Programmable Majority Gate for Compute-in-Memory Applications

In this study, a compact and novel ferroelectric (FE) programmable majority gate is proposed and its novel application in Binary Neural Network (BNNs) is investigated. We demonstrate: i) by integrating N metal-ferroelectric-metal (MFM) capacitors on the gate of a transistor (1T-N-MFM structure), a nonvolatile and programmable majority (MAJ) gate that performs MAJ of AND between the gate input and polarization is realized; ii) validation the functionality of our 3-input MAJ of AND gate through comprehensive theoretical and experimental investigations; iii) a compact implementation of 3-input MAJ of XNOR gate that leverages only five of our 3-input MAJ of AND gates connected in parallel; iv) application of MAJ of XNOR gates to replace the XNOR gates and the first layer of the adder tree in the BNNs for up to 21x area saving on top of eliminating the energy-hungry memory accesses due to the compute-in-memory nature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Defect‐Free Sound Insulator Using Single Metal‐Based Friction Stir Process Array

Metals are excellent conductors for phonon transportation such as vibration, sound, and heat. Generally, metal sound insulators require multimaterial structure or defects and unimetal sound insulators are challenging. Therefore, a design of a defect-free sound insulator made by single alloys with multiple friction stir processes (FSPs) is proposed. Periodic friction stir processing can induce superlattice-like local mechanical properties’ modifications. By experimental acoustic characterization, it is observed that FSP can introduce clear acoustic–elastic property contrast on an aluminum plate by the presence of stir zone and heat-affected zones. In numerical simulations, the signature FSP-induced property profile is periodically and parallelly arranged on a long aluminum plate. The transmission gap frequencies are present on the frequency spectrum with the sound propagation direction perpendicular to the FSP paths. Disorder offsets on FSP periodicity are further introduced. Anderson localization is found on a resonance frequency, which provides −11 dB sound reduction by an exponential decay. Due to the finite design length, the slight disorder can also enhance sound insulation in the periodic transmission gap frequency. In conclusion, with analysis and comparison with different configurations, the best performance in the models can achieve −30 dB sound insulation in the 350 mm-long aluminum alloy plate with 14 parallel FSPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A geo-visual analysis for exploring the socioeconomic benefits of the heating electrification using geothermal energy

In parallel to population growth and climate change, the rapid pace of urbanization worldwide has led to an enormous increase in energy demand and costs in urban areas. The subsequent energy burden has become an increasing concern for many households in the U.S. Previous studies have revealed that geothermal resources can effectively lower the electricity demand and carbon emissions in large cities. In this paper, we focus on the socioeconomic impacts of geothermal energy on urban systems by presenting an interactive visual analytics dashboard. The dashboard allows urban planners to spatially examine geothermal energy's practical benefits on energy affordability, urban livability, and resilience across the U.S. We compiled a list of socioeconomic metrics by integrating the simulation results from multiple geothermal and building models with multi-domain urban datasets (socioeconomic, demographic, and electricity utility). These metrics are created to characterize the benefits of the heating electrification of buildings using Geothermal Heat Pumps (GHP) for lowering the energy burden of middle-and low-income households nationwide. The visual dashboard employs a combination of multivariate, glyph-based, and geospatial visualization to reveal the variability and patterns in our metrics. We present a pilot study to demonstrate the GHPs' potential as a renewable and affordable solution for increasing the economic and energy grid resilience in U.S cities.

Xu, Haowen↗
DOE OSTI2021

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗
DOE OSTI2024

Green's function methods for multiphysics simulations (Final Report)

Green's functions are important tools for analyzing mathematical properties of partial differential equations (PDEs), and for numerically solving PDEs, especially when equations for the same operator but with multiple right hand sides need to be solved simultaneously. This proposal aims at developing efficient and accurate numerical methods for computing Green's functions, which can be used to tackle a challenging question in multiphysics simulation of DOE-mission science problems: how to couple quantum physics with classical physics. The key mathematical difficulty is properly formulate a "boundary condition" for the region described by quantum physics, and conventional approaches often model such boundary conditions in an empirical way. The proposed Green's function methods use the Dirichlet-to-Neumann map to formulate a boundary condition that is non-empirical and can couple the quantum and classical regions in an in principle exact way. The key ingredient of the new methods is to construct the Dirichlet-to-Neumann map in an efficient, accurate and versatile manner. The new methods have provably low complexity and are ideally suited for massively parallel and emerging many-core computational systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

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

Generalized Power Flow Model of an Extra High-Power Multi-Terminal HVdc Transmission Grid with Parallel-Connected and Voltage-Stacked Converters

To increase power transfer capacity of high-voltage direct current (HVdc) transmission, a new extra high-power HVdc architecture with multiple standard modular multilevel converters (MMCs) per substation has recently been introduced. This paper proposes a power flow model for a multi-terminal HVdc (MTdc) grid with this innovative substation architecture. The proposed MTdc model can be integrated seamlessly with existing ac-dc power flow algorithms with minimal modifications. The model is applicable to various MTdc grid types and topologies, different numbers of dc buses, dc lines, and MMCs per substation, along with diverse control parameters. In addition, it accurately captures both balanced and unbalanced operations of the MTdc grid. The proposed model is verified using a 5-terminal bipole MTdc grid that spans 4 areas in the Eastern Interconnection system of the USA. The numerical solutions obtained from unified and sequential ac-dc power flow algorithms under different operating conditions closely match the time-domain simulation results in PSCAD, validating the accuracy and versatility of the proposed MTdc power flow model.

Nguyen, Quan H.↗
DOE OSTI2020

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N -body mock challenge for the quasar sample

ABSTRACT The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock–Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of halo occupation distribution mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover fσ8 to within 3 per cent and α∥ and α⊥ to within 1 per cent. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on fσ8, α∥, and α⊥ at the level of $\sigma _{f\sigma _8}=0.013$, $\sigma _{\alpha _\parallel }=0.012$, and $\sigma _{\alpha _\bot }=0.008$. The systematic error on the combined consensus is $\sigma _{f\sigma _8}=0.011$, $\sigma _{\alpha _\parallel }=0.008$, and $\sigma _{\alpha _\bot }=0.005$, which is used in the final DR16 analysis. For baryon acoustic oscillation fits in configuration and Fourier space, we take conservative systematic errors of $\sigma _{\alpha _\parallel }=0.010$ and $\sigma _{\alpha _\bot }=0.007$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Magnetic electron collimation in three-dimensional semi-metals

While electrons moving perpendicular to a magnetic field are confined to cyclotron orbits, they can move freely parallel to the field. This simple fact leads to complex current flow in clean, low carrier density semi-metals, such as long-ranged current jets forming along the magnetic field when currents pass through point-like constrictions. Occurring accidentally at imperfect current injection contacts, the phenomenon of "current jetting" plagues the research of longitudinal magneto-resistance, which is particularly important in topological conductors. Here we demonstrate the controlled generation of tightly focused electron beams in a new class of micro-devices machined from crystals of the Dirac semi-metal Cd 3 As 2 . The current beams can be guided by tilting a magnetic field and their range tuned by the field strength. Finite element simulations quantitatively capture the voltage induced at faraway contacts when the beams are steered towards them, supporting the picture of controlled electron jets. These experiments demonstrate direct control over the highly non-local signal propagation unique to 3D semi-metals in the current jetting regime, and may lead to applications akin to electron optics in free space.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Collective mode Brownian dynamics: A method for fast relaxation of statistical ensembles

Sampling equilibrium configurations of correlated systems of particles with long relaxation times (e.g., polymeric solutions) using conventional molecular dynamics and Monte Carlo methods can be challenging. This is especially true for systems with complicated, extended bond network topologies and other interactions that make the use and design of specialized relaxation protocols infeasible. We introduce a method based on Brownian dynamics simulations that can reduce the computational time it takes to reach equilibrium and draw decorrelated samples. Importantly, the method is completely agnostic to the particle configuration and the specifics of interparticle forces. In particular, we develop a mobility matrix that excites non-local, collective motion of N particles and can be computed efficiently in $\mathcal{O}(N)$ time. Particle motion in this scheme is computed by integrating the overdamped Langevin equation with an Euler–Maruyama scheme, in which Brownian displacements are drawn efficiently using a low-rank representation of the mobility matrix in position and wave space. We demonstrate the efficacy of the method with various examples from the realm of soft condensed matter and release a massively parallel implementation of the code as a plugin for the open-source package HOOMD-blue [J. A. Anderson et al., J. Comput. Phys. 227, 5342 (2008) and J. Glaser et al., Comput. Phys. Commun. 192, 97 (2015)] which runs on graphics processing units

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

Optical spin hall effect in exciton–polariton condensates in lead halide perovskite microcavities

An exciton–polariton condensate is a hybrid light–matter state in the quantum fluid phase. The photonic component endows it with characters of spin, as represented by circular polarization. Spin-polarization can form stochastically for quasi-equilibrium exciton–polariton condensates at parallel momentum vector k|| ∼ 0 from bifurcation or deterministically for propagating condensates at k|| > 0 from the optical spin-Hall effect (OSHE). Here, we report deterministic spin-polarization in exciton–polariton condensates at k|| ∼ 0 in microcavities containing methylammonium lead bromide perovskite (CH3NH3PbBr3) single crystals under non-resonant and linearly polarized excitation. We observe two energetically split condensates with opposite circular polarizations and attribute this observation to the presence of strong birefringence, which introduces a large OSHE at k|| ∼ 0 and pins the condensates in a particular spin state. Such spin-polarized exciton–polariton condensates may serve not only as circularly polarized laser sources but also as effective alternatives to ultracold atom Bose–Einstein condensates in quantum simulators of many-body spin–orbit coupling processes.

Chemistry↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Glass-bonded ceramic waste forms for immobilization of radioiodine from caustic scrubber wastes

Glass-bonded sodalite composite waste forms have been developed for the immobilization of liquid radioactive wastes resulting from off-gas treatment during aqueous reprocessing of used nuclear fuel, with a particular focus on 129I. The proposed composite waste form is comprised of aluminosilicate ceramic phases containing volatile radionuclides bonded with a glassy matrix. In this work, a suite of ten candidate low-temperature glass binders (ZnO-Bi2O3-based glasses and a Na2O-B2O3-SiO2 glass) were examined. Six glasses were mixed with caustic scrubber waste simulant previously converted into a sodalite-rich material (to provide glass fractions of 10 and 20 wt.%), uniaxially pressed into pellets, and sintered at 350 °C or 550 °C for 8 h in air. Iodine retention after heat treatment was assessed by neutron activation analysis, showing retention of 67-100 % of expected iodine. The aqueous durabilities of the resulting materials were then determined, following the ASTM C1308 standard test, showing iodine releases of 1 to 23 g m-2 after 4 d. The cumulative iodine release for the best performing system (a zinc-bismuth-borate glass binder) was <1 g m-2, and its iodine retention from processing was 67 %. The iodine releases compared favorably with other waste forms. In parallel, this best-performing composition was also consolidated via hot isostatic pressing (HIP) in a stainless-steel canister at 550 °C for 2 h under 100 MPa pressure. The HIPed sample was produced at the ~20 g scale and showed improved densification and minimal reaction with the canister.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

High-Throughput Electric-Field-Assisted Sintering and Characterization Techniques for Materials Discovery

Despite improvements in computing and modeling capabilities, the performance of new materials, particularly those which deviate greatly in composition from well-studied materials (e.g., high-entropy alloys), can be difficult to simulate given the lack of available experimental property data. While some modeling techniques may attempt to predict the properties of these exotic materials, most are forced to make extrapolations from more traditional materials. To fulfill the need for accelerated material synthesis and property measurement, a high-throughput methodology has been developed. Utilizing electric-field-assisted sintering (EFAS), also known as spark plasma sintering (SPS), equipped with custom tooling, samples of differing alloy compositions can be produced simultaneously as a single alloy array. Several arrays have been produced with compositions spanning the Co-Cr-Fe-Mn-Ni alloy family, including many high-entropy alloys, while the novel array geometry has enabled the samples to be polished and characterized in parallel, using X-ray diffraction, scanning-electron microscopy, and laser-based thermal diffusivity measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Pulse width dependence of magnetic field generation using laser-powered capacitor coils

Megagauss magnetic fields were generated by a current flowing through a U-shaped coil connecting two parallel copper foils. Two kJ-class lasers at various pulse widths from 2 ns to 9.9 ns passed through holes in the front foil and were focused on the back foil with an intensity of ~1.7×10 16 W/cm 2 . The coil current and resulting magnetic fields were characterized using ultrafast proton radiography, timed at the end of the laser pulses. Additionally, the measurements show that magnetic field strength decays with increasing laser pulse width. A lumped-circuit model was developed and showed consistency with the experimental measurements, demonstrating an ion shorting effect: as the ion current neutralizes the electron current contribution to interplate voltage, the coil current peaks on a timescale close to the ion transit time t i = d/v ion . FLASH simulations of the coil current are performed, and the calculated resistance values are used to constrain ion speed as a function of hot electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗
DOE OSTI2023

GPU acceleration of Swendsen–Wang dynamics

When simulating a lattice system near its critical temperature, local algorithms for modeling the system’s evolution can introduce very large autocorrelation times into sampled data. Here, this critical slowing down places restrictions on the analysis that can be completed in a timely manner of the behavior of systems around the critical point. Because it is often desirable to study such systems around this point, a new algorithm must be introduced. Therefore, we turn to cluster algorithms, such as the Swendsen–Wang algorithm and the Wolff clustering algorithm. They incorporate global updates which generate new lattice configurations with little correlation to previous states, even near the critical point. We look to accelerate the rate at which these algorithm are capable of running by implementing and benchmarking a parallel implementation of each algorithm designed to run on GPUs under NVIDIA’s CUDA framework. A 17 and 90 fold increase in the computational rate was, respectively, experienced when measured against the equivalent algorithm implemented in serial code.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗