Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geometry meshing”

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 271 records · Page 15

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Expansion of epicyclic gear dynamic analysis program

The multiple mesh/single stage dynamics program is a gear tooth analysis program which determines detailed geometry, dynamic loads, stresses, and surface damage factors. The program can analyze a variety of both epicyclic and single mesh systems with spur or helical gear teeth including internal, external, and buttress tooth forms. The modifications refine the options for the flexible carrier and flexible ring gear rim and adds three options: a floating Sun gear option; a natural frequency option; and a finite element compliance formulation for helical gear teeth. The option for a floating Sun incorporates two additional degrees of freedom at the Sun center. The natural frequency option evaluates the frequencies of planetary, star, or differential systems as well as the effect of additional springs at the Sun center and those due to a flexible carrier and/or ring gear rim. The helical tooth pair finite element calculated compliance is obtained from an automated element breakup of the helical teeth and then is used with the basic gear dynamic solution and stress postprocessing routines. The flexible carrier or ring gear rim option for planetary and star spur gear systems allows the output torque per carrier and ring gear rim segment to vary based on the dynamic response of the entire system, while the total output torque remains constant.

Boyd, Linda Smith↗

Dynamics of a split torque helicopter transmission

Split torque designs, proposed as alternatives to traditional planetary designs for helicopter main rotor transmissions, can save weight and be more reliable than traditional designs. This report presents the results of an analytical study of the system dynamics and performance of a split torque gearbox that uses a balance beam mechanism for load sharing. The Lagrange method was applied to develop a system of equations of motion. The mathematical model includes time-varying gear mesh stiffness, friction, and manufacturing errors. Cornell's method for calculating the stiffness of spur gear teeth was extended and applied to helical gears. The phenomenon of sidebands spaced at shaft frequencies about gear mesh fundamental frequencies was simulated by modeling total composite gear errors as sinusoid functions. Although the gearbox has symmetric geometry, the loads and motions of the two power paths differ. Friction must be considered to properly evaluate the balance beam mechanism. For the design studied, the balance beam is not an effective device for load sharing unless the coefficient of friction is less than 0.003. The complete system stiffness as represented by the stiffness matrix used in this analysis must be considered to precisely determine the optimal tooth indexing position.

Krantz, Timothy L.↗

A Cartesian grid approach with hierarchical refinement for compressible flows

Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.

Quirk, James J.↗

Adaptive hybrid prismatic-tetrahedral grids for viscous flows

The paper presents generation of adaptive hybrid prismatic/tetrahedral grids for complex 3-D geometries including multi-body domains. The prisms cover the region close to each body's surface, while tetrahedra are created elsewhere. Two developments are presented for hybrid grid generation around complex 3-D geometries. The first is a new octree/advancing front type of method for generation of the tetrahedra of the hybrid mesh. The main feature of the present advancing front tetrahedra generator that is different from previous such methods is that it does not require the creation of a background mesh by the user for the determination of the grid-spacing and stretching parameters. These are determined via an automatically generated octree. The second development is an Automatic Receding Method (ARM) for treating the narrow gaps in between different bodies in a multiply-connected domain. This method is applied to a two-element wing case. A hybrid grid adaptation scheme that employs both h-refinement and redistribution strategies is developed to provide optimum meshes for viscous flow computations. Grid refinement is a dual adaptation scheme that couples division of tetrahedra, as well as 2-D directional division of prisms.

Kallinderis, Yannis↗

Algorithms for the automatic generation of 2-D structured multi-block grids

Two different approaches to the fully automatic generation of structured multi-block grids in two dimensions are presented. The work aims to simplify the user interactivity necessary for the definition of a multiple block grid topology. The first approach is based on an advancing front method commonly used for the generation of unstructured grids. The original algorithm has been modified toward the generation of large quadrilateral elements. The second method is based on the divide-and-conquer paradigm with the global domain recursively partitioned into sub-domains. For either method each of the resulting blocks is then meshed using transfinite interpolation and elliptic smoothing. The applicability of these methods to practical problems is demonstrated for typical geometries of fluid dynamics.

Schoenfeld, Thilo↗

Initial Experiments of High-Speed Drive System Windage Losses

High speed gearing performance is very important to the overall drive system efficiency. Certain losses such as gear meshing and bearing drag can be minimized by design changes such as pressure angle of the gears and the geometry and type of bearings being used. One component that can have a large effect on the overall performance of high-speed drive systems is the parasitic drag known as gear windage. This loss mechanism is not well understood and minimizing this component is usually accomplished through much trial and error. The results presented in this paper will document some of the design parameter effects on the amount of windage losses. A new test facility at NASA Glenn has been assembled to systematically study the design variables. Results from recent tests will be presented. The tests are for a single gear, with and without lubricants, and some initial studies using shrouds

Handschuh, Robert F.↗

Initial Experiments of High-Speed Drive System Windage Losses

High speed gearing performance is very important to the overall drive system efficiency. Certain losses such as gear meshing and bearing drag can be minimized by design changes such as pressure angle of the gears and the geometry and type of bearings being used. One component that can have a large effect on the overall performance of high-speed drive systems is the parasitic drag known as gear windage. This loss mechanism is not well understood and minimizing this component is usually accomplished through much trial and error. The results presented in this paper will document some of the design parameter effects on the amount of windage losses. A new test facility at NASA Glenn has been assembled to systematically study the design variables. Results from recent tests will be presented. The tests are for a single gear, with and without lubricants, and some initial studies using shrouds.

Handschuh, Robert F.↗

Automatic Mesh Generation of Hybrid Mesh on Valves in Multiple Positions in Feedline Systems

Fluid flow simulations through a valve often require evaluation of the valve in multiple opening positions. A mesh has to be generated for the valve for each position and compounding. The problem is the fact that the valve is typically part of a larger feedline system. In this paper, we propose to develop a system to create meshes for feedline systems with parametrically controlled valve openings. Herein we outline two approaches to generate the meshes for a valve in a feedline system at multiple positions. There are two issues that must be addressed. The first is the creation of the mesh on the valve for multiple positions. The second is the generation of the mesh for the total feedline system including the valve. For generation of the mesh on the valve, we will describe the use of topology matching and mesh generation parameter transfer. For generation of the total feedline system, we will describe two solutions that we have implemented. In both cases the valve is treated as a component in the feedline system. In the first method the geometry of the valve in the feedline system is replaced with a valve at a different opening position. Geometry is created to connect the valve to the feedline system. Then topology for the valve is created and the portion of the topology for the valve is topology matched to the standard valve in a different position. The mesh generation parameters are transferred and then the volume mesh for the whole feedline system is generated. The second method enables the user to generate the volume mesh on the valve in multiple open positions external to the feedline system, to insert it into the volume mesh of the feedline system, and to reduce the amount of computer time required for mesh generation because only two small volume meshes connecting the valve to the feedline mesh need to be updated.

Ross, Douglass H.↗

Effect of geometry and operating conditions on spur gear system power loss

The results of an analysis of the effects of spur gear size, pitch, width, and ratio on total mesh power loss for a wide range of speeds, torques, and oil viscosities are presented. The analysis uses simple algebraic expressions to determine gear sliding, rolling, and windage losses and also incorporates an approximate ball bearing power loss expression. The analysis shows good agreement with published data. Large diameter and fine pitched gears had higher peak efficiencies but low part load efficiency. Gear efficiencies were generally greater than 98 percent except at very low torque levels. Tare (no-load) losses are generally a significant percentage of the full load loss except at low speeds.

Anderson, N. E.↗

Effect of geometry and operating conditions on spur gear system power loss

The results of an analysis of the effects of spur gear size, pitch, width and ratio on total mesh power loss for a wide range of speeds, torques and oil viscosities are presented. The analysis uses simple algebraic expressions to determine gear sliding, rolling and windage losses and also incorporates an approximate ball bearing power loss expression. The analysis shows good agreement with published data. Large diameter and fine-pitched gears had higher peak efficiencies but lower part-load efficiency. Gear efficiencies were generally greater than 98 percent except at very low torque levels. Tare (no-load) losses are generally a significant percentage of the full-load loss except at low speeds.

Anderson, N. E.↗

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗

Aerodynamic Shape Optimization of a Dual-Stream Supersonic Plug Nozzle

Aerodynamic shape optimization was performed on an isolated axisymmetric plug nozzle sized for a supersonic business jet. The dual-stream concept was tailored to attenuate nearfield pressure disturbances without compromising nozzle performance. Adjoint-based anisotropic mesh refinement was applied to resolve nearfield compression and expansion features in the baseline viscous grid. Deformed versions of the adapted grid were used for subsequent adjoint-driven shape optimization. For design, a nonlinear gradient-based optimizer was coupled to the discrete adjoint formulation of the Reynolds-averaged Navier- Stokes equations. All nozzle surfaces were parameterized using 3rd order B-spline interpolants and perturbed axisymmetrically via free-form deformation. Geometry deformations were performed using 20 design variables shared between the outer cowl, shroud and centerbody nozzle surfaces. Interior volume grid deformation during design was accomplished using linear elastic mesh morphing. The nozzle optimization was performed at a design cruise speed of Mach 1.6, assuming core and bypass pressure ratios of 6.19 and 3.24, respectively. Ambient flight conditions at design were commensurate with 45,000-ft standard day atmosphere.

shape optimization↗

A Computational Study to Investigate the Effect of Defect Geometries on the Fatigue Crack Driving Forces in Powder-Bed AM Materials

Powder-bed additive manufacturing (AM) processes are associated with the formation of multiple types of process-specific pores, including but not limited to lack-of-fusion (LoF) and keyhole pores. The performance of an AM component is dependent on the type of pores, their density and their proximity to the free surface, and other heterogeneities in the microstructure. In order to characterize the influence of porosity on the mechanical behavior of AM materials, it is imperative to quantitatively analyze the heterogeneous strain accumulation in the vicinity of porosity. Process-specific microstructure models are generated using SPPARKS, an open-source process simulation code. Spherical keyhole or irregular LoF pores are embedded into the microstructure models, which are meshed and input into a finite element code, ScIFEN, to solve for the heterogeneous strain localization in the vicinity of the pores. Given the non-smooth geometries of LoF pores, they readily promote strain accumulation in their vicinity thereby increasing the propensity of initiating fatigue cracks.

Saikumar R Yeratapally↗

Dynamic Response of a Planetary Gear System Using a Finite Element/Contact Mechanics Model

The dynamic response of a helicopter planetary gear system is examined over a wide range of operating speeds and torques. The analysis tool is a unique, semianalytical finite element formulation that admits precise representation of the tooth geometry and contact forces that are crucial in gear dynamics. Importantly, no a priori specification of static transmission error excitation or mesh frequency variation is required; the dynamic contact forces are evaluated internally at each time step. The calculated response shows classical resonances when a harmonic of mesh frequency coincides with a natural frequency. However, peculiar behavior occurs where resonances expected to be excited at a given speed are absent. This absence of particular modes is explained by analytical relationships that depend on the planetary configuration and mesh frequency harmonic. The torque sensitivity of the dynamic response is examined and compared to static analyses. Rotation mode response is shown to be more sensitive to input torque than translational mode response.

Parker, Robert G.↗

3-D Deformable Grid Electomagnetic PIC for Parallel Computers

We have developed a three-dimensional non-orthogonal grid electromagnetic PIC code for parallel supercomputers. The numerical formulation is based on hexahedral cells which are logically connected cubic cells but distorted to body-fit complex geometries in physical space. An outstanding area of research is to improve the accuracy of the electromagnetic field solve on distorted meshes.

3D deformable non-orthogonal grid code PIC codes p↗