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 145 records · Page 8

Finite Element Aircraft Simulation of Turbulence

A turbulence model has been developed for realtime aircraft simulation that accommodates stochastic turbulence and distributed discrete gusts as a function of the terrain. This model is applicable to conventional aircraft, V/STOL aircraft, and disc rotor model helicopter simulations. Vehicle angular activity in response to turbulence is computed from geometrical and temporal relationships rather than by using the conventional continuum approximations that assume uniform gust immersion and low frequency responses. By using techniques similar to those recently developed for blade-element rotor models, the angular-rate filters of conventional turbulence models are not required. The model produces rotational rates as well as air mass translational velocities in response to both stochastic and deterministic disturbances, where the discrete gusts and turbulence magnitudes may be correlated with significant terrain features or ship models. Assuming isotropy, a two-dimensional vertical turbulence field is created. A novel Gaussian interpolation technique is used to distribute vertical turbulence on the wing span or lateral rotor disc, and this distribution is used to compute roll responses. Air mass velocities are applied at significant centers of pressure in the computation of the aircraft's pitch and roll responses.

McFarland, R. E.↗

Energy transfer by radiation in non-grey atomic gases in isothermal and non-isothermal slabs

A multiband model for the absorption coefficient of atomic hydrogen-helium plasmas is constructed which includes continuum and line contributions. Emission from 28 stronger lines of 106 that have been screened is considered, of which 21 are from hydrogen and 7 belong to helium, with reabsorption due to line-line, line-continuum overlap accurately accounted for. The model is utilized in the computation of intensities and fluxes from shock-heated slabs of 85% H2-15% He mixtures for slab thicknesses from 1 to 30 cm, temperature from 10,000 to 20,000 K, and for different densities. In conjunction with the multiband model, simple numerical schemes have been devised which provide a quick and comprehensive way of computing radiative energy transfer in nonisothermal and nongrey gases.

Poon, P. T. Y.↗

17-25 micrometer spectra of Jupiter and Saturn

Ground-based spectra of Jupiter and Saturn in the region from 17 to 25 microns are compared with intensities computed from current thermal-structure models. Good agreement with the continuum of Jupiter is obtained for models which incorporate a temperature inversion, but published models give disagreement with the continuum of Saturn. Upper-limit abundances for sulfur and phosphorus in the Jovian atmosphere, as thermodynamically stable species S8 and P4, are found to be 0.04 and 2.1 times the solar abundance, respectively.

Tokunaga, A.↗

Photometry of Comet Halley

Observations of Comet Giacobini-Zinner were obtained on one night and of Comet Halley of three nights. When the comet becomes high enough it is hoped to obtain further observations of Comet Halley. Considerable progress was made in developing a probabalistic model for the emission of neutral dust by the nucleus of the comet. A preliminary comparison between the model and the observed intensity profiles in continuum filters is made. The model can also be used to compute the expected impact rates of dust on spacecraft flying past the comet.

Neff, J. S.↗

Anisotropic line emission and the geometry of the broad-line region in active galactic nuclei

Attention is given to an escape probability-based formalism, including the interaction with a background opacity, which is used to predict line-beaming factors for broad-line region (BLR) clouds of AGN. Several computational details necessary to model gas exposed to an intense nonthermal continuum are outlined. It is found that all hydrogen lines are strongly inwardly beamed, and that C IV 1549 is moderately beamed. The implications of this for the geometry of the emission-line regions of AGN are discussed. Beaming is found to cause the source-cloud separation to be overestimated by variability studies (so that the radiation field at the BLR is even more intense than previously thought); it can even cause spherical geometries to have little instantaneous response to continuum changes.

Ferland, G. J.↗

Pre- and postperihelon abundances of gas and dust in comet Halley

Photometrically calibrated spectra of comet P/Halley (1986 III) were recorded between 1985 September 12-1986 June 10 using the Ohio State University Image Dissector Scanner on the Perkins 72 inch telescope at the Lowell Observatory. Column densities of CN, C3, CH, C2, and NH2 were calculated from measured fluxes in these spectra, and molecular scale lengths were deduced from the radial distribution of CN, C3, C2, and NH2. Production rates were computed using the new scale lengths and a Haser model analysis. Continuum emission at 4260 A was used to derive gas-to-dust ratios. The data indicate than comet Halley was approximately 2-5 times more abundant in gas and dust at postperihelion than preperihelion. On 1986 June 8 we observed the onset of a cometary ourburst which appeared very strong in dust production. The gas-to-dust ratios appeared to be subject to changes as a result of short-term outbursts but otherwise did not exhibit any systematic dependence on heliocentric distance. Reflectivity gradients of the continuum were also measured from the spectra. While most of the continua were red, blue continua were also observed which may be correlated with dust outbursts.

Womack, Maria↗

Calculation of Free-Free Opacities

Free-free absorption is an important contribution to the opacity for radiation transport through hot materials Temperatures can be as high as several keV, such that it becomes a computational challenge to solve the Schrodinger equation efficiently for rapidly oscillating continuum functions for high angular momenta. Several groups\footnots, including ours, have studied the phase amplitude solution (PAS) of the Schrodinger equation, in which one solves equations for the wave function amplitude and phase, which are: smooth functions of the electron energy. It is also important to have an accurate Schroudinger benchmark for the development of the PAS method. We present results for dipole matrix elements, Gaunt factors, and cross sections for the absorption of radiation at various energies for Cs XIX at temperature=100 eV and density=0.187 g/cc for our newly developed PAS and Schrodinger benchmark.

Bhatia, A. K.↗

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗

Direct integration transmittance model

A transmittance model was developed for the 200-2000/cm region for interpretation of high spectral resolution measurements of laboratory absorption and of planetary thermal emission. The high spectral resolution requires transmittances to be computed monochromatically by summing the contribution of individual molecular absorption lines. A magnetic tape atlas of H2O,O3, and CO2 molecular line parameters serves as input to the transmittance model with simple empirical representations used for continuum regions wherever suitable laboratory data exist. The theoretical formulation of the transmittance model and the computational procedures used for the evaluation of the transmittances are discussed. Application is demonstrated of the model to several homogenous path laboratory absorption examples.

Kunde, V. G.↗

Direct integration transmittance model

A transmittance model has been developed for interpretation of high spectral resolution measurements of laboratory absorption and of planetary thermal emission. The high spectral resolution requires transmittances to be computed monochromatically by summing the contribution of individual molecular absorption lines. A magnetic tape atlas of H2O, O3, and CO2 molecular line parameters serves as input to the transmittance model with simple empirical representations used for continuum regions wherever suitable laboratory data exist. The theoretical formulation of the transmittance model and the computational procedures used for the evaluation of the transmittances are discussed, and application of the model to several homogeneous-path laboratory absorption examples is demonstrated.

Kunde, V. G.↗

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗

Hypersonic flow - Physical issues and computational simulation

Computational results are presented on hypersonic flows, with emphasis on thermal and chemical nonequilibrium effects and the role played by grids. The continuum equations and their boundary conditions are presented, and the direct simulation Monte Carlo method is described. The importance of using the correct computational and physical models in describing hypersonic nonequilibrium flows is demonstrated.

Hassan, H. A.↗

Theoretical Studies of Starburst Infrared Emission: Luminosity Indicators in Dusty Photoionized Environments

The luminosity of the central source in ionizing radiation is an essential parameter in a photoionized environment and is one of the most fundamental physical quantities one can measure. We outline a method of determining the luminosity for any emission-line region using only infrared data. In dusty environments, grains compete with hydrogen in absorbing continuum radiation. Grains produce infrared emission, and hydrogen produces recombination lines. We have computed a very large variety of photoionization models, using ranges of abundances, grain mixtures, ionizing continua, densities, and ionization parameters. The conditions were appropriate for such diverse objects as H(II) regions, planetary nebulae, starburst galaxies, and the narrow- and broad-line regions of active nuclei. The ratio of the total thermal grain emission relative to H-Beta (IR/H-Beta) is the primary indicator of whether the cloud behaves as a classical Stroemgren sphere (a hydrogen-bounded nebula) or whether grains absorb most of the incident continuum (a dust-bounded nebula). We find two global limits: when IR/H-Beta < 100, infrared recombination lines determine the source luminosity in ionizing photons; when IR/H-Beta >> 100, the grains act as a bolometer to measure the luminosity.

Bottorff, Mark↗

Polarimetric observations of Comet Kohoutek (1973f)

Radio, optical, and visual observations were made of Comet Kohoutek (1973f) using facilities of the National Radio Astronomy Observatory (NRAO) and of the Hawaii Night Sky Observatory, State University of New York at Albany (SUNYA). Transit observations were made using the NRAO 300 ft radio telescope on six consecutive days beginning November 24, 1973. Observations were carried out in the continuum at 11 cm using the 4-feed system at different position angles. The data were recorded on magnetic tape and were computer reduced. Results were negative, with an rms error of 62 milli flux units, giving an upper limit of 180 mfu for the 11 cm continuum emission from the Comet. This relatively high detection limit resulted from a combination of interference and bad weather that affected these pre-perihelion observations. No observations could be scheduled near or after perihelion.

Weinberg, J. L.↗

Assessment of impingement effects in the isentropic core of a small satellite control thruster plume

Two computational techniques commonly employed in the calculation of rocket and thruster expansion plumes are assessed. These are the method of characteristics (MOC), which is derived from the continuum Euler equations, and the direct simulation Monte Carlo (DSMC) method, which adopts a discrete particle approach. These techniques vary both in the computational expense and in the accuracy and detail of the solutions that they provide, depending upon the regime of application. The assessment is made with reference to the plume expanding from a small monopropellant hydrazine thruster and concentrates on the isentropic core of the jet for the flow regime lying between the continuum and free molecular limits. It is found that the more numerically intensive DSMC method offers the better correspondence to the available experimental data. In addition, large differences in typical impingement effects such as drag force and heat transfer are found at the free molecular limit of the plume expansion for the two predictive techniques. It is concluded that accurate estimation of impingement potential may only be achieved through application of the discrete particle method.

Boyd, I. D.↗

Space-Time Conservation Element and Solution Element Method Being Developed

The engineering research and design requirements of today pose great computer-simulation challenges to engineers and scientists who are called on to analyze phenomena in continuum mechanics. The future will bring even more daunting challenges, when increasingly complex phenomena must be analyzed with increased accuracy. Traditionally used numerical simulation methods have evolved to their present state by repeated incremental extensions to broaden their scope. They are reaching the limits of their applicability and will need to be radically revised, at the very least, to meet future simulation challenges. At the NASA Lewis Research Center, researchers have been developing a new numerical framework for solving conservation laws in continuum mechanics, namely, the Space-Time Conservation Element and Solution Element Method, or the CE/SE method. This method has been built from fundamentals and is not a modification of any previously existing method. It has been designed with generality, simplicity, robustness, and accuracy as cornerstones. The CE/SE method has thus far been applied in the fields of computational fluid dynamics, computational aeroacoustics, and computational electromagnetics. Computer programs based on the CE/SE method have been developed for calculating flows in one, two, and three spatial dimensions. Results have been obtained for numerous problems and phenomena, including various shock-tube problems, ZND detonation waves, an implosion and explosion problem, shocks over a forward-facing step, a blast wave discharging from a nozzle, various acoustic waves, and shock/acoustic-wave interactions. The method can clearly resolve shock/acoustic-wave interactions, wherein the difference of the magnitude between the acoustic wave and shock could be up to six orders. In two-dimensional flows, the reflected shock is as crisp as the leading shock. CE/SE schemes are currently being used for advanced applications to jet and fan noise prediction and to chemically reacting flows.

Chang, Sin-Chung↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗