Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scalar fields”

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.

447 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

The stability analysis of magnetohydrodynamic equilibria - Comparing the thermodynamic approach with the energy principle

This paper is a contribution to the stability analysis of current-carrying plasmas, i.e., plasma systems that are forced by external mchanisms to carry a nonrelaxing electrical current. Under restriction to translationally invariant configurations, the thermodynamic stability criterion for a multicomponent plasma is rederived within the framework of nonideal MHD. The chosen dynamics neglects scalar resistivity, but allows for other types of dissipation effects both in Ohm's law and in the equation of motion. In the second section of the paper, the thermodynamic stability criterion is compared with the ideal MHD based energy principle of Bernstein et al. With the help of Schwarz's inequality, it is shown that the former criterion is always more 'pessimistic' than the latter, i.e., that thermodynamic stability implies stability according to the MHD principle, but not vice versa. This reuslt confirms the physical plausible idea that dissipational effects tend to weaken the stability properties of current-carrying plasma equilibria by breaking the constraints of ideal MHD and allowing for possibly destabilizing effects such as magnetic field line reconfiguration.

Brinkmann, R. P.↗
NASA NTRS2010

Micromechanics Modeling of Composites Subjected to Multiaxial Progressive Damage in the Constituents

The high-fidelity generalized method of cells composite micromechanics model is extended to include constituent-scale progressive damage via a proposed damage model. The damage model assumes that all material nonlinearity is due to damage in the form of reduced stiffness, and it uses six scalar damage variables (three for tension and three for compression) to track the damage. Damage strains are introduced that account for interaction among the strain components and that also allow the development of the damage evolution equations based on the constituent material uniaxial stress strain response. Local final-failure criteria are also proposed based on mode-specific strain energy release rates and total dissipated strain energy. The coupled micromechanics-damage model described herein is applied to a unidirectional E-glass/epoxy composite and a proprietary polymer matrix composite. Results illustrate the capability of the coupled model to capture the vastly different character of the monolithic (neat) resin matrix and the composite in response to far-field tension, compression, and shear loading.

Bednarcyk, Brett A.↗
NASA NTRS1977

Interactive calculation procedure for supersonic flows

An interactive procedure was developed for supersonic viscous flows that can be used for either two-dimensional or axisymmetric configurations. The procedure is directed to supersonic internal flows as well as those supersonic external flows that require consideration of mutual interaction between the outer flow and the boundary layer flow. The flow field is divided into two regions: an inner region which is highly viscous and mostly subsonic and an outer region where the flow is supersonic and in which viscous effects are small but not negligible. For the outer region a numerical solution is obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines. The inner region is treated by a system of equations of the boundary layer type that includes higher order effects such as longitudinal and transverse curvature and normal pressure gradients. These equations are coupled and solved simultaneously in the physical coordinates by using an implicit finite difference scheme. This system can also be used to calculate laminar and turbulent boundary layers using a scalar eddy viscosity concept.

Tassa, Y.↗
NASA NTRS1995

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗
NASA NTRS2020

Small Aperture Telescopes for the Simons Observatory

The Simons Observatory (SO) is an upcoming cosmic microwave background(CMB) experiment located on Cerro Toco, Chile, that will map the microwave sky in temperature and polarization in six frequency bands spanning 27 to 285 GHz. SO will consist of one 6-meter Large Aperture Telescope (LAT) fielding∼30,000 detectors and an array of three 0.42-meter Small Aperture Telescopes (SATs) fielding an additional 30,000 detectors. This synergy will allow for the extremely sensitive characterization of the CMB over an-gular scales ranging from an arcmin to tens of degrees, enabling a wide range of scientific output. Here we focus on the SATs targeting degree angular scales with successive dichroic instruments observing at Mid-Frequency (MF: 93/145 GHz), Ultra-High-Frequency (UHF:225/285 GHz), and Low-Frequency (LF: 27/39 GHz). The three SATs will be able to map∼10% of the sky to a noise level of∼2 μK-arcmin when combining 93 and 145 GHz. The multiple frequency bands will allow the CMB to be separated from galactic foregrounds (primarily synchrotron and dust), with the primary science goal of characterizing the primordial tensor-to-scalar ratio, r, at a target level ofσ(r)≈0.003.

Small aperature telescope↗
NASA NTRS2003

A Scintillation Counter System Design To Detect Antiproton Annihilation using the High Performance Antiproton Trap(HiPAT)

The High Performance Antiproton Trap (HiPAT), a system designed to hold up to l0(exp 12) charge particles with a storage half-life of approximately 18 days, is a tool to support basic antimatter research. NASA's interest stems from the energy density represented by the annihilation of matter with antimatter, 10(exp 2)MJ/g. The HiPAT is configured with a Penning-Malmberg style electromagnetic confinement region with field strengths up to 4 Tesla, and 20kV. To date a series of normal matter experiments, using positive and negative ions, have been performed evaluating the designs performance prior to operations with antiprotons. The primary methods of detecting and monitoring stored normal matter ions and antiprotons within the trap includes a destructive extraction technique that makes use of a micro channel plate (MCP) device and a non-destractive radio frequency scheme tuned to key particle frequencies. However, an independent means of detecting stored antiprotons is possible by making use of the actual annihilation products as a unique indicator. The immediate yield of the annihilation event includes photons and pie mesons, emanating spherically from the point of annihilation. To "count" these events, a hardware system of scintillators, discriminators, coincident meters and multi channel scalars (MCS) have been configured to surround much of the HiPAT. Signal coincidence with voting logic is an essential part of this system, necessary to weed out the single cosmic ray events from the multi-particle annihilation shower. This system can be operated in a variety of modes accommodating various conditions. The first is a low-speed sampling interval that monitors the background loss or "evaporation" rate of antiprotons held in the trap during long storage periods; provides an independent method of validating particle lifetimes. The second is a high-speed sample rate accumulating information on a microseconds time-scale; useful when trapped antiparticles are extracted against a target, providing an indication of quantity. This paper details the layout of this system, setup of the hardware components around HiPAT, and applicable checkouts using normal matter radioactive sources.

Martin, James J.↗
NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu↗
NASA NTRS1989

Marine stratocumulus structure

Thirty-three Landsat TM scenes of California stratocumulus cloud fields were acquired as part of the FIRE Marine Stratocumulus Intensive Field Observations in July 1987. They exhibit a wide variety of stratocumulus structures. Analysis has so far focused upon the July 7 scene, in which aircraft from NASA, NCAR, and the British Meteorological Office repeatedly gathered data across a stratocumulus-fair weather cumulus transition. The aircraft soundings validate the cloud base temperature threshold determined by spatial coherence analysis of the TM thermal band. Brightness variations in the stratocumulus region exhibit a -5/3 power-law decrease of the wavenumber spectra for scales larger than the cloud thickness, about 200 m, changing to a -3 power at smaller scales. Observations by an upward-looking three-channel microwave radiometer on San Nicolas Island also show the -5/3 power-law in total integrated liquid water, suggesting that the largest-scale TM brightness variations are primarily due to variations in the liquid water. The Kolmogorov 5/3 power suggests that for some purposes liquid water in turbulent stratocumulus clouds may be treated as a passive scalar, simply reflecting variations in vertical velocity. This may be tested using the velocities measured by the aircraft.

Cahalan, Robert F.↗
NASA NTRSDate not supplied

Distribution of SLS Integrated Load Uncertainty to Surface Pressures and Sectional Loads

Aerodynamic loads that are important to launch vehicle programs such as NASA’s Space Launch System (SLS) include both integrated loads such as the force & moment on the entire vehicle and distributed loads. In this work two examples of distributed loads are considered: the pressure field on the surface of the vehicle and sectional loads, which are one-dimensional distributions along the axis of the launch vehicle. In some modern flight programs, the integrated loads, such as lift and drag, used to design the guidance and control laws for the vehicle come from wind tunnel testing, while distributed loads are produced using Computational Fluid Dynamics(CFD). The first task that is addressed in this paper, then, is to provide a formal method to adjust the distributed loads so that integrating them matches the prescribed integrated load. In addition, the integrated loads in a launch vehicle typically include an uncertainty estimate. The second task is to distribute this prescribed integrated uncertainty to each point in a distributed load. Both tasks are addressed using the same technique, which is to create distributed load profiles that isolate adjustments to one integrated load while leaving the others unaffected. These adjustments are informed by Proper Orthogonal Decomposition (POD) of the entire CFD-based distributed load database. Once applied, the adjusted distributed loads can be used to evaluate any scalar quantity of interest that might be needed by downstream users such as structural analysis or trajectory modelers.

SLS↗
NASA NTRS1996

Evaluation of joint probability density function models for turbulent nonpremixed combustion with complex chemistry

Two types of mixing sub-models are evaluated in connection with a joint-scalar probability density function method for turbulent nonpremixed combustion. Model calculations are made and compared to simulation results for homogeneously distributed methane-air reaction zones mixing and reacting in decaying turbulence within a two-dimensional enclosed domain. The comparison is arranged to ensure that both the simulation and model calculations a) make use of exactly the same chemical mechanism, b) do not involve non-unity Lewis number transport of species, and c) are free from radiation loss. The modified Curl mixing sub-model was found to provide superior predictive accuracy over the simple relaxation-to-mean submodel in the case studied. Accuracy to within 10-20% was found for global means of major species and temperature; however, nitric oxide prediction accuracy was lower and highly dependent on the choice of mixing sub-model. Both mixing submodels were found to produce non-physical mixing behavior for mixture fractions removed from the immediate reaction zone. A suggestion for a further modified Curl mixing sub-model is made in connection with earlier work done in the field.

Smith, N. S. A.↗
NASA NTRS1995

A theoretical derivation of the dependence of the remotely sensed reflectance of the ocean on the inherent optical properties

An expression for the ratio of the upwelling nadir radiance L(pi, z) and the downwelling scalar irradiance E(sub od)(Z) is derived from the following equation of radiative transfer. This expression is given by RSR(z) = (L(pi, z))/E(sub od) = (f(sub b)(z)b(sub b)(z))/2 pi(k(pi, z) + c(z) - F(sub L)(z)b(sub f)(z)), where b(sub b)(z) is the backscattering coefficient, k(pi, z) is the vertical attenuation coefficient of the nadir radiance, c(z) is the beam attenuation coefficient, and f(sub b)(z) and f(sub L)(z) are shape parameters that depend on the shape of the volume scattering function and the radiance distribution. Successive approximations are subsequently applied to the above exact equation. These are f(sub b)(z) = (2 pi beta(pi - theta(sub m), z)/(b(sub b)(z))), where beta(pi - theta(sub m), z) is the volume scattering function at 180 deg minus the zenith angle of the maximum radiance, and k(pi, z) = am = c(1 - 0.52 b/c - 0.44 (b/c)(exp 2)), where m is a parameter that is numerically equal to the inverse of the average cosine of the asymptotic light field for a medium with the same inherent optical properties, a is the absorption coefficient, and b/c is the single scattering albedo. Together with f(sub L)(z) = 1.05 and application of Gershun's equation, it is shown that for nearly all oceanic cases RSR(z) identical to L(pi, z)/E(sub od)(z) = (Beta(pi - theta(sub m), z))/(a(z)(1 + m(z))).

Zaneveld, J. Ronald V.↗
NASA NTRS1998

Constitutive Theory Developed for Monolithic Ceramic Materials

With the increasing use of advanced ceramic materials in high-temperature structural applications such as advanced heat engine components, the need arises to accurately predict thermomechanical behavior that is inherently time-dependent and that is hereditary in the sense that the current behavior depends not only on current conditions but also on the material's thermomechanical history. Most current analytical life prediction methods for both subcritical crack growth and creep models use elastic stress fields to predict the time-dependent reliability response of components subjected to elevated service temperatures. Inelastic response at high temperatures has been well documented in the materials science literature for these material systems, but this issue has been ignored by the engineering design community. From a design engineer's perspective, it is imperative to emphasize that accurate predictions of time-dependent reliability demand accurate stress field information. Ceramic materials exhibit different time-dependent behavior in tension and compression. Thus, inelastic deformation models for ceramics must be constructed in a fashion that admits both sensitivity to hydrostatic stress and differing behavior in tension and compression. A number of constitutive theories for materials that exhibit sensitivity to the hydrostatic component of stress have been proposed that characterize deformation using time-independent classical plasticity as a foundation. However, none of these theories allow different behavior in tension and compression. In addition, these theories are somewhat lacking in that they are unable to capture the creep, relaxation, and rate-sensitive phenomena exhibited by ceramic materials at high temperatures. The objective of this effort at the NASA Lewis Research Center has been to formulate a macroscopic continuum theory that captures these time-dependent phenomena. Specifically, the effort has focused on inelastic deformation behavior associated with these service conditions by developing a multiaxial viscoplastic constitutive model that accounts for time-dependent hereditary material deformation (such as creep and stress relaxation) in monolithic structural ceramics. Using continuum principles of engineering mechanics, we derived the complete viscoplastic theory from a scalar dissipative potential function.

Janosik, Lesley A.↗
NASA NTRS1997

A Viscoplastic Constitutive Theory for Monolithic Ceramic Materials

With increasing use of ceramic materials in high temperature structural applications such as advanced heat engine components, the need arises to accurately predict thermomechanical behavior. This paper, which is the first of two in a series, will focus on inelastic deformation behavior associated with these service conditions by providing an overview of a viscoplastic constitutive model that accounts for time-dependent hereditary material deformation (e.g., creep, stress relaxation, etc.) in monolithic structural ceramics. Early work in the field of metal plasticity indicated that inelastic deformations are essentially unaffected by hydrostatic stress. This is not the case, however, for ceramic-based material systems, unless the ceramic is fully dense. The theory presented here allows for fully dense material behavior as a limiting case. In addition, ceramic materials exhibit different time-dependent behavior in tension and compression. Thus, inelastic deformation models for ceramics must be constructed in a fashion that admits both sensitivity to hydrostatic stress and differing behavior in tension and compression. A number of constitutive theories for materials that exhibit sensitivity to the hydrostatic component of stress have been proposed that characterize deformation using time-independent classical plasticity as a foundation. However, none of these theories allow different behavior in tension and compression. In addition, these theories are somewhat lacking in that they are unable to capture creep, relaxation, and rate-sensitive phenomena exhibited by ceramic materials at high temperature. When subjected to elevated service temperatures, ceramic materials exhibit complex thermomechanical behavior that is inherently time-dependent, and hereditary in the sense that current behavior depends not only on current conditions, but also on thermo-mechanical history. The objective of this work is to present the formulation of a macroscopic continuum theory that captures these time-dependent phenomena. Specifically, the overview contained in this paper focuses on the multiaxial derivation of the constitutive model, and examines the scalar threshold function and its attending geometrical implications.

Janosik, Lesley A.↗
NASA NTRS1993

LIF measurements of scalar mixing in turbulent shear layers

The structure of shear layer flows at high Reynolds numbers remains a very interesting problem. Straight mixing layers have been studied and yielded information on the probability density function (pdf) of a passive scalar across the layer. Konrad and Koochesfahani & Dimotakis measured the pdf of the mixture fraction for mixing layers of moderate Reynolds numbers, each about 25,000 (Re based on velocity difference and visual thickness). Their measurements showed a 'non-marching' pdf (central hump which is invariant from edge to edge across the layer), a result which is linked to the visualizations of the spanwise Kelvin-Helmholtz (K-H) instability mode, which is the primary instability for plane shear layer flows. A secondary instability mode, the Taylor-Gortler (T-G) instability, which is associated with streamwise vortical structures, has also been observed in shear layers. Image reconstruction by Jimenez et al. and volume renderings by Karasso & Mungal at low Re numbers have demonstrated that the K-H and the T-G instability modes occur simultaneously in a non-mutually destructive way, evidence that supports the quasi two-dimensional aspect of these flows and the non-marching character of the pdf at low Reynolds numbers. At higher Re numbers though, the interaction of these two instability modes is still unclear and may affect the mixing process. In this study, we perform measurements of the concentration pdf of plane mixing layers for different operating conditions. At a speed ratio of r = U(sub 1)/U(sub 2) = 4:1, we examine three Reynolds number cases: Re = 14,000, Re = 31,000, and Re = 62,000. Some other Re number cases' results, not presented in detail, are invoked to explain the behavior of the pdf of the concentration field. A case of r = 2.6:1 at Re = 20,000 is also considered. The planar laser-induced fluorescence technique is used to yield quantitative measurements. The different Re are obtained by changing the velocity magnitudes of the two streams. The question of resolution of these measurements is addressed. In order to investigate the effects of the initial conditions on the development and the structure of the mixing layer, the boundary layer on the high-speed side of the splitter plate is tripped. The average concentration and the average mixed fluid concentration are also calculated to further understand the changes in the shear layer for the different cases examined.

Karasso, Paris S.↗
NASA NTRS2011

A New Inversion Routine to Produce Vertical Electron-Density Profiles from Ionospheric Topside-Sounder Data

Two software applications have been produced specifically for the analysis of some million digital topside ionograms produced by a recent analog-to-digital conversion effort of selected analog telemetry tapes from the Alouette-2, ISIS-1 and ISIS-2 satellites. One, TOPIST (TOPside Ionogram Scalar with True-height algorithm) from the University of Massachusetts Lowell, is designed for the automatic identification of the topside-ionogram ionospheric-reflection traces and their inversion into vertical electron-density profiles Ne(h). TOPIST also has the capability of manual intervention. The other application, from the Goddard Space Flight Center based on the FORTRAN code of John E. Jackson from the 1960s, is designed as an IDL-based interactive program for the scaling of selected digital topside-sounder ionograms. The Jackson code has also been modified, with some effort, so as to run on modern computers. This modification was motivated by the need to scale selected ionograms from the millions of Alouette/ISIS topside-sounder ionograms that only exist on 35-mm film. During this modification, it became evident that it would be more efficient to design a new code, based on the capabilities of present-day computers, than to continue to modify the old code. Such a new code has been produced and here we will describe its capabilities and compare Ne(h) profiles produced from it with those produced by the Jackson code. The concept of the new code is to assume an initial Ne(h) and derive a final Ne(h) through an iteration process that makes the resulting apparent-height profile fir the scaled values within a certain error range. The new code can be used on the X-, O-, and Z-mode traces. It does not assume any predefined profile shape between two contiguous points, like the exponential rule used in Jackson s program. Instead, Monotone Piecewise Cubic Interpolation is applied in the global profile to keep the monotone nature of the profile, which also ensures better smoothness in the final profile than in Jackson s program. The new code uses the complete refractive index expression for a cold collisionless plasma and can accommodate the IGRF, T96, and other geomagnetic field models.

Wang, Yongli↗