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,225 records · Page 68

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗

Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak

The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψ N ) or large (>4% ψ N ) island width depending on the magnitude of the edge bootstrap current, where ψ N is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χ | /χ ⊥ ) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗

Life cycle of streaks in the buffer layer of wall-bounded turbulence

In this work, streaks in the buffer layer of wall-bounded turbulence are tracked in time to study their life cycle. Spatially and temporally resolved direct numerical simulation data are used to analyze the strong wall-parallel movements conditioned to low-speed streamwise flow. The analysis of the streaks shows that there is a clear distinction between wall-attached and detached streaks, and that the wall-attached streaks can be further categorized into streaks that are contained in the buffer layer and the ones that reach the outer region. The results reveal that streaks are born in the buffer layer, coalescing with each other to create larger streaks that are still attached to the wall. Once the streak becomes large enough, it starts to meander due to the large streamwise-to-wall-normal aspect ratio, and consequently the elongation in the streamwise direction, which makes it more difficult for the streak to be oriented strictly in the streamwise direction. While the continuous interaction of the streaks allows the superstructure to span extremely long temporal and length scales, individual streak components are relatively small and short-lived. Tall-attached streaks eventually split into wall-attached and wall-detached components. These wall-detached streaks have a strong wall-normal velocity away from the wall, similar to ejections or bursts observed in the literature. Conditionally averaging the flow fields to these split events show that the detached streak has not only a larger wall-normal velocity compared to the wall-attached counterpart, it also has a larger (less negative) streamwise velocity, similar to the velocity field at the tip of a vortex cluster.

42 ENGINEERING↗

A Decoupled Droop Control Strategy for Cascaded Multicell Inverter with Low-Frequency Modulation

Cascaded multicell inverter (CMI), featured by its merits of modularity and fault-tolerance, is suitable for medium-voltage applications without using bulky step-up transformers. However, little research has reported grid-forming control of the CMI for the application in the power grid with high penetrations of inverter-interfaced generation. To bridge this research gap, this paper proposes a novel droop control strategy for the CMI with a low-frequency modulation to offer inertia and reactive power support to the power grid. A decoupled Q/V droop control strategy is proposed to eliminate the coupling effect from the P/ω control loop caused by the phase-angle difference between the voltages of the CMI and the grid. The output of the Q/V control loop is set as the reference voltage for each cell of the CMI, which makes the CMI naturally a voltage source to provide a grid-forming capability. Finally, the proposed control strategy shows great active and reactive power regulation and sharing capability, which is validated by simulation on two 100 kW, 2.4 kV CMIs connected in parallel to a power grid.

14 SOLAR ENERGY↗

PARMOO

ParMOO is a Python library for solving multiobjective simulation optimization problems, while exploiting problem structure. ParMOO stands for "parallel multiobjective optimization".

WILD, STEFAN↗

Detector and Beamline Simulation for Next-Generation High Energy Physics Experiments

The success of high energy physics programs relies heavily on accurate detector simulations and beam interaction modeling. The increasingly complex detector geometries and beam dynamics require sophisticated techniques in order to meet the demands of current and future experiments. Common software tools used today are unable to fully utilize modern computational resources, while data-recording rates are often orders of magnitude larger than what can be produced via simulation. In this paper, we describe the state, current and future needs of high energy physics detector and beamline simulations and related challenges, and we propose a number of possible ways to address them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A steering law for three double-gimbal control moment gyro systems

A steering law for three double-gimbal control moment gyro (CMG) systems is proposed. This steering law is applicable to systems with almost any configuration of CMGs and the CMG-out operation needs no special modification. Examples of three double gimbal CMG systems in an orthogonal configuration and in a parallel configuration are shown along with the results of digital simulations. Simulation results show that any command torque can always be met except when the system is in a singular state and that, whenever the system is in, or close to, a singularity, the steering law drives the system state out of the vicinity of the singularity.

Yoshikawa, T.↗

The ASSESS program documenting and evaluating low cost techniques for conducting research in the flight environment

The Airborne Science Office (ASO) of the Ames Research Center has for 10 years operated an airborne scientific research program in infrared astronomy and other disciplines. The Lear Jet, CV-990, and C-141 flying laboratories are flown by ASO crews, while the major responsibility for defining, developing, and operating the experimental equipment is placed on individual researchers, who have included scientists from many countries. The ASSESS (Airborne Science/Shuttle Experiment Systems Simulation) program consists of two phases: Phase A documents the present management and operational practices of the ASO, and Phase B consists of airborne research missions constrained (for example, by crew confinement) to simulate certain aspects of experimental operations on Shuttle/Spacelab missions. Various parallels between the Airborne Science Program and Spacelab are pointed out and their applications to Spacelab planning are discussed.

Weaver, L. B.↗

Shock acceleration of diffuse ions at the earth's bow shock Acceleration efficiency and A/Z enhancement

Observations of particle spectra, intensity, and enhancement of alpha particles over protons at diffuse ion events at the quasi-parallel earth bow shock are compared to a Monte Carlo simulation of diffusive shock acceleration. The simulation includes the back reaction of accelerated particles on the shock structure, particle escape at an upstream free escape boundary, and a low energy per nucleon threshold for thermal leakage of downstream, shock-heated particles into the upstream region. The simulation assumes that the same scattering operator that gives rise to shock acceleration can also describe a viscous shock governed by hydrodynamic turbulence. This implies that accelerated ions can be drawn directly from the thermal solar wind with no separate superthermal seed population. Good agreement between the simulation and observations made during nearly radial magnetic field configurations lends support to thermal leakage of downstream, shock-heated ions as the mode of injection for diffusion ion events.

Ellison, D. C.↗

Multi-stage design of an optimal momentum management controller for the Space Station

This paper presents a multistage design scheme for determining an optimal control-moment-gyro momentum-management and attitude-control system for the Space Station Freedom. The Space Station equations of motion are linearized and block-decomposed into two block-decoupled subsystems using the matrix-sign algorithm. A sequential procedure is utilized for designing a linear-quadratic regulator for each subsystem, which optimally places the eigenvalues of the closed-loop subsystem in the region of an open sector, bounded by lines inclined at + or - pi/2k (for k = 2 or 3) from the negative real axis, and the left-hand side of a line parallel to the imaginary axis in the s-plane. Simulation results are presented to compare the resultant designs.

Sunkel, J. W.↗

Order (n) DISCOS for multibody systems with gear reduction

Recent developments in O(n) algorithms (where n is the number of bodies in the system) and parallel processing have drastically reduced the computer time needed to simulate systems involving many bodies. This paper presents a gear-reduction model for the O(n) version of DISCOS, a standard software package for simulation and analysis of flexible multibody systems. The gear-reduction model allows the accurate modeling of harmonic drives, which are commonly used in robot joints. The formulation has been implemented and validated with known results. The gear model can also be used for gear-train, rack-and-pinion, and screw joints.

Chun, Hon M.↗

The energetic ion signature of an O-type neutral line in the geomagnetic tail

An energetic ion signature is presented which has the potential for remote sensing of an O-type neutral line embedded in a current sheet. A source plasma with a tailward flowing Kappa distribution yields a strongly non-Kappa distribution after interacting with the neutral line: sharp jumps, or ridges, occur in the velocity space distribution function f(nu-perpendicular, nu-parallel) associated with both increases and decreases in f. The jumps occur when orbits are reversed in the x-direction: a reversal causing initially earthward particles (low probability in the source distribution) to be observed results in a decrease in f, while a reversal causing initially tailward particles to be observed produces an increase in f. The reversals, and hence the jumps, occur at approximately constant values of perpendicular velocity in both the positive nu parallel and negative nu parallel half planes. The results were obtained using single particle simulations in a fixed magnetic field model.

Martin, R. F., Jr.↗

Impact of Load Balancing on Unstructured Adaptive Grid Computations for Distributed-Memory Multiprocessors

The computational requirements for an adaptive solution of unsteady problems change as the simulation progresses. This causes workload imbalance among processors on a parallel machine which, in turn, requires significant data movement at runtime. We present a new dynamic load-balancing framework, called JOVE, that balances the workload across all processors with a global view. Whenever the computational mesh is adapted, JOVE is activated to eliminate the load imbalance. JOVE has been implemented on an IBM SP2 distributed-memory machine in MPI for portability. Experimental results for two model meshes demonstrate that mesh adaption with load balancing gives more than a sixfold improvement over one without load balancing. We also show that JOVE gives a 24-fold speedup on 64 processors compared to sequential execution.

Biswas, Rupak↗

Magnetic Fluxtube Tunneling

We present numerical simulations of the collision and subsequent interaction of two initially orthogonal, twisted, force free field magnetic fluxtubes. The simulations were carried out using a new three dimensional explicit parallelized Fourier collocation algorithm for solving the viscoresistive equations of compressible magnetohydrodynamics. It is found that, under a wide range of conditions, the fluxtubes can 'tunnel' through each other. Two key conditions must be satisfied for tunneling to occur: the magnetic field must be highly twisted with a field line pitch much greater than 1, and the magnetic Lundquist number must be somewhat large, greater than or equal to 2880. This tunneling behavior has not been seen previously in studies of either vortex tube or magnetic fluxtube interactions. An examination of magnetic field lines shows that tunneling is due to a double reconnection mechanism. Initially orthogonal field lines reconnect at two specific locations, exchange interacting sections and 'pass' through each other. The implications of these results for solar and space plasmas are discussed.

Dahlburg, Russell B.↗

Impact of Load Balancing on Unstructured Adaptive Grid Computations for Distributed-Memory Multiprocessors

The computational requirements for an adaptive solution of unsteady problems change as the simulation progresses. This causes workload imbalance among processors on a parallel machine which, in turn, requires significant data movement at runtime. We present a new dynamic load-balancing framework, called JOVE, that balances the workload across all processors with a global view. Whenever the computational mesh is adapted, JOVE is activated to eliminate the load imbalance. JOVE has been implemented on an IBM SP2 distributed-memory machine in MPI for portability. Experimental results for two model meshes demonstrate that mesh adaption with load balancing gives more than a sixfold improvement over one without load balancing. We also show that JOVE gives a 24-fold speedup on 64 processors compared to sequential execution.

Sohn, Andrew↗