Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
NASA NTRSDate not supplied

PM 2.5 Concentrations over Major Metropolitan Regions Inferred from Airborne High Spectral Resolution Lidar Measurements Using Machine Learning Regression

We use measurements of near-surface aerosol backscatter, extinction, and depolarization acquired by four NASA Langley Research Center airborne High Spectral Resolution Lidars (HSRLs) to develop a machine learning regression methodology to infer PM2.5 concentrations at the surface and aloft. These airborne HSRL measurements were acquired over major metropolitan regions in the United States and Asia during more than 170 flights since 2010. Hourly surface PM2.5 measurements from the EPA air quality system and similar networks in other countries acquired within 10 km and 15 minutes of these near-surface HSRL measurements are used to train models that compute PM2.5 concentrations from the HSRL measurements. We examine several regression methods and find that exponential Gaussian Process algorithms consistently give the best performance in terms of the lowest root-mean-square (RMS) errors and the highest correlations. Model performance varies significantly depending on various combinations of HSRL aerosol measurements (e.g., aerosol backscatter, extinction, depolarization, backscatter color ratios, lidar ratios, aerosol optical thickness) and retrievals (e.g., mixed layer height, aerosol type) used in the regressions. Models that use near-surface measurements of aerosol backscatter and aerosol intensive properties such as depolarization, backscatter color ratio, and lidar ratio typically give the best performance with RMS errors around 4 mg/m3 and correlation coefficients above 0.9. HSRL measurements were often acquired when the aircraft flew systematic “raster-scan” patterns for several hours over these cities. These flight patterns enabled measurements of the spatial, temporal, and vertical variabilities in the distributions of aerosol backscatter and aerosol intensive properties and allowed us to derive the corresponding variabilities in PM2.5 concentrations. We present examples of such variabilities over urban areas in the United States as well as Asia. We describe also how the distribution of surface PM2.5 varies with aerosol type and use these retrievals to examine model simulations of surface PM2.5 in these metropolitan regions. We also discuss how this methodology may be applied to measurements from satellite lidars such as CALIOP on CALIPSO and ATLID on EarthCARE.

lidar↗
NASA NTRS1984

Ground-state-energy theorem and the virial theorem of a many-particle system in d dimensions

The equivalence of Pauli's ground-state-energy theorem and the virial theorem is demonstrated for a many-particle system interacting with an interparticle potential in d dimensions at zero and finite temperatures. Pauli's theorem has an integral form in which the variable is the coupling constant e-squared, while the virial theorem has a differential form in which the variable has the number density n. The essence of the equivalence proof consists in changing the variable from n to e-squared by noting the dependence of the excess free energy on dimensionless quantities for zero-temperature and classical cases.

Iwamoto, N.↗
NASA NTRS1988

Mechanics of fatigue crack closure

Papers are presented on plasticity induced crack closure, crack closure in fatigue crack growth, the dependence of crack closure on fatigue loading variables, and a procedure for standardizing crack closure levels. Also considered are a statistical approach to crack closure determination, the crack closure behavior of surface cracks under pure bending, closure measurements on short fatigue cracks, and crack closure under plane strain conditions. Other topics include fatigue crack closure behavior at high stress ratios, the use of acoustic waves for the characterization of closed fatigue cracks, and the influence of fatigue crack wake length and state of stress on crack closure.

Newman, J. C., Jr.↗
NASA NTRS1993

Variable temperature pressure broadening of the 4(1,4)-3(2,1) transition of H2O by O2 and N2

The O2 and N2 pressure-broadening parameters of the 4(1,4)-3(2,1) rotational transition in the ground vibrational state of H2O have been measured in the temperature range between 100 and 520 K. Above 250 K the measurements were made in an equilibrium cell. Below 250 K a cell, which uses collisional cooling to circumvent the temperature limits imposed by the vapor pressure of the sample gas, was used. The data were fitted to an exponential temperature-dependence for data above 150 K with resultant n values of 0.81(3) for O2 and 0.70(3) for N2, where the entries in parentheses represent the error in the last significant figure of the parameter. Below 150 K the measured pressure-broadening parameters are smaller than those calculated using these values of n.

Goyette, Thomas M.↗
NASA NTRS1989

Internal state variable approach for predicting stiffness reductions in fibrous laminated composites with matrix cracks

A mathematical model utilizing the internal state variable concept is proposed for predicting the upper bound of the reduced axial stiffnesses in cross-ply laminates with matrix cracks. The axial crack opening displacement is explicitly expressed in terms of the observable axial strain and the undamaged material properties. A crack parameter representing the effect of matrix cracks on the observable axial Young's modulus is calculated for glass/epoxy and graphite/epoxy material systems. The results show that the matrix crack opening displacement and the effective Young's modulus depend not on the crack length, but on its ratio to the crack spacing.

Lee, Jong-Won↗
NASA NTRS1947

On the one-dimensional theory of steady compressible fluid flow in ducts with friction and heat addition

Steady, diabatic (nonadiabatic), frictional, variable-area flow of a compressible fluid is treated in differential form on the basis of the one-dimensional approximation. The basic equations are first stated in terms of pressure, temperature, density, and velocity of the fluid. Considerable simplification and unification of the equations are then achieved by choosing the square of the local Mach number as one of the variables to describe the flow. The transformed system of equations thus obtained is first examined with regard to the existence of a solution. It is shown that, in general, a solution exists whose calculation requires knowledge only of the variation with position of any three of the dependent variables of the system. The direction of change of the flow variables can be obtained directly from the transformed equations without integration. As examples of this application of the equations, the direction of change of the flow variables is determined for two special flows. In the particular case when the local Mach number m = 1, a special condition must be satisfied by the flow if a solution is to exist. This condition restricts the joint rate of variation of heating, friction, and area at m = 1. Further analysis indicates that when a solution exists at this point it is not necessarily unique. Finally it is shown that the physical phenomenon of choking, which is known to occur in certain simple flow situations, is related to restrictions imposed on the variables by the form of the transformed equations. The phenomenon of choking is thus given a more general significance in that the transformed equations apply to a more general type of flow than has hitherto been treated. (author)

HIGH-SPEED↗
NASA NTRS1992

Materials, structures, and devices for high-speed electronics

Advances in materials, devices, and instrumentation made under this grant began with ex-situ null ellipsometric measurements of simple dielectric films on bulk substrates. Today highly automated and rapid spectroscopic ellipsometers are used for ex-situ characterization of very complex multilayer epitaxial structures. Even more impressive is the in-situ capability, not only for characterization but also for the actual control of the growth and etching of epitaxial layers. Spectroscopic ellipsometry has expanded from the research lab to become an integral part of the production of materials and structures for state of the art high speed devices. Along the way, it has contributed much to our understanding of the growth characteristics and material properties. The following areas of research are summarized: Si3N4 on GaAs, null ellipsometry; diamondlike carbon films; variable angle spectroscopic ellipsometry (VASE) development; GaAs-AlGaAs heterostructures; Ta-Cu diffusion barrier films on GaAs; GaAs-AlGaAs superlattices and multiple quantum wells; superconductivity; in situ elevated temperature measurements of III-V's; optical constants of thermodynamically stable InGaAs; doping dependence of optical constants of GaAs; in situ ellipsometric studies of III-V epitaxial growth; photothermal spectroscopy; microellipsometry; and Si passivation and Si/SiGe strained-layer superlattices.

Woollam, John A.↗
NASA NTRS2009

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗
NASA NTRS1983

The evolution equations for Taylor vortices in the small gap limit

The centrifugal instability of the viscous fluid flow between concentric circular cylinders in the small gap limit is considered. The amplitude of the Taylor vortex is allowed to depend on a slow time variable, a slow axial variable, and the polar angle. It is shown that the amplitude of the vortex cannot in general be described by a single amplitude equation. However, if the axial variations are periodic a single amplitude equation can be derived. In the absence of any slow axial variations it is shown that a Taylor vortex remains stable to wavy vortex perturbations. Furthermore, in this situation, stable nonaxisymmetric modes can occur but do not bifurcate from the Taylor vortex state. The stability of these modes is shown to be governed by a modified form of the Eckhaus criterion.

Hall, P.↗
NASA NTRS1984

Evolution equations for Taylor vortices in the small-gap limit

The centrifugal instability of the viscous fluid flow between concentric circular cylinders in the small gap limit is considered. The amplitude of the Taylor vortex is allowed to depend on a slow time variable, a slow axial variable, and the polar angle. It is shown that the amplitude of the vortex cannot in general be described by a single amplitude equation. However, if the axial variations are periodic a single amplitude equation can be derived. In the absence of any slow axial variations it is shown that a Taylor vortex remains stable to wavy vortex perturbations. Furthermore, in this situation, stable nonaxisymmeric modes can occur but do not bifurcate from the Taylor vortex state. The stability of these modes is shown to be governed by a modified form of the Eckhaus criterion. Previously announced in STAR as N84-10501

Hall, P.↗
NASA NTRS1988

Developing fluid flow in a curved duct of square cross-section and its fully developed dual solutions

Developing fluid flow in a curved duct of square cross-section is studied numerically by a factored ADI finite-difference method on a staggered grid. A central-difference scheme with primitive variables is used inside the computational domain to reduce numerical diffusion. Two Reynolds numbers, 574 and 790, based upon a bulk velocity and hydraulic diameter are chosen for curvature ratios of 1/6.45 and 1/2.3, respectively. It is found that the secondary flow is far more complicated than expected, with the appearance of at least two pairs of vortices. Main-flow separation is also observed for the higher curvature ratio. Furthermore, it is observed that the flow develops into two quite different states downstream, depending upon the inlet conditions. Solutions of the fully developed Navier-Stokes equations is shown to be not unique beyond a certain critical Reynolds number. Developing flow seems to evolve into the fully developed state along a particular branch into which the fully developed solution bifurcates.

Soh, W. Y.↗
NASA NTRS2009

Adjoint-Based Design of Rotors using the Navier-Stokes Equations in a Noninertial Reference Frame

Optimization of rotorcraft flowfields using an adjoint method generally requires a time-dependent implementation of the equations. The current study examines an intermediate approach in which a subset of rotor flowfields are cast as steady problems in a noninertial reference frame. This technique permits the use of an existing steady-state adjoint formulation with minor modifications to perform sensitivity analyses. The formulation is valid for isolated rigid rotors in hover or where the freestream velocity is aligned with the axis of rotation. Discrete consistency of the implementation is demonstrated using comparisons with a complex-variable technique, and a number of single- and multi-point optimizations for the rotorcraft figure of merit function are shown for varying blade collective angles. Design trends are shown to remain consistent as the grid is refined.

Nielsen, Eric J.↗
NASA NTRS2014

Background Error Covariance Estimation Using Information from a Single Model Trajectory with Application to Ocean Data Assimilation

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory.SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Error Covariance↗
NASA NTRS2014

Background Error Covariance Estimation using Information from a Single Model Trajectory with Application to Ocean Data Assimilation into the GEOS-5 Coupled Model

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory. SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Data Assimilation↗
NASA NTRS1998

Self-Similar Compressible Free Vortices

Lie group methods are used to find both exact and numerical similarity solutions for compressible perturbations to all incompressible, two-dimensional, axisymmetric vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for which the solutions are a set of two-dimensional, self-similar, incompressible vortices. These solutions are augmented by deriving a conserved quantity for each eigenvalue, and identifying a Lie group which leaves the reference flow equations invariant. The partial differential equations governing the compressible perturbations to these reference flows are also invariant under the action of the same group. The similarity variables found with this group are used to determine the decay rates of the velocities and thermodynamic variables in the self-similar flows, and to reduce the governing partial differential equations to a set of ordinary differential equations. The ODE's are solved analytically and numerically for a Taylor vortex reference flow, and numerically for an Oseen vortex reference flow. The solutions are used to examine the dependencies of the temperature, density, entropy, dissipation and radial velocity on the Prandtl number. Also, experimental data on compressible free vortex flow are compared to the analytical results, the evolution of vortices from initial states which are not self-similar is discussed, and the energy transfer in a slightly-compressible vortex is considered.

vonEllenrieder, Karl↗
NASA NTRS2010

Adjoint-Based Design of Rotors Using the Navier-Stokes Equations in a Noninertial Reference Frame

Optimization of rotorcraft flowfields using an adjoint method generally requires a time-dependent implementation of the equations. The current study examines an intermediate approach in which a subset of rotor flowfields are cast as steady problems in a noninertial reference frame. This technique permits the use of an existing steady-state adjoint formulation with minor modifications to perform sensitivity analyses. The formulation is valid for isolated rigid rotors in hover or where the freestream velocity is aligned with the axis of rotation. Discrete consistency of the implementation is demonstrated by using comparisons with a complex-variable technique, and a number of single- and multipoint optimizations for the rotorcraft figure of merit function are shown for varying blade collective angles. Design trends are shown to remain consistent as the grid is refined.

Nielsen, Eric J.↗
NASA NTRS1977

Finite state modeling of aeroelastic systems

A general theory of finite state modeling of aerodynamic loads on thin airfoils and lifting surfaces performing completely arbitrary, small, time-dependent motions in an airstream is developed and presented. The nature of the behavior of the unsteady airloads in the frequency domain is explained, using as raw materials any of the unsteady linearized theories that have been mechanized for simple harmonic oscillations. Each desired aerodynamic transfer function is approximated by means of an appropriate Pade approximant, that is, a rational function of finite degree polynomials in the Laplace transform variable. The modeling technique is applied to several two dimensional and three dimensional airfoils. Circular, elliptic, rectangular and tapered planforms are considered as examples. Identical functions are also obtained for control surfaces for two and three dimensional airfoils.

Vepa, R.↗