Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D equilibrium”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

The VATMOS-SR Mission Concept: DSMC Studies of the Gas Sampling

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept that would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a complete picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O, and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of volatile atmospheric signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process to determine how measured isotopic ratios of noble gases in the sample compared to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open-source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed on Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters, as well as the modeling parameters, will be performed.

direct simulation Monte Carlo↗
NASA NTRSDate not supplied

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗
NASA NTRS1999

3D-Stereoscopic Analysis of Solar Active Region Loops: SoHo/EIT Observations at Temperatures of 1.5-2.5 MK - 2

In this paper we study the three-dimensional (3D) structure of hot (T(sub e) approximately equals 1.5 - 2.5 MK) loops in solar active region NOAA 7986, observed on 1996 August 30 with the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SoHO). This complements a first study on cooler (T(sub e) approximately equals 1.0 - 1.5 MK) loops of the same active region, using the same method of Dynamic Stereoscopy to reconstruct the 3D geometry. We reconstruct the 3D-coordinates x(s), y(s), z(s), the density n(sub e)(s), and temperature profile T(sub e)(s) of 35 individual loop segments (as function of the loop coordinate s) using EIT 195 A and 284 A images. The major findings are: (1) All loops are found to be in hydrostatic equilibrium, in the entire temperature regime of T(sub e) = 1.0 - 2.5 MK; (2) The analyzed loops have a height of 2-3 scale heights, and thus only segments extending over about one vertical scale height have sufficient emission measure contrast for detection; (3) The temperature gradient over the lowest scale height is of order dT/ds is approximately 1 - 4 K/km; (4) The radiative loss rate is found to exceed the conductive loss rate by about two orders or magnitude, making thermal conduction negligible to explain the temperature structure of the loops; (5) A steady-state can only be achieved when the heating rate E(sub H) matches the radiative loss rate in hydrostatic equilibrium, requiring a heat deposition length lambda(sub H) of the half density scale height lambda, predicting a scaling law with the loop base pressure, EH varies as p(sub 0 exp 2). This favors coronal heating mechanisms that operate near the loop footpoints; (6) We find a reciprocal correlation between the loop pressure p(sub 0) and loop length L, i.e. p(sub 0) varies as 1/L, implying a scaling law of the steady-state requirement with loop length, i.e. E(sub H ) varies as 1/L(exp 2). The heating rate shows no correlation with the loop-aligned magnetic field component B(sub z) at the footpoints, but is correlated with the azimuthal field B(sub phi) = Bz(RDelta Phi/L) of a twisted loop, and is thus consistent with heating mechanisms based on field-aligned currents.

Aschwanden, Markus J.↗
NASA NTRS1990

X-ray spectral signatures of photoionized plasmas

Plasma emission codes have become a standard tool for the analysis of spectroscopic data from cosmic X-ray sources. However, the assumption of collisional equilibrium, typically invoked in these codes, renders them inapplicable to many important astrophysical situations, particularly those involving X-ray photoionized nebulae. This point is illustrated by comparing model spectra which have been calculated under conditions appropriate to both coronal plasmas and X-ray photoionized plasmas. It is shown that the (3s-2p)/(3d-2p) line ratios in the Fe L-shell spectrum can be used to effectively discriminate between these two cases. This diagnostic will be especially useful for data analysis associated with AXAF and XMM, which will carry spectroscopic instrumentation with sufficient sensitivity and resolution to identify X-ray photoionized nebulae in a wide range of astrophysical environments.

Liedahl, Duane A.↗
NASA NTRS1992

A numerical study of the effects of geometry on the performance of a supersonic combustor

A numerical study is conducted to investigate the effects of two combustor scale parameters, the expansion angle and the length of the constant area combustor section, on the mixing and combustion within a supersonic combustor. This study uses the SPARK 3D Navier-Stokes code. Turbulence is modeled utilizing an algebraic eddy viscosity model and the chemical mechanism is modeled with a 7-reaction, 7-species finite-rate chemistry model. The calculations show that, even at the relatively low flight Mach numbers (5-7) associated with the conditions examined, the chemical constituents are far from equilibrium, and therefore that kinetics effects are important.

Eklund, Dean R.↗
NASA NTRS2004

Atomic Data and Spectral Line Intensities for Ne III

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for Ne III. The configurations used are 2s(sup 2)2p(sup 4), 2s2p(sup 5), 2p(sup 6), 2s(sup 2)2p(sup 3)3(sub s), 2s(sup 2)2p(sup 3)3p, and 2s(sup 2)2p(sup 3)3d giving rise to 86 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies, 5.2, 10, 15, 20 and 25 Ry. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8) - 10(exp 14)/cc at an electron temperature of logTe/K = 5.0, corresponding to maximum abundance of Ne III. Relative spectral line intensities are calculated.

Landi, E.↗
NASA NTRS2006

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗
NASA NTRS2012

Upper Boundary Extension of the NASA Ames Mars General Circulation Model

Extending the NASA Ames Mars General Circulation Model (MGCM) upper boundary will expand our understanding of the connection between the lower and upper atmosphere of Mars through the middle atmosphere. The extension's main requirements is incorporation of Non-local thermodynamic equilibrium (NLTE) heating (visible) and cooling (infrared). NLTE occurs when energy is exchanged more rapidly with the radiation field (or other energy sources) rather than collisions with other molecules. Without NLTE above approximately 80km/approximately 60km in Mars' atmosphere the IR/visible heating rates are overestimated. Currently NLTE has been applied successfully into the 1D RT code and is in progress for the 3D application.

Brecht, Amanda S.↗
NASA NTRS2017

Flow-Tube Reactor Experiments on the High Temperature Oxidation of Carbon Weaves

Under entry conditions carbon weaves used in thermal protection systems (TPS) decompose via oxidation. Modeling this phenomenon is challenging due to the different regimes encountered along a flight trajectory. Approaches using equilibrium chemistry may lead to over-estimated mass loss and recession at certain conditions. Concurrently, there is a shortcoming of experimental data on carbon weaves to enable development of improved models. In this work, a flow-tube test facility was used to measure the oxidation of carbon weaves at temperatures up to 1500 K. The material tested was the 3D carbon weave used for the heat shield of the NASA Adaptive Deployable Entry and Placement Technology, ADEPT. Oxidation was characterized by quantifying decomposition gases (CO and CO2), by mass measurements, and by microscale surface analysis. The current set of measurements contributes to the development of finite rate chemistry models for carbon fabrics used in woven TPS materials.

Woven TPS↗
NASA NTRS1993

Analysis and design of an ultrahigh temperature hydrogen-fueled MHD generator

A coupled gas dynamics/radiative heat transfer analysis of partially ionized hydrogen, in local thermodynamic equilibrium, flowing through an ultrahigh temperature (10,000-20,000 K) magnetohydrodynamic (MHD) generator is performed. Gas dynamics are modeled by a set of quasi-one-dimensional, nonlinear differential equations which account for friction, convective and radiative heat transfer, and the interaction between the ionized gas and applied magnetic field. Radiative heat transfer is modeled using nongray, absorbing-emitting 2D and 3D P-1 approximations which permit an arbitrary variation of the spectral absorption coefficient with frequency. Gas dynamics and radiative heat transfer are coupled through the energy equation and through the temperature- and density-dependent absorption coefficient. The resulting nonlinear elliptic problem is solved by iterative methods. Design of such MHD generators as onboard, open-cycle, electric power supplies for a particular advanced airbreathing propulsion concept produced an efficient and compact 128-MWe generator characterized by an extraction ratio of 35.5 percent, a power density of 10,500 MWe/cu m, and a specific (extracted) energy of 324 MJe/kg of hydrogen. The maximum wall heat flux and total wall heat load were 453 MW/sq m and 62 MW, respectively.

Moder, Jeffrey P.↗
NASA NTRS1991

Relaxation processes in a low-order three-dimensional magnetohydrodynamics model

The time asymptotic behavior of a Galerkin model of 3D magnetohydrodynamics (MHD) has been interpreted using the selective decay and dynamic alignment relaxation theories. A large number of simulations has been performed that scan a parameter space defined by the rugged ideal invariants, including energy, cross helicity, and magnetic helicity. It is concluded that time asymptotic state can be interpreted as a relaxation to minimum energy. A simple decay model, based on absolute equilibrium theory, is found to predict a mapping of initial onto time asymptotic states, and to accurately describe the long time behavior of the runs when magnetic helicity is present. Attention is also given to two processes, operating on time scales shorter than selective decay and dynamic alignment, in which the ratio of kinetic to magnetic energy relaxes to values 0(1). The faster of the two processes takes states initially dominant in magnetic energy to a state of near-equipartition between kinetic and magnetic energy through power law growth of kinetic energy. The other process takes states initially dominant in kinetic energy to the near-equipartitioned state through exponential growth of magnetic energy.

Stribling, Troy↗
NASA NTRS2001

Atomic Data and Spectral Line Intensities for Ne III

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ne III. The configurations used are 2s(sup 2) 2p(sup 4),2s2p(sup 5),2s(sup 2) 2p(sup 3)3s, and 2s(sup 2)3p(sup 3)3d giving rise to 57 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies, 5, 10, 15, 20, and 25 Ry. Excitation rate coefficients are calculated by assuming a Maxwellian electron velocity distribution at an electron temperature of logT,(K)=5.0, corresponding to maximum abundance of Ne III. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8)-10(exp 14) per cubic centimeter. Relative spectral line intensities are calculated. Proton excitation rates between the lowest three levels have been included in the statistical equilibrium equations. The predicted Ne III line intensities are compared with SERTS rocket measurements of a solar active region and of a laboratory EUV light source.

Bhatia, A. K.↗
NASA NTRS1989

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗
NASA NTRS2001

Atomic Data and Spectral Line Intensities for Ni XXI

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ni XXI. The configurations used are 2s(sup 2)2p(sup 4), 2s2p(sup 5), 2p(sup 6), 2s(sup 2)2p(sup 3)3s, and 2s(sup 2)3p(sup 3)3d giving rise to 58 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies, 85, 170, 255, 340, and 425 Ry. Excitation rate coefficients are calculated by assuming a Maxwellian electron velocity distribution at an electron temperature of log T(sub e)(K)=6.9, corresponding to maximum abundance of Ni XXI. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities 10(exp 8)-10(exp 14) per cubic centimeter. Relative spectral line intensities are calculated. Proton excitation rates between the lowest three levels have been included in the statistical equilibrium equations. The predicted intensity ratios are compared with available observations.

Bhatia, A. K.↗
NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗