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.

SourceResearch connectionFrom yearThrough year

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

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Application of a new finite difference algorithm for computational aeroacoustics

Acoustic problems have become extremely important in recent years because of research efforts such as the High Speed Civil Transport program. Computational aeroacoustics (CAA) requires a faithful representation of wave propagation over long distances, and needs algorithms that are accurate and boundary conditions that are unobtrusive. This paper applies a new finite difference method and boundary algorithm to the Linearized Euler Equations (LEE). The results demonstrate the ability of a new fourth order propagation algorithm to accurately simulate the genuinely multidimensional wave dynamics of acoustic propagation in two space dimensions with the LEE. The results also show the ability of a new outflow boundary condition and fourth order algorithm to pass the evolving solution from the computational domain with no perceptible degradation of the solution remaining within the domain.

Goodrich, John W.↗

NASA NTRS1994

The simulated Indian monsoon: A GCM sensitivity study

A series of sensitivity experiments are conducted in an attempt to understand and correct deficiencies in the simulation of the seasonal mean Indian monsoon with a global atmospheric general circulation model. The seasonal mean precipitation is less than half that observed. This poor simulation in seasonal integrations is independent of the choice of initial conditions and global sea surface temperature data used. Experiments are performed to test the sensitivity of the Indian monsoon simulation to changes in orography, vegetation, soil, wetness, and cloudiness. The authors find that the deficiency of the model precipitation simulation may be attributed to the use of an enhanced orography in the integrations. Replacement of this orography with a mean orography results in a much more realistic simulation of Indian monsoon circulation and rainfall. Experiments with a linear primitive equation model on the sphere suggest that this striking improvement is due to modulations of the orographically forced waves in the lower troposphere. This improvement in the monsoon simulation is due to the kinematic and dynamical effects of changing the topography, rather than the thermal effects, which were minimal. The magnitude of the impact on the Indian monsoon of the other sensitivity experiments varied considerably, but was consistently less than the impact of using the mean orography. However, results from the soil moisture sensitivity experiments suggest a possibly important role for soil moisture in simulating tropical precipitation, including that associated with the Indian monsoon.

Fennessy, M. J.↗

NASA NTRS1994

Planetary-scale circulations in the presence of climatological and wave-induced heating

Interaction between the large-scale circulation and the convective pattern is investigated in a coupled system governed by the linearized primitive equations. Convection is represented in terms of two components of heating: A 'climatological component' is prescribed stochastically to represent convection that is maintained by fixed distributions of land and sea and sea surface temperature (SST). An 'induced component' is defined in terms of the column-integrated moisture flux convergence to represent convection that is produced through feedback with the circulation. Each component describes the envelope organizing mesoscale convective activity. As SST on the equator is increased, induced heating amplifies in the gravest zonal wavenumbers at eastward frequencies, where positive feedback offsets dissipation. Under barotropic stratification, a critical SST of 29.5 C results in positive feedback exactly cancelling dissipation in wavenumber 1 for an eastward phase speed of 6 m/s. Sympathetic interaction between the circulation and the induced heating is the basis for 'frictional wave-Conditional Instability of the Second Kind (CISK)', which is distinguished from classical wave-CISK by rendering the gravest zonal dimensions most unstable. Under baroclinic stratification, the coupled system exhibits similar behavior. The critical SST is only 26.5 C for conditions representative of equinox, but in excess of 30 C for conditions representative of solstice. Having the form of an unsteady Walker circulation, the disturbance produced by frictional wave-CISK compares favorably with the observed life cycle of the Madden-Julian oscillation (MJO). SST above the critical value produces an amplifying disturbance in which enhanced convection coincides with upper-tropospheric westerlies and is positively correlated with temperature and surface convergence. Conversely, SST below the critical value produces a decaying disturbance in which enhanced convection coincides with upper-tropospheric easterlies and is nearly in quadrature with temperature and surface convergence. While sharing essential features with the MJO in the Eastern Hemisphere, frictional wave-CISK does not explain observed behavior in the Western Hemisphere, where the convective signal is largely absent. Comprised of Kelvin structure with the same frequency, observed behavior in the Western Hemisphere can be understood as a propagating response that is excited in and radiates away from the fluctuation of convection in the Eastern Hemisphere.

Salby, Murry L↗

NASA NTRS1994

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗

NASA NTRS1994

Practical formulation of a positively conservative scheme

Approximate Riemann solvers have been highly successful for computing the Euler/Navier-Stokes equations, but linearized Riemann solvers are known to fail occasionally by predicting non-physical states with negative density or internal energy. Positively conservative schemes, in contrast, guarantee physical solutions from realistic input. The Harten-Lax-van Leer-Einfeldt (HLLE) scheme is a typical example of a positively conservative scheme. However, the HLLE scheme is highly dissipative at contact discontinuities and shear layers and thus it is not applicable to practicle simulations. An existing modification to the HLLE scheme, known as HLLEM, enhances the resolution to that of the Roe scheme. However, this modification violates the positivity of density and internal energy. Precise derivation of the modification yields a quatratic inequality and thus requires a case-by-case treatment. This Note describes a new, modified HLLE scheme that satisfies the positively conservative condition approximately. Sample computationa are included to demonstrate the resolution and the robustness of the scheme.

Obayashi, Shigeru↗

NASA NTRS1995

Relaxation approximation in the theory of shear turbulence

Leslie's perturbative treatment of the direct interaction approximation for shear turbulence (Modern Developments in the Theory of Turbulence, 1972) is applied to derive a time dependent model for the Reynolds stresses. The stresses are decomposed into tensor components which satisfy coupled linear relaxation equations; the present theory therefore differs from phenomenological Reynolds stress closures in which the time derivatives of the stresses are expressed in terms of the stresses themselves. The theory accounts naturally for the time dependence of the Reynolds normal stress ratios in simple shear flow. The distortion of wavenumber space by the mean shear plays a crucial role in this theory.

Rubinstein, Robert↗

NASA NTRS1995

A far-field non-reflecting boundary condition for two-dimensional wake flows

Far-field boundary conditions for external flow problems have been developed based upon long-wave perturbations of linearized flow equations about a steady state far field solution. The boundary improves convergence to steady state in single-grid temporal integration schemes using both regular-time-stepping and local-time-stepping. The far-field boundary may be near the trailing edge of the body which significantly reduces the number of grid points, and therefore the computational time, in the numerical calculation. In addition the solution produced is smoother in the far-field than when using extrapolation conditions. The boundary condition maintains the convergence rate to steady state in schemes utilizing multigrid acceleration.

Danowitz, Jeffrey S.↗

NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

NASA NTRS1995

The relation between the period of oscillations and brightness in chromospheric bright points

The chromospheric bright points are the sites where intense heating occurs in three minute period waves. The bright points are grouped into three classes depending on the amount of intensity enhancement and the pattern of their dynamical evolution. A 35-minute time series of photographic spectra in the Ca(II) H line on a quiet region ofthe center of the solar disk was used to show that the period of intensity oscillations seen at sites of the bright points is independent of their intensity enhancements. The series was also used to show that the period may not depend on the strength of the magnetic fields with which they are associated. A linear regression equation was fitted to a curve representing the variation of the period of intensity oscillations with the peak value of I(sub H2V). The correlation coefficient was found to be 0.19.

Kariyappa, R.↗

NASA NTRS1996

A dynamic response model for pressure sensors in continuum and high Knudsen number flows with large temperature gradients

This paper develops a dynamic model for pressure sensors in continuum and rarefied flows with longitudinal temperature gradients. The model was developed from the unsteady Navier-Stokes momentum, energy, and continuity equations and was linearized using small perturbations. The energy equation was decoupled from momentum and continuity assuming a polytropic flow process. Rarefied flow conditions were accounted for using a slip flow boundary condition at the tubing wall. The equations were radially averaged and solved assuming gas properties remain constant along a small tubing element. This fundamental solution was used as a building block for arbitrary geometries where fluid properties may also vary longitudinally in the tube. The problem was solved recursively starting at the transducer and working upstream in the tube. Dynamic frequency response tests were performed for continuum flow conditions in the presence of temperature gradients. These tests validated the recursive formulation of the model. Model steady-state behavior was analyzed using the final value theorem. Tests were performed for rarefied flow conditions and compared to the model steady-state response to evaluate the regime of applicability. Model comparisons were excellent for Knudsen numbers up to 0.6. Beyond this point, molecular affects caused model analyses to become inaccurate.

Whitmore, Stephen A.↗

NASA NTRS1996

On buffer layers as non-reflecting computational boundaries

We examine an absorbing buffer layer technique for use as a non-reflecting boundary condition in the numerical simulation of flows. One such formulation was by Ta'asan and Nark for the linearized Euler equations. They modified the flow inside the buffer zone to artificially make it supersonic in the layer. We examine how this approach can be extended to the nonlinear Euler equations. We consider both a conservative and a non-conservative form modifying the governing equations in the buffer layer. We compare this with the case that the governing equations in the layer are the same as in the interior domain. We test the effectiveness of these buffer layers by a simulation of an excited axisymmetric jet based on a nonlinear compressible Navier-Stokes equations.

Hayder, M. Ehtesham↗

NASA NTRS1995

A kinetic study of solar wind electrons in the transition region from collision dominated to collisionless flow

We have studied the evolution of the velocity distribution function of a test population of electrons in the solar corona and inner solar wind region, using a recently developed kinetic model. The model solves the time dependent, linear transport equation, with a Fokker-Planck collision operator to describe Coulomb collisions between the 'test population' and a thermal background of charged particles, using a finite differencing scheme. The model provides information on how non-Maxwellian features develop in the distribution function in the transition region from collision dominated to collisionless flow. By taking moments of the distribution the evolution of higher order moments, such as the heat flow, can be studied.

Lie-Svendsen, O.↗

NASA NTRS1996

Characterization of the Dynamic Material Properties of Magnetostrictive Terfenol-D

A major limitation in use of electromagnetic and/or magnetomechanical models for design of Terfenol-D actuators is the lack of reliable material property data for Terfenol-D. In particular data on the performance of Terfenol-D as employed in a transducer, operating under real world dynamic conditions is needed. To provide this information, Terfenol-D rod properties need to be measured under as run prestressed and magnetically biased states. Using a Terfenol-D actuator, the following properties can be measured and/or calculated: mechanical quality factor, speed of sound in the material, the resonant frequency, the anti-resonant frequency, two magnetic permeabilities (one at constant stress and one at constant strain), two Young's moduli (one at constant amplitude applied magnetic field and one at constant amplitude magnetic flux density in the material), the magnetomechanical coupling, and the axial strain coefficient. The development of the material properties measurements and calculations is based on the model of low signal, linear, magnetostriction from Clark, the linear transduction equations for a transducer from Hunt, and a one degree of freedom mechanical model of the transducer. The electrical impedance and admittance mobility loops are used to determine the resonant, anti-resonant, and half power point frequencies. The rest of the material properties indicated above can then be calculated using these frequencies, acceleration from an accelerometer mounted on the actuator arm, and readily measurable transducer and Terfenol-D rod parameters.

Calkins, Frederick T.↗

NASA NTRS1996

X-29A Lateral-Directional Stability and Control Derivatives Extracted From High-Angle-of-Attack Flight Data

The lateral-directional stability and control derivatives of the X-29A number 2 are extracted from flight data over an angle-of-attack range of 4 degrees to 53 degrees using a parameter identification algorithm. The algorithm uses the linearized aircraft equations of motion and a maximum likelihood estimator in the presence of state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics over the aircraft at angles of attack above 15 degrees. The results supported the flight-envelope-expansion phase of the X-29A number 2 by helping to update the aerodynamic mathematical model, to improve the real-time simulator, and to revise flight control system laws. Effects of the aircraft high gain flight control system on maneuver quality and the estimated derivatives are also discussed. The derivatives are plotted as functions of angle of attack and compared with the predicted aerodynamic database. Agreement between predicted and flight values is quite good for some derivatives such as the lateral force due to sideslip, the lateral force due to rudder deflection, and the rolling moment due to roll rate. The results also show significant differences in several important derivatives such as the rolling moment due to sideslip, the yawing moment due to sideslip, the yawing moment due to aileron deflection, and the yawing moment due to rudder deflection.

Iliff, Kenneth W.↗

NASA NTRS1996

Dust Coagulation in Protoplanetary Accretion Disks

The time evolution of dust particles in circumstellar disk-like structures around protostars and young stellar objects is discussed. In particular, we consider the coagulation of grains due to collisional aggregation. The coagulation of the particles is calculated by solving numerically the non-linear Smoluchowski equation. The different physical processes leading to relative velocities between the grains are investigated. The relative velocities may be induced by Brownian motion, turbulence and drift motion. Starting from different regimes which can be identified during the grain growth we also discuss the evolution of dust opacities. These opacities are important for both the derivation of the circumstellar dust mass from submillimeter/millimeter continuum observations and the dynamical behavior of the disks. We present results of our numerical studies of the coagulation of dust grains in a turbulent protoplanetary accretion disk described by a time-dependent one-dimensional (radial) alpha-model. For several periods and disk radii, mass distributions of coagulated grains have been calculated. From these mass spectra, we determined the corresponding Rosseland mean dust opacities. The influence of grain opacity changes due to dust coagulation on the dynamical evolution of a protostellar disk is considered. Significant changes in the thermal structure of the protoplanetary nebula are observed. A 'gap' in the accretion disk forms at the very frontier of the coagulation, i.e., behind the sublimation boundary in the region between 1 and 5 AU.

Schmitt, W.↗

NASA NTRS1997

Extraction of Lateral-Directional Stability and Control Derivatives for the Basic F-18 Aircraft at High Angles of Attack

The results of parameter identification to determine the lateral-directional stability and control derivatives of an F-18 research aircraft in its basic hardware and software configuration are presented. The derivatives are estimated from dynamic flight data using a specialized identification program developed at NASA Dryden Flight Research Center. The formulation uses the linearized aircraft equations of motions in their continuous/discrete form and a maximum likelihood estimator that accounts for both state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics, such as separated and vortical flows, over the aircraft. The derivatives are plotted as functions of angle of attack between 3 deg and 47 deg and compared with wind-tunnel predictions. The quality of the derivative estimates obtained by parameter identification is somewhat degraded because the maneuvers were flown with the aircraft's control augmentation system engaged, which introduced relatively high correlations between the control variables and response variables as a result of control motions from the feedback control system.

Iliff, Kenneth W.↗

NASA NTRS1996

Inhomogeneous Radiation Boundary Conditions Simulating Incoming Acoustic Waves for Computational Aeroacoustics

A set of nonhomogeneous radiation and outflow conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as nondispersive ans nondissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is necessary The proposed nonhomogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes.

Tam, Christopher K. W.↗

NASA NTRS1997

Numerical Simulation of the Effect of Heating on Supersonic Jet Noise

The axisymmetric linearized Euler equations are used to simulate noise amplification and radiation from a supersonic jet. The effect of heating on the noise field of the jet is studied and compared to experimental results. Special attention was given to boundary treatment, and the resulting solution is stable and nearly free from boundary reflections.

Hixon, R.↗