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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Evaluation of Cloud-resolving and Limited Area Model Intercomparison Simulations using TWP-ICE Observations. Part 1: Deep Convective Updraft Properties

Ten 3D cloud-resolving model (CRM) simulations and four 3D limited area model (LAM) simulations of an intense mesoscale convective system observed on January 23-24, 2006 during the Tropical Warm Pool – International Cloud Experiment (TWP-ICE) are compared with each other and with observed radar reflectivity fields and dual-Doppler retrievals of vertical wind speeds in an attempt to explain published results showing a high bias in simulated convective radar reflectivity aloft. This high bias results from ice water content being large, which is a product of large, strong convective updrafts, although hydrometeor size distribution assumptions modulate the size of this bias. Snow reflectivity can exceed 40 dBZ in a two-moment scheme when a constant bulk density of 100 kg m-3 is used. Making snow mass more realistically proportional to area rather than volume should somewhat alleviate this problem. Graupel, unlike snow, produces high biased reflectivity in all simulations. This is associated with large amounts of liquid water above the freezing level in updraft cores. Peak vertical velocities in deep convective updrafts are greater than dual-Doppler retrieved values, especially in the upper troposphere. Freezing of large rainwater contents lofted above the freezing level in simulated updraft cores greatly contributes to these excessive upper tropospheric vertical velocities. Strong simulated updraft cores are nearly undiluted, with some showing supercell characteristics. Decreasing horizontal grid spacing from 900 meters to 100 meters weakens strong updrafts, but not enough to match observational retrievals. Therefore, overly intense simulated updrafts may partly be a product of interactions between convective dynamics, parameterized microphysics, and large-scale environmental biases that promote different convective modes and strengths than observed.

clould-resolving↗
DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗
DOE OSTI2016

ENERGY DISSIPATION AND LANDAU DAMPING IN TWO- AND THREE-DIMENSIONAL PLASMA TURBULENCE

Plasma turbulence is ubiquitous in space and astrophysical plasmas, playing an important role in plasma energization, but the physical mechanisms leading to dissipation of the turbulent energy remain to be definitively identified. Kinetic simulations in two dimensions (2D) have been extensively used to study the dissipation process. How the limitation to 2D affects energy dissipation remains unclear. This work provides a model of comparison between two- and three-dimensional (3D) plasma turbulence using gyrokinetic simulations; it also explores the dynamics of distribution functions during the dissipation process. It is found that both 2D and 3D nonlinear gyrokinetic simulations of a low-beta plasma generate electron velocity-space structures with the same characteristics as that of the linear Landau damping of Alfvén waves in a 3D linear simulation. The continual occurrence of the velocity-space structures throughout the turbulence simulations suggests that the action of Landau damping may be responsible for the turbulent energy transfer to electrons in both 2D and 3D, and makes possible the subsequent irreversible heating of the plasma through collisional smoothing of the velocity-space fluctuations. Although, in the 2D case where variation along the equilibrium magnetic field is absent, it may be expected that Landau damping is not possible, a common trigonometric factor appears in the 2D resonant denominator, leaving the resonance condition unchanged from the 3D case. The evolution of the 2D and 3D cases is qualitatively similar. However, quantitatively, the nonlinear energy cascade and subsequent dissipation is significantly slower in the 2D case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Extraction of effective solid-liquid interfacial free energies for full 3D solid crystallites from equilibrium MD simulations

Molecular dynamics simulations of an embedded atom copper system in the NPH ensemble are used to study the e ective solid-liquid interfacial free energy of quasispherical solid crystals within a liquid. This is within the larger context of MD simulations of this system undergoing solidi cation, where single individually-prepared crystallites of di erent sizes grow until they reach a thermodynamically stable nal state. The resulting equilibrium shapes possess the full structural details expected for solids with weakly anisotropic surface free energies (in these cases, ~5 % radial attening and rounded [111] octahedral faces). The simplifying assumption of sphericity and perfect isotropy leads to an e ective interfacial free energy as appearing in the Gibbs-Thomson equation, which we determine to be ~179 erg/cm 2 , roughly independent of crystal size for radii in the 50 - 250 A range. This quantity may be used in atomistically-informed models of solidi cation kinetics for this system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using photoelectron spectroscopy to measure resonant inelastic X-ray scattering: a computational investigation

Resonant inelastic X-ray scattering (RIXS) has become an important scientific tool. Nonetheless, conventional high-resolution (few hundred meV or less) RIXS measurements, especially in the soft X-ray range, require low-throughput grating spectrometers, which limits measurement accuracy. Here, the performance of a different method for measuring RIXS, i.e. photoelectron spectrometry for analysis of X-rays (PAX), is computationally investigated. This method transforms the X-ray measurement problem of RIXS to an electron measurement problem, enabling use of high-throughput, compact electron spectrometers. X-rays to be measured are incident on a converter material and the energy distribution of the resultant photoelectrons, the PAX spectrum, is measured with an electron spectrometer. A deconvolution algorithm for analysis of such PAX data is proposed. It is shown that the deconvolution algorithm works well on data recorded with ~0.5 eV resolution. Additional simulations show the potential of PAX for estimation of RIXS features with smaller widths. For simulations using the 3d levels of Ag as a converter material, and with 10 5 simulated detected electrons, it is estimated that features with a few hundred meV width can be accurately estimated in a model RIXS spectrum. For simulations using a sharp Fermi edge to encode RIXS spectra, it is estimated that one can accurately distinguish 100 meV FWHM peaks separated by 45 meV with 10 5 simulated detected electrons that were photoemitted from within 0.4 eV of the Fermi level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

The meteorological environment of the tropospheric ozone maximum over the tropical South Atlantic Ocean

Atmospheric flow patterns are examined over the South Atlantic Ocean where a maximum of tropospheric ozone has been observed just west of southern Africa. We investigate the flow climatology during October and perform a case study for six days during October 1989. Horizontal and vertical motions are examined and used to prepare 3D backward trajectories from the region of greatest ozone. An initially zonally symmetric distribution of ozone is treated as a passive tracer and advected by 3D flows forecast by the global model. Results from the passive tracer simulation indicate that 3D advection alone can produce a maximum of tropospheric ozone in the observed location. In addition, the trajectories suggest that by-products of biomass burning could be transported to the area of maximum ozone. Low-level flow from commonly observed regions of burning in Africa streams westward to the area of interest. Over Brazil, if the burning by-products are carried into the upper troposphere by convective process, they then could be transported eastward to the ozone feature in approximately five days. There is considerable subsidence over the tropical southern Atlantic, such that stratospheric influences also are a factor in producing the ozone maximum. Both planetary-scale and transient synoptic-scale circulation features play major roles in the various transport processes that influence the region. In summary, the observed tropospheric ozone maximum appears to be caused by a complex set of horizontal and vertical advections, transport from regions of biomass burning, and stratospheric influences.

Krishnamurti, T. N.↗
NASA NTRS2011

Calculation of Dose Deposition in 3D Voxels by Heavy Ions and Simulation of gamma-H2AX Experiments

The biological response to high-LET radiation is different from low-LET radiation due to several factors, notably difference in energy deposition and formation of radiolytic species. Of particular importance in radiobiology is the formation of double-strand breaks (DSB), which can be detected by -H2AX foci experiments. These experiments has revealed important differences in the spatial distribution of DSB induced by low- and high-LET radiations [1,2]. To simulate -H2AX experiments, models based on amorphous track with radial dose are often combined with random walk chromosome models [3,4]. In this work, a new approach using the Monte-Carlo track structure code RITRACKS [5] and chromosome models have been used to simulate DSB formation. At first, RITRACKS have been used to simulate the irradiation of a cubic volume of 5 m by 1) 450 1H+ ions of 300 MeV (LET 0.3 keV/ m) and 2) by 1 56Fe26+ ion of 1 GeV/amu (LET 150 keV/ m). All energy deposition events are recorded to calculate dose in voxels of 20 m. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. Many differences are found in the spatial distribution of dose voxels for the 56Fe26+ ion. The track structure can be distinguished, and voxels with very high dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and indicate clustered energy deposition, which may be responsible for complex DSB. In the second step, assuming that DSB will be found only in voxels where energy is deposited by the radiation, the intersection points between voxels with dose > 0 and simulated chromosomes were obtained. The spatial distribution of the intersection points is similar to -H2AX foci experiments. These preliminary results suggest that combining stochastic track structure and chromosome models could be a good approach to understand radiation-induced DSB and chromosome aberrations.

Plante, I.↗
NASA NTRS2013

3D Guided Wave Motion Analysis on Laminated Composites

Ultrasonic guided waves have proved useful for structural health monitoring (SHM) and nondestructive evaluation (NDE) due to their ability to propagate long distances with less energy loss compared to bulk waves and due to their sensitivity to small defects in the structure. Analysis of actively transmitted ultrasonic signals has long been used to detect and assess damage. However, there remain many challenging tasks for guided wave based SHM due to the complexity involved with propagating guided waves, especially in the case of composite materials. The multimodal nature of the ultrasonic guided waves complicates the related damage analysis. This paper presents results from parallel 3D elastodynamic finite integration technique (EFIT) simulations used to acquire 3D wave motion in the subject laminated carbon fiber reinforced polymer composites. The acquired 3D wave motion is then analyzed by frequency-wavenumber analysis to study the wave propagation and interaction in the composite laminate. The frequency-wavenumber analysis enables the study of individual modes and visualization of mode conversion. Delamination damage has been incorporated into the EFIT model to generate "damaged" data. The potential for damage detection in laminated composites is discussed in the end.

Tian, Zhenhua↗
NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗
NASA NTRS2020

A Jupyter Notebook Environment For Multibody Dynamics

DARTS is a rigid/flexible multibody dynamics toolkit for themodeling and simulation of aerospace and robotic vehicles forengineering applications. In this paper we describe an on-line,browser-based environment using Jupyter notebooks to supporttraining needs for the DARTS software. The suite of curated tutorial notebooks is organized into different topic areas, and intomultiple themes within each topic area. The notebooks within atheme use a progression of examples for users to expand theirunderstanding of the software. The topic areas include one onthe DARTS multibody dynamics software and another one on thetheory underlying the multibody dynamics formulation. We alsodescribe a number of Jupyter extensions that were used - andsome developed in house - to enhance the notebook interface foruse with the dynamics simulation software. One significant extension we implemented allows the embedding of live 3D visualizations within simulation notebooks.

Gaut, Aaron↗
DOE OSTI2014

Fracture Propagation, Fluid Flow, and Geomechanics of Water-Based Hydraulic Fracturing in Shale Gas Systems and Electromagnetic Geophysical Monitoring of Fluid Migration

We investigate fracture propagation induced by hydraulic fracturing with water injection, using numerical simulation. For rigorous, full 3D modeling, we employ a numerical method that can model failure resulting from tensile and shear stresses, dynamic nonlinear permeability, leak-off in all directions, and thermo-poro-mechanical effects with the double porosity approach. Our numerical results indicate that fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics, whereas water saturation is determined by fluid flow. At early times, the water saturation front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with injected water. However, at late times, advance of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, which is filled with reservoir gas. We also find considerable leak-off of water to the reservoir. The inconsistency between the fracture volume and the volume of injected water cannot properly calculate the fracture length, when it is estimated based on the simple assumption that the fracture is fully saturated with injected water. As an example of flow-geomechanical responses, we identify pressure fluctuation under constant water injection, because hydraulic fracturing is itself a set of many failure processes, in which pressure consistently drops when failure occurs, but fluctuation decreases as the fracture length grows. We also study application of electromagnetic (EM) geophysical methods, because these methods are highly sensitive to changes in porosity and pore-fluid properties due to water injection into gas reservoirs. Employing a 3D finite-element EM geophysical simulator, we evaluate the sensitivity of the crosswell EM method for monitoring fluid movements in shaly reservoirs. For this sensitivity evaluation, reservoir models are generated through the coupled flow-geomechanical simulator and are transformed via a rock-physics model into electrical conductivity models. It is shown that anomalous conductivity distribution in the resulting models is closely related to injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for distribution of injected fluids (i.e., migration pathways) during hydraulic fracturing operations

02 PETROLEUM↗