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 415 records · Page 23

Infrared cooling rate calculations in operational general circulation models - Comparisons with benchmark computations

The performance of several parameterized models is described with respect to numerical prediction and climate research at GFDL, NCAR, and GISS. The radiation codes of the models were compared to benchmark calculations and other codes for the intercomparison of radiation codes in climate models (ICRCCM). Cooling rates and fluxes calculated from the models are examined in terms of their application to established general circulation models (GCMs) from the three research institutions. The newest radiation parameterization techniques show the most significant agreement with the benchmark line-by-line (LBL) results. The LBL cooling rates correspond to cooling rate profiles from the models, but the parameterization of the water vapor continuum demonstrates uncertain results. These uncertainties affect the understanding of some lower tropospheric cooling, and therefore more accurate parameterization of the water vapor continuum, as well as the weaker absorption bands of CO2 and O3 is recommended.

Kiehl, J. T.↗

Electrical conductivity and velocity of highly ionized plasma flows - Theory and experiment.

Use of an immersible, three-coil, magnetic-induction probe, previously tested in a low-density supersonic argon jet, to measure electrical conductivity and velocity profiles of a highly ionized high-density nitrogen jet in the continuum flow regime where effects due to probe bow shocks and boundary layers might not be negligible. Measured centerline values of electrical conductivity and velocity were compared with predictions based on a theoretical analysis previously developed to study the gas as it expanded adiabatically and inviscidly from an equilibrium sonic state to the nozzle exit. The resulting numerical exit plane values for electron density and electron temperature were then substituted into the Spitzer-Haerm conductivity formula to compute a theoretical conductivity which agreed within 40% of the measured conductivity, while the calculated and experimental velocity values differed by as much as 50%. The lack of agreement was attributed to the possible use of invalid assumptions and boundary conditions in the computer analysis or to the unknown effects of shocks on the probe data.

Vendell, E. W.↗

X-ray bumps, iron K-alpha lines, and X-ray suppression by obscuring tori in Seyfert galaxies

We investigate the X-ray spectral properties of unobscured type 1 and obscured type 2 Seyferts as predicted by the unified Seyfert scheme. We consider the reprocessing of X-ray photons by photoelectric absorption, iron fluorescence, and Compton downscattering in the obscuring tori surrounding these active nuclei, and compute by Monte Carlo methods the reprocessed spectra as a function of the viewing angle. Depending on the optical depth and shape of the torus, and on the viewing angle, the X-ray flux can be suppressed by substantial factors when our line of sight is obscured. We show that an immediate consequence of the existence of an obscuring thick torus is the production in the spectra of type 1 Seyfert galaxies of a bump in the continuum above 10-20 keV and an Fe K-alpha line with significant equivalent width. In those type 2 Seyferts for which the hard X-ray spectrum has been substantially suppressed, the equivalent width of the Fe K-alpha line in the transmitted spectrum can be very large.

Krolik, Julian H.↗

Aerodynamic characteristics of a hypersonic viscous optimized waverider at high altitudes

The present paper addresses the applicability of the basic concept of waveriding at high altitudes, and the extent to which the large viscous forces degrade the aerodynamic performance of waveriders. The waverider under consideration was designed using a continuum flow methodology. It is shown that the lift-to-drag ratio of high-altitude/high-Knudsen-number waveriders can be expected to be significantly lower than their low altitude/low Knudsen number counterparts. The aerodynamic performance of a representative waverider which was optimized for a 90-km, Mach-25 application is studied for altitudes ranging from 97 km to 145 km and incidence angles of 0 to 30 deg. Typical values of the lift-to-drag ratio were computed to be in the range of 0 to 0.3. Friction forces are mostly responsible for this poor performance. Friction forces account for more than 93 percent of the drag and significantly reduce lift.

Rault, Didier F. G.↗

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the space shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modeling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the Space Shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modelling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗

DRF: A Software Architecture for a Data Marketplace to Support Advanced Air Mobility

Advanced Air Mobility is a new aviation vision, where unmanned aerial systems will trans- port passengers and cargo across urban and rural areas. Critical to the realization of this vision is the development of a digital marketplace, which allows service providers and consumers operating in the airspace ecosystem to securely exchange data and reasoning insights. In this paper, we present the architecture of a decentralized data marketplace that connects data and reasoning service providers to vehicles and other service consumers along the cloud-to-edge continuum. We also present two example use cases to demonstrate the value of our approach.

Autonomy↗

Interactive Spectral Analysis and Computation (ISAAC)

Isaac is a task in the NSO external package for IRAF. A descendant of a FORTRAN program written to analyze data from a Fourier transform spectrometer, the current implementation has been generalized sufficiently to make it useful for general spectral analysis and other one dimensional data analysis tasks. The user interface for Isaac is implemented as an interpreted mini-language containing a powerful, programmable vector calculator. Built-in commands provide much of the functionality needed to produce accurate line lists from input spectra. These built-in functions include automated spectral line finding, least squares fitting of Voigt profiles to spectral lines including equality constraints, various filters including an optimal filter construction tool, continuum fitting, and various I/O functions.

Lytle, D. M.↗

Combined structures-controls optimization of lattice trusses

The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.

Balakrishnan, A. V.↗

Study of Plume Impingement Effects in the Lunar Lander Environment

Plume impingement effects from the descent and ascent engine firings of the Lunar Lander were analyzed in support of the Lunar Architecture Team under the Constellation Program. The descent stage analysis was performed to obtain shear and pressure forces on the lunar surface as well as velocity and density profiles in the flow field in an effort to understand lunar soil erosion and ejected soil impact damage which was analyzed as part of a separate study. A CFD/DSMC decoupled methodology was used with the Bird continuum breakdown parameter to distinguish the continuum flow from the rarefied flow. The ascent stage analysis was performed to ascertain the forces and moments acting on the Lunar Lander Ascent Module due to the firing of the main engine on take-off. The Reacting and Multiphase Program (RAMP) method of characteristics (MOC) code was used to model the continuum region of the nozzle plume, and the Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) was used to model the impingement results in the rarefied region. The ascent module (AM) was analyzed for various pitch and yaw rotations and for various heights in relation to the descent module (DM). For the ascent stage analysis, the plume inflow boundary was located near the nozzle exit plane in a region where the flow number density was large enough to make the DSMC solution computationally expensive. Therefore, a scaling coefficient was used to make the DSMC solution more computationally manageable. An analysis of the effectiveness of this scaling technique was performed by investigating various scaling parameters for a single height and rotation of the AM. Because the inflow boundary was near the nozzle exit plane, another analysis was performed investigating three different inflow contours to determine the effects of the flow expansion around the nozzle lip on the final plume impingement results.

Marichalar, Jeremiah↗

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗

Computational Study of Electron-Molecule Collisions Related to Low-Temperature Plasmas

Computational study of electron-molecule collisions not only complements experimental measurements, but can also be used to investigate processes not readily accessible experimentally. A number of ab initio computational methods are available for these types of calculations. Here we describe a recently developed technique, the finite element Z-matrix method, Analogous to the R-matrix, method, it partitions the space into regions and employs real matrix elements. However, unlike the implementation of the R-matrix method commonly used in atomic and molecular physics, the Z-matrix method is fully variational. In the present implementation, a mixed basis of finite elements and Gaussians is used to represent the continuum electron, thus offering full flexibility without imposing fixed boundary conditions. Numerical examples include the electron-impact dissociation of N2 via the metastable A3Su+ state, a process which may be important in the lower thermosphere, and the dissociation of the CF radical, a process of interest to plasma etching. To understand the dissociation pathways, large scale quantum chemical calculations have been carried out for all target states which dissociate to the lowest five limits in the case of N2, and to the lowest two limits in the case of CF. For N2, the structural calculations clearly show the preference for predissociation if the initial state is the ground X1 Sg+ state, but direct dissociation appears to be preferable if the initial state is the A3Su+ state. Multi-configuration SCF target functions are used in the collisional calculation.

Huo, Winifred M.↗

Gamma radiation and photospheric white-light flare continuum

It is noted that recent gamma-ray observations of solar flares have provided a better means for estimating the heating of the solar atmosphere by energetic protons. This type of heating has been suggested as the explanation of the continuum emission of the white-light flare. The effects on the photosphere of high-energy particles capable of producing the intense gamma-ray emission observed in the flare of July 11, 1978, are analyzed. A simple energy-balance argument is used, and hydrogen ionization is taken into account. It is found that energy deposition increases with height for the inferred proton spectra and is not strongly dependent upon the assumed angle of incidence. At the top of the photosphere, the computed energy inputs fall in the range 10-100 ergs/cu cm-s.

Hudson, H. S.↗

Cloud-particle galactic gas dynamics and star formation

Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.

Roberts, W. W., Jr.↗

Free-free absorption of infrared radiation in collisions of electrons with neutral rare-gas atoms

A relationship between the inverse bremsstrahlung absorption cross section and the electron neutral momentum transfer cross section has been utilized to determine the infrared free-free continuum absorption coefficient for the negative ions of helium, neon, argon, krypton, and xenon. The values of the momentum transfer cross section for this calculation have been obtained from experimental measurements. Analytical expressions for the absorption coefficient have also been developed. From the results of this calculation, it is possible to determine the absorption coefficient per unit electron density per neutral atom for temperatures in the range from 2500 to 25,000 K. The results are compared with those from tabulations of previous calculations and those computed from theoretical values of the phase shifts for the elastic scattering of electrons by neutral atoms.

Stallcop, J. R.↗

Determinants of the mechanical properties of bones

The mechanical properties of bones are governed by the same principles as those of man-made load-bearing structures, but the organism is able to adapt its bone structure to changes in skeletal loading. In this overview of the determinants of the strength and stiffness of bone, a continuum approach has been taken, in which the behavior of a macroscopic structure depends on its shape and size, and on the mechanical properties of the material within. The latter are assumed to depend on the composition (porosity and mineralization) and organization (trabecular or cortical bone architecture, collagen fiber orientation, fatigue damage) of the bone. The effects of each of these factors are reviewed. Also, the possible means of non-invasively estimating the strength or other mechanical properties of a bone are reviewed, including quantitative computed tomography, photon absorptiometry, and ultrasonic measurements. The best estimates of strength have been obtained with photon absorptiometry and computed tomography, which at best are capable of accounting for 90% of the strength variability in a simple in vitro test, but results from different laboratories have been highly variable.

Non-NASA Center↗

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

Interpretation of the prominence differential emissions measure for 3 geometries

Researchers have used prominence extreme ultraviolet line intensities observed from Skylab to derive the differential emission measure Q(T) in the prominence-corona (PC) interface from 3 x 10,000 to 3 times 1 million K, including the effects of Lyman Continuum absorption. Using lines both shortward and longward of the Lyman limit, researchers have estimated the importance of absorption as function of temperature. The magnitude of the absorption, as well as its rate of increase as a function of temperature, place limits on the thread scales and the character of the interfilar medium. Researchers have calculated models based on three assumed geometries: (1) threads with hot sheaths and cool cores; (2) isothermal threads; and (3) threads with longitudinal temperature gradients along the magnetic field. Comparison of the absorption computed from these models with the observed absorption in prominences shows that none of the geometries is totally satisfactory.

Schmahl, E. J.↗