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 415 records · Page 23

Boeing 18-Inch Fan Rig Broadband Noise Test

The purposes of the subject test were to identify and quantify the mechanisms by which fan broadband noise is produced, and to assess the validity of such theoretical models of those mechanisms as may be available. The test was conducted with the Boeing 18-inch fan rig in the Boeing Low-Speed Aeroacoustic Facility (LSAF). The rig was designed to be particularly clean and geometrically simple to facilitate theoretical modeling and to minimize sources of interfering noise. The inlet is cylindrical and is equipped with a boundary layer suction system. The fan is typical of modern high-by-pass ratio designs but is capable of operating with or without fan exit guide vanes (stators), and there is only a single flow stream. Fan loading and tip clearance are adjustable. Instrumentation included measurements of fan performance, the unsteady flow field incident on the fan and stators, and far-field and in-duct acoustic fields. The acoustic results were manipulated to estimate the noise generated by different sources. Significant fan broadband noise was found to come from the rotor self-noise as measured with clean inflow and no boundary layer. The rotor tip clearance affected rotor self-noise somewhat. The interaction of the rotor with inlet boundary layer turbulence is also a significant source, and is strongly affected by rotor tip clearance. High level noise can be generated by a high-order nonuniform rotating at a fraction of the fan speed, at least when tip clearance and loading are both large. Stator-generated noise is the loudest of the significant sources, by a small margin, at least on this rig. Stator noise is significantly affected by propagation through the fan.

Ganz, Ulrich W.↗

Evaluation of Boundary Conditions for the Gust-Cascade Problem

Using a high-order accuracy finite-difference time-domain algorithm, the acoustic scattering from a flat-plate cascade is computed. Keeping the grid and time step fixed, the effect of four different boundary conditions on the accuracy and stability of the computed solution is compared.

Hixon, R.↗

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.↗

Parallelization of an Object-Oriented Unstructured Aeroacoustics Solver

A computational aeroacoustics code based on the discontinuous Galerkin method is ported to several parallel platforms using MPI. The discontinuous Galerkin method is a compact high-order method that retains its accuracy and robustness on non-smooth unstructured meshes. In its semi-discrete form, the discontinuous Galerkin method can be combined with explicit time marching methods making it well suited to time accurate computations. The compact nature of the discontinuous Galerkin method also makes it well suited for distributed memory parallel platforms. The original serial code was written using an object-oriented approach and was previously optimized for cache-based machines. The port to parallel platforms was achieved simply by treating partition boundaries as a type of boundary condition. Code modifications were minimal because boundary conditions were abstractions in the original program. Scalability results are presented for the SCI Origin, IBM SP2, and clusters of SGI and Sun workstations. Slightly superlinear speedup is achieved on a fixed-size problem on the Origin, due to cache effects.

Baggag, Abdelkader↗

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗

Implementation of Nonlinear Control Laws for an Optical Delay Line

This paper discusses the implementation of a globally stable nonlinear controller algorithm for the Real-Time Interferometer Control System Testbed (RICST) brassboard optical delay line (ODL) developed for the Interferometry Technology Program at the Jet Propulsion Laboratory. The control methodology essentially employs loop shaping to implement linear control laws. while utilizing nonlinear elements as means of ameliorating the effects of actuator saturation in its coarse, main, and vernier stages. The linear controllers were implemented as high-order digital filters and were designed using Bode integral techniques to determine the loop shape. The nonlinear techniques encompass the areas of exact linearization, anti-windup control, nonlinear rate limiting and modal control. Details of the design procedure are given as well as data from the actual mechanism.

Hench, John J.↗

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.↗

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗

Selected contribution: a three-dimensional model for assessment of in vitro toxicity in balaena mysticetus renal tissue

This study established two- and three-dimensional renal proximal tubular cell cultures of the endangered species bowhead whale (Balaena mysticetus), developed SV40-transfected cultures, and cloned the 61-amino acid open reading frame for the metallothionein protein, the primary binding site for heavy metal contamination in mammals. Microgravity research, modulations in mechanical culture conditions (modeled microgravity), and shear stress have spawned innovative approaches to understanding the dynamics of cellular interactions, gene expression, and differentiation in several cellular systems. These investigations have led to the creation of ex vivo tissue models capable of serving as physiological research analogs for three-dimensional cellular interactions. These models are enabling studies in immune function, tissue modeling for basic research, and neoplasia. Three-dimensional cellular models emulate aspects of in vivo cellular architecture and physiology and may facilitate environmental toxicological studies aimed at elucidating biological functions and responses at the cellular level. Marine mammals occupy a significant ecological niche (72% of the Earth's surface is water) in terms of the potential for information on bioaccumulation and transport of terrestrial and marine environmental toxins in high-order vertebrates. Few ex vivo models of marine mammal physiology exist in vitro to accomplish the aforementioned studies. Techniques developed in this investigation, based on previous tissue modeling successes, may serve to facilitate similar research in other marine mammals.

Whales/genetics/physiology↗

622-Mbps Orthogonal Frequency Division Multiplexing Modulator Developed

The Communications Technology Division at the NASA Lewis Research Center is developing advanced electronic technologies for the space communications and remote sensing systems of tomorrow. As part of the continuing effort to advance the state-of-the art in satellite communications and remote sensing systems, Lewis is developing a programmable Orthogonal Frequency Division Multiplexing (OFDM) modulator card for high-data-rate communication links. The OFDM modulator is particularly suited to high data-rate downlinks to ground terminals or direct data downlinks from near-Earth science platforms. It can support data rates up to 622 megabits per second (Mbps) and high-order modulation schemes such as 16-ary quadrature amplitude modulation (16-ary QAM) or 8- phase shift keying (8PSK). High order modulations can obtain the bandwidth efficiency over the traditional binary phase shift keying (BPSK) or quadrature phase shift keying (QPSK) modulator schemes. The OFDM modulator architecture can also be precompensated for channel disturbances and alleviate amplitude degradations caused by nonlinear transponder characteristics.

Nguyen, Na T.↗

Statistical Mechanics and Dynamics of the Outer Solar System.I. The Jupiter/Saturn Zone

We report on numerical simulations designed to understand how the solar system evolved through a winnowing of planetesimals accreeted from the early solar nebula. This sorting process is driven by the energy and angular momentum and continues to the present day. We reconsider the existence and importance of stable niches in the Jupiter/Saturn Zone using greatly improved numerical techniques based on high-order optimized multi-step integration schemes coupled to roundoff error minimizing methods.

numerical simulations planetesimals planet formati↗

Direct Numerical Simulations of Boundary Layer Transition on a Flat Plate

In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Rai, Man Mohan↗

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗