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,207 records · Page 67

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components

Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.

47 OTHER INSTRUMENTATION↗

Quantifying and Accounting for the Importance of Aerosol Mixing State. Final Report

The objective of this project was to develop a comprehensive framework of modeling and theory to rigorously investigate the impacts of aerosol mixing state on climate-relevant aerosol properties. The project led to several breakthroughs in the development of new methods as well as new science results. The significant accomplishments of this project are as follows: 1. We developed a conceptual framework to quantify mixing state based on rigorous information-theoretic entropy concepts, and used this to quantify climate-relevant aerosol impacts. Our new mixing state index has seen significant adoption by experimentalists and modelers. 2. We developed new numerical algorithms to enable the efficient simulation of particle-resolved aerosol microphysics within 3D atmospheric models, and implemented these in our coupled WRF-PartMC model. This makes WRF-PartMC the world’s most accurate model of aerosol microphysics coupled to atmospheric dynamics, and we performed the first ever 3D regional-scale simulations that fully resolve aerosol mixing state. 3. We integrated our modeling and conceptual framework with observations including ASR field campaigns (CARES) and ASR-supported laboratory campaigns (Boston College BC-3).

58 GEOSCIENCES↗

High Performance Interactive System Dynamics Visualization

This brochure describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

High Performance Interactive System Dynamics Visualization

This presentation describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

Accuracy of Length of Virtual Stents in Treatment of Intracranial Wide-Necked Aneurysms

Background and Purpose: Precise stent deployment is important for successful treatment of intracranial aneurysms by stent-assisted coiling (SAC). We evaluated the accuracy of virtual stents generated using commercial stent planning software by comparing the length of virtual and actually deployed intracranial laser cut stents on three-dimensional digital subtraction angiography (3D-DSA) images. Methods: We retrospectively analyzed the data of 75 consecutive cases of intracranial wide-necked aneurysms treated with the SAC technique using laser cut stents. Based on 3D-DSA images acquired by C-arm CT, stent sizing and placement were intraoperatively simulated by a commercial software application. The difference in length of the stents was estimated by measuring proximal discrepancies between the end points of the virtual and actually deployed stents on fused pre-procedural and post-procedural 3D-DSA images. Discrepancies between distal stent end points were manually minimized. The Kruskal–Wallis test was applied to test whether stent location, type, and length had an effect on difference in length between virtual and real stent. Results: The median difference in length between virtual and real stents was 1.58 mm with interquartile range 1.12–2.12 mm. There was no evidence for an effect of stent location (p = 0.23), stent type (p = 0.33), or stent length (p = 0.53) on difference in length between virtual and real stents. Conclusions: Stent planning software allows 3D simulation of laser cut stents overlain on 3D-DSA images of vessels and may thus be useful for stent selection and deployment of laser cut stents during stent-assisted coiling of intracranial aneurysms.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

IMPRINTS OF EXPANSION ON THE LOCAL ANISOTROPY OF SOLAR WIND TURBULENCE

We study the anisotropy of II-order structure functions (SFs) defined in a frame attached to the local mean field in three-dimensional (3D) direct numerical simulations of magnetohydrodynamic turbulence, with the solar wind expansion both included and not included. We simulate spacecraft flybys through the numerical domain by taking increments along the radial (wind) direction that form an angle of 45° with the ambient magnetic field. We find that only when expansion is taken into account do the synthetic observations match the 3D anisotropy observed in the solar wind, including the change of anisotropy with scale. Our simulations also show that the anisotropy changes dramatically when considering increments oblique to the radial directions. Both results can be understood by noting that expansion reduces the radial component of the magnetic field at all scales, thus confining fluctuations in the plane perpendicular to the radial. Expansion is thus shown to affect not only the (global) spectral anisotropy, but also the local anisotropy of second-order SF by influencing the distribution of the local mean field, which enters this higher-order statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum simulations of nuclei and nuclear pasta with the multiresolution adaptive numerical environment for scientific simulations

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation, and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle three-dimensional (3D) Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of 16 O, 208 Pb, and 238 U nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24 fm with an average density of about ρ = 0.05fm -3 , proton fraction of Yp = 0.3, and temperature of T = 0 MeV. Overall, our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. We compare calculations with and without spin-orbit forces. We find that while subtle differences are present, the pasta phase remains in the waffle geometry. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Monte Carlo Simulation with CAD Interface for Calculation of 3D Maps of Residual Dose (CRADA)

Objective: To develop an easy-to-use software application to predict and mitigate radiation effects in research environment, space instruments, nuclear plants and medical facilities and help nonproliferation and national security efforts. Tech-X will develop standalone software libraries and command-line tools for ( 1) translating CAD into tessellated surfaces and tetrahedral meshes in GDML (for Geant4 and MARS 15), ROOT (for MARS 15) and HDF5 (for compact representation and for the visualization) formats, (2) healing CAD geometries to make them suitable for Monte Carlo simulations; (3) creating uniform and variable Cartesian and cylindrical meshes for detailed scoring; and ( 4) efficient Monte Carlo navigation in CAD geometries. JLAB will finish automation of simulations of residual dose in CAD geometries and integrate Tech-X software into Geant4 and MARS15. Finally, Tech-X will develop a Graphical User Interface to set up and heal CAD geometries, create input files, run and visualize simulations for residual dose. This application will run on local desktops, local and remote clusters and supercomputers and will be made available through public clouds, such as Amazon Web Services.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Composite Structured Surfaces for Durable Dropwise Condensation

Condensation heat transfer is a vital process for a plethora of industrial application. Dropwise condensation of steam on hydrophobic polymers, typified by the formation of discrete water droplets which shed and clear the surface for re-nucleation, has a 10X higher heat transfer coefficient when compared to filmwise condensation on bare metallic surfaces. However, 8 decades of research have conclusively demonstrated a lack of durability of the hydrophobic coatings during dropwise condensation to quantify industrial implementation. They key bottleneck is the coating conundrum, which states that thick hydrophobic coatings (> 10 µm) are required for higher durability, while thin hydrophobic coatings (< 1 µm) are required for low parasitic thermal resistance due to the low intrinsic thermal conductivity of polymeric materials (k ~ 0.1 W/m·K). Here, we attempt to solve the conundrum by developing rationally-designed metal-polymer composite coatings to increase the effective thermal conductivity of the hydrophobic coating, thereby enabling thicker coatings for higher durability. We study nanowire, inverse opal, and sintered sphere metallic structures infilled with hydrophobic polymers. To quantify the heat transfer performance, we analyze the vapor-to-coolant and vapor-to-surface heat transfer coefficients of composite coated flat surfaces having varying wire/sphere diameters (D = 100 nm, 1 μm, 10 μm) and coating thickness (10 nm < t < 100 μm) using three-dimensional (3D), steady state, Comsol simulations. We show that analytical models such as Maxwell or effective medium theory fail to predict the composite thermal conductivity. By conducting simulation with two independent sets of boundary conditions; isothermal surfaces and convection, as well as by varying the interfacial metal fractions, we demonstrate the importance of 3D heat spreading near the coating-vapor and coating-substrate interface. To predict the wetting characteristics to ensure stable dropwise condensation, we coupled our 3D heat transfer simulations to a unified model for contact angle hysteresis on heterogeneous surfaces. Our wetting predictions identify the key tradeoff between low spreading resistance (high k, filmwise condensation) and low contact angle hysteresis (low k, dropwise condensation) on high and low surface metal fractions, respectively, and identify the optimum metallic structure, volume fraction, interfacial solid fraction, and thickness for enhanced durability. Our work provides a novel approach to utilize rational metallic nanostructuring to create composite hydrophobic coatings which promise enhanced durability.

Chang, Ho Chan↗

3D Printing Prototypes for a Reflector-Based Broadband Axion Search

This work sought to begin research and development on a horn antenna which can be used to transfer the energy of microwave photons produced by axions in a strong magnetic field to a coaxial cable for measurement. 3D prototypes of pyramidal horn antennas and a parabolic reflector were printed and then coated with conductive spray paint. The S 11 parameter was measured and simulated for both 3D printed horn antenna and the conical horn antenna both with and without the reflector present to determine how much of the radiated power was reflected back. Ultimately, for both the conical and pyramidal antennas, less than one percent of the energy emitted by the antennas was reflected back. When the pyramidal antenna was tilted by about 45 degrees, however, around 30-40 percent of the microwaves emitted by the antennas was reflected back.

43 PARTICLE ACCELERATORS↗

SU-F-BRD-11: A Virtual Simulator Designed for Collision Prevention in Proton Therapy

Purpose: In proton therapy, collisions between patient and nozzle potentially occur in attaining minimal air gap due to the large nozzle structure. Thus, we developed software predicting the collisions of the nozzle and patient by simulating treatments. Methods: 3D modeling of a gantry inner-floor, nozzle and robotic-couch was done by using the SolidWorks based on the manufacturer’s machine data. To obtain patient body information, a 3D-scanner was utilized to scan a patient right before CT scanning. From the acquired images, a 3D-image of the patient’s body contour was reconstructed. The accuracy of the image was confirmed against the CT image for a humanoid phantom. The machine components and the virtual patient were combined on the treatment-room coordinate system, resulting in a virtual simulator. The simulator simulated the motion of its components such as rotation and translation of gantry, nozzle and couch, in real scale. Collision, if any, was examined both in static mode and dynamic mode. The static mode checks only at fixed positions of the machine’s components while dynamic mode examines while one component is in motion. Collision was notified if any voxel of two components, for example a nozzle and a patient or couch, overlapped when calculating volume locations. The event and collision point are visualized and colliding volumes are reported. Results: All components were successfully assembled and the motions could be accurately controlled. The 3D-shape of a phantom agreed with CT images within a deviation of 2 mm. Collision situations can be simulated within minutes and the results are displayed and reported. Conclusion: The developed software will be useful in improving patient safety and clinical efficiency for proton therapy. This work was supported by the National Research Foundation of Korea funded by Ministry of Science, ICT & Future Planning (2012M3A9B6055201, 2013M2A2A7043507), and Samsung Medical Center grant (GFO1130081)

60 APPLIED LIFE SCIENCES↗

Heat loads to divertor nearby components from secondary radiation evolved during plasma instabilities

A fundamental issue in tokamak operation related to power exhaust during plasma instabilities is the understanding of heat and particle transport from the core plasma into the scrape-off layer and to plasma-facing materials. During abnormal and disruptive operation in tokamaks, radiation transport processes play a critical role in divertor/edge-generated plasma dynamics and are very important in determining overall lifetimes of the divertor and nearby components. This is equivalent to or greater than the effect of the direct impact of escaped core plasma on the divertor plate. We have developed and implemented comprehensive enhanced physical and numerical models in the upgraded HEIGHTS package for simulating detailed photon and particle transport in the evolved edge plasma during various instabilities. The paper describes details of a newly developed 3D Monte Carlo radiation transport model, including optimization methods of generated plasma opacities in the full range of expected photon spectra. Response of the ITER divertor's nearby surfaces due to radiation from the divertor-developed plasma was simulated by using actual full 3D reactor design and magnetic configurations. We analyzed in detail the radiation emission spectra and compared the emission of both carbon and tungsten as divertor plate materials. The integrated 3D simulation predicted unexpectedly high damage risk to the open stainless steel legs of the dome structure in the current ITER design from the intense radiation during a disruption on the tungsten divertor plate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aircraft observations and cloud-microwave radiative model simulations of a deep convective system

A 3D cloud model is used to simulate the storm structure of the June 29, 1986 storm system observed during the Cooperative Huntsville Meteorological Experiment (Dodge et al., 1986), and the results are linked to a microwave radiative transfer model for simulating aircraft observations. The results of cloud and radiative model simulations are then compared with radar data and aircraft microwave observations, respectively. It is shown that the degree of agreement between model simulations and the aircraft microwave observations depends on the model-derived microphysics of cloud and precipitation. It is found that the calculated values of the brightness temperature (Tb) at 92 GHz are close to the observed temperatures, while the calculated Tb values at 18 and 37 GHz are lower than the observed Tb.

Yeh, H.-Y. M.↗

Computational Simulation of Damage Propagation in Three-Dimensional Woven Composites

Three dimensional (3D) woven composites have demonstrated multi-directional properties and improved transverse strength, impact resistance, and shear characteristics. The objective of this research is to develop a new model for predicting the elastic constants, hygrothermal effects, thermomechanical response, and stress limits of 3D woven composites; and to develop a computational tool to facilitate the evaluation of 3D woven composite structures with regard to damage tolerance and durability. Fiber orientations of weave and braid patterns are defined with reference to composite structural coordinates. Orthotropic ply properties and stress limits computed via micromechanics are transformed to composite structural coordinates and integrated to obtain the 3D properties. The various stages of degradation, from damage initiation to collapse of structures, in the 3D woven structures are simulated for the first time. Three dimensional woven composite specimens with various woven patterns under different loading conditions, such as tension, compression, bending, and shear are simulated in the validation process of this research. Damage initiation, growth, accumulation, and propagation to fracture are included in these simulations.

Huang, Dade↗

Dynamics and Structure of Main-Sequence Stars with Shallow Convection Zones

A dramatic increase in observational data from NASA’s Kepler, K2, and TESS missions and supporting ground-based observatories has opened new opportunities to investigate the internal structure, dynamics, and evolution of stars and their atmospheres. We present 3D radiative MHD simulations for several main-sequence stars with masses from 1.4 to 1.5M sun. The simulations are performed using the "StellarBox" code developed for modeling stellar turbulent convection and atmospheres with a high degree of realism. This presentation discusses similarities and differences between 3D realistic-type and 1D mixing-length models with regard to structural, thermodynamic, and turbulent property variations from the radiative zone to the convection zone and photosphere.

SMD↗

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗