Solution-adaptive grid procedure for the parabolized Navier-Stokes equations
Explore the source record for details and available documents.
Harvey, Albert D., III↗
SEARCH · Search NASA
Results for “grid solutions”
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.
Application of a solution adaptive grid scheme to complex three-dimensional flows
Explore the source record for details and available documents.
Davies, Carol B.↗
Hyperbolic Prismatic Grid Generation and Solution of Euler Equations on Prismatic Grids
A hyperbolic grid generation method is used to generate prismatic grids and an approach using prismatic grids to solve the Euler equations is presented. The theory of the stability and feasibility of the hyperbolic grid generation method is presented. The hyperbolic grid generation method of Steger et al for structured grids is applied to a three dimensional triangularized surface definition to generate a grid that is unstructured on each successive layer. The grid, however, retains structure in the body-normal direction and has a computational cell shaped like a triangular prism. In order to take advantage of the structure in the normal direction, a finite-volume scheme that treats the unknowns along the normal direction implicitly is introduced and the flow over a sphere is simulated.
Pandya, S. A.↗
GRAPEVINE: Grids about anything by Poisson's equation in a visually interactive networking environment
A proven 3-D multiple-block elliptic grid generator, designed to run in 'batch mode' on a supercomputer, is improved by the creation of a modern graphical user interface (GUI) running on a workstation. The two parts are connected in real time by a network. The resultant system offers a significant speedup in the process of preparing and formatting input data and the ability to watch the grid solution converge by replotting the grid at each iteration step. The result is a reduction in user time and CPU time required to generate the grid and an enhanced understanding of the elliptic solution process. This software system, called GRAPEVINE, is described, and certain observations are made concerning the creation of such software.
Sorenson, Reese L.↗
Achieving geometric accuracy in FFT-based micromechanical models using conformal grid
Owing to its efficiency, simplicity and robustness, the FFT-based method has become the standard for computation of mechanical fields in a heterogeneous periodic unit cell. One of the main disadvantages of the FFT-based method is the inaccurate representation of the initial microstructure on a regular grid of voxels, which can be alleviated through the use of distorted initial grids. Here, in this paper, a method for generation of distorted initial grids conforming to the microstructural features (e.g. straight/curved boundaries) is proposed. The method determines the positions of the grid nodes in the initial configuration by solving a system of springs connecting the nodes. Microstructures consisting of layers, Voronoi tessellation and circular/spherical inclusions are considered, and mechanical fields simulated using the FFT-based method. It is found that distorted initial grids, conforming to the microstructural features, lead to more accurate mechanical fields in comparison to the corresponding non-distorted initial grid solution. The effect of initial grid distortion on the convergence of the FFT-based method is analyzed and discussed.
36 MATERIALS SCIENCE↗
A truncation error injection approach to viscous-inviscid interaction
An approach to viscous-inviscid interaction which is based on truncation error injection is presented in the context of solving flow over an airfoil. A two-dimensional interpolation scheme is used to restrict the fine grid solutions to the global coarse grid. Details on the current implementation of the approach are given, and the boundary conditions being used are discussed. Inviscid results from a NACA0012 airfoil test case and the viscous results are presented.
Goble, B. D.↗
An adaptive grid/Navier-Stokes methodology for the calculation of nozzle afterbody base flows with a supersonic freestream
This paper describes an adaptive grid method for base flows in a supersonic freestream. The method is based on the direct finite-difference statement of the equidistribution principle. The weighting factor is a combination of the Mach number, density, and velocity first-derivative gradients in the radial direction. Two key ideas of the method are to smooth the weighting factor by using a type of implicit smoothing and to allow boundary points to move in the grid adaptation process. An AGARD nozzle afterbody base flow configuration is used to demonstrate the performance of the adaptive grid methodology. Computed base pressures are compared to experimental data. The adapted grid solutions offer a dramatic improvement in base pressure prediction compared to solutions computed on a nonadapted grid. A total-variation-diminishing (TVD) Navier-Stokes scheme is used to solve the governing flow equations.
Williams, Morgan↗
Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids
We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.
Jeffery A. White↗
Unsteady airfoil flow solutions on moving zonal grids
Euler and Navier-Stokes solutions for airfoil flows on zonal grids are presented. The governing equations are solved with an implicit, iterative, factorized numerical scheme. The inviscid fluxes are examined with a third-order accurate upwind method. Zonal grid solutions are compared with experimental measurements for flows over airfoils at fixed angles of incidence. The computed unsteady solutions for rapidly pitching and oscillating airfoils are in good agreement with experiments.
Cricelli, Antonio S.↗
Computation of supersonic jet noise under imperfectly expanded conditions
The turbulent mixing noise of supersonic jet under imperfectly expanded conditions is calculated for convergent and convergent-divergent (CD) axisymmetric nozzle geometries. The noise prediction incorporates CFD solution of Navier-Stokes equations. The effect of grid resolution on shock structure computation is demonstrated. Mixing noise spectra predicted from fine and coarse grid solutions exhibit little sensitivity to the grid resolution. A proper grid resolution, however, results in a significant improvement in shock capturing capability and helps predictions agree favorably with experimental data. Good agreement between predicted noise spectra and data shows that the CFD-incorporated noise prediction scheme, which was demonstrated for shock-free conditions, works as well for shock-containing flow conditions.
Kim, Chan M.↗
Computation of supersonic jet noise under imperfectly expanded conditions
The turbulent mixing noise of supersonic jet under imperfectly expanded conditions is calculated for convergent and convergent-divergent (CD) axisymmetric nozzle geometries. The noise prediction incorporates CFD solution of Navier-Stokes equations. The effect of grid resolution on shock structure computation is demonstrated. Mixing noise spectra predicted from fine and coarse grid solutions exhibit little sensitivity to the grid resolution. A proper grid resolution, however, results in a significant improvement in shock capturing capability and helps predictions agree favorably with experimental data. Good agreement between predicted noise spectra and data shows that the CFD-incorporated noise prediction scheme, which was demonstrated for shock-free conditions, works as well for shock-containing flow conditions.
Kim, Chan M.↗
CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5
Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.
Lemieux, Jean-François (ORCID:0000000320845759)↗
Solution of heat conduction problems by the grid method
Solution of heat conduction problems by grid method
CONDUCTIVE HEAT TRANSFER↗
Adapting Grids For Computing Two-Dimensional Flows
SAGE2D is two-dimensional implementation of Self Adaptive Grid Evolution computer program that intelligently redistributes initial grid points on basis of initial flow-field solution. Grids modified according to initial computed flows enabling recomputation at greater accuracy. Written in FORTRAN 77.
Davies, Carol B.↗
ICASE/LaRC Workshop on Adaptive Grid Methods
Solution-adaptive grid techniques are essential to the attainment of practical, user friendly, computational fluid dynamics (CFD) applications. In this three-day workshop, experts gathered together to describe state-of-the-art methods in solution-adaptive grid refinement, analysis, and implementation; to assess the current practice; and to discuss future needs and directions for research. This was accomplished through a series of invited and contributed papers. The workshop focused on a set of two-dimensional test cases designed by the organizers to aid in assessing the current state of development of adaptive grid technology. In addition, a panel of experts from universities, industry, and government research laboratories discussed their views of needs and future directions in this field.
South, Jerry C., Jr.↗
Grid generation and flow solution method for Euler equations on unstructured grids
A grid generation and flow solution algorithm for the Euler equations on unstructured grids is presented. The grid generation scheme, which uses Delaunay triangulation, generates the field points for the mesh based on cell aspect ratios and allows clustering of grid points near solid surfaces. The flow solution method is an implicit algorithm in which the linear set of equations arising at each time step is solved using a Gauss-Seidel procedure that is completely vectorizable. Also, a study is conducted to examine the number of subiterations required for good convergence of the overall algorithm. Grid generation results are shown in two dimensions for an NACA 0012 airfoil as well as a two element configuration. Flow solution results are shown for a two dimensional flow over the NACA 0012 airfoil and for a two element configuration in which the solution was obtained through an adaptation procedure and compared with an exact solution. Preliminary three dimensional results also are shown in which the subsonic flow over a business jet is computed.
Anderson, W. Kyle↗
Nonreflecting boundary conditions for the complete unsteady transonic small-disturbance equation
Nonreflecting far-field boundary conditions that are consistent with the complete transonic small-disturbance (TSD) equations are derived. They are implemented in a new code for solving the complete TSD equation and are tested for a harmonically oscillating NACA 64A010 airfoil in transonic flow and for a flat-plate airfoil with a pulse in the angle of attack. Using the new boundary conditions on a relatively small grid, solutions for the airfoil that are obtained that agree with large-grid calculations, resulting in a 44 percent savings in computer time. Frequency responses for the flat plate show that most of the disturbances incident on the computational boundaries are absorbed by the boundary conditions.
Whitlow, W., Jr.↗
Three dimensional unstructured multigrid for the Euler equations
The three-dimensional Euler equations are solved on unstructured tetrahedral meshes using a multigrid strategy. The driving algorithm consists of an explicit vertex-based finite-element scheme, which employs an edge-based data-structure to assemble the residuals. The multigrid approach employs a sequence of independently generated coarse and fine meshes to accelerate the convergence to steady-state of the fine grid solution. Variables, residuals and corrections are passed back and forth between the various grids of the sequence using linear interpolation. The addresses and weights for interpolation are determined in a preprocessing stage using an efficient graph traversal algorithm. The preprocessing operation is shown to require a negligible fraction of the CPU time required by the overall solution procedure, while gains in overall solution efficiencies greater than an order of magnitude are demonstrated on meshes containing up to 350,000 vertices. Solutions using globally regenerated fine meshes as well as adaptively refined meshes are given.
Mavriplis, D. J.↗