Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equations”

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 1,315 records · Page 73

Pitch angle diffusion of newborn ions due to intrinsic turbulence in the solar wind

The objective of the present study is to understand the interaction of the solar wind with newborn ions in far upstream regions of a comet where the level of intrinsic turbulence is moderately low. Based on the assumption that quasi-linear theory is adequate and applicable, the pitch angle diffusion process and the time evolution of the newborn ion distribution function are investigated. Numerical solutions to the quasi-linear diffusion equation, including the effect of resonance broadening and that of continuous creation of newborn ions due to the ionization process, are obtained under several assumptions and approximations. It is found that theoretical results are consistent with the Giotto observations recently reported by Neugebauer et al. (1989).

Ziebell, L. F.↗

Normal mode projections of GCM response to tropical forcing

The transient response of a general circulation model (GCM) is examined for the first 10 days of integration after a tropical heating source is switched on. Model output from 10 control and 10 experiment GCM runs are projected onto the normal modes of an adiabatic primitive equation model linearized about a basic state at rest to quantify the contributions from different modes to the 200 mb wind response. The response in the vicinity of the source is similar to that predicted by models linearized about uniform basic states, with important contributions from the internal Rossby and Kelvin modes in the tropics. The heating also excites internal gravity modes with contributions to the 200 mb winds about 40 percent of those of the internal Rossby modes through the 10-day integration time. The response in the external Rossby modes in the extratropics becomes statistically significant after five days.

Zhong, Weili↗

Galerkin approximations for dissipative magnetohydrodynamics

A Galerkin approximation scheme is proposed for voltage-driven, dissipative magnetohydrodynamics. The trial functions are exact eigenfunctions of the linearized continuum equations and represent helical deformations of the axisymmetric, zero-flow, driven steady state. The lowest nontrivial truncation is explored: one axisymmetric trial function and one helical trial function each for the magnetic and velocity fields. The system resembles the Lorenz approximation to Benard convection, but in the region of believed applicability, its dynamical behavior is rather different, including relaxation to a helically deformed state similar to those that have emerged in the much higher resolution computations of Dahlburg et al.

Chen, Hudong↗

An approach to optimal guidance of an advanced launch vehicle concept

An approximate solution for the maximum payload trajectory of a two-stage launch vehicle using a regular perturbation technique is presented. A zero-order solution for a two-stage vehicle based on a flat-earth approximation and negligible atmospheric effects is obtained in closed form. High-order correction terms are obtained from the solution of nonhomogeneous, first-order linear differential equations by quadrature. This promises the capability for an onboard optimal guidance law implementation.

Leung, Martin S.↗

Convective instability in superposed fluid and anisotropic porous layers

The onset of thermal convection due to heating from below in a system consisting of a fluid layer overlying a porous layer with anisotropic permeability and thermal diffusivity is studied. The linear perturbation equations are solved numerically. It is found that the effects of anisotropy on the onset of thermal convection are most profound for small values of the depth ratio xi (ratio of fluid layer thickness to porous layer thickness), since in that case, the onset of convection corresponds to significant motion in both layers. For fixed vlues of the vertical permeability in the porous medium, decreasing the value of zeta (ratio of horizontal to vertical permeability) leads to stabilization of the superposed layer configuration because of increased resistance to motion in the porous medium. For larger values of xi, the onset of motion is increasingly confined to the fluid layer, with the transport of heat through the porous layer occurring primarily by conduction. Accordingly, the influence of zeta on the stability characteristics for larger xi is less significant than the effects of an anisotropic thermal conductivity.

Chen, Falin↗

Unsteady supersonic flow around delta wings with symmetric and asymmetric flaps oscillation

A parametric study is presented to investigate the effect of reduced frequency of the leading-edge flaps on the locally-conical, unsteady, supersonic flow around a delta wing. This study covers symmetric and antisymmetric forced oscillation of the leading-edge flaps. The effects of the freestream Mach number and angle of attack are also presented. The problem is solved using time-accurate integration of the unsteady, compressible, thin-layer Navier-Stokes equations and the unsteady, linearized, Navier-displacement equations. The delta wing is of aspect ratio of 1.5 and its leading-edge flaps are hinged at 65 percent of the local-half span length. The reduced frequency is varied between 2 pi and pi/2. Two supersonic flow conditions have been investigated; the first is for a freestream Mach number of 2.4 and an angle of attack of 19 deg and the second is for a freestream Mach number of 1.5 and an angle attack of 15 deg.

Kandil, Osama A.↗

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗

Probabilistic methods for rotordynamics analysis

This paper summarizes the development of the methods and a computer program to compute the probability of instability of dynamic systems that can be represented by a system of second-order ordinary linear differential equations. Two instability criteria based upon the eigenvalues or Routh-Hurwitz test functions are investigated. Computational methods based on a fast probability integration concept and an efficient adaptive importance sampling method are proposed to perform efficient probabilistic analysis. A numerical example is provided to demonstrate the methods.

Wu, Y.-T.↗

Analysis of radiative and phase-change phenomena with application to space-based thermal energy storage

The simplified geometry for the analysis is an infinite, axis symmetric annulus with a specified solar flux at the outer radius. The inner radius is either adiabatic (modeling Flight Experiment conditions), or convective (modeling Solar Dynamic conditions). Liquid LiF either contacts the outer wall (modeling ground based testing), or faces a void gap at the outer wall (modeling possible space based conditions). The analysis is presented in three parts: Part 3 considers and adiabatic inner wall and linearized radiation equations; part 2 adds effects of convection at the inner wall; and part 1 includes the effect of the void gap, as well as previous effects, and develops the radiation model further. The main results are the differences in melting behavior which can occur between ground based 1 g experiments and the microgravity flight experiments. Under 1 gravity, melted PCM will always contact the outer wall having the heat flux source, thus providing conductance from this source to the phase change front. In space based tests where a void gap may likely form during solidification, the situation is reversed; radiation is now the only mode of heat transfer and the majority of melting takes place from the inner wall.

Lund, Kurt O.↗

Field quantization and squeezed states generation in resonators with time-dependent parameters

The problem of electromagnetic field quantization is usually considered in textbooks under the assumption that the field occupies some empty box. The case when a nonuniform time-dependent dielectric medium is confined in some space region with time-dependent boundaries is studied. The basis of the subsequent consideration is the system of Maxwell's equations in linear passive time-dependent dielectric and magnetic medium without sources.

Dodonov, V. V.↗

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 2: User's manual (version 3.0)

A comprehensive description of user problem definition for the PAN AIR (Panel Aerodynamics) system is given. PAN AIR solves the 3-D linear integral equations of subsonic and supersonic flow. Influence coefficient methods are used which employ source and doublet panels as boundary surfaces. Both analysis and design boundary conditions can be used. This User's Manual describes the information needed to use the PAN AIR system. The structure and organization of PAN AIR are described, including the job control and module execution control languages for execution of the program system. The engineering input data are described, including the mathematical and physical modeling requirements. Version 3.0 strictly applies only to PAN AIR version 3.0. The major revisions include: (1) inputs and guidelines for the new FDP module (which calculates streamlines and offbody points); (2) nine new class 1 and class 2 boundary conditions to cover commonly used modeling practices, in particular the vorticity matching Kutta condition; (3) use of the CRAY solid state Storage Device (SSD); and (4) incorporation of errata and typo's together with additional explanation and guidelines.

Sidwell, Kenneth W.↗

Application of optical distributed sensing and computation to control of large space structures

A real time holographic sensing technique is introduced and its advantages are investigated from the filtering and control point of view. A feature of holographic sensing is its capability to make distributed measurements of the position and velocity of moving objects, such as a vibrating flexible space structure. This work is based upon the distributed parameter models of linear time invariant systems, particularly including the linear oscillator equations describing the vibration of large flexible space structures. The general conclusion is that application of optical distributed sensors bring gains in the situation where Kalman filtering is necessary for state estimation. In this case, both steady state and transient filtering error covariance become smaller. This in turn results in smaller cost in the LQG problem.

Balakrishnan, A. V.↗

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗

The mirror and ion cyclotron anisotropy instabilities

The linear dispersion equation for fully electromagnetic waves and instabilities at arbitrary directions of propagation relative to a background magnetic field B(0) in a homogeneous Vlasov plasma is solved numerically for bi-Maxwellian particle distributions. For isotropic plasmas the dispersion and damping of the three modes below the proton cyclotron frequency are studied as functions of Beta(i) and T(e)/T(i). The transport ratios of helicity, cross-helicity, Alfven ratio, compressibility, and parallel compressibility are defined. Under the condition that the proton temperature perpendicular to B(0) is greater than the parallel temperature, the growth rates and transport ratios of the mirror instability and the ion cyclotron anisotropy instability are examined and compared. Both the proton parallel compressibility and the proton Alfven ratio are significantly different for the two growing modes.

Gary, S. P.↗

Acoustic modes in fluid networks

Pressure and flow rate eigenvalue problems for one-dimensional flow of a fluid in a network of pipes are derived from the familiar transmission line equations. These equations are linearized by assuming small velocity and pressure oscillations about mean flow conditions. It is shown that the flow rate eigenvalues are the same as the pressure eigenvalues and the relationship between line pressure modes and flow rate modes is established. A volume at the end of each branch is employed which allows any combination of boundary conditions, from open to closed, to be used. The Jacobi iterative method is used to compute undamped natural frequencies and associated pressure/flow modes. Several numerical examples are presented which include acoustic modes for the Helium Supply System of the Space Shuttle Orbiter Main Propulsion System. It should be noted that the method presented herein can be applied to any one-dimensional acoustic system involving an arbitrary number of branches.

Michalopoulos, C. D.↗

Planetary circulations in the presence of transient and self-induced heating

The research program focuses on large-scale circulations and their interaction with the global convective pattern. An 11-year record of global cloud imagery and contemporaneous fields of motion and temperature have been used to investigate organized convection and coherent variability of the tropical circulation operating on intraseasonal time scales. This study provides a detailed portrait of tropical variability associated with the so-called Madden-Julian Oscillation (MJO). It reveals the nature, geographical distribution, and seasonality of discrete convective signal, which is a measure of feedback between the circulation and the convective pattern. That discrete spectral behavior has been evaluated in light of natural variability of the ITCZ associated with climatological convection. A composite signature of the MJO, based on cross-covariance statistics of cloud cover, motion, and temperature, has been constructed to characterize the lifecycle of the disturbance in terms of these properties. The composite behavior has also been used to investigate the influence the MJO exerts on the zonal-mean circulation and the involvement of the MJO in transfers of momentum between the atmosphere and the solid Earth. The aforementioned observational studies have led to the production of two animations. One reveals the convective signal in band-pass filtered OLR and compares it to climatological convection. The other is a 3-dimensional visualization of the composite lifecycle of the MJO. With a clear picture of the MJO in hand, feedback between the circulation and the convective pattern can be diagnosed meaningfully in numerical simulations. This process is being explored in calculations with the linearized primitive equations on the sphere in the presence of realistic stability and shear. The numerical framework represents climatological convection as a space-time stochastic process and wave-induced convection in terms of the vertically-integrated moisture flux convergence. In these calculations, frictional convergence near the equator emerges as a key to feedback between the circulation and the convective pattern. At low latitudes, nearly geostrophic balance in the boundary layer gives way to frictional balance. This shifts the wave-induced convection into phase with the temperature anomaly and allows the attending heating to feed back positively onto the circulation. The calculations successfully reproduce the salient features of the MJO. They are being used to understand the growth and decay phases of the composite lifecycle and the conditions that favor amplification of the MJO.

Salby, Murry L.↗

Optimal fixed-finite-dimensional compensator for Burgers' equation with unbounded input/output operators

The problem of using reduced order dynamic compensators to control a class of nonlinear parabolic distributed parameter systems was considered. Concentration was on a system with unbounded input and output operators governed by Burgers' equation. A linearized model was used to compute low-order-finite-dimensional control laws by minimizing certain energy functionals. Then these laws were applied to the nonlinear model. Standard approaches to this problem employ model/controller reduction techniques in conjunction with linear quadratic Gaussian (LQG) theory. The approach used is based on the finite dimensional Bernstein/Hyland optimal projection theory which yields a fixed-finite-order controller.

Burns, John A.↗

On the establishment of stationary waves in the Northern Hemisphere winter

The establishment of stationary waves in the Northern Hemisphere winter is investigated using stationary and time-dependent linear primitive equation models. Confirming the results of Nigam and Lindzen, we find that small displacements of the subtropical jet can cause significant changes in the stationary-wave response. The time scale for stationary establishment is found to be on the order of 5 days, both in the troposphere and in the lower stratosphere. The exception is for a northward shift of the subtropical jet, in which case the establishment of the new stationary solution in the stratosphere occurs on a longer time scale, which is mainly determined by dissipation. Implications for low-frequency atmospheric variability and mid- and long-range weather forecasting are discussed.

Da Silva, Arlindo M.↗