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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ORNL/grit

Grit is a software library designed specifically for conducting particle-based Lagrangian simulations with CPU/GPUperformance portability. This library allows researchers and engineers to perform a wide range of simulations, including multiphase flow, and more. Grit employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs.

Ge, Wenjun [Oak Ridge National Laboratory (ORNL), ↗
DOE OSTI2020

Hydration Mechanisms in Nanoparticle Interaction and Surface Energetics

Work on the project advanced molecular level understanding, prediction and control of nanoscale hydration in salt solutions under electric stimuli or ionic patterning. Permeation of nanoporous electrodes at preset voltage underlies the function of ultracapacitors. Transitory regulation of wetting in nanoporous media by electric field spans an array of applications in materials, energy storage, and separation sciences. The physically related modulation of nanoparticle solubility by surface charges can significantly extend the range of the nanomaterial applications, improve processing techniques, and potentially alleviate environmental concerns. Pore/solution equilibria, and activated kinetics of liquid gating by nanoconfined electrolytes are challenging problems at the forefront of experimental and theoretical research. Addressing these problems from a molecular perspective required the development of state-of-the-art simulation algorithms in statistical mechanics to capture complex processes in open systems under electric control. Parallel studies of wetting and dispersibility of polar and ionizing particles aim to uncover predictive relations between electrowetting, chemical functionalization, and geometry of nanomaterial particles. By nonequilibrium dynamic modeling, we elucidated electrolyte flow in nanochanels and associated electrokinetic energy conversion. Research on the program provided training opportunities for the next generation of scientists in computational chemistry. Insights, and methods from the project implicate broad segment of researchers in energy and nanosciences, materials and surface chemistry.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

SST-GPU: A Scalable SST GPU Component for Performance Modeling and Profiling

Programmable accelerators have become commonplace in modern computing systems. Advances in programming models and the availability of unprecedented amounts of data have created a space for massively parallel accelerators capable of maintaining context for thousands of concurrent threads resident on-chip. These threads are grouped and interleaved on a cycle-by-cycle basis among several massively parallel computing cores. One path for the design of future supercomputers relies on an ability to model the performance of these massively parallel cores at scale. The SST framework has been proven to scale up to run simulations containing tens of thousands of nodes. A previous report described the initial integration of the open-source, execution-driven GPU simulator, GPGPU-Sim, into the SST framework. This report discusses the results of the integration and how to use the new GPU component in SST. It also provides examples of what it can be used to analyze and a correlation study showing how closely the execution matches that of a Nvidia V100 GPU when running kernels and mini-apps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗
DOE OSTI2022

GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multistart gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 s. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of ${ \mathcal O }({10}^{5})$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1978

Measurement of transient strain and surface temperature on simulated turbine blades using noncontacting techniques

Noncontacting techniques were used to measure strain and temperature in thermally cycled simulated turbine blades. An electro-optical extensometer was used to measure the displacement between parallel targets mounted on the leading edge of the blades throughout a complete heating and cooling cycle. An infrared photographic pyrometry method was used to measure blade steady state surface temperature. The blade was cyclically heated and cooled by moving it into and out of a Mach 1 hot-gas stream. Transient leading edge strain and steady state surface temperature distributions are presented for blades of three different configurations.

Calfo, F. D.↗
NASA NTRS1979

Calculation of transonic aileron buzz

An implicit finite-difference computer code that uses a two-layer algebraic eddy viscosity model and exact geometric specification of the airfoil has been used to simulate transonic aileron buzz. The calculated results, which were performed on both the Illiac IV parallel computer processor and the Control Data 7600 computer, are in essential agreement with the original expository wind-tunnel data taken in the Ames 16-Foot Wind Tunnel just after World War II. These results and a description of the pertinent numerical techniques are included.

Steger, J. L.↗
NASA NTRS1982

Pockels-effect cell for gas-flow simulation

A Pockels effect cell using a 75 cu cm DK*P crystal was developed and used as a gas flow simulator. Index of refraction gradients were produced in the cell by the fringing fields of parallel plate electrodes. Calibration curves for the device were obtained for index of refraction gradients in excess of .00025 m.

Weimer, D.↗
NASA NTRS1984

Two dimensional blade-vortex interaction flow visualization study

Blade-vortex interaction occurs when a rotor blade encounters the tip vortex from a previous rotor blade. To obtain details of the close encounter process, the results from a flow visualization study of an airfoil representing a rotor blade in the wake of an oscillating airfoil serving as a vortex generator are described. A distinguishing feature of this study is that the vortex filament is oriented parallel to the blade span, orthogonal to the test section free stream velocity. This orientation simulates the case of two-dimensional blade-vortex interaction, which is known to produce the most impulsive and most intensive BVI noise. Photographic data are examined to deduce qualitative and quantitative details of the close encounter interaction process with emphasis on structural changes in the vortex filament and its trajectory.

Booth, E. R., Jr.↗
NASA NTRS1986

Ice banding as a response of the coupled ice-ocean system to temporally varying winds

This study models formation of ice bands in the marginal ice zones. A one-dimensional coupled ice-ocean model is used in which the ice model is coupled to a reduced gravity ocean model through interfacial stresses. The internal ice stresses are important only at high ice concentrations (90-100 percent); otherwise, the main balance for the ice motion is between the air-ice stress and the ice-water stress, i.e., free drift. The drag coefficients were chosen so that the air-ice momentum flux is 3 times greater than the air-ocean momentum flux. Thus the Ekman transport is larger under the ice than in the open water, so that winds parallel to the ice edge, with the ice on the right, produce upwelling. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4-day period. This forcing resembles successive cyclone passings perpendicular to the ice edge. When the oceanic upper layer was thin, which means that the dynamics are strongly nonlinear, the ice bands were formed. The up/downwelling signals do not disappear in wind reversals because of nonlinear advection. This leads to convergences and divergences in oceanic and ice velocities that manifest themselves as ice banding. At least one wind reversal is needed to produce one ice band.

Hakkinen, S.↗
NASA NTRS1988

Thermostructural analysis of simulated cowl lips

Three dimensional finite element analyses using MSC/NASTRAN and MARC are performed to predict the thermal and structural response of various cooling schemes under high heat loads. Steady state heat transfer analyses and elastic stress analyses are performed using MSC/NASTRAN. Elastic/plastic analyses are done using MARC. To help verify these analyses experimentally, a hydrogen-oxygen rocket engine was modified to use the exhaust stream as a high enthalpy, high heat flux source to evaluate various actively cooled, simulated cowl lip (leading edges) segments as well as flat structural segments. Cross flow and parallel flow cooling configurations were tested and analyzed using cooling fluids of water and gaseous hydrogen. In addition, various material types, including high conductivity copper, nickel, and a copper and graphite metal matrix composite were tested and compared.

Melis, Matthew E.↗
NASA NTRS1988

Evolution of nonlinear polarization in localized and finite amplitude Alfven waves

Theoretical and computational study for a polarization change of localized and finite amplitude Alfven waves propagating parallel to an applied magnetic field is presented using both reductive perturbation theory and numerical simulations. In the magnetohydrodynamic limit, where right-hand and left-hand circularly polarized waves are degenerate, one of the transverse components of the circularly polarized Alfven wave is stable as it propagates, but the other component is unstable to either self-focusing or diffraction effects. Consequently, the wave changes its polarization from circular to linear. In the high frequency regime, where two circularly polarized waves (right-hand and left-hand circularly polarized waves) are not degenerate, two transverse components of the circularly polarized Alfven wave are strongly coupled to each other, and there is almost no polarization change.

Hoshino, M.↗
NASA NTRS1990

Multi-level Hierarchical Poly Tree computer architectures

Based on the concept of hierarchical substructuring, this paper develops an optimal multi-level Hierarchical Poly Tree (HPT) parallel computer architecture scheme which is applicable to the solution of finite element and difference simulations. Emphasis is given to minimizing computational effort, in-core/out-of-core memory requirements, and the data transfer between processors. In addition, a simplified communications network that reduces the number of I/O channels between processors is presented. HPT configurations that yield optimal superlinearities are also demonstrated. Moreover, to generalize the scope of applicability, special attention is given to developing: (1) multi-level reduction trees which provide an orderly/optimal procedure by which model densification/simplification can be achieved, as well as (2) methodologies enabling processor grading that yields architectures with varying types of multi-level granularity.

Padovan, Joe↗
NASA NTRS1990

Prominence formation in a coronal loop

A model is presented which depends on the preferential deposition of heating in the legs of a coronal loop and which produces a stable prominence-scale condensation at the loop top. Dynamic stability is attained by the subsequent adjustment of local parallel gravity by a magnetic inversion at the loop (or arcade) apex. A nonlinear numerical simulation of this process, which includes a deep chromosphere, a heating rate with a fixed dissipation length, and full solar gravity is described.

Mok, Y.↗
NASA NTRS1992

Electron acceleration by low frequency electric field fluctuations - Electron conics

Electron conics can be generated by fluctuations at roughly one Hz in a parallel electric field accelerating auroral electrons at altitudes of several thousand kilometers. The results from simulations including such time-varying electric fields are in good agreement with observed electron distributions. The simulation parameters are consistent with the ion beams and broadband waves associated with electron conics. Calculations are used to predict for example oscillations in the low altitude electron fluxes. These results can be used to compare different possible mechanisms for the generation of electron conics.

Andre, Mats↗
NASA NTRS1992

Three-dimensional structure of coronal mass ejections

The 3D structure of coronal mass ejections (CMEs) are examined on the basis of numerical solutions of the MHD equations. The numerical solution in a meridional plane centered over the region of newly emerged flux has features similar to those computed in previous 2D simulations. The CME is found to be shaped more like an arcade with an axis parallel to the Equator than like a bubble. This characteristic arcade shape occurs both in the density distribution in meridional and azimuthal planes and in the integrated white-light brightness in meridional and polar views. Examination of the brightness change due to the CME shows that the CME in a meridional projection as viewed from the earth has the same loop shape even if the CME originates almost directly in front of or behind the sun.

Steinolfson, R. S.↗
NASA NTRS1994

A Procedure for 3-D Contact Stress Analysis of Spiral Bevel Gears

Contact stress distribution of spiral bevel gears using nonlinear finite element static analysis is presented. Procedures have been developed to solve the nonlinear equations that identify the gear and pinion surface coordinates based on the kinematics of the cutting process and orientate the pinion and the gear in space to mesh with each other. Contact is simulated by connecting GAP elements along the intersection of a line from each pinion point (parallel to the normal at the contact point) with the gear surface. A three dimensional model with four gear teeth and three pinion teeth is used to determine the contact stresses at two different contact positions in a spiral bevel gearset. A summary of the elliptical contact stress distribution is given. This information will be helpful to helicopter and aircraft transmission designers who need to minimize weight of the transmission and maximize reliability.

Kumar, A.↗