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 199 records · Page 11

Planetary Nebulae with Supporting IR Data

We present new HST/WFPC2 imagery for the planetary nebula (PN) NGC 6818. Observations were made in line filters F437N, F487N, F502N, and F656N plus continuum filter F547M. The primary goal was to develop a high spatial resolution (approx. 0.1 in.) map of the intrinsic line ratio [O III] 4363/5007 and thereby evaluate the electron temperature (T(sub e)) and the mean-square T(sub e) variation (t(sup 2) across the PN. In this process we developed an extinction map from the F487N (H(beta)) and F656N (H(alpha)) images by comparing the observed line ratios in each pixel to the theoretical ratio and computing a c(H(beta)) map which was used to correct the observed 4363/5007 ratios for reddening. We also adjust for the continuum contribution to the line filter data. We present color-coded pictures of the reddening (c(H(beta))) map, the [O III] T(sub e) map, as well as our determinations of t(sup 2). The T(sub e) map shows a decline from approx. 14000 K in the inner regions to approx. 11000 K at the outer edge. Such a radial T(sub e) gradient is expected for a high-excitation nebula with a prominent He(++) zone such as NGC 6818. A composite of images taken in 3 filters (F656N, red; F487N, blue; and F502N, [O III] 5007, green) shows a roughly spherical outer envelope as well as a brighter vase-shaped interior "bubble". There is a prominent orifice to the North and a smaller one to the South, along the major axis, possibly caused by a blow-out from a fast wind. This nebula has an appearance remarkably similar to that of the PN NGC 3918 previously imaged with HST by H. Bond. We note from the continuum images (F547M) two stars in the nebular field that are fainter than the prominent central star; these are roughly 2-4 sec. N and NE of the central star. Further study is needed to establish whether or not there may be a physical association of either star with the central star.

Harrington, J. Patrick↗

Ablation and radiation coupled viscous hypersonic shock layers, volume 1

The results for a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are presented. The analysis includes flow field coupling with the ablator surface, equilibrium chemistry, a step-function diffusion model and a coupled line and continuum radiation calculation. This report serves as the documentation, i e. users manual and operating instructions for the computer programs listed in the report.

Engel, C. D.↗

Optical observables in stars with non-stationary atmospheres

Experience gained by use of Cepheid modeling codes to predict the dimensional and photometric behavior of nuclear fireballs is used as a means of validating various computational techniques used in the Cepheid codes. Predicted results from Cepheid models are compared with observations of the continuum and lines in an effort to demonstrate that the atmospheric phenomena in Cepheids are quite complex but that they can be quantitatively modeled.

Hillendahl, R. W.↗

Numerical simulation of solidification in a horizontal cylindrical annulus charged with an aqueous salt solution

A computational study of solidification of a binary Na2CO3 solution in a horizontal cylindrical annulus is performed using a continuum formulation with a control-volume based, finite-difference scheme. The initial conditions were selected to facilitate the study of counter thermal and solutal convection, accompanied by extensive mushy region growth. Numerical results are compared with experimental data with mixed success. Qualitative agreement is obtained for the overall solidification process and associated physical phenomena. However, the plume thickness calculated for the solutally-driven convective upflow is substantially smaller than the observed value. Evolution of double-diffusive layers is predicted, but over a time scale much smaller than that observed experimentally. Good agreement is obtained between predicted and measured results for solid growth, but the mushy region thickness is significantly overpredicted.

Neilson, D. G.↗

The wavelength dependence of broadband linear polarization in cool stars

A model is developed to investigate the effect of wavelength dependence of magnetically generated broadband linear polarization (BLP) in late-type stars as a function of line blanketing, observation filters, and other variables. Linear polarization is assumed to be a single 'average' line and the line opacity, magnetic area filling factor, and continuum intensity are incorporated to establish an equation for disk-generated BLP. The method is investigated by computing special models for the flare star BY Dra and the sun which are compared to previous calculations. The model is shown to agree with the calculations of BLP, and the magnetic BLP dependent on wavelengths varies significantly from the simple power law dependences on wavelength. The model sets lower limits to the area occupied by magnetic regions, and an upper limit to the expected BLP also results from the formulation.

Saar, Steven H.↗

Accretion dynamics and polarized X-ray emission of magnetized neutron stars

The basic ideas of accretion onto magnetized neutron stars are outlined. These are applied to a simple model of the structure of the plasma mound sitting at the magnetic poles of such a star, in which upward diffusion of photons is balanced by their downward advection. This steady flow model of the plasma's dynamical state is used to compute the emission of polarized X-raysfrom the optically thick, birefringent medium. The linear polarization of the continuum radiation emerging from the quasi-static mound is found to be as much as 40 percent at some rotation phases, but is insensitive to the geometry of the accretion flow. The role of the accretion shock, whose detailed polarimetric and spectral characteristics have yet to be calculated, is emphasized as the final determinant of the properties of the emerging X-rays. Some results describing the fully time dependent dynamics of the flow are also presented. In particular, steady flow onto a neutron star is shown to exhibit formation of 'photon bubbles', regions of greatly reduced plasma density filled with radiation which form and rise on millisecond time scale. The possible role of these complex structures in the flow for the formation of the emergent spectrum is briefly outlined.

Arons, Jonathan↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗

Comparison of shock structure solutions using independent continuum and kinetic theory approaches

A vehicle traversing the atmosphere will experience flight regimes at high altitudes in which the thickness of a hypersonic shock wave is not small compared to the shock standoff distance from the hard body. When this occurs, it is essential to compute accurate flow field solutions within the shock structure. In this paper, one-dimensional shock structure is investigated for various monatomic gases from Mach 1.4 to Mach 35. Kinetic theory solutions are computed using the Direct Simulation Monte Carlo method. Steady-state solutions of the Navier-Stokes equations and of a slightly truncated form of the Burnett equations are determined by relaxation to a steady state of the time-dependent continuum equations. Monte Carlo results are in excellent agreement with published experimental data and are used as bases of comparison for continuum solutions. For a Maxwellian gas, the truncated Burnett equations are shown to produce far more accurate solutions of shock structure than the Navier-Stokes equations.

Fiscko, Kurt A.↗

The ultraviolet spectrum of Herbig-Haro object 2H

IUE spectra of Herbig-Haro object 2H are presented. The spectra show a strong 'excess' UV continuum and prominent emission lines of C, N, O, Si, Mg, and possibly Al. The continuum, F(lambda), exhibits a turnover shortward of about 1450 A, confirming for the first time the H0 two-photon nature of the emission source. A possible absorption feature near 1680 A, which could result from a new grain or molecular constituent in these protostellar objects is also noted. Recently computed models of steady shocks into partially ionized gas reproduce the two-photon spectral shape, but its observed intensity relative to H-beta and the Balmer continuum is anomalously high. It is suggested that a range of shock velocities, 70-100 km/s, or nonsteady, 'truncated' shocks may be responsible. Future high-sensitivity UV observations of HH objects may be used to probe grain extinction curves in star-forming regions.

Brugel, E. W.↗

Input-output characterization of fiber composites by SH waves

Input-output characterization of fiber composites is studied theoretically by tracing SH waves in the media. A fiberglass epoxy composite is modeled as a homogeneous transversely isotropic continuum plate. The reflection of an SH wave at a stress-free plane boundary in a semi-infinite transversely isotropic medium is considered first. It is found that an incident SH wave reflects only a similar SH wave back into the medium. It is also established that the angle of reflection of the reflected wave is equal to the angle of incidence of the incident wave. The phase velocity of the SH waves and the delay time of the SH waves in reaching the receiving transducer are computed as functions of a reflection index, defined as the number of reflections of the SH waves from the bottom face of the continuum plate. The directivity function corresponding to the shear stress associated with the SH waves in the continuum plate is also derived as a function of the reflection index. A theoretical output voltage from the receiving transducer is calculated for a tone burst (a periodic input voltage of finite duration). The output voltage is shown for tone bursts of duration 60 microseconds and center frequencies of 0.75, 1.00, and 1.25 MHz. The study enhances the quantitative and qualitative understanding of the nondestructive evaluation (NDE) of fiber composites which can be modeled as transversely isotropic media.

Renneisen, John D.↗

Monte Carlo study of vibrational relaxation processes

A new model is proposed for the computation of vibrational nonequilibrium in the direct simulation Monte Carlo method (DSMC). This model permits level to level vibrational transitions for the first time in a Monte Carlo flowfield simulation. The model follows the Landau-Teller theory for a harmonic oscillator in which the rates of transition are related to an experimental correlation for the vibrational relaxation time. The usual method for simulating such processes in the DSMC technique applies a constant exchange probability to each collision and the vibrational energy is treated as a continuum. A comparison of these two methods is made for the flow of nitrogen over a wedge. Significant differences exist for the vibrational temperatures computed. These arise as a consequence of the incorrect application of a constant exchange probability in the old method. It is found that the numerical performances of the two vibrational relaxation models are equal.

Boyd, Iain D.↗

Computational strategy for the solution of large strain nonlinear problems using the Wilkins explicit finite-difference approach

The STEALTH code system, which solves large strain, nonlinear continuum mechanics problems, was rigorously structured in both overall design and programming standards. The design is based on the theoretical elements of analysis while the programming standards attempt to establish a parallelism between physical theory, programming structure, and documentation. These features have made it easy to maintain, modify, and transport the codes. It has also guaranteed users a high level of quality control and quality assurance.

Hofmann, R.↗

Development of a moderately sized finite element program for nonlinear structural analysis

AGGIE 1 is a computer program for predicting the linear and nonlinear, static and dynamic structural response of two- and three-dimensional continuum solids. The program is based on isoparametric finite elements and allows for 2-D plane stress, plane strain, and axisymmetric analyses and general 3-D analyses. Large strain kinematics is based on the total Lagrangian formulation. Materially nonlinear models include several elastic-plastic work-hardening models as well as an incompressible Mooney-Rivlin model. Included in this report is a brief description of the theoretical bases of the program, the material models used, the element library and the overall program organization. Instructions for data input preparation are given in detail. Several sample problems are given along with the required program input and program generated solutions.

Haisler, W. E.↗

Accurate Determination of Tranverse Normal Stresses in Sandwich Panels Subjected to Thermomechanical Loadings

A two-stage computational procedure is presented for the accurate determination of transverse normal stresses in sandwich panels subjected to thermomechanical loadings. The procedure involves the use of a first-order shear deformation model in the first stage, and an iterational process for successive improvement of the accuracy of the displacement and stress fields in the second stage. The effectiveness of the procedure is demonstrated by means of numerical studies (if thin and moderately thick flat rectangular sandwich panels. Two sets of boundary conditions are considered; namely, one with all edges simply supported. and the other with two opposite edges simply supported and the remaining two clamped. The displacement components and the transverse shear and normal stresses obtained by the proposed computational procedure are found to be in close agreement with the Solutions of the three-dimensional (3-D) continuum models.

Noor, Ahmed K.↗

Micro-Mechanical Modeling of Ductile Fracture in Welded Aluminum-Lithium Alloys

This study describes a 3-D computational framework to model stable extension of a macroscopic crack under mode I conditions in ductile metals. Material separation occurs through a local fracture mechanism based on the growth and coalescence of micro voids dispersed in the material. The Gurson-Tvergaard plasticity model for void materials describes the eventual loss of material stress capacity under sustained loading. Fixed-sized computational cell elements defined on a thin layer adjacent to the crack plane provide an implicit length scale for continuum damage model. Outside this layer, the material remains undamaged by void growth, consistent with metallurgical observations. The key micro-mechanics parameters are D, the thickness of cell layer, and f 0 , the initial cell porosity.

Ahmed Ibrahim↗

Integrated Force Method for Indeterminate Structures

Two methods of solving indeterminate structural-mechanics problems have been developed as products of research on the theory of strain compatibility. In these methods, stresses are considered to be the primary unknowns (in contrast to strains and displacements being considered as the primary unknowns in some prior methods). One of these methods, denoted the integrated force method (IFM), makes it possible to compute stresses, strains, and displacements with high fidelity by use of modest finite-element models that entail relatively small amounts of computation. The other method, denoted the completed Beltrami Mitchell formulation (CBMF), enables direct determination of stresses in an elastic continuum with general boundary conditions, without the need to first calculate displacements as in traditional methods. The equilibrium equation, the compatibility condition, and the material law are the three fundamental concepts of the theory of structures. For almost 150 years, it has been commonly supposed that the theory is complete. However, until now, the understanding of the compatibility condition remained incomplete, and the compatibility condition was confused with the continuity condition. Furthermore, the compatibility condition as applied to structures in its previous incomplete form was inconsistent with the strain formulation in elasticity.

Hopkins, Dale A.↗

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields, volume 3

The computer programs developed to calculate the shock wave precursor and the method of using them are described. This method calculated the precursor flow field in a nitrogen gas including the effects of emission and absorption of radiation on the energy and composition of gas. The radiative transfer is calculated including the effects of absorption and emission through the line as well as the continuum process in the shock layer and through the continuum processes only in the precursor. The effects of local thermodynamic nonequilibrium in the shock layer and precursor regions are also included in the radiative transfer calculations. Three computer programs utilized by this computational scheme to calculate the precursor flow field solution for a given shock layer flow field are discussed.

Carlson, Leland A.↗