SEARCH · Search NASA

Results for “Mesh optimization”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Quasi-Optimal Schwarz Methods for the Conforming Spectral Element Discretization

Fast methods are proposed for solving the system K(sub N)x = b resulting from the discretization of self-adjoint elliptic equations in three dimensional domains by the spectral element method. The domain is decomposed into hexahedral elements, and in each of these elements the discretization space is formed by polynomials of degree N in each variable. Gauss-Lobatto-Legendre (GLL) quadrature rules replace the integrals in the Galerkin formulation. This system is solved by the preconditioned conjugate gradients method. The conforming finite element space on the GLL mesh consisting of piecewise Q(sub 1) elements produces a stiffness matrix K(sub h) that is spectrally equivalent to the spectral element stiffness matrix K(sub N). The action of the inverse of K(sub h) is expensive for large problems, and is therefore replaced by a Schwarz preconditioner B(sub h) of this finite element stiffness matrix. The preconditioned operator then becomes B(sub h)(exp -l)K(sub N). The technical difficulties stem from the nonregularity of the mesh. Tools to estimate the convergence of a large class of new iterative substructuring and overlapping Schwarz preconditioners are developed. This technique also provides a new analysis for an iterative substructuring method proposed by Pavarino and Widlund for the spectral element discretization.

Casarin, Mario↗

NASA NTRS1991

A geometric nonlinear degenerated shell element using a mixed formulation with independently assumed strain fields

A mixed formulation is chosen to overcome deficiencies of the standard displacement-based shell model. Element development is traced from the incremental variational principle on through to the final set of equilibrium equations. Particular attention is paid to developing specific guidelines for selecting the optimal set of strain parameters. A discussion of constraint index concepts and their predictive capability related to locking is included. Performance characteristics of the elements are assessed in a wide variety of linear and nonlinear plate/shell problems. Despite limiting the study to geometric nonlinear analysis, a substantial amount of additional insight concerning the finite element modeling of thin plate/shell structures is provided. For example, in nonlinear analysis, given the same mesh and load step size, mixed elements converge in fewer iterations than equivalent displacement-based models. It is also demonstrated that, in mixed formulations, lower order elements are preferred. Additionally, meshes used to obtain accurate linear solutions do not necessarily converge to the correct nonlinear solution. Finally, a new form of locking was identified associated with employing elements designed for biaxial bending in uniaxial bending applications.

Graf, Wiley E.↗

NASA NTRS1996

Computational Issues Associated with Temporally Deforming Geometries Such as Thrust Vectoring Nozzles

During the past decade, computational simulation of fluid flow around complex configurations has progressed significantly and many notable successes have been reported, however, unsteady time-dependent solutions are not easily obtainable. The present effort involves unsteady time dependent simulation of temporally deforming geometries. Grid generation for a complex configuration can be a time consuming process and temporally varying geometries necessitate the regeneration of such grids for every time step. Traditional grid generation techniques have been tried and demonstrated to be inadequate to such simulations. Non-Uniform Rational B-splines (NURBS) based techniques provide a compact and accurate representation of the geometry. This definition can be coupled with a distribution mesh for a user defined spacing. The present method greatly reduces cpu requirements for time dependent remeshing, facilitating the simulation of more complex unsteady problems. A thrust vectoring nozzle has been chosen to demonstrate the capability as it is of current interest in the aerospace industry for better maneuverability of fighter aircraft in close combat and in post stall regimes. This current effort is the first step towards multidisciplinary design optimization which involves coupling the aerodynamic heat transfer and structural analysis techniques. Applications include simulation of temporally deforming bodies and aeroelastic problems.

Boyalakuntla, Kishore↗

NASA NTRS1990

Vectorization of a particle code used in the simulation of rarefied hypersonic flow

A limitation of the direct simulation Monte Carlo (DSMC) method is that it does not allow efficient use of vector architectures that predominate in current supercomputers. Consequently, the problems that can be handled are limited to those of one- and two-dimensional flows. This work focuses on a reformulation of the DSMC method with the objective of designing a procedure that is optimized to the vector architectures found on machines such as the Cray-2. In addition, it focuses on finding a better balance between algorithmic complexity and the total number of particles employed in a simulation so that the overall performance of a particle simulation scheme can be greatly improved. Simulations of the flow about a 3D blunt body are performed with 10 to the 7th particles and 4 x 10 to the 5th mesh cells. Good statistics are obtained with time averaging over 800 time steps using 4.5 h of Cray-2 single-processor CPU time.

Baganoff, D.↗

NASA NTRSDate not supplied

Toward Fully 3D-Printed Two Degree of Freedom Acoustic Liners

An acoustic liner optimization tool is developed for designing two degree of freedom (2DOF) liners with high absorption over a wide range of frequencies and sound pressure levels (SPLs). Two additively manufactured 2DOF liners (one constant and one variable chamber depth) are designed and printed with an embedded perforate layer as the septum. A normal incidence impedance tube study is performed to directly compare impedance and absorption spectra to a more traditionally manufactured 2DOF liner with embedded mesh caps in phenolic honeycomb. Comparisons of test data to predictions are shown as well as SPL sensitivities for each acoustic liner. Broadband absorption is achieved with the 3D printed 2DOF liner containing constant chamber depths, although due to current limitations in printed embedded perforate hole size, not as broad as the traditional mesh cap liner. However, results also show that the 3D printed variable chamber depth sample is more comparable to the mesh cap liner absorption, demonstrating the viability of printed embedded perforates in novel concepts.

acoustic liner 2DOF two degree of freedom 3D print↗

NASA NTRS2012

Inviscid Analysis of Extended Formation Flight

Flying airplanes in extended formations, with separation distances of tens of wingspans, significantly improves safety while maintaining most of the fuel savings achieved in close formations. The present study investigates the impact of roll trim and compressibility at fixed lift coefficient on the benefits of extended formation flight. An Euler solver with adjoint-based mesh refinement combined with a wake propagation model is used to analyze a two-body echelon formation at a separation distance of 30 spans. Two geometries are examined: a simple wing and a wing-body geometry. Energy savings, quantified by both formation drag fraction and span efficiency factor, are investigated at subsonic and transonic speeds for a matrix of vortex locations. The results show that at fixed lift and trimmed for roll, the optimal location of vortex impingement is about 10% inboard of the trailing airplane s wing-tip. Interestingly, early results show the variation in drag fraction reduction is small in the neighborhood of the optimal position. Over 90% of energy benefits can be obtained with a 5% variation in transverse and 10% variation in crossflow directions. Early results suggest control surface deflections required to achieve trim reduce the benefits of formation flight by 3-5% at subsonic speeds. The final paper will include transonic effects and trim on extended formation flight drag benefits.

Kless, James↗

NASA NTRS1995

Supersonic wing and wing-body shape optimization using an adjoint formulation

This paper describes the implementation of optimization techniques based on control theory for wing and wing-body design of supersonic configurations. The work represents an extension of our earlier research in which control theory is used to devise a design procedure that significantly reduces the computational cost by employing an adjoint equation. In previous studies it was shown that control theory could be used to~eviseransonic design methods for airfoils and wings in which the shape and the surrounding body-fitted mesh are both generated analytically, and the control is the mapping function. The method has also been implemented for both transonic potential flows and transonic flows governed by the Euler equations using an alternative formulation which employs numerically generated grids, so that it can treat more general configurations. Here results are presented for three-dimensional design cases subject to supersonic flows governed by the Euler equation.

Reuther, James↗

NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗

NASA NTRS1980

A Navier-Stokes fast solver for turbulence modeling applications

A computer code for the evaluation and/or optimization of the predicative potential of second order turbulent closure models in simple two dimensional flow configurations is discussed. A procedure for the numerical solution of the steady constant property Navier-Stokes equations are described together with algebraic, one dimensional and two dimensional equations of turbulence closure models. Four turbulence models are compared with several sets of experimental data. The effects of initial conditions and boundary conditions are also described. The effects of purely numerical parameters, such as mesh size, boundary locations, and convergence criteria are presented.

Murphy, J. D.↗

NASA NTRS2015

Enabling Rapid and Robust Structural Analysis During Conceptual Design

This paper describes a multi-year effort to add a structural analysis subprocess to a supersonic aircraft conceptual design process. The desired capabilities include parametric geometry, automatic finite element mesh generation, static and aeroelastic analysis, and structural sizing. The paper discusses implementation details of the new subprocess, captures lessons learned, and suggests future improvements. The subprocess quickly compares concepts and robustly handles large changes in wing or fuselage geometry. The subprocess can rank concepts with regard to their structural feasibility and can identify promising regions of the design space. The automated structural analysis subprocess is deemed robust and rapid enough to be included in multidisciplinary conceptual design and optimization studies.

Eldred, Lloyd B.↗

NASA NTRS2004

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Dadone, A.↗

NASA NTRS1997

Aerodynamic Design Optimization on Unstructured Grids with a Continuous Adjoint Formulation

A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analyzed. The derivation of the costate equations is presented, and a second-order accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finite-difference gradients and several design examples are presented.

Anderson, W. Kyle↗

NASA NTRS1994

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.↗

NASA NTRS1995

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗

NASA NTRS1996

3D Structured Grid Adaptation

Grid adaptation for structured meshes is the art of using information from an existing, but poorly resolved, solution to automatically redistribute the grid points in such a way as to improve the resolution in regions of high error, and thus the quality of the solution. This involves: (1) generate a grid vis some standard algorithm, (2) calculate a solution on this grid, (3) adapt the grid to this solution, (4) recalculate the solution on this adapted grid, and (5) repeat steps 3 and 4 to satisfaction. Steps 3 and 4 can be repeated until some 'optimal' grid is converged to but typically this is not worth the effort and just two or three repeat calculations are necessary. They also may be repeated every 5-10 time steps for unsteady calculations.

Banks, D. W.↗

NASA NTRS2000

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

NASA NTRS2001

Hydrogen/Air Fuel Nozzle Emissions Experiments

The use of hydrogen combustion for aircraft gas turbine engines provides significant opportunities to reduce harmful exhaust emissions. Hydrogen has many advantages (no CO2 production, high reaction rates, high heating value, and future availability), along with some disadvantages (high current cost of production and storage, high volume per BTU, and an unknown safety profile when in wide use). One of the primary reasons for switching to hydrogen is the elimination of CO2 emissions. Also, with hydrogen, design challenges such as fuel coking in the fuel nozzle and particulate emissions are no longer an issue. However, because it takes place at high temperatures, hydrogen-air combustion can still produce significant levels of NOx emissions. Much of the current research into conventional hydrocarbon-fueled aircraft gas turbine combustors is focused on NOx reduction methods. The Zero CO2 Emission Technology (ZCET) hydrogen combustion project will focus on meeting the Office of Aerospace Technology goal 2 within pillar one for Global Civil Aviation reducing the emissions of future aircraft by a factor of 3 within 10 years and by a factor of 5 within 25 years. Recent advances in hydrocarbon-based gas turbine combustion components have expanded the horizons for fuel nozzle development. Both new fluid designs and manufacturing technologies have led to the development of fuel nozzles that significantly reduce aircraft emissions. The goal of the ZCET program is to mesh the current technology of Lean Direct Injection and rocket injectors to provide quick mixing, low emissions, and high-performance fuel nozzle designs. An experimental program is planned to investigate the fuel nozzle concepts in a flametube test rig. Currently, a hydrogen system is being installed in cell 23 at NASA Glenn Research Center's Research Combustion Laboratory. Testing will be conducted on a variety of fuel nozzle concepts up to combustion pressures of 350 psia and inlet air temperatures of 1200 F. Computational fluid dynamics calculations, with the Glenn developed National Combustor Code, are being performed to optimize the fuel nozzle designs.

Smith, Timothy D.↗

NASA NTRS2009

Microstrip Antenna for Remote Sensing of Soil Moisture and Sea Surface Salinity

This compact, lightweight, dual-frequency antenna feed developed for future soil moisture and sea surface salinity (SSS) missions can benefit future soil and ocean studies by lowering mass, volume, and cost of the antenna system. It also allows for airborne soil moisture and salinity remote sensors operating on small aircraft. While microstrip antenna technology has been developed for radio communications, it has yet to be applied to combined radar and radiometer for Earth remote sensing. The antenna feed provides a key instrument element enabling high-resolution radiometric observations with large, deployable antennas. The design is based on the microstrip stacked-patch array (MSPA) used to feed a large, lightweight, deployable, rotating mesh antenna for spaceborne L-band (approximately equal to 1 GHz) passive and active sensing systems. The array consists of stacked patches to provide dual-frequency capability and suitable radiation patterns. The stacked-patch microstrip element was designed to cover the required L-band center frequencies at 1.26 GHz (lower patch) and 1.413 GHz (upper patch), with dual-linear polarization capabilities. The dimension of patches produces the required frequencies. To achieve excellent polarization isolation and control of antenna sidelobes for the MSPA, the orientation of each stacked-patch element within the array is optimized to reduce the cross-polarization. A specialized feed-distribution network was designed to achieve the required excitation amplitude and phase for each stacked-patch element.

Ramhat-Samii, Yahya↗