Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “direct solver”

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 325 records · Page 18

The Prediction of Broadband Shock-Associated Noise Including Propagation Effects

An acoustic analogy is developed based on the Euler equations for broadband shock- associated noise (BBSAN) that directly incorporates the vector Green's function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) as the mean flow. The vector Green's function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation. An adjoint vector Green's function solver is implemented to determine the vector Green's function based on a locally parallel mean flow at streamwise locations of the SRANS solution. However, the developed acoustic analogy could easily be based on any adjoint vector Green's function solver, such as one that makes no assumptions about the mean flow. The newly developed acoustic analogy can be simplified to one that uses the Green's function associated with the Helmholtz equation, which is consistent with the formulation of Morris and Miller (AIAAJ 2010). A large number of predictions are generated using three different nozzles over a wide range of fully expanded Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise labs. In addition, two models for the so-called 'fine-scale' mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include the propagation effects, especially in the upstream direction of the jet.

Miller, Steven↗

An O(Nm(sup 2)) Plane Solver for the Compressible Navier-Stokes Equations

A hierarchical multigrid algorithm for efficient steady solutions to the two-dimensional compressible Navier-Stokes equations is developed and demonstrated. The algorithm applies multigrid in two ways: a Full Approximation Scheme (FAS) for a nonlinear residual equation and a Correction Scheme (CS) for a linearized defect correction implicit equation. Multigrid analyses which include the effect of boundary conditions in one direction are used to estimate the convergence rate of the algorithm for a model convection equation. Three alternating-line- implicit algorithms are compared in terms of efficiency. The analyses indicate that full multigrid efficiency is not attained in the general case; the number of cycles to attain convergence is dependent on the mesh density for high-frequency cross-stream variations. However, the dependence is reasonably small and fast convergence is eventually attained for any given frequency with either the FAS or the CS scheme alone. The paper summarizes numerical computations for which convergence has been attained to within truncation error in a few multigrid cycles for both inviscid and viscous ow simulations on highly stretched meshes.

Thomas, J. L.↗

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗

Time-Varying Loads of Co-Axial Rotor Blade Crossings

The blade crossing event of a coaxial counter-rotating rotor is a potential source of noise and impulsive blade loads. Blade crossings occur many times during each rotor revolution. In previous research by the authors, this phenomenon was analyzed by simulating two airfoils passing each other at specified speeds and vertical separation distances, using the compressible Navier-Stokes solver OVERFLOW. The simulations explored mutual aerodynamic interactions associated with thickness, circulation, and compressibility effects. Results revealed the complex nature of the aerodynamic impulses generated by upperlower airfoil interactions. In this paper, the coaxial rotor system is simulated using two trains of airfoils, vertically offset, and traveling in opposite directions. The simulation represents multiple blade crossings in a rotor revolution by specifying horizontal distances between each airfoil in the train based on the circumferential distance between blade tips. The shed vorticity from prior crossing events will affect each pair of upperlower airfoils. The aerodynamic loads on the airfoil and flow field characteristics are computed before, at, and after each airfoil crossing. Results from the multiple-airfoil simulation show noticeable changes in the airfoil aerodynamics by introducing additional fluctuation in the aerodynamic time history.

Time-Varying Loads↗

Progress in NEXT Ion Optics Modeling

Results are presented from an ion optics simulation code applied to the NEXT ion thruster geometry. The error in the potential field solver of the code is characterized, and methods and requirements for reducing this error are given. Results from a study on electron backstreaming using the improved field solver are given and shown to compare much better to experimental results than previous studies. Results are also presented on a study of the beamlet behavior in the outer radial apertures of the NEXT thruster. The low beamlet currents in this region allow over-focusing of the beam, causing direct impingement of ions on the accelerator grid aperture wall. Different possibilities for reducing this direct impingement are analyzed, with the conclusion that, of the methods studied, decreasing the screen grid aperture diameter eliminates direct impingement most effectively.

Emhoff, Jerold W.↗

AEROELASTIC SIMULATION TOOL FOR INFLATABLE BALLUTE AEROCAPTURE

A multidisciplinary analysis tool is under development for predicting the impact of aeroelastic effects on the functionality of inflatable ballute aeroassist vehicles in both the continuum and rarefied flow regimes. High-fidelity modules for continuum and rarefied aerodynamics, structural dynamics, heat transfer, and computational grid deformation are coupled in an integrated multi-physics, multi-disciplinary computing environment. This flexible and extensible approach allows the integration of state-of-the-art, stand-alone NASA and industry leading continuum and rarefied flow solvers and structural analysis codes into a computing environment in which the modules can run concurrently with synchronized data transfer. Coupled fluid-structure continuum flow demonstrations were conducted on a clamped ballute configuration. The feasibility of implementing a DSMC flow solver in the simulation framework was demonstrated, and loosely coupled rarefied flow aeroelastic demonstrations were performed. A NASA and industry technology survey identified CFD, DSMC and structural analysis codes capable of modeling non-linear shape and material response of thin-film inflated aeroshells. The simulation technology will find direct and immediate applications with NASA and industry in ongoing aerocapture technology development programs.

Liever, P. A.↗

An efficient three-dimensional Poisson solver for SIMD high-performance-computing architectures

We present an algorithm that solves the three-dimensional Poisson equation on a cylindrical grid. The technique uses a finite-difference scheme with operator splitting. This splitting maps the banded structure of the operator matrix into a two-dimensional set of tridiagonal matrices, which are then solved in parallel. Our algorithm couples FFT techniques with the well-known ADI (Alternating Direction Implicit) method for solving Elliptic PDE's, and the implementation is extremely well suited for a massively parallel environment like the SIMD architecture of the MasPar MP-1. Due to the highly recursive nature of our problem, we believe that our method is highly efficient, as it avoids excessive interprocessor communication.

Cohl, H.↗

Application of the computational aeroacoustics method to an advanced counterrotating propfan configuration

The acoustical near field of a counterrotating propeller configuration of unducted fan (UDF) type is presently predicted by means of a numerical methodology employing the transonic flow field generated by the three-dimensional numerical flow solver designated SSTAGE. This numerical approach simulates the time-averaged, three-dimensional inviscid flow field of the UDF by means of a four-stage Runge-Kutta integration scheme which marches the equations forward in time toward an asymptotic limit. The near field pressure time history is obtained directly from the flow field solution; by converting these histories with a Fourier series representation, specific data on the acoustic signal frequency spectra can be obtained.

Kim, Jinhan↗

Convergence acceleration of an aeroelastic Navier-Stokes solver

New capabilities have been added to a Navier-Stokes solver to perform steady-state simulations more efficiently. The flow solver for solving the Navier-Stokes equations is completely rewritten with a combination of the LU-SGS (Lower-Upper factored Symmetric Gauss-Seidel) implicit method and the modified HLLE (Harten-Lax-van Leer-Einfeldt) upwind scheme. A pseudo-time marching method is used for the directly coupled structural equations to improve overall convergence rates for static aeroelastic analysis. Results are demonstrated for transonic flows over rigid and flexible wings.

Obayashi, S.↗

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

Comparisons Between NO PLIF Imaging and CFD Simulations of Mixing Flowfields for High-Speed Fuel Injectors

The current work compares experimentally and computationally obtained nitric oxide (NO) planar laser-induced fluorescence (PLIF) images of the mixing flowfields for three types of high-speed fuel injectors: a strut, a ramp, and a rectangular flush-wall. These injection devices, which exhibited promising mixing performance at lower flight Mach numbers, are currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, and improve the understanding of underlying physical processes relevant to flight Mach numbers greater than eight. In the experiments, conducted in the NASA Langley Arc-Heated Scramjet Test Facility (AHSTF), the injectors are placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number air flow at the entrance of a scramjet combustor. Helium is used as an inert substitute for hydrogen fuel. The PLIF is obtained by using a tunable laser to excite the NO, which is present in the AHSTF air as a direct result of arc-heating. Consequently, the absence of signal is an indication of pure helium (fuel). The PLIF images computed from the computational fluid dynamics (CFD) simulations are obtained by combining a fluorescence model for NO with the Reynolds-Averaged Simulation results carried out using the VULCAN-CFD solver to obtain a computational equivalent of the experimentally measured PLIF signal. The measured NO PLIF signal is mainly a function of NO concentration allowing for semi-quantitative comparisons between the CFD and the experiments. The PLIF signal intensity is also sensitive to pressure and temperature variations in the flow, allowing additional flow features to be identified and compared with the CFD. Good agreement between the PLIF and the CFD results provides increased confidence in the CFD simulations for investigations of injector performance.

Drozda, Tomasz, G.↗

A Study of Coaxial Rotor Performance and Flow Field Characteristics

As a precursor to studying the acoustics of a coaxial rotor system, the aerodynamics and flow field of a coaxial rotor were simulated to better understand the interaction between the two rotors. RotUNS, an unsteady Navier-Stokes solver that uses a simplified blade aerodynamics model, was used to predict coaxial rotor performance in hover and forward flight. RotUNS steady hover calculations showed improved performance compared to blade element momentum theory. Prior to examining the complex 3D flow field of a coaxial rotor in detail, two airfoils traveling in opposite directions with a vertical separation distance equivalent to the separation between the upper and lower rotor of the coaxial system were simulated. The pressure field generated by the two airfoils aided our interpretation of the more complex coaxial rotor system flow field. The pressure fields above, between, and below the coaxial rotor system were then examined for different azimuth positions of the upper and lower rotor blades.

Barbely, Natasha L.↗

Development of Unsteady Aerodynamic State-Space Models from CFD-Based Pulse Responses

A method for computing discrete-time state-space models of linearized unsteady aerodynamic behavior directly from aeroelastic CFD codes is presented. The method involves the treatment of CFD-based pulse responses as Markov parameters for use in a system identification /realization algorithm. Results are presented for the AGARD 445.6 Aeroelastic Wing with four aeroelastic modes at a Mach number of 0.96 using the EZNSS Euler/Navier-Stokes flow solver with aeroelastic capability. The System/Observer/Controller Identification Toolbox (SOCIT) algorithm, based on the Ho-Kalman realization algorithm, is used to generate 15th- and 32nd-order discrete-time state-space models of the unsteady aerodynamic response of the wing over the entire frequency range of interest.

Silva, Walter A.↗

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

Tandem Cylinder Noise Predictions

In an effort to better understand landing-gear noise sources, we have been examining a simplified configuration that still maintains some of the salient features of landing-gear flow fields. In particular, tandem cylinders have been studied because they model a variety of component level interactions. The present effort is directed at the case of two identical cylinders spatially separated in the streamwise direction by 3.7 diameters. Experimental measurements from the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF) at NASA Langley Research Center (LaRC) have provided steady surface pressures, detailed off-surface measurements of the flow field using Particle Image Velocimetry (PIV), hot-wire measurements in the wake of the rear cylinder, unsteady surface pressure data, and the radiated noise. The experiments were conducted at a Reynolds number of 166 105 based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent shedding process and simulate the effects of a high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The current calculations further explore the influence of the grid resolution and spanwise extent on the flow and associated radiated noise. Extensive comparisons with the experimental data are used to assess the ability of the computations to simulate the details of the flow. The results show that the pressure fluctuations on the upstream cylinder, caused by vortex shedding, are smaller than those generated on the downstream cylinder by wake interaction. Consequently, the downstream cylinder dominates the noise radiation, producing an overall directivity pattern that is similar to that of an isolated cylinder. Only calculations based on the full length of the model span were able to capture the complete decay in the spanwise correlation, thereby producing reasonable noise radiation levels.

Lockard, David P.↗

A perspective on unstructured grid flow solvers

This survey paper assesses the status of compressible Euler and Navier-Stokes solvers on unstructured grids. Different spatial and temporal discretization options for steady and unsteady flows are discussed. The integration of these components into an overall framework to solve practical problems is addressed. Issues such as grid adaptation, higher order methods, hybrid discretizations and parallel computing are briefly discussed. Finally, some outstanding issues and future research directions are presented.

Venkatakrishnan, V.↗