SEARCH · Engineering Papers

Results for “Mesh Refinement”

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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Adaptive finite element strategies for shell structures

The present paper extends existing finite element adaptive refinement (AR) techniques to shell structures, which have heretofore been neglected in the AR literature. Specific challenges in applying AR to shell structures include: (1) physical discontinuities (e.g., stiffener intersections); (2) boundary layers; (3) sensitivity to geometric imperfections; (4) the sensitivity of most shell elements to mesh distortion, constraint definition and/or thinness; and (5) intrinsic geometric nonlinearity. All of these challenges but (5) are addressed here.

Stanley, G.↗

NASA NTRS1988

Accurate multigrid solution of the Euler equations on unstructured and adaptive meshes

A method for accurately solving inviscid compressible flow in the subcritical and supercritical regimes about complex configurations is presented. The method is based on the use of unstructured triangular meshes in two dimensions, and special emphasis is placed on the accuracy and efficiency of the solutions. High accuracy is achieved by careful scaling of the artificial dissipation terms, and by reformulating the inner and outer boundary conditions for both the convective and dissipative operators. An adaptive grid refinement strategy is presented which enhances the solution accuracy for complex flows. When coupled with an unstructured multigrid algorithm, this method is shown to produce an efficient solver for flows about arbitrary configurations.

Mavriplis, Dimitri J.↗

NASA NTRS1988

Accurate multigrid solution of the Euler equations on unstructured and adaptive meshes

Mavriplis, Dimitri J.↗

NASA NTRS2012

Repair Concepts as Design Constraints of a Stiffened Composite PRSEUS Panel

A design and analysis of a repair concept applicable to a stiffened thin-skin composite panel based on the Pultruded Rod Stitched Efficient Unitized Structure is presented. The concept is a bolted repair using metal components, so that it can easily be applied in the operational environment. The damage scenario considered is a midbay-to-midbay saw-cut with a severed stiffener, flange and skin. In a previous study several repair configurations were explored and their feasibility confirmed but refinement was needed. The present study revisits the problem under recently revised design requirements and broadens the suite of loading conditions considered. The repair assembly design is based on the critical tension loading condition and subsequently its robustness is verified for a pressure loading case. High fidelity modeling techniques such as mesh-independent definition of compliant fasteners, elastic-plastic material properties for metal parts and geometrically nonlinear solutions are utilized in the finite element analysis. The best repair design is introduced, its analysis results are presented and factors influencing the design are assessed and discussed.

Przekop, Adam↗

NASA NTRS1993

A solution-adaptive hybrid-grid method for the unsteady analysis of turbomachinery

A solution-adaptive method for the time-accurate analysis of two-dimensional flows in turbomachinery is described. The method employs a hybrid structured-unstructured zonal grid topology in conjunction with appropriate modeling equations and solution techniques in each zone. The viscous flow region in the immediate vicinity of the airfoils is resolved on structured O-type grids while the rest of the domain is discretized using an unstructured mesh of triangular cells. Implicit, third-order accurate, upwind solutions of the Navier-Stokes equations are obtained in the inner regions. In the outer regions, the Euler equations are solved using an explicit upwind scheme that incorporates a second-order reconstruction procedure. An efficient and robust grid adaptation strategy, including both grid refinement and coarsening capabilities, is developed for the unstructured grid regions. Grid adaptation is also employed to facilitate information transfer at the interfaces between unstructured grids in relative motion. Results for grid adaptation to various features pertinent to turbomachinery flows are presented. Good comparisons between the present results and experimental measurements and earlier structured-grid results are obtained.

Mathur, Sanjay R.↗

NASA NTRS1994

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

NASA NTRS1994

Effects of a Forward-swept Front Rotor on the Flowfield of a Counterrotation Propeller

The effects of a forward-swept front rotor on the flowfield of a counterrotation model propeller at takeoff conditions at zero degree angle of attack are studied by solving the unsteady three-dimensional Euler equations. The configuration considered is an uneven blade count counterrotation model with twelve forward-swept blades on the fore rotor and ten aft-swept blades on the aft rotor. The flowfield is compared with that of a reference aft-swept counterrotation geometry and Laser Doppler Velocimeter (LDV) measurements. At the operating conditions considered, the forward-swept blade experiences a higher tip loading and produces a stronger tip vortex compared to the aft-swept blade, consistent with the LDV and acoustic measurements. Neither the solution nor the LDV data indicated the formation of a leading edge vortex. The predicted radial distribution of the circumferentially averaged axial velocity at the measurement station agreed very closely with LDV data, while crossflow velocities showed poor agreement. The discrepancy between prediction and LDV data of tangential and radial velocities is due in part to the insufficient mesh resolution in the region between the rotors and in the tip region to track the tip vortex. The vortex is diffused by the time it arrives at the measurement station. The uneven blade count configuration requires the solution to be carried out for six blade passages of the fore rotor and five passages of the aft rotor, thus making grid refinement prohibitive.

Nallasamy, M.↗

NASA NTRS1988

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

NASA NTRS1988

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗

NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗

NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme-based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗

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 NTRSDate not supplied

High Fidelity Adaptively Refined CFD and Reduced Order Models of a High Aspect Ratio Aeroelastic Wing Wind-Tunnel Model

The NASA Advanced Air Transport Technology (AATT) goal of reduced fuel burn for transport aircraft has led to the NASA N+3 High Aspect Ratio Wing (HARW) subproject. This project requires identifying, developing, and demonstrating key technologies and integrated multidisciplinary solutions to enable a safe, high performance, aeroelastic wing. Since this aircraft will have a high aspect ratio wing, aeroelasticity is expected to be a major issue in the design. In this paper high fidelity computational fluid dynamics (CFD) is performed with flow adapted meshes. A system identification of the aerodynamics is developed using both a multi-modal multi-sine time-marching and a multi-mode linear frequency domain method. GLA, MLA and flutter suppression simulations will be performed.

Robert Bartels↗

NASA NTRS2014

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

NASA NTRS2002

On Accuracy of Adaptive Grid Methods for Captured Shocks

The accuracy of two grid adaptation strategies, grid redistribution and local grid refinement, is examined by solving the 2-D Euler equations for the supersonic steady flow around a cylinder. Second- and fourth-order linear finite difference shock-capturing schemes, based on the Lax-Friedrichs flux splitting, are used to discretize the governing equations. The grid refinement study shows that for the second-order scheme, neither grid adaptation strategy improves the numerical solution accuracy compared to that calculated on a uniform grid with the same number of grid points. For the fourth-order scheme, the dominant first-order error component is reduced by the grid adaptation, while the design-order error component drastically increases because of the grid nonuniformity. As a result, both grid adaptation techniques improve the numerical solution accuracy only on the coarsest mesh or on very fine grids that are seldom found in practical applications because of the computational cost involved. Similar error behavior has been obtained for the pressure integral across the shock. A simple analysis shows that both grid adaptation strategies are not without penalties in the numerical solution accuracy. Based on these results, a new grid adaptation criterion for captured shocks is proposed.

Yamaleev, Nail K.↗

NASA NTRSDate not supplied

Transition prediction in hypersonic regime on complex geometries with RANS-based models

In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.

CFD↗

NASA NTRSDate not supplied

Transition Prediction in Hypersonic Regime on Complex Geometries With Rans-Based Models

CFD↗

NASA NTRSDate not supplied

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗