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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Time-Domain Impedance Boundary Conditions for Computational Aeroacoustics

It is an accepted practice in aeroacoustics to characterize the properties of an acoustically treated surface by a quantity known as impedance. Impedance is a complex quantity. As such, it is designed primarily for frequency-domain analysis. Time-domain boundary conditions that are the equivalent of the frequency-domain impedance boundary condition are proposed. Both single frequency and model broadband time-domain impedance boundary conditions are provided. It is shown that the proposed boundary conditions, together with the linearized Euler equations, form well-posed initial boundary value problems. Unlike ill-posed problems, they are free from spurious instabilities that would render time-marching computational solutions impossible.

Tam, Christopher K. W.↗
NASA NTRS1997

Flight-Determined Subsonic Longitudinal Stability and Control Derivatives of the F-18 High Angle of Attack Research Vehicle (HARV) with Thrust Vectoring

The subsonic longitudinal stability and control derivatives of the F-18 High Angle of Attack Research Vehicle (HARV) are extracted from dynamic flight data using a maximum likelihood parameter identification technique. The technique uses the linearized aircraft equations of motion in their continuous/discrete form and accounts for state and measurement noise as well as thrust-vectoring effects. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics over the aircraft, particularly at high angles of attack. Thrust vectoring was implemented using electrohydraulically-actuated nozzle postexit vanes and a specialized research flight control system. During maneuvers, a control system feature provided independent aerodynamic control surface inputs and independent thrust-vectoring vane inputs, thereby eliminating correlations between the aircraft states and controls. Substantial variations in control excitation and dynamic response were exhibited for maneuvers conducted at different angles of attack. Opposing vane interactions caused most thrust-vectoring inputs to experience some exhaust plume interference and thus reduced effectiveness. The estimated stability and control derivatives are plotted, and a discussion relates them to predicted values and maneuver quality.

Iliff, Kenneth W.↗
NASA NTRS1998

Stability of Capillary Surfaces in Rectangular Containers: The Right Square Cylinder

The linearized governing equations for an ideal fluid are presented for numerical analysis for the stability of free capillary surfaces in rectangular containers against unfavorable disturbances (accelerations,i.e. Rayleigh-Taylor instability). The equations are solved for the case of the right square cylinder. The results are expressed graphically in term of a critical Bond number as a function of system contact angle. A critical wetting phenomena in the corners is shown to significantly alter the region of stability for such containers in contrast to simpler geometries such as the right circular cylinder or the infinite rectangular slot. Such computational results provide additional constraints for the design of fluids systems for space-based applications.

Weislogel, M. M.↗
NASA NTRS1997

Development of Improved Surface Integral Methods for Jet Aeroacoustic Predictions

The accurate prediction of aerodynamically generated noise has become an important goal over the past decade. Aeroacoustics must now be an integral part of the aircraft design process. The direct calculation of aerodynamically generated noise with CFD-like algorithms is plausible. However, large computer time and memory requirements often make these predictions impractical. It is therefore necessary to separate the aeroacoustics problem into two parts, one in which aerodynamic sound sources are determined, and another in which the propagating sound is calculated. This idea is applied in acoustic analogy methods. However, in the acoustic analogy, the determination of far-field sound requires the solution of a volume integral. This volume integration again leads to impractical computer requirements. An alternative to the volume integrations can be found in the Kirchhoff method. In this method, Green's theorem for the linear wave equation is used to determine sound propagation based on quantities on a surface surrounding the source region. The change from volume to surface integrals represents a tremendous savings in the computer resources required for an accurate prediction. This work is concerned with the development of enhancements of the Kirchhoff method for use in a wide variety of aeroacoustics problems. This enhanced method, the modified Kirchhoff method, is shown to be a Green's function solution of Lighthill's equation. It is also shown rigorously to be identical to the methods of Ffowcs Williams and Hawkings. This allows for development of versatile computer codes which can easily alternate between the different Kirchhoff and Ffowcs Williams-Hawkings formulations, using the most appropriate method for the problem at hand. The modified Kirchhoff method is developed primarily for use in jet aeroacoustics predictions. Applications of the method are shown for two dimensional and three dimensional jet flows. Additionally, the enhancements are generalized so that they may be used in any aeroacoustics problem.

Pilon, Anthony R.↗
NASA NTRS1998

Supersonic Coaxial Jets: Noise Predictions and Measurements

The noise from perfectly expanded coaxial jets was measured in an anechoic chamber for different operating conditions with the same total thrust, mass flow, and exit area. The shape of the measured noise spectrum at different angles to the jet axis was found to agree with spectral shapes for single, axisymmetric jets. Based on these spectra, the sound was characterized as being generated by large turbulent structures or fine-scale turbulence. Modeling the large scale structures as instability waves, a stability analysis was conducted for the coaxial jets to identify the growing and decaying instability waves in each shear layer and predict their noise radiation pattern outside the jet. When compared to measured directivity, the analysis identified the region downstream of the outer potential core, where the two shear layers were merging, as the source of the peak radiated noise where instability waves, with their origin in the inner shear layer, reach their maximum amplitude. Numerical computations were also performed using a linearized Euler equation solver. Those results were compared to both the results from the instability wave analysis and to measured data.

Dahl, Milo D.↗
NASA NTRS1997

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗
NASA NTRS1961

A Method for Longitudinal and Lateral Range Control for a High-Drag Low-Lift Vehicle Entering the Atmosphere of a Rotating Earth

A study has been made of a method for controlling the trajectory of a high-drag low-lift entry vehicle to a desired longitude and latitude on the surface of a rotating earth. By use of this control technique the vehicle can be guided to the desired point when the present position and heading of the vehicle are known and the desired longitude and latitude are specified. The present study makes use of a single reference trajectory and an estimate of the lift and side-force capabilities of the vehicle. This information is stored in a control-logic system and used with linear control equations to guide the vehicle to the desired destination. Results are presented of a number of trajectory studies which describe the operation of the control system and illustrate its ability to control the vehicle trajectory to the desired landing area.

Young, John W.↗
NASA NTRS1999

Real-Time Simulation of the X-33 Aerospace Engine

This paper discusses the development and performance of the X-33 Aerospike Engine RealTime Model. This model was developed for the purposes of control law development, six degree-of-freedom trajectory analysis, vehicle system integration testing, and hardware-in-the loop controller verification. The Real-Time Model uses time-step marching solution of non-linear differential equations representing the physical processes involved in the operation of a liquid propellant rocket engine, albeit in a simplified form. These processes include heat transfer, fluid dynamics, combustion, and turbomachine performance. Two engine models are typically employed in order to accurately model maneuvering and the powerpack-out condition where the power section of one engine is used to supply propellants to both engines if one engine malfunctions. The X-33 Real-Time Model is compared to actual hot fire test data and is been found to be in good agreement.

Aguilar, Robert↗
NASA NTRS1999

Martian Aerocapture Terminal Point Guidance: A Reference Path Optimization Study

An effective method of terminal point guidance is to employ influence coefficients, which are solved from a set of differential equations adjoint to the linearized perturbations of the equations of motion about a reference trajectory. Hence, to optimize this type of guidance, one must first optimize the reference trajectory that the guidance is based upon. This study concentrates on various methods to optimize a reference trajectory for a Martian aerocapture maneuver, including a parametric analysis and first order gradient method. Resulting reference trajectories were tested in separate 2000 6-DOF Monte Carlo runs, using the Atmospheric Guidance Algorithm Testbed for the Mars Surveyor Program 2001 (MSP '01) Orbiter. These results were compared to an August 1998 study using the same terminal point control guidance algorithm and simulation testbed. Satisfactory improvements over the 1998 study are amply demonstrated.

Ro, Theodore U.↗
NASA NTRS2000

A Reduced Order Model of the Linearized Incompressible Navier-Strokes Equations for the Sensor/Actuator Placement Problem

A reduced order modeling approach of the Navier-Stokes equations is presented for the design of a distributed optimal feedback kernel. This approach is based oil a Krylov subspace method where significant modes of the flow are captured in the model This model is then used in all optimal feedback control design where sensing and actuation is performed oil tile entire flow field. This control design approach yields all optimal feedback kernel which provides insight into the placement of sensors and actuators in the flow field. As all evaluation of this approach, a two-dimensional shear layer and driven cavity flow are investigated.

Allan, Brian G.↗
NASA NTRS2000

A Very High Order, Adaptable MESA Implementation for Aeroacoustic Computations

Since computational efficiency and wave resolution scale with accuracy, the ideal would be infinitely high accuracy for problems with widely varying wavelength scales. Currently, many of the computational aeroacoustics methods are limited to 4th order accurate Runge-Kutta methods in time which limits their resolution and efficiency. However, a new procedure for implementing the Modified Expansion Solution Approximation (MESA) schemes, based upon Hermitian divided differences, is presented which extends the effective accuracy of the MESA schemes to 57th order in space and time when using 128 bit floating point precision. This new approach has the advantages of reducing round-off error, being easy to program. and is more computationally efficient when compared to previous approaches. Its accuracy is limited only by the floating point hardware. The advantages of this new approach are demonstrated by solving the linearized Euler equations in an open bi-periodic domain. A 500th order MESA scheme can now be created in seconds, making these schemes ideally suited for the next generation of high performance 256-bit (double quadruple) or higher precision computers. This ease of creation makes it possible to adapt the algorithm to the mesh in time instead of its converse: this is ideal for resolving varying wavelength scales which occur in noise generation simulations. And finally, the sources of round-off error which effect the very high order methods are examined and remedies provided that effectively increase the accuracy of the MESA schemes while using current computer technology.

Dydson, Roger W.↗
NASA NTRS2000

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗
NASA NTRS2000

An Automated Approach to Very High Order Aeroacoustic Computations in Complex Geometries

Computational aeroacoustics requires efficient, high-resolution simulation tools. And for smooth problems, this is best accomplished with very high order in space and time methods on small stencils. But the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewslci recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that are located near wall boundaries. These procedures are used to automatically develop and implement very high order methods (>15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗
NASA NTRS2001

An Experimental Evaluation of Generalized Predictive Control for Tiltrotor Aeroelastic Stability Augmentation in Airplane Mode of Flight

The results of a joint NASA/Army/Bell Helicopter Textron wind-tunnel test to assess the potential of Generalized Predictive Control (GPC) for actively controlling the swashplate of tiltrotor aircraft to enhance aeroelastic stability in the airplane mode of flight are presented. GPC is an adaptive time-domain predictive control method that uses a linear difference equation to describe the input-output relationship of the system and to design the controller. The test was conducted in the Langley Transonic Dynamics Tunnel using an unpowered 1/5-scale semispan aeroelastic model of the V-22 that was modified to incorporate a GPC-based multi-input multi-output control algorithm to individually control each of the three swashplate actuators. Wing responses were used for feedback. The GPC-based control system was highly effective in increasing the stability of the critical wing mode for all of the conditions tested, without measurable degradation of the damping in the other modes. The algorithm was also robust with respect to its performance in adjusting to rapid changes in both the rotor speed and the tunnel airspeed.

Kvaternik, Raymond G.↗
NASA NTRS2001

Automated Approach to Very High-Order Aeroacoustic Computations

Computational aeroacoustics requires efficient, high-resolution simulation tools. For smooth problems, this is best accomplished with very high-order in space and time methods on small stencils. However, the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewski recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that am located near wall boundaries. These procedures are used to develop automatically and to implement very high-order methods (> 15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗
NASA NTRS2000

A Low Order and a High Order Solution for a Converging-Diverging Nozzle Problem

The Computational Aeroacoustics Workshops on Benchmark problems are conducted in order to generate solutions with a variety of methods for problems that are chosen both to be representative of significant and relevant acoustic phenomena. It is generally recognized that CAA applications require the resolution in time and in space of solution details over a significantly broader range of scales than typical computational fluid dynamics applications. These requirements are forcing the development of new methodologies for CAA applications. This paper presents numerical results from a third and a seventh order algorithm for the propagation of an acoustic signal through a converging-diverging nozzle. These two algorithms are members of a new family of high accuracy methods that have the same order accuracy in both space and time, and are an extension of previous work for linearized Euler equations to fully nonlinear time dependent problems. The simulations are all with the fully nonlinear quasi-1D Euler equations for the total solution, which includes both the fluid dynamics and the acoustics. The acoustic solution is obtained from the time dependent nonlinear solution by subtracting the steady solution.

Goodrich, John W.↗
NASA NTRS2000

Comparison with Analytical Solution: Generation and Radiation of Acoustic Waves from a 2-D Shear Layer

An acoustic source inside of a 2-D jet excites an instability wave in the shear layer resulting in sound radiating away from the shear layer. Solve the linearized Euler equations to predict the sound radiation outside of the jet. The jet static pressure is assumed to be constant. The jet flow is parallel and symmetric about the x-axis. Use a symmetry boundary condition along the x-axis.

Dahl, Milo D.↗
NASA NTRS1994

Steady Secondary Momentum and Enthalpy Streaming in the Pulse Tube Refrigerator

Our study investigates the steady secondary streaming that occurs in the pulse tube refrigerator through application of the conservation laws and support by flow visualization experiments. A series expansion of the mass, momentum and energy conservation laws for an axisymmetric system yields a set of linear differential equations. The series solution is expanded in the small parameter epsilon = 1/Str where Str is the Strouhal number. The O(1) basic flow problem is solved for the case of a compressible ideal gas oscillating within a cylindrical tube. The phase and amplitude of the boundary conditions on axial velocity are treated as independent parameters. The basic solution for the gas domain is solved for the case of temperature and heat flux coupling to the tube wall. The temperature coupled solution shows how the Fourier number of the tube wall affects the gas temperature, phase angle, and time-averaged enthalpy flow. The basic solution is also examined for the time-averaged Reynolds stresses which are shown to be the driving force that produces steady secondary momentum streaming. Calculated solutions to the steady momentum streaming are shown to be in qualitative agreement with flow visualization experiments.

Lee, J. M.↗