SEARCH · Engineering Papers

Results for “spatial discretization”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Revisiting and Extending Interface Penalties for Multi-Domain Summation-by-Parts Operators

General interface coupling conditions are presented for multi-domain collocation methods, which satisfy the summation-by-parts (SBP) spatial discretization convention. The combined interior/interface operators are proven to be L2 stable, pointwise stable, and conservative, while maintaining the underlying accuracy of the interior SBP operator. The new interface conditions resemble (and were motivated by) those used in the discontinuous Galerkin finite element community, and maintain many of the same properties. Extensive validation studies are presented using two classes of high-order SBP operators: 1) central finite difference, and 2) Legendre spectral collocation.

Carpenter, Mark H.↗

NASA NTRS2008

Grid Quality and Resolution Issues from the Drag Prediction Workshop Series

The drag prediction workshop series (DPW), held over the last six years, and sponsored by the AIAA Applied Aerodynamics Committee, has been extremely useful in providing an assessment of the state-of-the-art in computationally based aerodynamic drag prediction. An emerging consensus from the three workshop series has been the identification of spatial discretization errors as a dominant error source in absolute as well as incremental drag prediction. This paper provides an overview of the collective experience from the workshop series regarding the effect of grid-related issues on overall drag prediction accuracy. Examples based on workshop results are used to illustrate the effect of grid resolution and grid quality on drag prediction, and grid convergence behavior is examined in detail. For fully attached flows, various accurate and successful workshop results are demonstrated, while anomalous behavior is identified for a number of cases involving substantial regions of separated flow. Based on collective workshop experiences, recommendations for improvements in mesh generation technology which have the potential to impact the state-of-the-art of aerodynamic drag prediction are given.

Mavriplis, Dimitri J.↗

NASA NTRS2008

Boundary Layer Receptivity due to Roughness and Freestream Sound for Supersonic Flows over Axisymmetric Cones

The receptivity of supersonic boundary layers due to the interaction of a two-dimensional acoustic wave with a three-dimensional roughness is numerically investigated over a 5-degree sharp tipped straight cone at a free stream Mach number of 3.5 and a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. The simulations are performed with different roughness heights and roughness locations. The steady computations with the roughness reveal that the perturbations induced by the roughness remain very close to the wall and do not penetrate to the outer part of the boundary layer. The effects of the roughness persist for about 10 boundary layer thicknesses downstream of each element. The unsteady simulation shows the generation of unstable three-dimensional first mode disturbances within a short distance from the roughness. However, the amplitudes of the instability waves are very small near the neutral point, on the order of 10(exp -4) compared to the free stream disturbances. The receptivity coefficients due to the roughness elements located very close to the tip are one order of magnitude smaller than those due to the roughness elements located closer to the neutral point.

Balakumar, P.↗

NASA NTRS2009

Stability of Supersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of three-dimensional supersonic boundary layers over a 7deg sharp tipped straight cone at an angle of attack of 4.2deg is numerically investigated at a free stream Mach number of 3.5 and at two high Reynolds numbers, 0.25 and 0.50x10(exp 6)/inch. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier-Stokes equations in cylindrical coordinates using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using the third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations reveal that the azimuthal wavenumbers are in the range of m approx. 25-50 for the most amplified traveling disturbances and in the range of m approx. 40-70 for the stationary disturbances. The N-Factor computations predicted that transition would occur further forward in the middle of the cone compared to the transition fronts near the windward and the leeward planes. The simulations revealed that the crossflow vortices originating from the nose region propagate towards the leeward plane. No perturbations were observed in the lower part of the cone.

Balakumar, Ponnampalam↗

NASA NTRS2009

Adiabatic Shock Capturing in Perfect Gas Hypersonic Flows

This paper considers the streamline-upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier-Stokes equations in conservation-variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock capturing operator, which is required to regularize the scheme in the presence of strong, shock-induced gradients. A standard shock capturing operator which has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator which are designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator.

Kirk, Benjamin S.↗

NASA NTRS2010

Stability of Hypersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of a three-dimensional hypersonic boundary layer over a 7deg half-angle straight cone at an angle of attack of 6deg is numerically investigated at a freestream Mach number of 6.0 and a Reynolds number of 10.4x10(exp 6)/m. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier- Stokes equations in cylindrical coordinates using a fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations produced azimuthal wavenumbers in the range of m approx. 20-50 for the most amplified traveling disturbances and in the range of m approx.30-70 for the stationary disturbances. The frequency of the unstable second-mode ranges from 400 kHz to 900 kHz along the windward ray. The N-Factor computations predicted transition would occur more forward on the sides of the cone as compared to the transition fronts near the windward and the leeward rays. The simulations also show the crossflow vortices originating from the nose region propagate towards the leeward ray. No perturbations were observed toward the windward half of the cone.

Balakumar, P.↗

NASA NTRS2011

The Detector Subsystem for the SXS Instrument on the Astro-H Observatory

The Soft X-ray Spectrometer (SXS) instrument on the Astro-H observatory is based on a 36 pixel x-ray calorimeter array cooled to 50 mK in a sophisticated spaceflight cryostat. The SXS is a true spatial-spectral instrument, where each spatially discrete pixel functions as a high-resolution spectrometer. Here we discuss the SXS detector subsystem that includes the detector array, the anticoincidence detector, the first stage amplifiers, the thermal and mechanical staging of the detector, and the cryogenic bias electronics. The design of the SXS detector subsystem has significant heritage from the Suzaku/XRS instrument but has some important modifications that increase performance margins and simplify the focal plane assembly. Notable improvements include x-ray absorbers with significantly lower heat capacity, improved load resistors, improved thermometry, and a decreased sensitivity to thermal radiation. These modifications have yielded an energy resolution of 3.5-4.0 eV FWHM at 6 keV for representative devices in the laboratory, giving considerable margin against the 7 eV instrument requirement. We expect similar performance in flight

Porter, Frederick↗

NASA NTRS2012

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

NASA NTRS2014

Grid Sensitivity Study for Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations in conjunction with a Ffowcs Williams-Hawkings acoustics solver. Many previous simulations have been performed for the configuration, and the case was introduced as a new category for the Second AIAA workshop on Benchmark problems for Airframe Noise Configurations (BANC-II). However, the cost of the simulations has restricted the study of grid resolution effects to a baseline grid and coarser meshes. In the present study, two different approaches are being used to investigate the effect of finer resolution of near-field unsteady structures. First, a standard grid refinement by a factor of two is used, and the calculations are performed by using the same CFL3D solver employed in the majority of the previous simulations. Second, the OVERFLOW code is applied to the baseline grid, but with a 5th-order upwind spatial discretization as compared with the second-order discretization used in the CFL3D simulations. In general, the fine grid CFL3D simulation and OVERFLOW calculation are in very good agreement and exhibit the lowest levels of both surface pressure fluctuations and radiated noise. Although the smaller scales resolved by these simulations increase the velocity fluctuation levels, they appear to mitigate the influence of the larger scales on the surface pressure. These new simulations are used to investigate the influence of the grid on unsteady high-lift simulations and to gain a better understanding of the physics responsible for the noise generation and radiation.

Lockard, David P.↗

NASA NTRS2016

An Optimized Multicolor Point-Implicit Solver for Unstructured Grid Applications on Graphics Processing Units

In the field of computational fluid dynamics, the Navier-Stokes equations are often solved using an unstructuredgrid approach to accommodate geometric complexity. Implicit solution methodologies for such spatial discretizations generally require frequent solution of large tightly-coupled systems of block-sparse linear equations. The multicolor point-implicit solver used in the current work typically requires a significant fraction of the overall application run time. In this work, an efficient implementation of the solver for graphics processing units is proposed. Several factors present unique challenges to achieving an efficient implementation in this environment. These include the variable amount of parallelism available in different kernel calls, indirect memory access patterns, low arithmetic intensity, and the requirement to support variable block sizes. In this work, the solver is reformulated to use standard sparse and dense Basic Linear Algebra Subprograms (BLAS) functions. However, numerical experiments show that the performance of the BLAS functions available in existing CUDA libraries is suboptimal for matrices representative of those encountered in actual simulations. Instead, optimized versions of these functions are developed. Depending on block size, the new implementations show performance gains of up to 7x over the existing CUDA library functions.

Zubair, Mohammad↗

NASA NTRS2019

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

NASA NTRSDate not supplied

TPSAS-NF1676L-10767-DND

Version 2.1ae of the NASA CFD code OVERFLOW has been used to simulate the flow past tandem cylinders. The multiblock grid system consisted of 15 overset meshes with a total grid count of 15.7 million points. The grid extended 3 cylinder diameters in the spanwise direction with periodic boundary conditions at the ends. Riemann conditions were used at the outer boundaries, and a no-slip condition on the cylinder surfaces. The 3rd order HLLC upwind spatial discretization scheme was used with a dual-time stepping algorithm. A modified version of the Delayed Detached Eddy Simulation hybrid RANS/LES turbulence model proposed by Spalart was used in a fully turbulent mode. The code was run on 53 Core 2 Duo nodes using MPI across nodes and OPENMP within each node. Simulation results will be compared with the available experimental data.

David P. Lockard↗

NASA NTRSDate not supplied

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

NASA NTRSDate not supplied

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

acoustics↗

NASA NTRSDate not supplied

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗

NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗

NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗