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

Usability and Optimization Improvements in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to MOOSE’s broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. Such improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities were implemented in the same order as they are covered in this report: multiple nonlinear systems in the same input file, implementation of generic field transfers and other transfer system enhancements, support for stateful material property redistribution for adaptivity and distributed meshes, and dynamic linking and loading of individually compiled applications. These four additions are fundamental capabilities that will be leveraged by many NEAMS applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Solar panel acceptance testing using a pulsed solar simulator

Utilizing specific parameters as area of an individual cell, number in series and parallel, and established coefficient of current and voltage temperature dependence, a solar array irradiated with one solar constant at AMO and at ambient temperature can be characterized by a current-voltage curve for different intensities, temperatures, and even different configurations. Calibration techniques include: uniformity in area, depth and time, absolute and transfer irradiance standards, dynamic and functional check out procedures. Typical data are given for individual cell (2x2 cm) to complete flat solar array (5x5 feet) with 2660 cells and on cylindrical test items with up to 10,000 cells. The time and energy saving of such testing techniques are emphasized.

Hershey, T. L.↗
DOE OSTI2021

Investigation of turbulent inflow specification in Euler–Lagrange simulations of mid-field spray

The process of atomization of a liquid jet by a parallel high-speed gas stream results in a spray, whose downstream development is of considerable interest to several applications. The round jet spray can be spatially divided into (i) a near-field (near-nozzle) region of liquid atomization and (ii) a downstream mid-field region of fully-dispersed droplets. In order to accurately model mid-field droplet dispersion, this work aims at developing a rigorous and robust injection model for Euler–Lagrange spray simulations. Results from experiments are used to obtain the relevant droplet number density, size distribution, and mean and standard deviation velocity distributions of the injection model, systematically in a step-by-step process. Two-phase large eddy simulations are performed by stochastically generating the Lagrangian droplets at the inlet of the mid-field region. Number flux, diameter distribution, mean velocity, and other time-averaged statistics at several downstream locations are shown to agree well with the corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗
NASA NTRS2016

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗
NASA NTRS1995

High Performance Parallel Computational Nanotechnology

At a recent press conference, NASA Administrator Dan Goldin encouraged NASA Ames Research Center to take a lead role in promoting research and development of advanced, high-performance computer technology, including nanotechnology. Manufacturers of leading-edge microprocessors currently perform large-scale simulations in the design and verification of semiconductor devices and microprocessors. Recently, the need for this intensive simulation and modeling analysis has greatly increased, due in part to the ever-increasing complexity of these devices, as well as the lessons of experiences such as the Pentium fiasco. Simulation, modeling, testing, and validation will be even more important for designing molecular computers because of the complex specification of millions of atoms, thousands of assembly steps, as well as the simulation and modeling needed to ensure reliable, robust and efficient fabrication of the molecular devices. The software for this capacity does not exist today, but it can be extrapolated from the software currently used in molecular modeling for other applications: semi-empirical methods, ab initio methods, self-consistent field methods, Hartree-Fock methods, molecular mechanics; and simulation methods for diamondoid structures. In as much as it seems clear that the application of such methods in nanotechnology will require powerful, highly powerful systems, this talk will discuss techniques and issues for performing these types of computations on parallel systems. We will describe system design issues (memory, I/O, mass storage, operating system requirements, special user interface issues, interconnects, bandwidths, and programming languages) involved in parallel methods for scalable classical, semiclassical, quantum, molecular mechanics, and continuum models; molecular nanotechnology computer-aided designs (NanoCAD) techniques; visualization using virtual reality techniques of structural models and assembly sequences; software required to control mini robotic manipulators for positional control; scalable numerical algorithms for reliability, verifications and testability. There appears no fundamental obstacle to simulating molecular compilers and molecular computers on high performance parallel computers, just as the Boeing 777 was simulated on a computer before manufacturing it.

Saini, Subhash↗
NASA NTRS1998

A Domain-Decomposed Multi-Level Method for Adaptively Refined Cartesian Grids with Embedded Boundaries

The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.

Aftosmis, M. J.↗
NASA NTRS1993

Visual information processing II; Proceedings of the Meeting, Orlando, FL, Apr. 14-16, 1993

Various papers on visual information processing are presented. Individual topics addressed include: aliasing as noise, satellite image processing using a hammering neural network, edge-detetion method using visual perception, adaptive vector median filters, design of a reading test for low-vision image warping, spatial transformation architectures, automatic image-enhancement method, redundancy reduction in image coding, lossless gray-scale image compression by predictive GDF, information efficiency in visual communication, optimizing JPEG quantization matrices for different applications, use of forward error correction to maintain image fidelity, effect of peanoscanning on image compression. Also discussed are: computer vision for autonomous robotics in space, optical processor for zero-crossing edge detection, fractal-based image edge detection, simulation of the neon spreading effect by bandpass filtering, wavelet transform (WT) on parallel SIMD architectures, nonseparable 2D wavelet image representation, adaptive image halftoning based on WT, wavelet analysis of global warming, use of the WT for signal detection, perfect reconstruction two-channel rational filter banks, N-wavelet coding for pattern classification, simulation of image of natural objects, number-theoretic coding for iconic systems.

Huck, Friedrich O.↗
DOE OSTI2020

Dynamics of Laterally Propagating Flames in X-Ray Bursts. I. Burning Front Structure

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Finally, our studies complement other multidimensional studies of burning in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
NASA NTRS1990

Ion Beam Heating in the Auroral Zone

Recent satellite observations at high altitudes (greater than 5000 km) in the auroral zone have shown the existence of hybrid or bimodal ion beam distributions that are evidence of both parallel and perpendicular ion acceleration. Acceleration parallel to the magnetic field is most likely due to quasi-static electric fields (double layers) which can create outflowing ion beams; since ions of different mass will have different drift speeds due to this acceleration, a plasma configuration unstable to the ion-ion two-stream acoustic mode develops. When the net drift velocity (U) between the two ion species is greater than the sound speed (C(sub 0)), the ion-ion instability has maximum growth at oblique wave propagation. To study the nonlinear effects of the ion-ion instability in terms of plasma heating, a numerical simulation parametric study has been performed. It was found that the parallel acceleration that forms the ion beams occurs on a time scale faster than ion- ion wave growth at low drifts; thus ion-ion wave growth is expected to occur primarily for higher drift speeds (U greater than C(sub 0)) which results in strong oblique heating of the ions (both hydrogen and oxygen) forming elevated ion conics (sometimes called 'bowl' distributions). Also, strong parallel electron heating in the direction of the ion beams can occur, and electrons near the top of the acceleration region may attain a net upward drift along with the elevated ion conics. Variation of the oxygen density greatly affects the ion heating due to the ion-ion instability; as the oxygen density decreases, oxygen heating increases, in agreement with observations (Collin et al., 1987). Ion-ion electrostatic wave properties and the plasma heating that results over a wide range of auroral zone parameters are included.

Schriver, David↗
DOE OSTI2022

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Motion Planning Algorithms for Safety and Quantum Computing Efficiency

Motion planning remains a fundamental problem in robotics. Sampling-based algorithms use randomization to allow efficient solutions to this complex problem. As mobile robots and autonomous vehicles become more prevalent in everyday life, motion planning must be applied to increasingly challenging scenarios. Safety has become a paramount concern in motion planning for ensuring robotic applications enrich human lives. To date, many motion planning techniques to increase safety in the face of uncertain and dynamic environments have been developed. This dissertation first addresses distributional safety of Rapidly-Exploring Random Trees (RRT) through our algorithm W-Safe RRT. To acknowledge distributional uncertainty and poor modeling, W-Safe RRT uses the Wasserstein metric to provide a probabilistic bound on the distributional distance between a robot and obstacles. Human-interpretable environmental agent classification allows online safety margin adaptation. We propose and analyze an integrating region method for online classification that increases actor labeling accuracy based on behavioral feature values when compared to state of the art methods. The method performs class assignments based on local maximum likelihood in a created behavioral feature-space, allowing a notion of classification uncertainty. Model-based methods with safety guarantees can quickly become computationally in tractable, especially with multiple agents, higher dimensions, and plentiful unknowns. Sampling based algorithms have been parallelized for computation with multi-core computers and GPUs. We consider the use of quantum algorithms and computers for sampling-based motion planning for the first time. Quantum computing performs operations on superpositions of states and can solve certain problems much more efficiently than classical computers, but introduces previously unseen challenges. With Quantum-RRT, we recast the motion planning problem into a database-search structure and use Quantum Amplitude Amplification to find reachable states in the database with a quadratic performance increase over classical methods. We address two error sources with this method: quantum measurement and quantum oracle errors. We then extend this method to Parallel Quantum-RRT, which uses a manager-worker architecture with multiple parallel quantum workers to increase database search efficiency. We compare algorithm architectures and characterize probabilities of multiple workers finding solutions. Lastly, we test in simulation the quantum algorithms against classical versions in a wide variety of scenarios, concluding that a similar parallelization improvement is to be found in the quantum case as was found in the parallelization of classical RRT.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

A Multi-Level Parallelization Concept for High-Fidelity Multi-Block Solvers

The integration of high-fidelity Computational Fluid Dynamics (CFD) analysis tools with the industrial design process benefits greatly from the robust implementations that are transportable across a wide range of computer architectures. In the present work, a hybrid domain-decomposition and parallelization concept was developed and implemented into the widely-used NASA multi-block Computational Fluid Dynamics (CFD) packages implemented in ENSAERO and OVERFLOW. The new parallel solver concept, PENS (Parallel Euler Navier-Stokes Solver), employs both fine and coarse granularity in data partitioning as well as data coalescing to obtain the desired load-balance characteristics on the available computer platforms. This multi-level parallelism implementation itself introduces no changes to the numerical results, hence the original fidelity of the packages are identically preserved. The present implementation uses the Message Passing Interface (MPI) library for interprocessor message passing and memory accessing. By choosing an appropriate combination of the available partitioning and coalescing capabilities only during the execution stage, the PENS solver becomes adaptable to different computer architectures from shared-memory to distributed-memory platforms with varying degrees of parallelism. The PENS implementation on the IBM SP2 distributed memory environment at the NASA Ames Research Center obtains 85 percent scalable parallel performance using fine-grain partitioning of single-block CFD domains using up to 128 wide computational nodes. Multi-block CFD simulations of complete aircraft simulations achieve 75 percent perfect load-balanced executions using data coalescing and the two levels of parallelism. SGI PowerChallenge, SGI Origin 2000, and a cluster of workstations are the other platforms where the robustness of the implementation is tested. The performance behavior on the other computer platforms with a variety of realistic problems will be included as this on-going study progresses.

Hatay, Ferhat F.↗
NASA NTRS2018

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗
DOE OSTI2022

Towards Low-Overhead Resilience for Data Parallel Deep Learning

Data parallel techniques have been widely adopted both in academia and industry as a tool to enable scalable training of deep learning models. At scale, DL training jobs can fail due to software or hardware bugs, may need to be preempted or terminated due to unexpected events, or may perform suboptimally because they were misconfigured. Under such circumstances, there is a need to recover and/or reconfigure data-parallel DL training jobs on-the-fly, while minimizing the impact on the accuracy of the DNN model and the runtime overhead. In this regard, state-of-art techniques adopted by the HPC community mostly rely on checkpoint-restart, which inevitably leads to loss of progress, thus increasing the runtime overhead. In this paper we explore alternative techniques that exploit the properties of modern deep learning frameworks (overlapping of gradient averaging and weight updates with local gradient computations through pipeline parallelism) to reduce the overhead of resilience/elasticity. To this end we introduce a failure simulation framework and two resilience strategies (immediate mini-batch rollback and lossy forward recovery), which we study compared with checkpoint-restart approaches in a variety of settings in order to understand the trade-offs between the accuracy loss of the DNN model and the runtime overhead.

data-parallel training↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗