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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1992

New optimality criteria methods - Forcing uniqueness of the adjoint strains by corner-rounding at constraint intersections

In new, iterative continuum-based optimality criteria (COC) methods, the strain in the adjoint structure becomes non-unique if the number of active local constraints is greater than the number of design variables for an element. This brief note discusses the use of smooth envelope functions (SEFs) in overcoming economically computational problems caused by the above non-uniqueness.

Rozvany, G. I. N.↗

NASA NTRS2003

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗

NASA NTRS2006

Cross Sections for Electron Impact Excitation of Astrophysically Abundant Atoms and Ions

Electron collisional excitation rates and transition probabilities are important for computing electron temperatures and densities, ionization equilibria, and for deriving elemental abundances from emission lines formed in the collisional and photoionized astrophysical plasmas. Accurate representation of target wave functions that properly account for the important correlation and relaxation effects and inclusion of coupling effects including coupling to the continuum are essential components of a reliable collision calculation. Non-orthogonal orbitals technique in multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities. The effect of coupling to the continuum spectrum is included through the use of pseudostates which are chosen to account for most of the dipole polarizabilities of target states. The B-spline basis is used in the R-matrix approach to calculate electron excitation collision strengths and rates. Results for oscillator strengths and electron excitation collision strengths for transitions in N I, O I, O II, O IV, S X and Fe XIV have been produced

Tayal, S. S.↗

NASA NTRS1993

Viewpoint 9--molecular structure of aqueous interfaces

In this review we summarize recent progress in our understanding of the structure of aqueous interfaces emerging from molecular level computer simulations. It is emphasized that the presence of the interface induces specific structural effects which, in turn, influence a wide variety of phenomena occurring near the phase boundaries. At the liquid-vapor interface, the most probable orientations of a water molecule is such that its dipole moment lies parallel to the interface, one O-H bond points toward the vapor and the other O-H bond is directed toward the liquid. The orientational distributions are broad and slightly asymmetric, resulting in an excess dipole moment pointing toward the liquid. These structural preferences persist at interfaces between water and nonpolar liquids, indicating that the interactions between the two liquids in contact are weak. It was found that liquid-liquid interfaces are locally sharp but broadened by capillary waves. One consequence of anisotropic orientations of interfacial water molecules is asymmetric interactions, with respect to the sign of the charge, of ions with the water surface. It was found that even very close to the surface ions retain their hydration shells. New features of aqueous interfaces have been revealed in studies of water-membrane and water-monolayer systems. In particular, water molecules are strongly oriented by the polar head groups of the amphiphilic phase, and they penetrate the hydrophilic head-group region, but not the hydrophobic core. At infinite dilution near interfaces, amphiphilic molecules exhibit behavior different from that in the gas phase or in bulk water. This result sheds new light on the nature of hydrophobic effect in the interfacial regions. The presence of interfaces was also shown to affect both equilibrium and dynamic components of rates of chemical reactions. Applications of continuum models to interfacial problems have been, so far, unsuccessful. This, again, underscores the importance of molecular-level information about interfaces.

Review, Tutorial↗

NASA NTRS2008

Growth Control and Disease Mechanisms in Computational Embryogeny

This paper presents novel approach to applying growth control and diseases mechanisms in computational embryogeny. Our method, which mimics fundamental processes from biology, enables individuals to reach maturity in a controlled process through a stochastic environment. Three different mechanisms were implemented; disease mechanisms, gene suppression, and thermodynamic balancing. This approach was integrated as part of a structural evolutionary model. The model evolved continuum 3-D structures which support an external load. By using these mechanisms we were able to evolve individuals that reached a fixed size limit through the growth process. The growth process was an integral part of the complete development process. The size of the individuals was determined purely by the evolutionary process where different individuals matured to different sizes. Individuals which evolved with these characteristics have been found to be very robust for supporting a wide range of external loads.

Genetic algorithm↗

NASA NTRS1972

Production cross sections for Lee-Wick massive electromagnetic bosons and for spin-zero and spin-one W bosons at high energies.

Production cross sections for three types of hypothetical particles are calculated in the presented paper. Several (Z, Z') cases were studied corresponding to elastic scattering off protons and neutrons (either free or embedded within a Fermi sea), coherent scattering off a nucleus, and inelastic scattering off a proton (in which case Z' denotes a nucleon resonance or hadronic system in the continuum). Detailed structure-function data are used to improve the accuracy of the inelastic scattering calculation. Results of calculations are given for beam energies between 50 and 10,000 GeV, and masses between 5 and 40 GeV for the massive Lee-Wick spin-1 boson. Cross sections were computed for resonant and semiweak processes. The production cross section of spin-zero weak intermediate bosons was found to be at least one order of magnitude smaller than for spin-1 weak bosons in nearly all regions of interest. The production cross section of spin-zero weak intermediate bosons for inelastic scattering off protons compares with that for elastic scattering in the regions of interest. In the case of massive spin-1 bosons and spin-1 weak intermediates, the main contribution to total production cross section off protons is elastic.

Linsker, R.↗

NASA NTRS1971

Ablation and radiation coupled viscous hypersonic shock layers, volume 2

A computer program can be used to determine the radiative flux and flux divergence through a nonisothermal planar slab of gas is described. The program considers species typical of air at high temperatures and nylon or carbon phenolic ablation products. The equations solved are for a variable optical depth line and continuum gas model using species cross sections and line widths as the basic data.

Engel, C. D.↗

NASA NTRS1990

Energy balance in the solar transition region. I - Hydrostatic thermal models with ambipolar diffusion

The energy balance in the lower transition region is analyzed by constructing theoretical models which satisfy the energy balance constraint. The energy balance is achieved by balancing the radiative losses and the energy flowing downward from the corona. This energy flow is mainly in two forms: conductive heat flow and hydrogen ionization energy flow due to ambipolar diffusion. Hydrostatic equilibrium is assumed, and, in a first calculation, local mechanical heating and Joule heating are ignored. In a second model, some mechanical heating compatible with chromospheric energy-balance calculations is introduced. The models are computed for a partial non-LTE approach in which radiation departs strongly from LTE but particles depart from Maxwellian distributions only to first order. The results, which apply to cases where the magnetic field is either absent, or uniform and vertical, are compared with the observed Lyman lines and continuum from the average quiet sun. The approximate agreement suggests that this type of model can roughly explain the observed intensities in a physically meaningful way, assuming only a few free parameters specified as chromospheric boundary conditions.

Fontenla, J. M.↗

NASA NTRS1999

Laminar Soot Processes

Soot formation within hydrocarbon-fueled flames is an important unresolved problem of combustion science for several reasons: soot emissions are responsible for more deaths than any other combustion pollutant, thermal loads due to continuum radiation from soot limit the durability of combustors, thermal radiation from soot is mainly responsible for the growth and spread of unwanted fires, carbon monoxide associated with soot emissions is responsible for most fire deaths, and limited understanding of soot processes is a major impediment to the development of computational combustion. Thus, soot processes within laminar nonpremixed (diffusion) flames are being studied, emphasizing space-based experiments at microgravity. The study is limited to laminar flames due to their experimental and computational tractability, noting the relevance of these results to practical flames through laminar flamelet concepts. The microgravity environment is emphasized because buoyancy affects soot processes in laminar diffusion flames whereas effects of buoyancy are small for most practical flames. Results discussed here were obtained from experiments carried out on two flights of the Space Shuttle Columbia. After a brief discussion of experimental methods, results found thus far are described, including soot concentration measurements, laminar flame shapes, laminar smoke points and flame structure. The present discussion is brief.

Lin, K. -C.↗

NASA NTRS1999

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗

NASA NTRS2005

Physics-Based Spectra of Accretion Disks around Black Holes

The purpose of this grant was to begin the process of deriving the light output of accretion disks around black holes directly from the actual processes that inject heat into the accreting matter, rather than from guessed dependences of heating rate on physical parameters. At JHU, the effort has focussed so far on models of accretion onto "intermediate mass black holes", a possible class of black holes, examples of which may have recently been discovered in nearby galaxies. There, Krolik and his student (Yawei Hui) have computed stellar atmospheres for uniformly-heated disks around this class of black holes. Their models serve two purposes: they are the very first serious attempts to compute the spectrum from accreting black holes in this mass range; and a library of such models can be used later in this program as contrasts for those computed on the basis of real disk dynamics. The output from these local disk calculations has also been successfully coupled to a program that applies the appropriate relativistic transformations and computes photon trajectories in order to predict the spectrum received by observers located at different polar angles. The principal new result of these calculations is the discovery of potentially observable ionization edges of H-like C and O at frequencies near the peak in flux from these objects. Most of the grant money at UCSB was spent on supporting graduate student Shane Davis. In addition. some money was spent on supporting two other students: Ari Socrates (now a Hubble Fellow at Princeton), and Laura Melling. Davis spent the year constructing stellar atmosphere models of accretion disks appropriate for the high/soft (thermal) state of black hole X-ray binaries. As with AGN models published previously by our collaboration with NASA support. our models include a complete general relativistic treatment of both the disk structure and the propagation of photons from the disk to a distant observer. They also include all important continuum opacity sources, including Compton scattering and bound-free opacity from abundant metal species. The principal new result is that bound-free opacity is very significant in altering the continuum spectral shape, resulting for example in quite different "color correction factors" compared to those predicted previously. In addition, the models predict a relationship between luminosity and inner disk temperature that is, for the first time, in accord with that observed. The primary purpose of the grant was to incorporate more realistic accretion disk physics, learned largely from simulations, into such spectral models. The Davis et al. paper includes consideration of a vertical dissipation profile computed from radiation magneto-hydrodynamic simulations of MRI turbulence by N. J. Turner (2004). So long as the disk is effectively thick, such dissipation profiles do not affect the predicted spectrum significantly. (More work needs to be done on these simulations, however.) A potentially more serious issue is that MRI turbulence produces substantial inhomogeneities, as do photon bubble instabilities. These inhomogeneities can affect the spectra by enhancing the effects of absorption opacity over scattering opacity. We have done some preliminary Monte Carlo calculations to explore these effects.

Krolik, Julian H.↗

NASA NTRS1970

Code CUGEL: A code to unfold Ge(Li) spectrometer polyenergetic gamma photon experimental distributions

A FORTRAN code was developed for the Univac 1108 digital computer to unfold lithium-drifted germanium semiconductor spectrometers, polyenergetic gamma photon experimental distributions. It was designed to analyze the combination continuous and monoenergetic gamma radiation field of radioisotope volumetric sources. The code generates the detector system response matrix function and applies it to monoenergetic spectral components discretely and to the continuum iteratively. It corrects for system drift, source decay, background, and detection efficiency. Results are presented in digital form for differential and integrated photon number and energy distributions, and for exposure dose.

Steyn, J. J.↗

NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗

NASA NTRS2012

Atomic Data for the CHIANTI Database

The CHIANTI spectral code consists of an atomic database and a suite of computer programs to calculate the optically thin spectrum of astrophysical objects and to carry out spectroscopic plasma diagnostics. The database includes atomic energy levels, wavelengths, radiative transition rates, collisional excitation, ionization and recombination rate coefficients, as well as data to calculate free-free, free-bound and two-photon continuum emission. In recent years, we have been pursuing a program to calculate atomic data for ions whose lines have been observed in astrophysical spectra but have been neglected in the literature, and to provide CHIANTI with all the data necessary to predict line intensities. There are two types of such ions: those for which calculations are available for low-energy configurations but not for high-energy configurations (i.e., C-like, N-like, O-like systems), and ions that have never or only seldom been studied. This poster will summarize the current status of this project and indicate the future activities .

Bhatia, Anand K.↗

NASA NTRS1988

Shape optimization of 3D continuum structures via force approximation techniques

The existing need to develop methods whereby the shape design efficiency can be improved through the use of high quality approximation methods is addressed. An efficient approximation method for stress constraints in 3D shape design problems is proposed based on expanding the nodal forces in Taylor series with respect to shape variations. The significance of this new method is shown through elementary beam theory calculations and via numerical computations using 3D solid finite elements. Numerical examples including the classical cantilever beam structure and realistic automotive parts like the engine connecting rod are designed for optimum shape using the proposed method. The numerical results obtained from these methods are compared with other published results, to assess the efficiency and the convergence rate of the proposed method.

Vanderplaats, Garret N.↗

NASA NTRS1986

Thermodynamics and combustion modeling

Modeling fluid phase phenomena blends the conservation equations of continuum mechanics with the property equations of thermodynamics. The thermodynamic contribution becomes especially important when the phenomena involve chemical reactions as they do in combustion systems. The successful study of combustion processes requires (1) the availability of accurate thermodynamic properties for both the reactants and the products of reaction and (2) the computational capabilities to use the properties. A discussion is given of some aspects of the problem of estimating accurate thermodynamic properties both for reactants and products of reaction. Also, some examples of the use of thermodynamic properties for modeling chemically reacting systems are presented. These examples include one-dimensional flow systems and the internal combustion engine.

Zeleznik, Frank J.↗

NASA NTRS1992

A far wing line shape theory and its application to the foreign-broadened water continuum absorption. III

The far wing line shape theory developed previously and applied to the calculation of the continuum absorption of pure water vapor is extended to foreign-broadened continua. Explicit results are presented for H2O-N2 and H2O-CO2 in the frequency range from 0 to 10,000/cm. For H2O-N2 the positive and negative resonant frequency average line shape functions and absorption coefficients are computed for a number of temperatures between 296 and 430 K for comparison with available laboratory data. In general the agreement is very good.

Ma, Q.↗

NASA NTRS1993

Efficient Three-Dimensional Direct Simulation Monte Carlo for Complex Geometry Problems

The simulation of flowfields in the transition flow regime is notoriously difficult with high demands on computer resources (CPU time and storage) and user expertise/labor. This paper describes a new, efficient code which has been developed to simulate high Knudsen number flowfields in three dimensions about bodies of arbitrarily complex geometry. The algorithm has been tested over a wide range of conditions, from free molecular to near-continuum flow regimes, for slender and blunt bodies, for re-entry vehicles and spacecraft. A series of validation tests have been conducted using both wind-tunnel measurements and flight data.

Rault, Didier F. G.↗