SEARCH · Engineering Papers

Results for “grid codes”

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 847 records · Page 47

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Analytical calculation of a single jet in cross-flow and comparison with experiment

A series of progressively finer grid systems to calculate the single jet in cross-flow experimentally measured were employed. The experimental measurements provide a collection of velocities, turbulence intensities, and jet concentration profiles with measurements of the inlet field. The use of a series of progressively finer grid systems allows a differentiation between numerical errors and the hydrodynamic modeling assumptions embodied in the 3D combustor code.

Claus, R. W.↗

NASA NTRS1984

Calculating C-grids with fine and embedded mesh regions

A program for calculating a C-type mesh around airfoil like shapes is described. The Jameson/Caughey approach is used: a parabolic transformation coupled with a shearing transformation. The algebraic algorithm is capable of efficiently generating nearly orthogonal grids. A high degree of grid control is possible. The user may specify grid boundaries, number of grid lines, and location of (and spacing in) trailing edge and leading edge fine mesh areas. The capability of embedding fine mesh regions, for use with new adaptive grid techniques, is being developed. Grids generated by the program were used in Euler flow flow calculatons by W. Usab. Results superior to results calculated on previous O-type grids were obtained. Specifically, calculations converged faster using C-grids rather than 0-grids, total pressure loss spikes at the trailing edge of the airfoil were eliminated, and the Ni method converged with zero artificial smoothing for a subcritical case (resulting overall total pressure loss was then nearly zero). These improvements were attributed to higher grid orthogonality, especially at the trailing edge. The program itself is fairly straightforward. Roughly half of the 800 code lines are comment lines.

Loyd, B.↗

NASA NTRS1989

Computation of hypersonic flow fields

The present discussion of the theory, solution methods, and development status of chemically reactive flow CFD codes pertaining to the flow field around a hypersonic vehicle encompasses the formulation of multicomponent viscous equations, grid-generation techniques, and finite-difference algorithms. The Navier-Stokes equations presented focus on such particularities of high speed vehicles as their aerodynamic configuration, lee and wake flow, bow shock, and chemistry and low density effects. Issues for which further development is desirable encompass boundary treatments, grid quality, shock capturing, and the decoupling of chemistry from flow variables.

Li, Chien-Peng↗

NASA NTRS1991

The NASA aircraft noise prediction program improved propeller analysis system

The improvements and the modifications of the NASA Aircraft Noise Prediction Program (ANOPP) and the Propeller Analysis System (PAS) are described. Comparisons of the predictions and the test data are included in the case studies for the flat plate model in the Boundary Layer Module, for the effects of applying compressibility corrections to the lift and pressure coefficients, for the use of different weight factors in the Propeller Performance Module, for the use of the improved retarded time equation solution, and for the effect of the number grids in the Transonic Propeller Noise Module. The DNW tunnel test data of a propeller at different angles of attack and the Dowty Rotol data are compared with ANOPP predictions. The effect of the number of grids on the Transonic Propeller Noise Module predictions and the comparison of ANOPP TPN and DFP-ATP codes are studied. In addition to the above impact studies, the transonic propeller noise predictions for the SR-7, the UDF front rotor, and the support of the enroute noise test program are included.

Nguyen, L. Cathy↗

NASA NTRS2002

Time-Dependent Simulations of Turbopump Flows

Unsteady flow simulations for RLV (Reusable Launch Vehicles) 2nd Generation baseline turbopump for one and half impeller rotations have been completed by using a 34.3 Million grid points model. MLP (Multi-Level Parallelism) shared memory parallelism has been implemented in INS3D, and benchmarked. Code optimization for cash based platforms will be completed by the end of September 2001. Moving boundary capability is obtained by using DCF module. Scripting capability from CAD (computer aided design) geometry to solution has been developed. Data compression is applied to reduce data size in post processing. Fluid/Structure coupling has been initiated.

Kiris, Cetin↗

NASA NTRS2011

Current Status on Radiation Modeling for the Hayabusa Re-entry

On June 13, 2010 the Japanese Hayabusa capsule performed its reentry into the Earths atmosphere over Australia after a seven year journey to the asteroid Itokawa. The reentry was studied by numerous imaging and spectroscopic instruments onboard NASA's DC-8 Airborne Laboratory and from three sites on the ground, in order to measure surface and plasma radiation generated by the Hayabusa Sample Return Capsule (SRC). Post flight, the flow solutions were recomputed to include the whole flow field around the capsule at 11 points along the reentry trajectory using updated trajectory information. Again, material response was taken into account to obtain most reliable surface temperature information. These data will be used to compute thermal radiation of the glowing heat shield and plasma radiation by the shock/post-shock layer system to support analysis of the experimental observation data. For this purpose, lines of sight data are being extracted from the flow field volume grids and plasma radiation will be computed using NEQAIR [4] which is a line-by-line spectroscopic code with one-dimensional transport of radiation intensity. The procedures being used were already successfully applied to the analysis of the observation of the Stardust reentry [5].

Winter, Michael W.↗

NASA NTRS2007

An Application of Overset Grids to Payload/Fairing Three-Dimensional Internal Flow CFD Analysis

The application of overset grids to the computational fluid dynamics analysis of three-dimensional internal flow in the payload/fairing of an expendable launch vehicle is described. In conjunction with the overset grid system, the flowfield in the payload/fairing configuration is obtained with the aid of OVERFLOW Navier-Stokes code. The solution exhibits a highly three dimensional complex flowfield with swirl, separation, and vortices. Some of the computed flow features are compared with the measured Laser-Doppler Velocimetry (LDV) data on a 1/5th scale model of the payload/fairing configuration. The counter-rotating vortex structures and the location of the saddle point predicted by the CFD analysis are in general agreement with the LDV data. Comparisons of the computed (CFD) velocity profiles on horizontal and vertical lines in the LDV measurement plane in the faring nose region show reasonable agreement with the LDV data.

Kandula, Max↗

NASA NTRS2016

An Optimal-Estimation-Based Aerosol Retrieval Algorithm Using OMI Near-UV Observations

An optimal-estimation(OE)-based aerosol retrieval algorithm using the OMI (Ozone Monitoring Instrument) near-ultraviolet observation was developed in this study. The OE-based algorithm has the merit of providing useful estimates of errors simultaneously with the inversion products. Furthermore, instead of using the traditional lookup tables for inversion, it performs online radiative transfer calculations with the VLIDORT (linearized pseudo-spherical vector discrete ordinate radiative transfer code) to eliminate interpolation errors and improve stability. The measurements and inversion products of the Distributed Regional Aerosol Gridded Observation Network campaign in northeast Asia (DRAGON NE-Asia 2012) were used to validate the retrieved aerosol optical thickness (AOT) and single scattering albedo (SSA). The retrieved AOT and SSA at 388 nm have a correlation with the Aerosol Robotic Network (AERONET) products that is comparable to or better than the correlation with the operational product during the campaign. The OEbased estimated error represented the variance of actual biases of AOT at 388 nm between the retrieval and AERONET measurements better than the operational error estimates. The forward model parameter errors were analyzed separately for both AOT and SSA retrievals. The surface reflectance at 388 nm, the imaginary part of the refractive index at 354 nm, and the number fine-mode fraction (FMF) were found to be the most important parameters affecting the retrieval accuracy of AOT, while FMF was the most important parameter for the SSA retrieval. The additional information provided with the retrievals, including the estimated error and degrees of freedom, is expected to be valuable for relevant studies. Detailed advantages of using the OE method were described and discussed in this paper.

Air Pollution↗

NASA NTRS2015

Influence of Spanwise Boundary Conditions on Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations with the OVERFLOW code in conjunction with a Ffowcs Williams-Hawkings acoustics solver. In the present study, two different spanwise grids are being used to investigate the effect of the spanwise extent and periodicity on the near-field unsteady structures and radiated noise. The baseline grid with periodic boundary conditions has a short span equal to 1/9th of the stowed chord, whereas the other, longer span grid adds stretched grids on both sides of the core, baseline grid to allow inviscid surface boundary conditions at both ends. The results indicate that the near-field mean statistics obtained using the two grids are similar to each other, as are the directivity and spectral shapes of the radiated noise. However, periodicity forces all acoustic waves with less than one wavelength across the span to be two-dimensional, without any variation in the span. The spanwise coherence of the acoustic waves is what is needed to make estimates of the noise that would be radiated from realistic span lengths. Simulations with periodic conditions need spans of at least six slat chords to allow spanwise variation in the low-frequencies associated with the peak of broadband slat noise. Even then, the full influence of the periodicity is unclear, so employing grids with a fine, central region and highly stretched meshes that go to slip walls may be a more efficient means of capturing the spanwise decorrelation of low-frequency acoustic phenomena.

Lockard, David P.↗

NASA NTRS2016

Preliminary Computational Study for Future Tests in the NASA Ames 9 foot' x 7 foot Wind Tunnel

The NASA Advanced Air Vehicles Program, Commercial Supersonics Technology Project seeks to advance tools and techniques to make over-land supersonic flight feasible. In this study, preliminary computational results are presented for future tests in the NASA Ames 9 foot x 7 foot supersonic wind tunnel to be conducted in early 2016. Shock-plume interactions and their effect on pressure signature are examined for six model geometries. Near- field pressure signatures are assessed using the CFD code USM3D to model the proposed test geometries in free-air. Additionally, results obtained using the commercial grid generation software Pointwise Reigistered Trademark are compared to results using VGRID, the NASA Langley Research Center in-house mesh generation program.

Pearl, Jason M.↗

NASA NTRS2015

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗

NASA NTRS1988

Application of advanced grid generation techniques for flow field computations about complex configurations

In the computation of flowfields about complex configurations, it is very difficult to construct a boundary-fitted coordinate system. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach; its applications are investigated. The method conservative providing conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Steady state solutions of the Euler equations are presented and discussed. The solutions include: low speed flow over a sphere, high speed flow over a slender body, supersonic flow through a duct, and supersonic internal/external flow interaction for an aircraft configuration at various angles of attack. The results demonstrate that the multiple grids approach along with the conservative interfacing is capable of computing the flows about the complex configurations where the use of a single grid system is not possible.

Kathong, Monchai↗

NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗

NASA NTRS1981

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗

NASA NTRS1990

Computational analysis of flow in 3D propulsive transition ducts

A numerical analysis of fully three dimensional, statistically steady flows in propulsive transition ducts being considered for use in future aircraft of higher maneuverability is investigated. The purpose of the transition duct is to convert axisymmetric flow from conventional propulsion systems to that of a rectangular geometry of high aspect ratio. In an optimal design, the transition duct would be of minimal length in order to reduce the weight penalty, while the geometrical change would be gradual enough to avoid detrimental flow perturbations. Recent experiments conducted at the Propulsion Aerodynamics Branch have indicated that thrust losses in ducts of superelliptic cross-section can be surprisingly low, even if flow separation occurs near the divergent walls. In order to address the objective of developing a rational design procedure for optimal transition ducts, it is necessary to have available a reliable computational tool for the analysis of flows achieved in a sequence of configurations. Current CFD efforts involving complicated geometries usually must contend with two separate but interactive aspects: namely, grid generation and flow solution. The first two avenues of the present investigation were comprised of suitable grid generation for a class of transition ducts of superelliptic cross-section, and the subsequent application of the flow solver PAB3D to this geometry. The code, PAB3D, was developed as a comprehensive tool for the solution of both internal and external high speed flows. The third avenue of investigation has involved analytical formulations to aid in the understanding of the nature of duct flows, and also to provide a basis of comparison for subsequent numerical solutions. Numerical results to date include the generation of two preliminary grid systems for duct flows, and the initial application of PAB3D to the corresponding geometries, which are of the class tested experimentally.

Sepri, Paavo↗

NASA NTRS1994

Incremental triangulation by way of edge swapping and local optimization

This document is intended to serve as an installation, usage, and basic theory guide for the two dimensional triangulation software 'HARLEY' written for the Silicon Graphics IRIS workstation. This code consists of an incremental triangulation algorithm based on point insertion and local edge swapping. Using this basic strategy, several types of triangulations can be produced depending on user selected options. For example, local edge swapping criteria can be chosen which minimizes the maximum interior angle (a MinMax triangulation) or which maximizes the minimum interior angle (a MaxMin or Delaunay triangulation). It should be noted that the MinMax triangulation is generally only locally optical (not globally optimal) in this measure. The MaxMin triangulation, however, is both locally and globally optical. In addition, Steiner triangulations can be constructed by inserting new sites at triangle circumcenters followed by edge swapping based on the MaxMin criteria. Incremental insertion of sites also provides flexibility in choosing cell refinement criteria. A dynamic heap structure has been implemented in the code so that once a refinement measure is specified (i.e., maximum aspect ratio or some measure of a solution gradient for the solution adaptive grid generation) the cell with the largest value of this measure is continually removed from the top of the heap and refined. The heap refinement strategy allows the user to specify either the number of cells desired or refine the mesh until all cell refinement measures satisfy a user specified tolerance level. Since the dynamic heap structure is constantly updated, the algorithm always refines the particular cell in the mesh with the largest refinement criteria value. The code allows the user to: triangulate a cloud of prespecified points (sites), triangulate a set of prespecified interior points constrained by prespecified boundary curve(s), Steiner triangulate the interior/exterior of prespecified boundary curve(s), refine existing triangulations based on solution error measures, and partition meshes based on the Cuthill-McKee, spectral, and coordinate bisection strategies.

Wiltberger, N. Lyn↗

NASA NTRS2004

Single Airfoil Gust Response Problem: Category 3, Problem 1

An unsteady aerodynamic code, called GUST3D (ref. 3), has been developed to solve equation (8) for flows with periodic vortical disturbances. The code uses a frequency-domain approach with second-order central differences and a pressure radiation condition in the far field. GUST3D requires as input certain mean flow quantities which are calculated separately by a potential flow solver. The solver calculates the mean ow using a Gothert's Rule approximation (ref. 3). On the airfoil surface, it uses the solution calculated by the potential code FLO36 (ref. 4). Figures 1-2 show the mean pressure along the airfoil surface for the two airfoil geometries. In Figures 3 - 8, we present the RMS pressure on the airfoil surface. Each figure shows three GUST3D solutions (calculated on grids with different far-field boundary locations). Three solutions are shown to provide some indication of the numerical uncertainty in the results. Figures 9 - 13 present the acoustic intensity. We again show three solutions per case. Note that no results are presented for the k1 = k2 = 2.0 loaded airfoil case, as an acceptable solution could not be obtained. A few comments need to be made about the results shown. First, since the last Workshop, the GUST3D code has been substantially upgraded. This includes implementing a more accurate far-field boundary condition (ref. 5) and developing improved gridding capabilities. This is the reason for any differences that may exist between the present results and results from the last Workshop. Second, the intensity results on the circle R = 4C were obtained using a Kirchoff method (ref. 6). The Kirchoff surface was the circle R = 2C. Finally, the GUST3D code is most accurate for low reduced frequencies. A new domain decomposition approach (ref. 7) has been developed to improve accuracy. Both the single domain and domain decomposition approaches were used in generating the present results.

Scott, James R.↗

NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

Context. Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims. Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods. We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results. The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

hydrodynamics↗