Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Generation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗
NASA NTRS2000

Handbook on Face Gear Drives with a Spur Involute Pinion

The use of face gears in power transmission and drive systems has a significant number of benefits. Face gears allow a variety of new transmission arrangements as well as high reduction ratio capability. This leads to drive system weight reduction and improvements in performance. In this work, basic information about the design and analysis of face gear drives is presented. The work considers face gears in mesh with spur involute pinions for both intersecting axes and offset drives. Tooth geometry, kinematics, generation of face gears with localized bearing contact by cutting and grinding, avoidance of tooth undercutting, avoidance of tooth pointing, tooth contact analysis, and algorithms for the simulation of meshing and contact arc all topics which are discussed. In addition, applications of face gear drives are presented. Included are design uses in aerospace applications such as helicopter transmissions, split-torque face gear arrangements, comparisons of face gears with bevel gears, and general design considerations.

Litvin, F. L.↗
NASA NTRS1990

An overlapped grid method for multigrid, finite volume/difference flow solvers: MaGGiE

The objective is to develop a domain decomposition method via overlapping/embedding the component grids, which is to be used by upwind, multi-grid, finite volume solution algorithms. A computer code, given the name MaGGiE (Multi-Geometry Grid Embedder) is developed to meet this objective. MaGGiE takes independently generated component grids as input, and automatically constructs the composite mesh and interpolation data, which can be used by the finite volume solution methods with or without multigrid convergence acceleration. Six demonstrative examples showing various aspects of the overlap technique are presented and discussed. These cases are used for developing the procedure for overlapping grids of different topologies, and to evaluate the grid connection and interpolation data for finite volume calculations on a composite mesh. Time fluxes are transferred between mesh interfaces using a trilinear interpolation procedure. Conservation losses are minimal at the interfaces using this method. The multi-grid solution algorithm, using the coaser grid connections, improves the convergence time history as compared to the solution on composite mesh without multi-gridding.

Baysal, Oktay↗
NASA NTRS1991

A new technique for simulating composite material

This project dealt with the development on new methodologies and algorithms for the multi-spectrum electromagnetic characterization of large scale nonmetallic airborne vehicles and structures. A robust, low memory, and accurate methodology was developed which is particularly suited for modern machine architectures. This is a hybrid finite element method that combines two well known numerical solution approaches. That of the finite element method for modeling volumes and the boundary integral method which yields exact boundary conditions for terminating the finite element mesh. In addition, a variety of high frequency results were generated (such as diffraction coefficients for impedance surfaces and material layers) and a class of boundary conditions were developed which hold promise for more efficient simulations. During the course of this project, nearly 25 detailed research reports were generated along with an equal number of journal papers. The reports, papers, and journal articles are listed in the appendices along with their abstracts.

Volakis, John L.↗
NASA NTRS2012

Integrated Design Engineering Analysis (IDEA) Environment Automated Generation of Structured CFD Grids using Topology Methods

This report documents the work performed from March 2010 to March 2012. The Integrated Design and Engineering Analysis (IDEA) environment is a collaborative environment based on an object-oriented, multidisciplinary, distributed framework using the Adaptive Modeling Language (AML) as a framework and supporting the configuration design and parametric CFD grid generation. This report will focus on describing the work in the area of parametric CFD grid generation using novel concepts for defining the interaction between the mesh topology and the geometry in such a way as to separate the mesh topology from the geometric topology while maintaining the link between the mesh topology and the actual geometry.

Kamhawi, Hilmi N.↗
NASA NTRS2004

Development of a Flexible Framework for Hypersonic Navier-Stoke Space Shuttle Orbiter Meshes

A flexible framework constructing block structured volume grids for hypersonic Navier-Strokes flow simulations was developed for the analysis of the Shuttle Orbiter Columbia. The development of the framework, which was partially basedon the requirements of the primary flow solvers used resulted in an ability to directly correlate solutions contributed by participating groups on a common surface mesh. A foundation was built through the assessment of differences between differnt solvers, which provided confidence for independent assessment of other damage scenarios by team members. The framework draws on the experience of NASA Langley and NASA Ames Research Centers in structured grid generation, and consists of a grid generation, and consist of a grid generation process implemented through a division of responsibilities. The nominal division of labor consisted of NASA Johnson Space Center coordinating the damage scenarios to be analyzed by the Aerothermodynamics Columbia Accident Investigation (ACAI) team, Ames developing the surface grids that described the computational volume about the Orbiter, and Langley improving grid quality of Ames generated data and constructing the final computational volume grids. Distributing the work among the participant in th ACAI team resulted in significantl less time required to construct complete meshes than possible by any individual participant. The approach demonstrated that the One-NASA grid generation team could sustain the demand of for five new meshes to explore new damage scenarios within an aggressive time-line.

Alter, Stephen J.↗
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 NTRS1989

Three dimensional unstructured grids for the solution of the Euler equations

The advancing front technique is being used to develop a code to generate grids around complex 3-D configurations for use in computing the invisid flow solutions by the Euler equations. By the advancing front technique points are introduced concurrently with the connectivity information so that a separate library is not required. The generation of a 3-D grid is accomplished in several steps. First the boundaries of the domain to be gridded must be described by two-, three- or four-sided surface patches. Next, a background mesh is required to control the grid spacing and stretching throughout the domain. This coarse tetrahedral grid is not required to conform to any of the boundaries. Next, each of the patches is mapped to 2-D, triangulated by the advancing front technique and mapped back to 3-D. These triangles form the initial front for the generation of the final tetrahedral mesh.

Gumbert, Clyde↗
NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit scheme can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗
NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit schemes can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗
NASA NTRS2000

A Very High Order, Adaptable MESA Implementation for Aeroacoustic Computations

Since computational efficiency and wave resolution scale with accuracy, the ideal would be infinitely high accuracy for problems with widely varying wavelength scales. Currently, many of the computational aeroacoustics methods are limited to 4th order accurate Runge-Kutta methods in time which limits their resolution and efficiency. However, a new procedure for implementing the Modified Expansion Solution Approximation (MESA) schemes, based upon Hermitian divided differences, is presented which extends the effective accuracy of the MESA schemes to 57th order in space and time when using 128 bit floating point precision. This new approach has the advantages of reducing round-off error, being easy to program. and is more computationally efficient when compared to previous approaches. Its accuracy is limited only by the floating point hardware. The advantages of this new approach are demonstrated by solving the linearized Euler equations in an open bi-periodic domain. A 500th order MESA scheme can now be created in seconds, making these schemes ideally suited for the next generation of high performance 256-bit (double quadruple) or higher precision computers. This ease of creation makes it possible to adapt the algorithm to the mesh in time instead of its converse: this is ideal for resolving varying wavelength scales which occur in noise generation simulations. And finally, the sources of round-off error which effect the very high order methods are examined and remedies provided that effectively increase the accuracy of the MESA schemes while using current computer technology.

Dydson, Roger W.↗
NASA NTRS2013

Method for Enhancing a Three Dimensional Image from a Plurality of Frames of Flash LIDAR Data

A method for enhancing a three dimensional image from frames of flash LIDAR data includes generating a first distance R(sub i) from a first detector i to a first point on a surface S(sub i). After defining a map with a mesh theta having cells k, a first array S(k), a second array M(k), and a third array D(k) are initialized. The first array corresponds to the surface, the second array corresponds to the elevation map, and the third array D(k) receives an output for the DEM. The surface is projected onto the mesh theta, so that a second distance R(sub k) from a second point on the mesh theta to the detector can be found. From this, a height may be calculated, which permits the generation of a digital elevation map. Also, using sequential frames of flash LIDAR data, vehicle control is possible using an offset between successive frames.

Vanek, Michael D.↗
NASA NTRS2014

Methods of Real Time Image Enhancement of Flash LIDAR Data and Navigating a Vehicle Using Flash LIDAR Data

A method for creating a digital elevation map ("DEM") from frames of flash LIDAR data includes generating a first distance R(sub i) from a first detector i to a first point on a surface S(sub i). After defining a map with a mesh THETA having cells k, a first array S(k), a second array M(k), and a third array D(k) are initialized. The first array corresponds to the surface, the second array corresponds to the elevation map, and the third array D(k) receives an output for the DEM. The surface is projected onto the mesh THETA, so that a second distance R(sub k) from a second point on the mesh THETA to the detector can be found. From this, a height may be calculated, which permits the generation of a digital elevation map. Also, using sequential frames of flash LIDAR data, vehicle control is possible using an offset between successive frames.

Vanek, Michael D.↗
NASA NTRS1980

An implicit finite-difference code for inviscid and viscous cascade flow

An implicit finite-difference code is developed to solve either inviscid or viscous flow about two-dimensional cascade blade elements. General coordinate transformations are used so that boundaries can coincide with coordinate lines, and an automatic grid generation routine based on elliptic partial differential equations is employed to mesh arbitrary cascade elements. Characteristic combinations of the differential equations are used at inflow and outflow boundaries. Computed results for both inviscid and viscous flow are compared with other existing cascade solutions and experimental data.

Steger, J. L.↗
NASA NTRS1983

Theoretical investigations of high lift aerodynamics

A program which generates a coordinate system for a two element airfoil with the mesh points concentrated in areas of significant vorticity, i.e., boundary layer and wake is operational. The 'imbedded' grid method developed allows a transition from the scale of the main airfoil to the scale of the flap. This requirement is essential for the modeling of viscous flows over the flap and slat of a multielement airfoil. An airfoil mounted in a 2-D wind tunnel was formulated. The program is ready for a fine grid and a large number of planes to explore the characteristics of a Navier-Stokes solver in a quasi-3D case. The program was converted to a form suitable for the STAR computer. Runs were made to map a three dimensional flow field for a wall airfoil intersection with and without lift.

Bennett, G.↗
NASA NTRS1985

An adaptive finite element procedure for compressible high speed flows

A practical finite element-based solution procedure for high-speed inviscid compressible flow problems is described. The method provides time-accurate solutions to the coompressible Euler equations, and is computationally more efficient than the one-step Taylor-Galerkin approach and better suited for implementation on the modern generation of vector computers. The method is coupled to an adaptive mesh refinement process that enables steady state solutions of improved quality to be obtained.

Loehner, R.↗
NASA NTRS1990

Euler analysis comparison with LDV data for an advanced counter-rotation propfan at cruise

A fine mesh Euler solution of the F4/A4 unducted fan (UDF) model flowfield is compared with laser Doppler velocimeter (LDV) data taken in the NASA Lewis 8- by 6-Foot Supersonic Wind Tunnel. The comparison is made primarily at one axial plane downstream of the front rotor where the LDV particle lag errors are reduced. The agreement between measured and predicted velocities in this axial plane is good. The results show that a dense mesh is needed in the centerbody stagnation region to minimize entropy generation that weakens the aft row passage shock. The predicted radial location of the tip vortex downstream of the front rotor agrees well with the experimental results but the strength is overpredicted. With 40 points per chord line, the integrated performance quantities are nearly converged, but more points are needed to resolve passage shocks and flow field details.

Miller, Christopher J.↗