SEARCH · Engineering Papers

Results for “finite density”

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 NTRS1974

The effects of spline interpolation on power spectral density

This paper discusses the power spectral effects of spline interpolators. A general technique is given for finding the steady-state spectral effects of splines of all orders, when applied following uniform sampling of the input function. The following observations are made: (1) the even order splines that were examined (second and fourth order) possessed divergent steady-state frequency transfer functions, (2) the degree of preservation of the power spectral density of the input process increased with the order of the (odd order) spline used for interpolation, and (3) the reconstruction of a stationary random process over a finite record length will, on the average, have less power than indicated by the steady-state transfer function.

Horowitz, L. L.↗

NASA NTRS1979

Influence of interface on composite failure

The influence of the interface on the composite system behavior is investigated by analytical modeling. The stress analysis is based on the two-dimensional finite element procedure in which twelve-node isoparametric elements with cubic shape functions are used. The location of possible failure sites is predicted by the strain energy density (SED) failure theory which assumes failure to coincide with locations of minimum SED while the locations of maximum SED correspond to regions of excessive distortion or yielding. The results of the analysis show that the way in which the modulus of elasticity varies within the interface is as important in modeling as the average interface modulus.

Sih, G. C.↗

NASA NTRS1983

MHD thermal instabilities in cool inhomogeneous atmospheres

The formation of a coronal state in a stellar atmosphere is investigated. A numerical code is used to study the effects of atmospheric gradients and finite loop dimension on the scale of unstable perturbations, solving for oscillatory perturbations as eigenfunctions of a boundary value problem. The atmosphere is considered as initially isothermal, with density and pressure having scale heights fixed by the hydrostatic equations. Joule mode instability is found to be an efficient mechanism for current filamentation and subsequent heating in initially cool atmospheres. This instability is mainly effective at the top of magnetic loops and is not suppressed by thermal conduction.

Bodo, G.↗

NASA NTRS1988

Probabilistic structural analysis computer code (NESSUS)

Probabilistic structural analysis has been developed to analyze the effects of fluctuating loads, variable material properties, and uncertain analytical models especially for high performance structures such as SSME turbopump blades. The computer code NESSUS (Numerical Evaluation of Stochastic Structure Under Stress) was developed to serve as a primary computation tool for the characterization of the probabilistic structural response due to the stochastic environments by statistical description. The code consists of three major modules NESSUS/PRE, NESSUS/FEM, and NESSUS/FPI. NESSUS/PRE is a preprocessor which decomposes the spatially correlated random variables into a set of uncorrelated random variables using a modal analysis method. NESSUS/FEM is a finite element module which provides structural sensitivities to all the random variables considered. NESSUS/FPI is Fast Probability Integration method by which a cumulative distribution function or a probability density function is calculated.

Shiao, Michael C.↗

NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

NASA NTRS2005

Second-order multiple-scattering theory associated with backscattering enhancement for a millimeter wavelength weather radar with a finite beam width

Effects of multiple scattering on reflectivity are studied for millimeter wavelength weather radars. A time-independent vector theory, including up to second-order scattering, is derived for a single layer of hydrometeors of a uniform density and a uniform diameter. In this theory, spherical waves with a Gaussian antenna pattern are used to calculate ladder and cross terms in the analytical scattering theory. The former terms represent the conventional multiple scattering, while the latter terms cause backscattering enhancement in both the copolarized and cross-polarized components. As the optical thickness of the hydrometeor layer increases, the differences from the conventional plane wave theory become more significant, and essentially, the reflectivity of multiple scattering depends on the ratio of mean free path to radar footprint radius. These results must be taken into account when analyzing radar reflectivity for use in remote sensing.

radar reflectivity↗

NASA NTRS2007

Role of Ionospheric Plasmas in Earth's Magnetotail

This tutorial will summarize observations and theories indicating a prominent role of ionospheric plasma in the Earth's magnetotail. At the Global scale, I will argue that it is ionospheric plasma momentum and dynamic pressure that are responsible for the production of plasmoids, through the action of a transient near-Earth neutral or X-line, which serves to release excessive plasma pressure from the magnetotail field. Ionospheric plasma gains the momentum and energy to produce plasmoids and their related effects through its interaction with the solar wind, beginning at the dayside reconnection region and extending across the polar caps through the magnetotail lobes. This distant neutral line can be depicted as a feature much like that found in cometary magnetospheres, where disconnection limits the amount of IMF hung up on the cometary coma. On the other hand, the near-Earth neutral one can be seen as a feature unique to planets with an intrinsic magnetic field and internal source of plasma, the heating of which produces pressures too large to be restrained. Ionospheric plasmas also have other more local roles to play in the magnetotail. The circulation influences the composition of the plasma sheet, and the resultant wave environment, giving rise to reduced wave propagation speeds. Important heavy ion cyclotron resonances, and enhanced finite gyro-radius effects including non-adiabatic particle acceleration. At minimum, the presence of ionospheric plasma must influence the rate of reconnection via its enhanced mass density. Other non-MHD effects of ionospheric plasma presence are likely to be important but need much more investigation to be well understood. The MMS mission is designed to penetrate the subtle diffusion region physics that is involved, and its ability to observe ionospheric plasma involvement in reconnection will contribute significantly toward that goal.

Moore, Thomas E.↗

NASA NTRS2016

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗

NASA NTRS1982

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

This paper describes a computer code for calculating the flow dynamics of constant density flow in the second stage trumpet shaped nozzle section of a two stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite difference computer code has been developed to allow the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for efficient solution of the equations. Turbulence simulation is by way of a two-equation Kappa-epsilon model. The code produces as output the flowfield map of the nondimensional stream function, axial, and swirl velocity. Good agreement was obtained between the theoretical predictions and the qualitative experimental results. The best seed injector location for uniform seed distribution at combustor exit is with injector located centrally on the combustor axis at entrance to the second stage combustor.

Gupta, A. K.↗

NASA NTRS1987

Comparison of measured temperatures, thermal stresses and creep residues with predictions on a built-up titanium structure

Temperature, thermal stresses, and residual creep stresses were studied by comparing laboratory values measured on a built-up titanium structure with values calculated from finite-element models. Several such models were used to examine the relationship between computational thermal stresses and thermal stresses measured on a built-up structure. Element suitability, element density, and computational temperature discrepancies were studied to determine their impact on measured and calculated thermal stress. The optimum number of elements is established from a balance between element density and suitable safety margins, such that the answer is acceptably safe yet is economical from a computational viewpoint. It is noted that situations exist where relatively small excursions of calculated temperatures from measured values result in far more than proportional increases in thermal stress values. Measured residual stresses due to creep significantly exceeded the values computed by the piecewise linear elastic strain analogy approach. The most important element in the computation is the correct definition of the creep law. Computational methodology advances in predicting residual stresses due to creep require significantly more viscoelastic material characterization.

Jenkins, Jerald M.↗

NASA NTRS1992

An iodine hypersonic wind tunnel for the study of nonequilibrium reacting flows

A pilot scale hypersonic wind tunnel operating on pure iodine vapor has been designed and tested. The wind tunnel operates intermittently with a run phase lasting approximately 20 minutes. Successful recirculation of the iodine used during the run phase has been achieved but can be improved. Relevant issues regarding the full scale facility's design and operation, and the use of iodine as a working gas are discussed. Continuous wave laser induced fluorescence was used to monitor number densities within the plume flowfield, while pulsed laser induced fluorescence was used in an initial attempt to measure vibrational energy state population distributions. Preliminary nozzle flow calculations based on finite rate chemistry are presented.

Pham-Van-diep, G. C.↗

NASA NTRS1995

Computations of Axisymmetric Flows in Hypersonic Shock Tubes

A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar flow, such as test times, and to predict density and velocity profiles. Boundary-layer parameters such as bar-delta(sub u), bar-delta(sup *), and bar-tau(sub w), and test time parameters such as bar-tau and particle time of flight t(sub f), are computed and compared with those evaluated by using Mirels' correlations. This article then discusses in detail the effects of flow nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis of flowfield parameters, using a computer code such as used in this study, must be performed.

Sharma, Surendra P.↗

NASA NTRS2003

Thermal Response of Cooled Silicon Nitride Plate Due to Thermal Conductivity Effects Analyzed

Lightweight, strong, tough high-temperature materials are required to complement efficiency improvements for next-generation gas turbine engines that can operate with minimum cooling. Because of their low density, high-temperature strength, and high thermal conductivity, ceramics are being investigated as materials to replace the nickelbase superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass (ref. 1). To complement the effectiveness of the ceramics and their applicability for turbine engine applications, a parametric study using the finite element method is being carried out. The NASA Glenn Research Center remains very active in conducting and supporting a variety of research activities related to ceramic matrix composites through both experimental and analytical efforts (ref. 1). The objectives of this work are to develop manufacturing technology, develop a thermal and environmental barrier coating (TBC/EBC), develop an analytical modeling capability to predict thermomechanical stresses, and perform a minimal burner rig test on silicon nitride (Si3N4) and SiC/SiC turbine nozzle vanes under simulated engine conditions. Moreover, we intend to generate a detailed database of the material s property characteristics and their effects on structural response. We expect to offer a wide range of data since the modeling will account for other variables, such as cooling channel geometry and spacing. Comprehensive analyses have begun on a plate specimen with Si3N4 cooling holes.

Baaklini, George Y.↗

NASA NTRS1996

Probabilistic Fatigue Life Analysis of High Density Electronics Packaging

The fatigue of thin film metal interconnections in high density electronics packaging subjected to thermal cycling has been evaluated using a probabilistic fracture mechanics methodology. This probabilistic methodology includes characterization of thin film stress using an experimentally calibrated finite element model and simulation of flaw growth in the thin films using a stochastic crack growth model.

thermal cycling fatigue of thin film metal interco↗

NASA NTRS2010

Broken Ergodicity in MHD Turbulence

Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.

Shebalin, John V.↗

NASA NTRS1979

Generalized lower-hybrid-drift instability

The theory of lower-hybrid-drift instability is extended to include a finite value of the component of wave vector parallel to the ambient magnetic field so that the analysis bridges the usual lower-hybrid-drift instability of flute modes and the modified-two-stream instability. The present theory also includes electromagnetic and ambient magnetic field-gradient effects. It is found that in the cold-electron limit the density and magnetic gradients can qualitatively modify the conclusion obtained in the early theory of the modified-two-stream instability. For example, even if the relative drift far exceeds the Alfven speed of the plasma, the instability may still persist. This result is in contrast to that established in the literature. When the electron temperature is finite, the problem is complicated. Numerical solutions are obtained for a number of cases.

Hsia, J. B.↗

NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

NASA NTRS1975

High-beta turbulence in two-dimensional magnetohydrodynamics

Incompressible turbulent flows were investigated in the framework of ideal magnetohydrodynamics. Equilibrium canonical distributions are determined in a phase whose coordinates are the real and imaginary parts of the Fourier coefficients for the field variables. The magnetic field and fluid velocity have variable x and y components, and all field quantities are independent of z. Three constants of the motion are found which survive the truncation in Fourier space and permit the construction of canonical distributions with three independent temperatures. Spectral densities are calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long wavelength, self-generated, magnetic fields ("magnetic islands"). In the presence of finite dissipation, energy cascades to higher wave numbers can be accompanied by vector potential cascades to lower wave numbers, in much the same way that in the fluid dynamic case, energy cascades to lower wave numbers accompany entropy cascades to higher wave numbers.

Fyfe, D.↗