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,063 records · Page 59

Understanding the Flow Physics of Shock Boundary-Layer Interactions Using CFD and Numerical Analyses

Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.

Friedlander, David J.↗

Detection of the Compressed Primary Stellar Wind in eta Carinae

A series of three HST/STIS spectroscopic mappings, spaced approximately one year apart, reveal three partial arcs in [Fe II] and [Ni II] emissions moving outward from eta Carinae. We identify these arcs with the shell-like structures, seen in the 3D hydrodynamical simulations, formed by compression of the primary wind by the secondary wind during periastron passages.

eta Carinae↗

Detection of the Compressed Primary Stellar Wind in eta Carinae*

A series of three Hubble Space Telescope Space Telescope Imaging Spectrograph (HST/STIS) spectroscopic mappings, spaced approximately one year apart, reveal three partial arcs in [Fe II] and [Ni II] emissions moving outward from η Carinae. We identify these arcs with the shell-like structures, seen in the 3D hydrodynamical simulations, formed by compression of the primary wind by the secondary wind during periastron passages.

spectroscopic mappings↗

Contamination Analysis Tools

This talk presents 3 different tools developed recently for contamination analysis:HTML QCM analyzer: runs in a web browser, and allows for data analysis of QCM log filesJava RGA extractor: can load in multiple SRS.ana files and extract pressure vs. time dataC++ Contamination Simulation code: 3D particle tracing code for modeling transport of dust particulates and molecules. Uses residence time to determine if molecules stick. Particulates can be sampled from IEST-STD-1246 and be accelerated by aerodynamic forces.

Climate of Earth-Like Planets With and Without Ocean Heat Transport Orbiting a Range of M and K Stars

The mean surface temperature of a planet is now acknowledged as insufficient to surmise its full potential habitability. Advancing our understanding requires exploration with 3D general circulation models (GCMs), which can take into account how gradients and fluxes across a planet's surface influence the distribution of heat, clouds, and the potential for heterogeneous distribution of liquid water. Here we present 3D GCM simulations of the effects of alternative stellar spectra, instellation, model resolution, and ocean heat transport, on the simulated distribution of heat and moisture of an Earth-like planet (ELP).

Planets↗

Effect of Dietary Countermeasures and Impact of Gravity on Renal Calculi Size Distributions Predicted by PBE-System and PBE-CFD Models

Renal stone disease is not only a concern on earth but can conceivably pose a serious risk to the astronauts health and safety in Space. In this work, two different deterministic models based on a Population Balance Equation (PBE) analysis of renal stone formation are developed to assess the risks of critical renal stone incidence for astronauts during space travel. In the first model, the nephron is treated as a continuous mixed suspension mixed product removal crystallizer and the PBE for the nucleating, growing and agglomerating renal calculi is coupled to speciation calculations performed by JESS. Predictions of stone size distributions in the kidney using this model indicate that the astronaut in microgravity is at noticeably greater but still subcritical risk and recommend administration of citrate and augmented hydration as effective means of minimizing and containing this risk. In the second model, the PBE analysis is coupled to a Computational Fluid Dynamics (CFD) model for flow of urine and transport of Calcium and Oxalate in the nephron to predict the impact of gravity on the stone size distributions. Results presented for realistic 3D tubule and collecting duct geometries, clearly indicate that agglomeration is the primary mode of size enhancement in both 1g and microgravity. 3D numerical simulations seem to further indicate that there will be an increased number of smaller stones developed in microgravity that will likely pass through the nephron in the absence of wall adhesion. However, upon reentry to a 1g (Earth) or 38g (Mars) partial gravitational fields, the renal calculi can lag behind the urinary flow in tubules that are adversely oriented with respect to the gravitational field and grow agglomerate to large sizes that are sedimented near the wall with increased propensity for wall adhesion, plaque formation, and risk to the astronauts.

Renal Stones↗

A Breakout Model for Solar Coronal Jets with Filaments

Recent observations have revealed that many solar coronal jets involve the eruption of miniature versions of large-scale filaments. Such ''mini-filaments'' are observed to form along the polarity inversion lines of strong, magnetically bipolar regions embedded in open (or distantly closing) unipolar field. During the generation of the jet, the filament becomes unstable and erupts. Recently we described a model for these mini-filament jets, in which the well-known magnetic-breakout mechanism for large-scale coronal mass ejections is extended to these smaller events. In this work we use 3D magnetohydrodynamic simulations to study in detail three realizations of the model. We show that the breakout-jet generation mechanism is robust and that different realizations of the model can explain different observational features. The results are discussed in relation to recent observations and previous jet models.

Sun: magnetic fields↗

Climate Modeling of a Potential Exovenus

The planetary mass and radius sensitivity of exoplanet discovery capabilities has reached into the terrestrial regime. The focus of such investigations is to search within the Habitable Zone where a modern Earth-like atmosphere may be a viable comparison. However, the detection bias of the transit and radial velocity methods lies close to the host star where the received flux at the planet may push the atmosphere into a runaway greenhouse state. One such exoplanet discovery, Kepler-1649b, receives a similar flux from its star as modern Venus does from the Sun, and so was categorized as a possible exoVenus. Here we discuss the planetary parameters of Kepler-1649b in relation to Venus to establish its potential as a Venus analog. We utilize the general circulation model ROCKE-3D to simulate the evolution of the surface temperature of Kepler-1649b under various assumptions, including relative atmospheric abundances. We show that in all our simulations the atmospheric model rapidly diverges from temperate surface conditions toward a runaway greenhouse with rapidly escalating surface temperatures. We calculate transmission spectra for the evolved atmosphere and discuss these spectra within the context of the James Webb Space Telescope Near-Infrared Spectrograph capabilities. We thus demonstrate the detectability of the key atmospheric signatures of possible runaway greenhouse transition states and outline the future prospects of characterizing potential Venus analogs.

Astrobiology↗

Massive Galaxies Impede Massive Outflows

A set of 66 3D hydrodynamical simulations explores how galactic stellar mass affects three-phase, starburst-driven outflows. Simulated velocities are compared to two basic analytic models: with and without a gravitational potential. For stellar mass <1010 M ⊙, simulated velocities match those of both analytical models and are unaffected by the potential; above they reduce significantly as expected from the analytic model with gravity. Gravity also affects total outflow mass and each of the three phases differently. Outflow masses in the hot, warm, and cold phases each scale with stellar mass as $\mathrm{log}{M}_{* }\,=$ −0.25, −0.97, and −1.70, respectively. Thus, the commonly used Chevalier & Clegg analytic model should be modified to include gravity when applied to higher-mass galaxies. In particular, using M82 as the canonical galaxy to interpret hydrodynamical simulations of starburst-driven outflows from higher-mass galaxies will underestimate the retarding effect of gravity. Using the analytic model of Johnson & Axford with realistic thermalization efficiency and mass loading, I find that only galaxy masses lesssim10(exp 11.5) M(⊙) can outflow.

Ryan Tanner↗

3D Representation of UAV-obstacle Collision Risk under off-nominal conditions

Safe operations of autonomous unmanned aerial vehicles (UAVs) in low-altitude airspace with beyond visual line-of-sight (BVLOS) flights demand robust risk monitoring of airspace as well as of people and property on ground. One of the safety critical factors for UAV flights is the risk of collision with static and dynamic obstacles in proximity to its flight path. This paper presents a detailed formulation of risk likelihood of obstacle collision incorporating the effects of off-nominal conditions introduced by component failures, degraded controllability and environmental disturbances such as wind gusts. The deviation in the planned trajectory caused due to wind is computed utilizing a point-mass 3D kinematic simulation model of the vehicle. Likelihood of risk for the flight plan is then analyzed based on generating the probability of collision for each point in the trajectory. The proposed risk factor is demonstrated on real flight data from experimental flights of an octocopter at NASA Langley Research Center in presence of simulated obstacles and wind conditions. Effect of varying wind conditions, distance from obstacles, level of controllability and obstacle measurement noise on the risk factor is demonstrated. The proposed approach enables risk-informed decision making for timely mitigation of current and future unsafe events in autonomous systems.

Portia Banerjee↗

Advances and Challenges in Observations and Modeling of the Global-Sun Dynamics and Dynamo

Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature of the solar dynamics, covering a vast range of scales, which can-not be solved with the current computational resources. In such a situation, significant progress has been achieved by the mean-field approach, based on the separation of small-scale turbulence and large-scale dynamics. The mean-field simulations can reproduce solar observations, qualitatively and quantitatively, and uncover new phenomena. However, they do not reveal the complex physics of large-scale convection, solar magnetic cycles, and the magnetic self-organization that causes sunspots and solar eruptions. Thus, developing a synergy of these approaches seems to be a necessary but very challenging task.

SMD↗

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

Magnetic Reconnection in a Sheared Magnetic Flux Tube: Slippage Versus Tearing

The process of magnetic reconnection in a flux tube can occur in a time-stationary fashion as slippage reconnection or in a time-dependent manner based on, for example, the tearing instability. However, it is not well known under which conditions a system can sustain slippage reconnection. Likewise, it is unclear whether systems can exhibit slippage reconnection and time-dependent reconnection simultaneously. In order to investigate these questions, we employ a set of 3D magnetohydrodynamic simulation. Using these simulations, we model a twisted flux tube with a spatially localized resistive region in the center of the simulation box, and an applied driver on one end plane of the simulation box. As a result, the Poynting flux injected at the boundary propagates to the resistive region and dissipates there. If the driver is nearly circular, we find a predominance of slippage reconnection, which decouples the applied driver entirely from one half of the flux tube. Second, a tearing-like instability occurs with slippage reconnection in the resistive region when the resistivity is strong enough and the cross section of the flux tube, as defined by the shape of the applied velocity perturbation, is elliptical enough. In this case we find the surprising feature that magnetic field perturbation generated by the tearing-like instability radiates away from the resistive region to influence the entire flux tube shape.

Hidetaka Kuniyoshi↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

Compressible Turbulence in the Near-Sun Solar Wind Parker Solar Probe's First Eight Perihelia

Many questions remain about the compressibility of solar wind turbulence with respect to its origins and properties. Low plasma beta (ratio of thermal to magnetic pressure) environments allow for the easier generation of compressible turbulence, enabling study of the relationship between density fluctuations and turbulent Mach number. Utilizing Parker Solar Probe plasma data, we examine the normalized proton density fluctuations as a function of turbulent Mach number M t conditioned on plasma beta and cross helicity. With consideration of statistical error in the parameters computed from in situ data, we find a general result consistent with both linear-wave theory and nearly incompressible turbulence in an inhomogeneous background field. We compare observational results conditioned on plasma beta and cross helicity with 3D magnetohydrodynamic simulations and observe rather significant similarities with respect to how those parameters affect the proportionality between density fluctuations and turbulent Mach number. This study further investigates the complexity of compressible turbulence as viewed by the density scaling relationship and may help better understand the compressible environment of the near-Sun solar wind.

Manuel Enrique Cuesta↗

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗