SEARCH · Engineering Papers

Results for “Triangular meshes”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Development of an unstructured mesh/Navier-Stokes method for aerodynamics of aircraft with ice accretions

An advanced flowfield prediction method for airfoils with leading edge ice accretions was developed and applied. The method is intended to be eventually used in an aircraft icing analysis. The flowfield is obtained by solving the Euler or Navier-Stokes equations on an unstructured triangular mesh. A new method has been developed to permit efficient mesh generation while providing high grid resolution near complicated airfoil ice accumulations. Unstructured-mesh Euler calculations are presented for clean and iced airfoils; the results are compared to similar structured-mesh calculations to demonstrate the new method's accuracy and efficiency.

Caruso, Steven C.↗

NASA NTRS1990

An implicit finite element method for high speed flows

A fast algorithm is presented for constructing continuous lines, consisting of element sides, on general unstructured two-dimensional triangular meshes. The lines must pass through each node of the mesh once and only once. The discussion focuses on the use of these lines in a relaxation method for the solution of the equation system arising from an implicit algorithm for the solution of two-dimensional Euler and Navier-Stokes equations on general unstructured grids. It is also shown that the method can be used for solving three-dimensional Navier-Stokes equations on a grid composed of both structured and unstructured regions.

Hassan, O.↗

NASA NTRS1992

Two-dimensional unstructured triangular grid generation

The capability of generating 2-D unstructured triangular meshes about arbitrary geometries is demonstrated. This work uses a distribution of boundary points and triangulates the computational domain using a Delaunay triangulation algorithm. Typically, initial cells are added based on cell aspect ratios or cell areas. A resulting mesh can then be used along with the connectivity of the cells to solve either a Euler or Navier-Stokes flow problem.

Jorgenson, Philip C. E.↗

NASA NTRS1993

An upwind multigrid algorithm for calculating flows on unstructured grids

An algorithm is described that calculates inviscid, laminar, and turbulent flows on triangular meshes with an upwind discretization. A brief description of the base solver and the multigrid implementation is given, followed by results that consist mainly of convergence rates for inviscid and viscous flows over a NACA four-digit airfoil section. The results show that multigrid does accelerate convergence when the same relaxation parameters that yield good single-grid performance are used; however, larger gains in performance can be realized by doing less work in the relaxation scheme.

Bonhaus, Daryl L.↗

NASA NTRS1993

A hybrid finite element-boundary integral for the analysis of cavity-backed antennas of arbitrary shape

This is the final report on this project which was concerned with the analysis of cavity-backed antennas and more specifically spiral antennas. The project was a continuation of a previous analysis, which employed rectangular brick elements, and was, thus, restricted to planar rectangular patch antennas. A total of five reports were submitted under this project and we expect that at least four journal papers will result from the research described in these reports. The abstracts of the four previous reports are included. The first of the reports (028918-1-T) is over 75 pages and describes the general formulation using tetrahedral elements and the computer program. Report 028918-2-T was written after the completion of the computer program and reviews the capability of the analysis and associated software for planar circular rectangular patches and for a rectangular planar spiral. Measurements were also done at the University of Michigan and at Mission Research Corp. for the purpose of validating the software. We are pleased to acknowledge a partial support from Mission Research Corp. in carrying out the work described in this report. The third report (028918-3-T) describes the formulation and partial validation (using 2D data) for patch antennas on a circular platform. The 3D validation and development of the formulation for patch antennas on circular platforms is still in progress. The fourth report (028918-4-T) is basically an invited journal paper which will appear in the 'J. Electromagnetic Waves and Applications' in early 1994. It describes the application of the finite element method in electromagnetics and is primarily based on our work here at U-M. This final report describes the culmination of our efforts in characterizing complex cavity-backed antennas on planar platforms. The report describes for the first time the analysis of non-planar spirals and non-rectangular slot antennas as well as traditional planar patch antennas. The comparisons between measurements and calculations are truly impressive. Another unique aspect of this work is the incorporation of the FFT as part of the BiCG solver by overlaying a structured triangular mesh over the unstructured mesh. The implementation of this BiCG-FFT solution algorithm is important in minimizing the CPU and storage requirements. This final report will be submitted for publication in a refereed journal.

Gong, Jian↗

NASA NTRS1995

Essentially nonoscillatory (ENO) reconstructions via extrapolation

In this paper, the algorithm for determining the stencil of a one-dimensional Essentially Nonoscillatory (ENO) reconstruction scheme on a uniform grid is reinterpreted as being based on extrapolation. This view leads to another extension of ENO reconstruction schemes to two-dimensional unstructured triangular meshes. The key idea here is to select several cells of the stencil in one step based on extrapolation rather than one cell at a time. Numerical experiments confirm that the new scheme yields sharp nonoscillatory reconstructions and that it is about five times faster than previous schemes.

Suresh, Ambady↗

NASA NTRS1995

Numerical study to assess sulfur hexafluoride as a medium for testing multielement airfoils

A methodology is described for computing viscous flows of air and sulfur hexafluoride (SF6). The basis is an existing flow solver that calculates turbulent flows in two dimensions on unstructured triangular meshes. The solver has been modified to incorporate the thermodynamic model for SF6 and used to calculate the viscous flow over two multielement airfoils that have been tested in a wind tunnel with air as the test medium. Flows of both air and SF6 at a free-stream Mach number of 0.2 and a Reynolds number of 9 x 10(exp 6) are computed for a range of angles of attack corresponding to the wind-tunnel test. The computations are used to investigate the suitability of SF6 as a test medium in wind tunnels and are a follow-on to previous computations for single-element airfoils. Surface-pressure, lift, and drag coefficients are compared with experimental data. The effects of heavy gas on the details of the flow are investigated based on computed boundary-layer and skin-friction data. In general, the predictions in SF6 vary little from those in air. Within the limitations of the computational method, the results presented are sufficiently encouraging to warrant further experiments.

Bonhaus, Daryl L.↗

NASA NTRS1995

Moving and adaptive grid methods for compressible flows

This paper describes adaptive grid methods developed specifically for compressible flow computations. The basic flow solver is a finite-volume implementation of Roe's flux difference splitting scheme or arbitrarily moving unstructured triangular meshes. The grid adaptation is performed according to geometric and flow requirements. Some results are included to illustrate the potential of the methodology.

Trepanier, Jean-Yves↗

NASA NTRS1993

An Upwind Multigrid Algorithm for Calculating Flows on Unstructured Grids

Bonhaus, Daryl L.↗

NASA NTRS2004

Computation of a Single Airfoil Gust Response and Gust-Cascade Interaction Using the CE/SE Method

The problems 1 and 2 in Category 3 are solved using the space-time conservation element and solution element (CE/SE) method. Problem 1 concerns the acoustic field generated by the interaction of a harmonic vertical gust with a single isolated airfoil. Problem 2 models rotor-stator interaction in a 2D cascade. Both problems involve complex geometries and flow physics including vortex shielding and acoustic radiation. An unstructured triangular mesh is used to solve both problems. For problem 2, the Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Numerical solution of both near and far fields of problem 1 are presented and compared with a frequency-domain solver GUST3D and a time-domain high-order Discontinuous Spectra Element Method (DSEM) solutions. For problem 2, numerical solutions on the blade surface, inlet and outlet planes are presented.

Wang, Xiao-Yen↗

NASA NTRS2004

Computation of Gust-Cascade Interaction Using the CE/SE Method

The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.

Wang, X.-Y.↗

NASA NTRS2008

Output-Adaptive Tetrahedral Cut-Cell Validation for Sonic Boom Prediction

A cut-cell approach to Computational Fluid Dynamics (CFD) that utilizes the median dual of a tetrahedral background grid is described. The discrete adjoint is also calculated, which permits adaptation based on improving the calculation of a specified output (off-body pressure signature) in supersonic inviscid flow. These predicted signatures are compared to wind tunnel measurements on and off the configuration centerline 10 body lengths below the model to validate the method for sonic boom prediction. Accurate mid-field sonic boom pressure signatures are calculated with the Euler equations without the use of hybrid grid or signature propagation methods. Highly-refined, shock-aligned anisotropic grids were produced by this method from coarse isotropic grids created without prior knowledge of shock locations. A heuristic reconstruction limiter provided stable flow and adjoint solution schemes while producing similar signatures to Barth-Jespersen and Venkatakrishnan limiters. The use of cut-cells with an output-based adaptive scheme completely automated this accurate prediction capability after a triangular mesh is generated for the cut surface. This automation drastically reduces the manual intervention required by existing methods.

Park, Michael A.↗

NASA NTRS2021

Using GPUs and the Parameterization Method for Rapid Search and Refinement of Connections between Tori in Periodically Perturbed Planar Circular Restricted 3-Body Problems

When the planar circular restricted 3-body problem is periodically perturbed, most unstable periodic orbits become invariant tori. However, 2D Poincare ́ sections no longer work to find their manifolds’ intersections; new methods are needed. In this study, we first review a method of restricting the intersection search to only certain manifold subsets. We then implement this search using Julia and OpenCL, representing the manifolds as triangular meshes and gaining a 30x speedup using GPUs. We finally show how to use manifold parametrizations to refine the ap- proximate connections found in the mesh search. We demonstrate the tools on the planar elliptic RTBP.

de la Llave, Rafael↗

NASA NTRS2022

Air Traffic Management TestBed: Weather Visualization using Map Tiles

Running realistic simulations and developing algorithms for decision support tools usually require weather data access. This document describes the weather data service capability implemented in the Air Traffic Management (ATM) TestBed. The capability provides ways for querying weather data from the weather products provided by the National Aeronautics and Space Administration’s (NASA) Sherlock ATM Data Warehouse. Depending on use cases, weather data can be queried using a Java application programming interface, web data access, and map tile image access. Map tile images are generated using triangular meshes obtained using the Matching Squares contour algorithm. An analysis determines a grid size that will produce both high image quality and fast image rendering. In addition, the visualization of weather data as map tiles on a TestBed visualization tool called Traffic Viewer is presented.

Chok Fung Lai↗

NASA NTRS1993

Dynamic mesh adaption for triangular and tetrahedral grids

The following topics are discussed: requirements for dynamic mesh adaption; linked-list data structure; edge-based data structure; adaptive-grid data structure; three types of element subdivision; mesh refinement; mesh coarsening; additional constraints for coarsening; anisotropic error indicator for edges; unstructured-grid Euler solver; inviscid 3-D wing; and mesh quality for solution-adaptive grids. The discussion is presented in viewgraph form.

Rupak Biswas↗

NASA NTRS2013

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗

NASA NTRS2003

Solution-Adaptive Program for Computing 2D/Axi Viscous Flow

A computer program solves the Navier- Stokes equations governing the flow of a viscous, compressible fluid in an axisymmetric or two-dimensional (2D) setting. To obtain solutions more accurate than those generated by prior such programs that utilize regular and/or fixed computational meshes, this program utilizes unstructured (that is, irregular triangular) computational meshes that are automatically adapted to solutions. The adaptation can refine to regions of high change in gradient or can be driven by a novel residual minimization technique. Starting from an initial mesh and a corresponding data structure, the adaptation of the mesh is controlled by use of minimization functional. Other improvements over prior such programs include the following: (1) Boundary conditions are imposed weakly; that is, following initial specification of solution values at boundary nodes, these values are relaxed in time by means of the same formulations as those used for interior nodes. (2) Eigenvalues are limited in order to suppress expansion shocks. (3) An upwind fluctuation-splitting distribution scheme applied to inviscid flux requires fewer operations and produces less artificial dissipation than does a finite-volume scheme, leading to greater accuracy of solutions.

Wood, William A.↗

NASA NTRS1994

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on an MIMD (multiple instruction/multiple data stream) parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed, one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the iPSC/860.

Venkatakrishnan, V.↗