Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-order”

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

Correlations in cosmic density fields

A method is proposed to place constraints on the functional form of the high-order correlation functions zeta(sub n) that arise in cosmic density fields at large scales. This technique is based on a mass-in-cell statistic and a difference of mass in partitions of a cell. The relationship between these measures is sensitive to the formal structure of the zeta(sub n) as well as their amplitudes. This relationship is quantified in several theoretical models of structure, based on the hierarchical clustering paradigm. The results lead to a test for specific types of hierarchical clustering that is sensitive to correlations of all orders. The method is applied to examples of simulated large-scaled structure dominated by cold dark matter. In the preliminary study, the hierarchical paradigm appears to be a realistic approximation over a broad range of the scales. Furthermore, there is evidence that graphs of low-order vertices are dominant. On the basis of simulated data a phenomological model is specified that gives a good representation of clustering from linear scales to the strongly clustered regime (zeta(sub 2) approximately 500).

Bromley, B. C.↗

Structure of supersonic jet flow and its radiated sound

The present paper explores the use of large-eddy simulations as a tool for predicting noise from first principles. A high-order numerical scheme is used to perform large-eddy simulations of a supersonic jet flow with emphasis on capturing the time-dependent flow structure representating the sound source. The wavelike nature of this structure under random inflow disturbances is demonstrated. This wavelike structure is then enhanced by taking the inflow disturbances to be purely harmonic. Application of Lighthill's theory to calculate the far-field noise, with the sound source obtained from the calculated time-dependent near field, is demonstrated. Alternative approaches to coupling the near-field sound source to the far-field sound are discussed.

Mankbadi, Reda R.↗

Reduction of blade-vortex interaction noise through porous leading edge

The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Lee, Soogab↗

Head-on parallel blade-vortex interaction

An experimental and computational study was carried out to investigate the parallel head-on blade-vortex interaction (BVI) and its noise generation mechanism. A shock tube, with an enlarged test section, was used to generate a compressible starting vortex which interacted with a target airfoil. The dual-pulsed holographic interferometry (DPHI) technique and airfoil surface pressure measurements were employed to obtain quantitative flow data during the BVI. A thin-layer Navier-Stokes code (BV12D), with a high-order upwind-biased scheme and a multizonal grid, was also used to simulate numerically the phenomena occurring in the head-on BVI. The detailed structure of a convecting vortex was studied through independent measurements of density and pressure distributions across the vortex center. Results indicate that, in a strong head-on BVI, the opposite pressure peaks are generated on both sides of the leading edge as the vortex approaches. Then, as soon as the vortex passes by the leading edge, the high-pressure peak suddenly moves toward the low-peak-reducing in magnitude as it moves--simultaneously giving rise to the initial sound wave. In both experiment and computation, it is shown that the viscous effect plays a significant role in head-on BVIs.

Lee, Soogab↗

Algorithmic trends in computational fluid dynamics; The Institute for Computer Applications in Science and Engineering (ICASE)/LaRC Workshop, NASA Langley Research Center, Hampton, VA, US, Sep. 15-17, 1991

The purpose here is to assess the state of the art in the areas of numerical analysis that are particularly relevant to computational fluid dynamics (CFD), to identify promising new developments in various areas of numerical analysis that will impact CFD, and to establish a long-term perspective focusing on opportunities and needs. Overviews are given of discretization schemes, computational fluid dynamics, algorithmic trends in CFD for aerospace flow field calculations, simulation of compressible viscous flow, and massively parallel computation. Also discussed are accerelation methods, spectral and high-order methods, multi-resolution and subcell resolution schemes, and inherently multidimensional schemes.

Hussaini, M. Y.↗

Universal Formulation For Symmetries In Computed Flows

Universal formulation for high-order symmetries in boundary conditions on flows devised. Eliminates need for special procedures to incorporate symmetries and corresponding boundary conditions into computer codes solving Navier-Stokes and Euler equations of flow.

Pao, S. Paul↗

System identification methods for aircraft flight control development and validation

System-identification methods compose a mathematical model, or series of models, from measurements of inputs and outputs of dynamic systems. The extracted models allow the characterization of the response of the overall aircraft or component subsystem behavior, such as actuators and on-board signal processing algorithms. This paper discusses the use of frequency-domain system-identification methods for the development and integration of aircraft flight-control systems. The extraction and analysis of models of varying complexity from nonparametric frequency-responses to transfer-functions and high-order state-space representations is illustrated using the Comprehensive Identification from FrEquency Responses (CIFER) system-identification facility. Results are presented for test data of numerous flight and simulation programs at the Ames Research Center including rotorcraft, fixed-wing aircraft, advanced short takeoff and vertical landing (ASTOVL), vertical/short takeoff and landing (V/STOL), tiltrotor aircraft, and rotor experiments in the wind tunnel. Excellent system characterization and dynamic response prediction is achieved for this wide class of systems. Examples illustrate the role of system-identification technology in providing an integrated flow of dynamic response data around the entire life-cycle of aircraft development from initial specifications, through simulation and bench testing, and into flight-test optimization.

Tischler, Mark B.↗

A High Order Finite Difference Scheme with Sharp Shock Resolution for the Euler Equations

We derive a high-order finite difference scheme for the Euler equations that satisfies a semi-discrete energy estimate, and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semi-discrete energy estimate is based on a symmetrization of the Euler equations that preserves the homogeneity of the flux vector, a canonical splitting of the flux derivative vector, and the use of difference operators that satisfy a discrete analogue to the integration by parts procedure used in the continuous energy estimate. Around discontinuities or sharp gradients, refined grids are created on which the discrete equations are solved after adding a newly constructed artificial viscosity. The positioning of the sub-grids and computation of the viscosity are aided by a detection algorithm which is based on a multi-scale wavelet analysis of the pressure grid function. The wavelet theory provides easy to implement mathematical criteria to detect discontinuities, sharp gradients and spurious oscillations quickly and efficiently.

Gerritsen, Margot↗

SCISEAL: A CFD Code for Analysis of Fluid Dynamic Forces in Seals

A 3D CFD code, SCISEAL, has been developed and validated. Its capabilities include cylindrical seals, and it is employed on labyrinth seals, rim seals, and disc cavities. State-of-the-art numerical methods include colocated grids, high-order differencing, and turbulence models which account for wall roughness. SCISEAL computes efficient solutions for complicated flow geometries and seal-specific capabilities (rotor loads, torques, etc.).

Althavale, Mahesh M.↗

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.↗

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

Exhausted Plume Flow Field Prediction Near the Afterbody of Hypersonic Flight Vehicles in High Altitudes

A two-dimensional computer code to solve the Burnett equations has been developed which computes the flow interaction between an exhausted plume and hypersonic external flow near the afterbody of a flight vehicle. This Burnett-2D code extends the capability of Navier-Stokes solver (RPLUS2D code) to include high-order Burnett source terms and slip-wall conditions for velocity and temperature. Higher-order Burnett viscous stress and heat flux terms are discretized using central-differencing and treated as source terms. Blocking logic is adopted in order to overcome the difficulty of grid generation. The computation of exhaust plume flow field is divided into two steps. In the first step, the thruster nozzle exit conditions are computed which generates inflow conditions in the base area near the afterbody. Results demonstrated that at high altitudes, the computations of nozzle exit conditions must include the effects of base flow since significant expansion exists in the base region. In the second step, Burnett equations were solved for exhaust plume flow field near the afterbody. The free stream conditions are set at an altitude equal to 80km and the Mach number is equal to 5.0. The preliminary results show that the plume expansion, as altitude increases, will eventually cause upstream flow separation.

Chou, Lynn Chen↗

The Quest for Pi

This article gives a brief history of the analysis and computation of the mathematical constant Pi=3.14159 ..., including a number of the formulas that have been used to compute Pi through the ages. Recent developments in this area are then discussed in some detail, including the recent computation of Pi to over six billion decimal digits using high-order convergent algorithms, and a newly discovered scheme that permits arbitrary individual hexadecimal digits of Pi to be computed.

Bailey, David H.↗

Closed-Loop System Identification Experience for Flight Control Law and Flying Qualities Evaluation of a High Performance Fighter Aircraft

This paper highlights some of the results and issues associated with estimating models to evaluate control law design methods and design criteria for advanced high performance aircraft. Experimental fighter aircraft such as the NASA-High Alpha Research Vehicle (HARV) have the capability to maneuver at very high angles of attack where nonlinear aerodynamics often predominate. HARV is an experimental F/A-18, configured with thrust vectoring and conformal actuated nose strakes. Identifying closed-loop models for this type of aircraft can be made difficult by nonlinearities and high order characteristics of the system. In this paper, only lateral-directional axes are considered since the lateral-directional control law was specifically designed to produce classical airplane responses normally expected with low-order, rigid-body systems. Evaluation of the control design methodology was made using low-order equivalent systems determined from flight and simulation. This allowed comparison of the closed-loop rigid-body dynamics achieved in flight with that designed in simulation. In flight, the On Board Excitation System was used to apply optimal inputs to lateral stick and pedals at five angles at attack : 5, 20, 30, 45, and 60 degrees. Data analysis and closed-loop model identification were done using frequency domain maximum likelihood. The structure of identified models was a linear state-space model reflecting classical 4th-order airplane dynamics. Input time delays associated with the high-order controller and aircraft system were accounted for in data preprocessing. A comparison of flight estimated models with small perturbation linear design models highlighted nonlinearities in the system and indicated that the closed-loop rigid-body dynamics were sensitive to input amplitudes at 20 and 30 degrees angle of attack.

Murphy, Patrick C.↗

Coupled thermal/chemical/mechanical modeling of energetic materials in ALE3D

We must improve our ability to model the response of energetic materials to thermal stimuli and the processes involved in the energetic response. We have developed and used a time step option to efficiently and accurately compute the hours that the energetic material can take to react. Since on these longer film scales, materials can be expected to have significant motion, it is even more important to provide high-order advection for all components, including the chemical species. We show an example cook-off problem to illustrate these capabilities.

WEAPONS↗

Towards Distributed Fluids/Controls Simulations

Computational prototyping of aircraft in nonlinear regimes will require the coupling of high-order disciplines on high performance computers. This effort investigates the coupling of the Reynolds-averaged Navier-Stokes equations with six degree-of-freedom nonlinear body dynamics and controls on distributed architectures. Towards this goal, solutions obtained on a wide-node IBM SP-2, an Ethernet-linked cluster of workstations, and a Cray C90 are compared in terms of accuracy and throughput.

Atwood, Christopher A.↗

Chaotic Rotation of Nereid

The shape and spin of Neptune's outermost satellite Nereid are still unknown. Ground-based photometry indicates large brightness variations, but different observers report very different lightcurve amplitudes and periods. On the contrary, Voyager 2 images spanning 12 days show no evidence of variations greater than 0.1 mag. The latter suggest either that Nereid is nearly spherical, or that it is rotating slowly. We propose that tides have already despun Nereid's rotation to a period of a few weeks, during the time before the capture of Triton when Nereid was closer to Neptune. Since Nereid reached its present orbit, tides have further despun Nereid to a period on the order of a month. For Nereid's orbital eccentricity of 0.75, tidal evolution ceases when the spin period is still approximately 1/8 of the orbital period. Furthermore, the synchronous resonance becomes quite weak for such high eccentricities, along with other low-order spin orbit commensurabilities. In contrast, high-order resonances become very strong particularly the 6:1, 6.5:1, 7:1, 7.5:1, and 8:1 spin states. If Nereid departs by more than approximately 1% from a sphere, however, these resonances overlap, generating chaos. Our simulations show that Nereid is likely to be in chaotic rotation for any spin period longer than about 2 weeks.

Dobrovolskis, Anthony R.↗

State Estimation of Main Rotor Flap and Lead-Lag Using Accelerometers and Laser Transducers on the RASCAL UH-60 Helicopter

Modern rotorcraft flight control system designs which promise to yield high vehicle response bandwidth and good gust rejection can benefit from the use of rotor-state feedbacks. The measurement of main rotor blade motions is also desirable to validate and improve rotorcraft simulation models, to identify high-order linear flight dynamics models, to provide rotor system health monitoring; during flight test, and to provide for correlation with acoustic measurements from wind tunnel and flight tests. However, few attempts have been made to instrument a flight vehicle in this manner, and no previous system has had the robustness and accuracy required for these diverse applications. A rotor blade motion measurement and estimation system has been developed by NASA and the U.S. Army for use on the Rotorcraft Aircrew Systems Concepts Airborne Laboratory (RASCAL) helicopter. RASCAL is a UH-60 Blackhawk which is being modified at Ames Research Center in a phased development program for use in flight dynamics and controls, navigation, airspace management, and rotorcraft human factors research. The aircraft will feature a full-authority, digital, fly-by-wire research flight control system; a coupled ring laser gyro, differential GPS based navigation system; a stereoscopic color wide field of view helmet, mounted display; programmable panel mounted displays; and advanced navigation sensors. The rotor blade motion system is currently installed for data acquisition only, but will be integrated with the research flight control system when it is installed later this year.

Fletcher, Jay W.↗