Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D 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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theory of Plasma-Cascade Instability

In this paper we present the theory of a novel micro-bunching instability occurring in charged particle beams propagating along a straight trajectory: a Plasma-Cascade Instability (PCI). This instability was confirmed by 3D numerical simulations and observed experimentally. It can by driven by variation of beam’s density and/or particle’s mobility. The PCI can strongly amplify noise in the beam and drastically reduce its quality. Conversely, such instability can drive novel high-power sources of radiation or can be used as a broadband amplifier.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Plasma Mirrors For High Power Lasers: A new approach for high repetition rates combined with realistic PIC simulations

We describe a research program to develop high quality plasma mirrors for use with high repetition rate, high power laser facilities to greatly facilitate experiments in high energy density physics, particularly when using solid density targets. A system delivering ultrathin, free standing liquid crystal films that could handle PW class laser pulses was developed that preserved the laser beam image quality and improved pulse contrast by over two orders of magnitude. The pulse contrast improvement was measured using a high dynamic range cross-correlation. The system operated at a repetition rate of several PMs formed in situ per minute. Extension to 1 Hz repetition rates should be possible. As a practical application, a PM system was used to separate laser accelerated electrons from the ultrashort pulse drive laser so that the electron beam emittance could be directly measured. The 20 nm thick PMs did not contribute significantly to emittance degradation. It was also demonstrated that the films could operate in the difficult plasma and laser environment present for GeV laser based electron acceleration. PM operation is still not well understood and much work with PMs is done empirically. Accordingly, we also demonstrated a particle-in-cell based simulation technique that would model PM operation, tracking the PM operation from cold, neutral atoms to a hot dense plasma by accounting for multiphoton ionization, plasma collisionality, and the dynamically changing index of refraction. Simulations in 3D were compared to a 2D based approach that was comparable. Finally, we developed a simpler model that can predict the intensity at which a plasma mirror will have its peak reflectivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Particle Acceleration during Magnetic Reconnection in Solar Flares and Its Parameter Dependence [Slide]

We performed 3D hybrid simulations (particle ions and fluid electrons) of magnetic reconnection and we have achieved, for the first time, nonthermal acceleration of all ion species (up to Iron) into power-law energy spectra. The spectra have a scaling of maximum energy per nucleon as a function of species consistent with satellite observations by NASA’s Parker Solar Probe mission.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Multi-fidelity modeling to predict the rheological properties of fiber suspensions

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications like biofuel production. The 3D numerical simulations of the suspension of fibers are often computationally expensive and time-consuming. Machine learning methods such as neural networks can simplify the prediction of rheological behavior; however, they require a relatively large training data set. Multi -fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified physical equations, can lead to optimized predictions. Here, we focus on a neural network with two levels of fidelity, i.e., high and low fidelity networks. To produce high-fidelity data, we perform direct numerical simulations to model the fibers as one-dimensional inextensible slender bodies that obey the Euler- Bernoulli beam equation. The Navier-Stokes equations govern the suspended fluid, and an immersed boundary method is used to couple the fluid and solid motion. The low-fidelity data is produced by using constitutive equations. Noticeable improvements have been observed in the accuracy of predicting the rheological behavior when a multi-fidelity network is used compared to the single-fidelity network.

Boodaghidizaji, Miad↗
DOE OSTI2024

Z-Target Radiography Postprocessing With A Deep Convolution Neural Network

Analyzing X-ray radiographs is crucial for understanding target behavior in Inertial Confinement Fusion (ICF) and High Energy Density (HED) platforms. However, the density of Magneto Raleigh Taylor (MRT) bands and limitations of target materials often obscure relevant spike growth and density information. To address this issue, machine learning postprocessing techniques can be applied to remove darkened regions in radiography images. In this study, a novel method is presented for removing MRT darkened regions from z-target radiographs using a convolutional neural network (CNN). The CNN, consisting of six layers, treats the darkened regions as noise and employs a mixed loss function and end-to-end frameworks to suppress them while preserving sharpness. The six-layer architecture is designed to effectively learn features when provided with a larger volume of learning space. Each layer is optimized using a mixed loss function that combines a standard loss pixel approach with a multi-scaled structural similarity index loss, which considers luminance, contrast, and structure in local neighborhoods. This approach is particularly beneficial for capturing the stochastic structure of MRT limbs. Due to the limited availability of experimental data, training is conducted using synthetic target radiography from 3D Alegra simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electric Drive Technologies Consortium (EDTC)/ Cost competitive, high-Performance, highly Reliable (CPR) Power Devices on 4H-SiC (Final Report)

4H-Silicon carbide (4H-SiC) is a wide bandgap semiconductor that offers superior material properties over silicon, including higher critical electric field, thermal conductivity, and electron saturation velocity. These advantages make 4H-SiC highly attractive for high-voltage, high-efficiency power electronics. However, realizing the full potential of SiC requires device technologies that are not only high-performing but also manufacturable and reliable under real-world operating conditions. This report summarizes the outcomes of a five-year R&D effort funded by the U.S. Department of Energy (DOE) under the Electric Drive Technologies Consortium (EDTC), focused on developing cost-competitive, high-performance, and highly reliable (CPR) power devices on 4H-SiC substrates. The program targeted scalable and manufacturable 1.2 kV-class SiC MOSFETs optimized for next-generation electric vehicles, renewable energy systems, and industrial power conversion. The project delivered transformative advancements in SiC power device performance and ruggedness. Particularly, Specific on-resistance (R on,sp ) was reduced by up to 37%, from ~4.0 m$\Omega \cdot$cm 2 in earlier designs to an industry-leading 2.40 m$\Omega \cdot$cm 2 , driven by optimized doping, refined JFET widths, and layout engineering. Breakdown voltages (BV) exceeded 1600 V, marking improvement over legacy baselines, and demonstrating the robustness of newly implemented junction profiles and edge terminations. Short-circuit withstand time (SCWT) saw a remarkable 4$\times$ increase, from ~2 $\mu$s to over 8 $\mu$s, achieved through the successful deployment of deep P-well structures (~1.8–2.0 $\mu$m) via channeling implantation. This innovative process breakthrough enabled precise junction formation without MeV-class implantation tools, reduced leakage under high field stress, and allowed even the shortest-channel devices (down to 0.3 $\mu$m) to achieve both high BV and excellent ruggedness—breaking the traditional trade-off between conduction efficiency and blocking capability. Several novel architectures pushed the performance envelope further. JBSFETs—featuring embedded Schottky portions—eliminated bipolar degradation and drastically reduced third-quadrant leakage, while Ladder MOSFETs introduced a clever orthogonal conduction path that achieved a 15.4% reduction in R on,sp over standard linear designs. Switching performance reached new benchmarks: short-channel devices showed a 31% reduction in total switching energy compared to 0.5 $\mu$m counterparts, while maintaining manageable gate drive requirements. Layout-optimized structures not only improved transconductance but also accelerated switching transitions, pointing to real-world benefits in converter-level efficiency. The devices also passed rigorous reliability validation. Stress-tested across TDDB, HTGB, HTRB, HVP, and burn-in, the devices screened under 30 V/10 hr and 43 V/1 s protocols consistently exhibited tighter lifetime distributions and long-term oxide robustness. These screening techniques proved effective in identifying latent defects and ensuring deployment-grade reliability. Meanwhile, advanced 3D TCAD simulations revealed and resolved electric field hotspots—particularly in HEXFET corners—where fields exceeding 4.8 MV/cm were mitigated through geometry-aware layout corrections. Overall, the results of this project demonstrate a manufacturable and scalable SiC power device platform that addresses key DOE performance targets for efficient, robust, and reliable 1.2kV 4H-SiC Power Devices. The developed technologies represent a meaningful step forward in the commercial readiness of high-voltage SiC solutions and provide a strong foundation for continued advancement in wide bandgap power electronics.

42 ENGINEERING↗
DOE OSTI2025

Validation of Photovoltaic Modeling Tool Solargraf Against Measured Data

Solargraf is a cloud-based 3D design tool by Enphase Energy that allows users to design solar and storage systems with a variety of elements. Through a Cooperative Research and Development Agreement (CRADA), Enphase Energy collaborated with the National Renewable Energy Laboratory (NREL) to validate Solargraf's 3D design simulation against measured PV system performance. This study follows the same methodology of similar validation studies completed at NREL. The predicted performance results from simulations in both Solargraf and NREL's System Advisor Model (SAM) tool were compared with measured data to evaluate performance predictions.

14 SOLAR ENERGY↗
DOE OSTI2019

Design of Comb Fabricated Halbach Undulators

An approach to fabricating Halbach array undulators using “combs” machined from single magnets is introduced. This technique is especially relevant to the fabrication of short period micro-undulators with period lengths considerably less than the few-centimeter-scale typical of current undulators. Manual, magnet-by-magnet assembly of micro-undulators would require the manipulation and alignment of thousands of magnets smaller than a grain of rice: comb fabrication dramatically increases the size of the basic unit cell of assembly with no increase in undulator period by creating many periods from a single piece, in a single machining modality. Further, as these comb teeth are intrinsically indexed to each other, tolerances are dictated by a single manufacturing step rather than accumulating errors by assembling many tiny magnets relative to each other. Different Halbach geometries, including M ′ = 2 , M ′ = 4 , isosceles triangle, and hybrid, are examined both from a theoretical perspective and with 3D magnetostatic simulations.

Majernik, Nathan↗
DOE OSTI2016

AN MHD AVALANCHE IN A MULTI-THREADED CORONAL LOOP

For the first time, we demonstrate how an MHD avalanche might occur in a multithreaded coronal loop. Considering 23 non-potential magnetic threads within a loop, we use 3D MHD simulations to show that only one thread needs to be unstable in order to start an avalanche even when the others are below marginal stability. This has significant implications for coronal heating in that it provides for energy dissipation with a trigger mechanism. The instability of the unstable thread follows the evolution determined in many earlier investigations. However, once one stable thread is disrupted, it coalesces with a neighboring thread and this process disrupts other nearby threads. Coalescence with these disrupted threads then occurs leading to the disruption of yet more threads as the avalanche develops. Magnetic energy is released in discrete bursts as the surrounding stable threads are disrupted. The volume integrated heating, as a function of time, shows short spikes suggesting that the temporal form of the heating is more like that of nanoflares than of constant heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

INTERPLAY AMONG COOLING, AGN FEEDBACK, AND ANISOTROPIC CONDUCTION IN THE COOL CORES OF GALAXY CLUSTERS

Feedback from the active galactic nuclei (AGNs) is one of the most promising heating mechanisms to circumvent the cooling-flow problem in galaxy clusters. However, the role of thermal conduction remains unclear. Previous studies have shown that anisotropic thermal conduction in cluster cool cores (CCs) could drive the heat-flux-driven buoyancy instabilities (HBIs) that reorient the field lines in the azimuthal directions and isolate the cores from conductive heating from the outskirts. However, how the AGN interacts with the HBI is still unknown. To understand these interwined processes, we perform the first 3D magnetohydrodynamic simulations of isolated CC clusters that include anisotropic conduction, radiative cooling, and AGN feedback. We find the following: (1) For realistic magnetic field strengths in clusters, magnetic tension can suppress a significant portion of HBI-unstable modes, and thus the HBI is either completely inhibited or significantly impaired, depending on the unknown magnetic field coherence length. (2) Turbulence driven by AGN jets can effectively randomize magnetic field lines and sustain conductivity at ∼1/3 of the Spitzer value; however, the AGN-driven turbulence is not volume filling. (3) Conductive heating within the cores could contribute to ∼10% of the radiative losses in Perseus-like clusters and up to ∼50% for clusters twice the mass of Perseus. (4) Thermal conduction has various impacts on the AGN activity and intracluster medium properties for the hottest clusters, which may be searched by future observations to constrain the level of conductivity in clusters. The distribution of cold gas and the implications are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MAGNETIC FIELDS IN EARLY PROTOSTELLAR DISK FORMATION

We consider formation of accretion disks from a realistically turbulent molecular gas using 3D MHD simulations. In particular, we analyze the effect of the fast turbulent reconnection described by the Lazarian and Vishniac model for the removal of magnetic flux from a disk. With our numerical simulations we demonstrate how the fast reconnection enables protostellar disk formation resolving the so-called “magnetic braking catastrophe.” In particular, we provide a detailed study of the dynamics of a 0.5 M{sub ⊙} protostar and the formation of its disk for up to several thousands years. We measure the evolution of the mass, angular momentum, magnetic field, and turbulence around the star. We consider effects of two processes that strongly affect the magnetic transfer of angular momentum, both of which are based on turbulent reconnection: the first, “reconnection diffusion,” removes the magnetic flux from the disk; the other involves the change of the magnetic field's topology, but does not change the absolute value of the magnetic flux through the disk. We demonstrate that for the first mechanism, turbulence causes a magnetic flux transport outward from the inner disk to the ambient medium, thus decreasing the coupling of the disk to the ambient material. A similar effect is achieved through the change of the magnetic field's topology from a split monopole configuration to a dipole configuration. We explore how both mechanisms prevent the catastrophic loss of disk angular momentum and compare both above turbulent reconnection mechanisms with alternative mechanisms from the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A CIRCULAR-CYLINDRICAL FLUX-ROPE ANALYTICAL MODEL FOR MAGNETIC CLOUDS

We present an analytical model to describe magnetic flux-rope topologies. When these structures are observed embedded in Interplanetary Coronal Mass Ejections (ICMEs) with a depressed proton temperature, they are called Magnetic Clouds (MCs). Our model extends the circular-cylindrical concept of Hidalgo et al. by introducing a general form for the radial dependence of the current density. This generalization provides information on the force distribution inside the flux rope in addition to the usual parameters of MC geometrical information and orientation. The generalized model provides flexibility for implementation in 3D MHD simulations. Here, we evaluate its performance in the reconstruction of MCs in in situ observations. Four Earth-directed ICME events, observed by the Wind spacecraft, are used to validate the technique. The events are selected from the ICME Wind list with the magnetic obstacle boundaries chosen consistently with the magnetic field and plasma in situ observations and with a new parameter (EPP, the Electron Pitch angle distribution Parameter) which quantifies the bidirectionally of the plasma electrons. The goodness of the fit is evaluated with a single correlation parameter to enable comparative analysis of the events. In general, at first glance, the model fits the selected events very well. However, a detailed analysis of events with signatures of significant compression indicates the need to explore geometries other than the circular-cylindrical. An extension of our current modeling framework to account for such non-circular CMEs will be presented in a forthcoming publication.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MODELING THE INITIATION OF THE 2006 DECEMBER 13 CORONAL MASS EJECTION IN AR 10930: THE STRUCTURE AND DYNAMICS OF THE ERUPTING FLUX ROPE

We carry out a 3D magnetohydrodynamic simulation to model the initiation of the coronal mass ejection (CME) on 2006 December 13 in the emerging δ -sunspot active region NOAA 10930. The setup of the simulation is similar to a previous simulation by Fan, but with a significantly widened simulation domain to accommodate the wide CME. The simulation shows that the CME can result from the emergence of a east–west oriented twisted flux rope whose positive, following emerging pole corresponds to the observed positive rotating sunspot emerging against the southern edge of the dominant pre-existing negative sunspot. The erupting flux rope in the simulation accelerates to a terminal speed that exceeds 1500 km s{sup −1} and undergoes a counter-clockwise rotation of nearly 180° such that its front and flanks all exhibit southward directed magnetic fields, explaining the observed southward magnetic field in the magnetic cloud impacting the Earth. With continued driving of flux emergence, the source region coronal magnetic field also shows the reformation of a coronal flux rope underlying the flare current sheet of the erupting flux rope, ready for a second eruption. This may explain the build up for another X-class eruptive flare that occurred the following day from the same region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

The magnetic furnace: intense core dynamos in B stars

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through three-dimensional (3D) global simulations of convective core dynamos operating within a young 10 M{sub ⊙} B-type star, using the anelastic spherical harmonic code. These simulations capture the inner 65% of this star by radius, encompassing the convective nuclear-burning core (about 23% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05% to 16% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

( Chandra ) view of magnetically confined wind in HD 191612: theory versus observations

High-resolution spectra of the magnetic star HD 191612 were acquired using the Chandra X-ray Observatory at both maximum and minimum emission phases. We confirm the flux and hardness variations previously reported with XMM-Newton, demonstrating the high repeatability of the behavior of HD 191612 over a decade. The line profiles appear typical for magnetic massive stars: no significant line shift, relatively narrow lines for high-Z elements, and formation radius at about 2 R∗. Line ratios confirm the softening of the X-ray spectrum at the minimum emission phase. Shift or width variations appear of limited amplitude at most (slightly lower velocity and slightly increased broadening at minimum emission phase, but within 1–2σ of values at maximum). In addition, a fully self-consistent 3D magnetohydrodynamic simulation of the confined wind in HD 191612 was performed. The simulation results were directly fitted to the data, leading to a remarkable agreement overall between them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Structure, dynamics, and deuterium fractionation of massive pre-stellar cores

High levels of deuterium fraction in N{sub 2}H{sup +} are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values (D{sub frac}{sup N{sub 2}H{sup +}}≡N{sub 2}D{sup +}/N{sub 2}H{sup +}≳0.1) is longer than the free-fall time, possibly 10 times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N{sub 2}H{sup +} and N{sub 2}D{sup +}. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N{sub 2}D{sup +} appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee and Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H{sub 2}, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios (≲0.01) and/or multiple free-fall times (≳3) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗