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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗
DOE OSTI2025

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗
NASA NTRS2000

A Shell/3D Modeling Technique for the Analysis of Delaminated Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with shell finite elements. Double Cantilever Beam, End Notched Flexure, and Single Leg Bending specimens were analyzed first using full 3D finite element models to obtain reference solutions. Mixed mode strain energy release rate distributions were computed using the virtual crack closure technique. The analyses were repeated using the shell/3D technique to study the feasibility for pure mode I, mode II and mixed mode I/II cases. Specimens with a unidirectional layup and with a multidirectional layup were simulated. For a local 3D model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗
DOE OSTI2022

Development of a synthetic phase contrast imaging diagnostic for turbulence studies at Wendelstein 7-X

Abstract We present a synthetic phase contrast imaging (PCI) diagnostic for studying turbulence at the Wendelstein 7-X (W7-X) stellarator. We first describe the implemented instrument response model, which captures diffraction effects, detector noise, and the long-wavelength cutoff due to the phase plate of the PCI system. To verify the instrument response model, we show that it is capable of reproducing the PCI signal generated by the sound wave speaker used for calibration at W7-X. Next, we discuss the calculation of synthetic PCI signals based on the global, nonlinear gyrokinetic codes GENE-3D and EUTERPE, including results from some of the first stellarator simulations of this type with kinetic electrons (KEs) in GENE-3D. While the simulations used in this work lack a neoclassical radial electric field, which is crucial for reproducing experimental PCI signals, they do indicate that the dominant rotation direction and velocities of the turbulent fluctuations can be inferred from the wave number-frequency spectra of the PCI signals, as expected. The synthetic PCI wave number spectra are further shown to be similar to those of the line-integrated fluctuating electron density, with distinct differences between adiabatic and KE simulations, explainable by previously published turbulence models. For example, the wave number spectra of all adiabatic electron simulations analyzed here follow a power law with an exponent close to −5 for sufficiently large wave numbers. This indicates that universal features of electron density turbulence at W7-X may be studied using the PCI system.

EUTERPE↗
DOE OSTI2020

Development of a predictive capability of short-pulse laser-driven broadband x-ray radiography

High intensity, short-pulse laser interaction with a solid metal target produces broadband hard x-rays potentially for various applications of x-ray radiography. In this work, experimental benchmarking of numerical modelling for short-pulse laser-driven broadband x-ray radiography is presented. Angular dependent x-ray spectra are first calculated with a hybrid particle-in-cell code, Large Scale Plasma (LSP), using fast electron parameters inferred from an analysis of measured bremsstrahlung signals. Subsequently, a calculated x-ray spectrum in the direction of radiography is used in photon transport calculations using a Monte Carlo code, Particle and Heavy Ion Transport code System (PHITS), to simulate a radiographic image including a modelled 3D test object, an x-ray attenuation filter and an image plate detector. Simulated radiographic images are compared with measurements obtained in an experiment using a 50-TW Leopard short-pulse laser at the University of Nevada Reno. Results show that simulations reproduce the experimental images well for three different attenuation filters (plastic, aluminium, and brass), while 1D transmission profiles for the plastic and aluminium filters are quantitatively in good agreement. The modelling approach established in this work could be used as a predictive tool to simulate radiographic images of complex 3D solid objects at any arbitrary angular position or to optimize experimental components such as the source spectrum, x-ray attenuation filters and a detector type depending on a radiographic object without carrying out radiographic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Verification of High-Fidelity Neutronics Code PROTEUS for C5G7 Benchmark Problems

Under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the 3D whole core transport code PROTEUS has been developed at Argonne National Laboratory to perform the detailed 3D transport calculation for various reactor types including irregular or complex geometric reactors. PROTEUS includes two distinct transport solvers: SN and MOC. Among them, the MOC solver based on the extruded geometry, thus named MOCEX, can produce the faithful 3-D transport solutions by employing the rigorous and robust 2D/3D MOC formulation. This formulation overcomes the drawbacks of a full 3D MOC method requiring enormous memory and computational resources and resolves the instability issues of the 2D/1D method as well. The solution of MOCEX converges to that of a Monte Carlo code as the discretization of angle, space and energy is refined, unlike the 2D/1D method. As a verification and validation (V and V) effort of PROTEUS-MOCEX, this paper presents verification test results using the well-known 3D C5G7 benchmark problems by solving and comparing eigenvalues and pin powers with MCNP solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2013

Damage Detection in Composite Structures with Wavenumber Array Data Processing

Guided ultrasonic waves (GUW) have the potential to be an efficient and cost-effective method for rapid damage detection and quantification of large structures. Attractive features include sensitivity to a variety of damage types and the capability of traveling relatively long distances. They have proven to be an efficient approach for crack detection and localization in isotropic materials. However, techniques must be pushed beyond isotropic materials in order to be valid for composite aircraft components. This paper presents our study on GUW propagation and interaction with delamination damage in composite structures using wavenumber array data processing, together with advanced wave propagation simulations. Parallel elastodynamic finite integration technique (EFIT) is used for the example simulations. Multi-dimensional Fourier transform is used to convert time-space wavefield data into frequency-wavenumber domain. Wave propagation in the wavenumber-frequency domain shows clear distinction among the guided wave modes that are present. This allows for extracting a guided wave mode through filtering and reconstruction techniques. Presence of delamination causes spectral change accordingly. Results from 3D CFRP guided wave simulations with delamination damage in flat-plate specimens are used for wave interaction with structural defect study.

Tian, Zhenhua↗
NASA NTRS2010

Monte-Carlo Simulation of Radiation Track Structure and Calculation of Dose Deposition in Nanovolumes

INTRODUCTION: The radiation track structure is of crucial importance to understand radiation damage to molecules and subsequent biological effects. Of a particular importance in radiobiology is the induction of double-strand breaks (DSBs) by ionizing radiation, which are caused by clusters of lesions in DNA, and oxidative damage to cellular constituents leading to aberrant signaling cascades. DSB can be visualized within cell nuclei with gamma-H2AX experiments. MATERIAL AND METHODS: In DSB induction models, the DSB probability is usually calculated by the local dose obtained from a radial dose profile of HZE tracks. In this work, the local dose imparted by HZE ions is calculated directly from the 3D Monte-Carlo simulation code RITRACKS. A cubic volume of 5 micron edge (Figure 1) is irradiated by a (Fe26+)-56 ion of 1 GeV/amu (LET approx.150 keV/micron) and by a fluence of 450 H+ ions, 300 MeV/amu (LET approx. 0.3 keV/micron). In both cases, the dose deposited in the volume is approx.1 Gy. The dose is then calculated into each 3D pixels (voxels) of 20 nm edge and visualized in 3D. RESULTS AND DISCUSSION: The dose is deposited uniformly in the volume by the H+ ions. The voxels which receive a high dose (orange) corresponds to electron track ends. The dose is deposited differently by the 56Fe26+ ion. Very high dose (red) is deposited in voxels with direct ion traversal. Voxels with electron track ends (orange) are also found distributed around the path of the track. In both cases, the appearance of the dose distribution looks very similar to DSBs seen in gammaH2AX experiments, particularly when the visualization threshold is applied. CONCLUSION: The refinement of the dose calculation to the nanometer scale has revealed important differences in the energy deposition between high- and low-LET ions. Voxels of very high dose are only found in the path of high-LET ions. Interestingly, experiments have shown that DSB induced by high-LET radiation are more difficult to repair. Therefore, this new approach may be useful to understand the nature of DSB and oxidative damage induced by ionizing radiation.

Plante, I.↗
DOE OSTI2025

3D Printing of Cement-Based Materials Using Seawater for Simulated Marine Environments

Global demand for adaptable and rapidly deployable construction solutions in offshore, coastal, and fluvial environments continues to rise, driven by pressing needs to develop energy platforms, improve coastal resilience, and support emergency response in the face of natural disasters. Increased investment in human-made coastal infrastructure, such as piers, support structures for power lines, offshore wind farms, and seawall protection systems, further underscores this trend. This study investigates the development of printable concrete mixtures for underwater environments using seawater as a replacement for freshwater, using a 3D printing syringe-based extrusion system. The effect of seawater addition and the printing medium (in air vs. underwater) was assessed via rheological and mechanical performance characterization. The results indicate rheological properties are favorable for seawater adoption by producing mixtures with higher yield stress and viscosity with the same levels of admixtures used for freshwater. Seawater-based mixtures demonstrated superior dimensional stability compared to freshwater counterparts, maintaining cross-sectional geometry, while compressive strength results showed no statistical differences between in-air and underwater samples. However, flexural strength was significantly influenced by geometry and printing medium. These findings establish critical rheological parameters for printable underwater mixtures and highlight the need for optimized curing strategies and layer bonding techniques to improve interfacial strength in underwater 3D printing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The return of the merging galaxy subclusters of El Gordo?

Merging galaxy clusters with radio relics provide rare insights to the merger dynamics as the relics are created by the violent merger process. We demonstrate one of the first uses of the properties of the radio relic to reduce the uncertainties of the dynamical variables and determine the three-dimensional (3D) configuration of a cluster merger, ACT-CL J0102-4915, nicknamed El Gordo. From the double radio relic observation and the X-ray observation of a comet-like gas morphology induced by motion of the cool core, it is widely believed that El Gordo is observed shortly after the first core passage of the subclusters. Here, we employ a Monte Carlo simulation to investigate the 3D configuration and dynamics of El Gordo. Using the polarization fraction of the radio relic, we constrain the estimate of the angle between the plane of the sky and the merger axis to be α=21°± $9\atop{11}$. We find the relative 3D merger speed of El Gordo to be 2400 ± $400\atop{200}$ km s -1 at pericentre. The two possible estimates of the time since pericentre (TSP) are 0.46 ± $0.09\atop{0.16}$ and 0.91± $0.22\atop{0.39}$ Gyr for the outgoing and returning scenario, respectively. We put our estimates of the TSP into context by showing that if the time-averaged shock velocity is approximately equal to or smaller than the pericentre velocity of the corresponding subcluster in the centre-of-mass frame, the two subclusters are more likely to be moving towards, rather than away, from each other, post apocentre. Finally, we compare and contrast the merger scenario of El Gordo with that of the Bullet Cluster, and show that this late-stage merging scenario explains why the south-east (SE) dark matter lensing peak of El Gordo is closer to the merger centre than the SE cool core.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Multiscale Failure Analysis of a 3D Woven Unit Cell Containing Defects

Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.

3D woven↗
DOE OSTI2015

Type I planet migration in a magnetized disk. I. Effect of large-scale vertical and azimuthal field components

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a combination of numerical simulations in 2D and 3D and a linear perturbation analysis. Steady-state models of such disks require the inclusion of magnetic diffusivity. To make progress using ideal MHD, we focus on simplified field configurations, involving purely vertical (B{sub z}) and azimuthal (B{sub φ}) field components and a combination of the two. For each of the models we calculate the locations of the relevant resonances and of the turning points, which delineate the propagation regions of the MHD waves that transport angular momentum from the planet to the disk. We use both numerical and semianalytic methods to evaluate the cumulative back torque acting on the planet, and explore the effect of spatial gradients in the disk’s physical variables on the results. We conclude that, under realistic (3D) circumstances, a large-scale magnetic field can slow down the inward migration that characterizes the underlying unmagnetized disk—by up to a factor of ∼2 when the magnetic pressure approaches the thermal pressure—but it cannot reverse it. A previous inference that a pure-B{sub ϕ} field whose amplitude decreases fast enough with radius leads to outward migration applies only in 2D. In fact, we find that, in 3D, a pure-B{sub ϕ} disk undergoes a rapid transition to turbulence on account of a magnetorotational instability that is triggered by the planet-induced appearance of a weak B{sub z} component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Application and Evaluation of Surrogate Models for Radiation Source Search

Surrogate models are increasingly required for applications in which first-principles simulation models are prohibitively expensive to employ for uncertainty analysis, design, or control. They can also be used to approximate models whose discontinuous derivatives preclude the use of gradient-based optimization or data assimilation algorithms. We consider the problem of inferring the 2D location and intensity of a radiation source in an urban environment using a ray-tracing model based on Boltzmann transport theory. Whereas the code implementing this model is relatively efficient, extension to 3D Monte Carlo transport simulations precludes subsequent Bayesian inference to infer source locations, which typically requires thousands to millions of simulations. Additionally, the resulting likelihood exhibits discontinuous derivatives due to the presence of buildings. To address these issues, we discuss the construction of surrogate models for optimization, Bayesian inference, and uncertainty propagation. Specifically, we consider surrogate models based on Legendre polynomials, multivariate adaptive regression splines, radial basis functions, Gaussian processes, and neural networks. We detail strategies for computing training points and discuss the merits and deficits of each method.

Cook, Jared A.↗
DOE OSTI2012

Non-linear magnetohydrodynamic simulations of density evolution in Tore Supra sawtoothing plasmas

The plasma density evolution in sawtooth regime on the Tore Supra tokamak is analyzed. The density is measured using fast-sweeping X-mode reflectometry which allows tomographic reconstructions. There is evidence that density is governed by the perpendicular electric flows, while temperature evolution is dominated by parallel diffusion. Postcursor oscillations sometimes lead to the formation of a density plateau, which is explained in terms of convection cells associated with the kink mode. A crescent-shaped density structure located inside q = 1 is often visible just after the crash and indicates that some part of the density withstands the crash. 3D full MHD nonlinear simulations with the code XTOR-2F recover this structure and show that it arises from the perpendicular flows emerging from the reconnection layer. The proportion of density reinjected inside the q = 1 surface is determined, and the implications in terms of helium ash transport are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Attosecond Probing of Coherent Vibrational Dynamics in CBr 4

A coherent vibrational wavepacket is launched and manipulated in the symmetric stretch (a1) mode of CBr 4 , by impulsive stimulated Raman scattering (ISRS) from nonresonant 400 nm laser pump pulses with various peak intensities on the order of tens of 10 12 W/cm 2 . Extreme ultraviolet (XUV) attosecond transient absorption spectroscopy (ATAS) records the wavepacket dynamics as temporal oscillations in XUV absorption energy at the bromine M 4,5 3d 3/2,5/2 edges around 70 eV. The results are augmented by nuclear timedependent Schrödinger equation simulations. Slopes of the (Br 3d 3/2,5/2 ) –1 10a 1 * core-excited state potential energy surface (PES) along the a 1 mode are calculated to be –9.4 eV/Å from restricted open-shell Kohn–Sham calculations. Using analytical relations derived for the small-displacement limit and the calculated slopes of the core-excited state PES, a deeper insight into the vibrational dynamics is obtained by retrieving the experimental excursion amplitude of the vibrational wavepacket and the amount of population transferred to the vibrational first-excited state as a function of pump-pulse peak intensity. Experimentally, the results show that XUV ATAS is capable of resolving oscillations in the XUV absorption energy on the order of a few to tens of meV with tens of femtosecond time precision. This corresponds to change in C–Br bond length on the order of 10 –4 to 10 –3 Å. The results and the analytic relationships offer a clear physical picture, on multiple levels of understanding, of how the pump-pulse peak intensity controls the vibrational dynamics launched by nonresonant ISRS in the small-displacement limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗
DOE OSTI2019

Numerical Simulation and Validation of Gas and Molten Metal Flows in Close-Coupled Gas Atomization

Molten metal atomization in close-coupled gas atomization dies can operate between two limiting conditions, jetting and filming, together with several complex mechanisms: liquid-gas drafting, downward/upward shearing, melt bouncing, etc. Liquid jet deformation depends on flow and geometric parameters, such as liquid Reynolds, liquid Weber, and gas Mach numbers, as well as gas jet apex angle and melt tube tip extension and aspect ratio, among others. Understanding their effect is of importance for the metal powder making industry. Numerical gas atomization studies can provide approximated flow information and consider a wide range of conditions, beyond experimental reach. Here, 3D high-resolution simulations employing a 5-equation compressible flow model coupled with the volume-of-fluid method are compared with experiments, for liquid Weber number in the range of 1–30 and liquid Reynolds number below 10,000. This validation explores the predicting capabilities of the numerical model.

Hernandez, F.↗
DOE OSTI2024

Improved Three-Dimensional Reconstructions in Electron Ptychography through Defocus Series Measurements

A detailed analysis of ptychography for three-dimensional (3D) phase reconstructions of thick specimens is performed. Here, we introduce multi-focus ptychography, which incorporates a 4D-STEM defocus series to enhance the quality of 3D reconstructions along the beam direction through a higher overdetermination ratio. This method is compared with established multi-slice ptychography techniques, such as conventional ptychography, regularized ptychography, and multi-mode ptychography. Additionally, we contrast multi-focus ptychography with an alternative method that uses virtual optical sectioning through a reconstructed scattering matrix (S-matrix), which offers more precise 3D structure information compared to conventional ptychography. Our findings from multiple 3D reconstructions based on simulated and experimental data demonstrate that multi-focus ptychography surpasses other techniques, particularly in accurately reconstructing the surfaces and interface regions of thick specimens.

3D reconstruction↗