Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Aerodynamics for the Mars Phoenix Entry Capsule

Pre-flight aerodynamics data for the Mars Phoenix entry capsule are presented. The aerodynamic coefficients were generated as a function of total angle-of-attack and either Knudsen number, velocity, or Mach number, depending on the flight regime. The database was constructed using continuum flowfield computations and data from the Mars Exploration Rover and Viking programs. Hypersonic and supersonic static coefficients were derived from Navier-Stokes solutions on a pre-flight design trajectory. High-altitude data (free-molecular and transitional regimes) and dynamic pitch damping characteristics were taken from Mars Exploration Rover analysis and testing. Transonic static coefficients from Viking wind tunnel tests were used for capsule aerodynamics under the parachute. Static instabilities were predicted at two points along the reference trajectory and were verified by reconstructed flight data. During the hypersonic instability, the capsule was predicted to trim at angles as high as 2.5 deg with an on-axis center-of-gravity. Trim angles were predicted for off-nominal pitching moment (4.2 deg peak) and a 5 mm off-axis center-ofgravity (4.8 deg peak). Finally, hypersonic static coefficient sensitivities to atmospheric density were predicted to be within uncertainty bounds.

Edquist, Karl T.↗
NASA NTRS2012

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗
NASA NTRS2000

Numerical Study of Rarefied Hypersonic Flow Interacting with a Continuum Jet

An uncoupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC) technique is developed and applied to provide solutions for continuum jets interacting with rarefied external flows. The technique is based on a correlation of the appropriate Bird breakdown parameter for a transitional-rarefied condition that defines a surface within which the continuum solution is unaffected by the external flow-jet interaction. The method is applied to two problems to assess and demonstrate its validity; one of a jet interaction in the transitional-rarefied flow regime and the other in the moderately rarefied regime. Results show that the appropriate Bird breakdown surface for uncoupling the continuum and non-continuum solutions is a function of a non-dimensional parameter relating the momentum flux and collisionality between the two interacting flows. The correlation is exploited for the simulation of a jet interaction modeled for an experimental condition in the transitional-rarefied flow regime and the validity of the correlation is demonstrated. The uncoupled technique is also applied to an aerobraking flight condition for the Mars Global Surveyor spacecraft with attitude control system jet interaction. Aerodynamic yawing moment coefficients for cases without and with jet interaction at various angles-of-attack were predicted, and results from the present method compare well with values published previously. The flow field and surface properties are analyzed in some detail to describe the mechanism by which the jet interaction affects the aerodynamics.

Glass, Christopher E.↗
NASA NTRS2002

Comparison of a 3-D CFD-DSMC Solution Methodology With a Wind Tunnel Experiment

A solution method for problems that contain both continuum and rarefied flow regions is presented. The methodology is applied to flow about the 3-D Mars Sample Return Orbiter (MSRO) that has a highly compressed forebody flow, a shear layer where the flow separates from a forebody lip, and a low density wake. Because blunt body flow fields contain such disparate regions, employing a single numerical technique to solve the entire 3-D flow field is often impractical, or the technique does not apply. Direct simulation Monte Carlo (DSMC) could be employed to solve the entire flow field; however, the technique requires inordinate computational resources for continuum and near-continuum regions, and is best suited for the wake region. Computational fluid dynamics (CFD) will solve the high-density forebody flow, but continuum assumptions do not apply in the rarefied wake region. The CFD-DSMC approach presented herein may be a suitable way to obtain a higher fidelity solution.

Glass, Christopher E.↗
NASA NTRS2009

A Model for Predicting Thermoelectric Properties of Bi2Te3

A parameterized orthogonal tight-binding mathematical model of the quantum electronic structure of the bismuth telluride molecule has been devised for use in conjunction with a semiclassical transport model in predicting the thermoelectric properties of doped bismuth telluride. This model is expected to be useful in designing and analyzing Bi2Te3 thermoelectric devices, including ones that contain such nano - structures as quantum wells and wires. In addition, the understanding gained in the use of this model can be expected to lead to the development of better models that could be useful for developing other thermoelectric materials and devices having enhanced thermoelectric properties. Bi2Te3 is one of the best bulk thermoelectric materials and is widely used in commercial thermoelectric devices. Most prior theoretical studies of the thermoelectric properties of Bi2Te3 have involved either continuum models or ab-initio models. Continuum models are computationally very efficient, but do not account for atomic-level effects. Ab-initio models are atomistic by definition, but do not scale well in that computation times increase excessively with increasing numbers of atoms. The present tight-binding model bridges the gap between the well-scalable but non-atomistic continuum models and the atomistic but poorly scalable ab-initio models: The present tight-binding model is atomistic, yet also computationally efficient because of the reduced (relative to an ab-initio model) number of basis orbitals and flexible parameterization of the Hamiltonian.

Lee, Seungwon↗
NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗
NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗
NASA NTRS1987

Unified viscoplastic constitutive equations and their applications

Unified constitutive equations for time- and temperature-dependent metallic plastic deformation have been applied in FEM simulations of forming processes; increasingly powerful computational tools and physical models are being used to numerically model complex engineering problems. Once confidence has been gained through adequate verification, these numerical models will increasingly replace experimental models. Attention is presently given to the contributions made by physical metallurgy, continuum mechanics, and computational mechanics.

Lindholm, U. S.↗
NASA NTRS1986

Composition of faint comets

The study uses an emission line, differential imaging camera built by the Science Operations Branch. This instrument allows photometric data to be obtained over a large area of a comet in a large number of resolution elements. The detector is a 100x100 Reticon array which with interchangeable optics can give resolutions from 2" to 30" over a field of 1' to 15'. The camera through its controlling computer can simultaneously take images in on-line and continuum filters and through computer subtraction and calibration present a photometric image of the comet produced by only the emission of the molecule under study. Initial work has shown two significant problems. First the auxiliary equipment of the telescope has not allowed the unambiguous location of faint comets so that systematic observations could be made, and secondly initial data has not shown much molecular emission from the faint comets which were located. Work last year on a software and hardware display system and this year on additional guide motors on the 36-inch telescope has allowed the differential camera to act as its own finder and guide scope. Comet IRAS was observed in C2 and CO+, as well as an occultation by the comet of SAO029103. The perodic comet Giacobini-Zinner was also observed in C2.

Brown, L. W.↗
NASA NTRS1979

Continuum Absorption Coefficient of Atoms and Ions

The rate of heat transfer to the heat shield of a Jupiter probe has been estimated to be one order of magnitude higher than any previously experienced in an outer space exploration program. More than one-third of this heat load is due to an emission of continuum radiation from atoms and ions. The existing computer code for calculating the continuum contribution to the total load utilizes a modified version of Biberman's approximate method. The continuum radiation absorption cross sections of a C - H - O - N ablation system were examined in detail. The present computer code was evaluated and updated by being compared with available exact and approximate calculations and correlations of experimental data. A detailed calculation procedure, which can be applied to other atomic species, is presented. The approximate correlations can be made to agree with the available exact and experimental data.

Armaly, B. F.↗
NASA NTRS1990

Bipropellant engine plume analyses for the Mars Observer scientific spacecraft

This paper discusses the methodology utilized to characterize exhaust plume expansion of the 490 N bipropellant engine and provides an analysis of exhaust plune impingement effects on the Mars Observer oxidizer tank structure. Nozzle and plume continuum flowfields are computed with a method of characteristics solution while transitional and rarefied regions of the exhaust gas flow were defined by a direct simulation Monte Carlo method. Results of this study provide sufficient data for an assessment of spacecraft thermal impacts due to plume heating effects and thrust losses from exhaust plume impingement.

Rattenni, L., Jr.↗
NASA NTRS1991

Infrared emission from isolated dust clouds in the presence of very small dust grains

Models of the effects of small grain-generated temperature fluctuations on the IR spectrum and surface brightness of externally heated interstellar dust clouds are presently constructed on the basis of a continuum radiation transport computer code which encompasses the transient heating of small dust grains. The models assume a constant fractional abundance of large and small grains throughout the given cloud. A comparison of model results with IRAS observations indicates that the observed 12-25 micron band emissions are associated with about 10-A radius grains, while the 60-100 micron emission is primarily due to large grains which are heated under the equilibrium conditions.

Lis, Dariusz C.↗
NASA NTRS1994

Thermal continua of AGN accretion disks

We have computed the thermal continuum energy distribution of thermal radiation from the atmospheres of supermassive accretion disks around supermassive black holes. Non-LTE radiative transfer is combined with a model of the vertical structure at each radius appropriate to the low effective gravities of these disks. Locally, the Lyman edge of H can be in emission or absorption. When the emission is summed over the disk with Doppler and gravitational redshifts taken into account, the observed continuum typically shows little sign of a discontinuity near the Lyman edge. For relatively cool disks, the Lyman edge is in absorption, but it appears as a slope change extending over several hundred angstroms, rather than an abrupt discontinuity. Disks around Kerr black holes can explain the observed range of soft X-ray luminosities of AGN, but disks around Schwarzschild holes are much too faint in soft X-rays.

Shields, G. A.↗
NASA NTRS1999

Does UV CETI Suffer from the MAD Syndrome?

Photometric data have been analysed and searched for events of flaring and other variability. Some flaring has been detected, though probably not at a level that will hinder our continuing spectral analysis. X-ray diagnostics for the very hot coronal emission measure are under investigation in order to determine whether or not the very hot coronal plasma contributes significantly to the observed X-ray flux in the EUV. The key test of the MAD syndrome lies in whether or not the coronal lines indicate a depletion in met- als in the corona relative to the underlying photosphere. While some progress has been made in this direction, not as much work has been accomplished as expected due to the increasing commitments of the PI to the CXO project as it nears launch. A no-cost extension has been granted in order to further the analysis and carry out the next stages of the investigation: to construct an emission measure distribution with which to compute a synthetic continuum to compare with the observed continuum.

Drake, Jeremy↗
NASA NTRS2005

Multi-Scale Modeling of a Graphite-Epoxy-Nanotube System

A multi-scale method is utilized to determine some of the constitutive properties of a three component graphite-epoxy-nanotube system. This system is of interest because carbon nanotubes have been proposed as stiffening and toughening agents in the interlaminar regions of carbon fiber/epoxy laminates. The multi-scale method uses molecular dynamics simulation and equivalent-continuum modeling to compute three of the elastic constants of the graphite-epoxy-nanotube system: C11, C22, and C33. The 1-direction is along the nanotube axis, and the graphene sheets lie in the 1-2 plane. It was found that the C11 is only 4% larger than the C22. The nanotube therefore does have a small, but positive effect on the constitutive properties in the interlaminar region.

Frankland, S. J. V.↗
NASA NTRS2018

Methane on Mars: New Insights into the Sensitivity of CH4 with the NOMAD/ExoMars Spectrometer Through Its First In-Flight Calibration

The Nadir and Occultation for MArs Discovery instrument (NOMAD), onboard the ExoMars Trace Gas Orbiter (TGO) spacecraft was conceived to observe Mars in solar occultation, nadir, and limb geometries, and will be able to produce an outstanding amount of diverse data, mostly focused on properties of the atmosphere. The infrared channels of the instrument operate by combining an echelle grating spectrometer with an Acousto- Optical Tunable Filter (AOTF). Using in-flight data, we characterized the instrument performance and parameterized its calibration. In particular: an accurate frequency calibration was achieved, together with its variability due to thermal effects on the grating. The AOTF properties and transfer function were also quantified, and we developed and tested a realistic method to compute the spectral continuum transmitted through the coupled grating and AOTF system. The calibration results enabled unprecedented insights into the important problem of the sensitivity of NOMAD to methane abundances in the atmosphere. We also deeply characterized its performance under realistic conditions of varying aerosol abundances, diverse albedos and changing illumination conditions as foreseen over the nominal mission. The results show that, in low aerosol conditions, NOMAD single spectrum, 1σ sensitivity to CH4 is around 0.33 ppbv at 20 km of altitude when performing solar occultation's, and better than 1 ppbv below 30 km. In dusty conditions, we show that the sensitivity drops to 0 below 10 km. In Nadir geometry, results demonstrate that NOMAD will be able to produce seasonal maps of CH4 with a sensitivity around 5 ppbv over most of planet's surface with spatial integration over 5×5° bins. Results show also that such numbers can be improved by a factor of ~10 to ~30 by data binning. Overall, our results quantify NOMAD's capability to address the variable aspects of Martian climate.

Liuzzi, Giuliano↗
NASA NTRS2018

Comparison of CFD and DSMC Using Calibrated Transport Parameters

Hypersonic re-entry flows span a wide range of length scales where regions of both rarefied and continuum flow exist. Traditional computational fluid dynamics (CFD) techniques do not provide an accurate solution for the rarefied regions of such ‘mixed’ flow fields. Although direct simulation Monte Carlo (DSMC) can be used to accurately capture both the continuum and rarefied features of ‘mixed’ flow fields, they are computationally expensive when employed to simulate the low Knudsen number continuum regimes. Thus, a hybrid framework for seamlessly combining the two methodologies, CFD and DSMC, continues to be a topic of significant research effort. Ensuring consistency in the reaction kinetics and transport models employed within CFD and DSMC is a crucial requirement for obtaining a reliable solution from a hybrid framework for combined continuum/rarefied high speed flows. This paper represents one of the first studies to utilize the calibrated transport parameters developed to ensure consistency between CFD and DSMC solvers. The new variable soft sphere (VSS) parameters are compared to both previous “standard” variable hard sphere (VHS) parameters and also to solutions from the CFD transport properties that the new parameters were developed to reproduce.

Liechty, D. S.↗