SEARCH · Engineering Papers

Results for “Galerkin method”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Application of p-Multigrid to Discontinuous Galerkin Formulations of the Poisson Equation

We investigate p-multigrid as a solution method for several different discontinuous Galerkin (DG) formulations of the Poisson equation. Different combinations of relaxation schemes and basis sets have been combined with the DG formulations to find the best performing combination. The damping factors of the schemes have been determined using Fourier analysis for both one and two-dimensional problems. One important finding is that when using DG formulations, the standard approach of forming the coarse p matrices separately for each level of multigrid is often unstable. To ensure stability the coarse p matrices must be constructed from the fine grid matrices using algebraic multigrid techniques. Of the relaxation schemes, we find that the combination of Jacobi relaxation with the spectral element basis is fairly effective. The results using this combination are p sensitive in both one and two dimensions, but reasonable convergence rates can still be achieved for moderate values of p and isotropic meshes. A competitive alternative is a block Gauss-Seidel relaxation. This actually out performs a more expensive line relaxation when the mesh is isotropic. When the mesh becomes highly anisotropic, the implicit line method and the Gauss-Seidel implicit line method are the only effective schemes. Adding the Gauss-Seidel terms to the implicit line method gives a significant improvement over the line relaxation method.

Helenbrook, B. T.↗

NASA NTRS1989

Model-free simulations of turbulent reactive flows

The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.

Givi, Peyman↗

NASA NTRS1980

Response of nonlinear panels to random loads

Lightweight aircraft structures exposed to a high intensity noise environment can fatigue prematurely if adequate consideration is not given to the problem. Design methods and design criteria for sonic fatigue prevention were developed based on analytical and experimental techniques. Most of the analytical work was based upon small deflection or linear structural theory which did not agree with the experimental results. A large deflection geometrical nonlinearity was incorporated into the analysis methods for determining the structural response to high intensity noise. The Karman-Herrmann large deflection equations with a single mode Galerkin approximation, and the method of equivalent linearization were used to predict mean square amplitude, mean square stresses, and nonlinear frequency at various acoustic loadings for rectangular panels. Both simply supported and clamped support conditions with immovable or movable inplane edges are considered. Comparisons with experimental results are presented.

Mei, C.↗

NASA NTRS2005

ON THE ROLE OF INVOLUTIONS IN THE DISCONTINUOUS GALERKIN DISCRETIZATION OF MAXWELL AND MAGNETOHYDRODYNAMIC SYSTEMS

The role of involutions in energy stability of the discontinuous Galerkin (DG) discretization of Maxwell and magnetohydrodynamic (MHD) systems is examined. Important differences are identified in the symmetrization of the Maxwell and MHD systems that impact the construction of energy stable discretizations using the DG method. Specifically, general sufficient conditions to be imposed on the DG numerical flux and approximation space are given so that energy stability is retained These sufficient conditions reveal the favorable energy consequence of imposing continuity in the normal component of the magnetic induction field at interelement boundaries for MHD discretizations. Counterintuitively, this condition is not required for stability of Maxwell discretizations using the discontinuous Galerkin method.

Barth, Timothy↗

NASA NTRS1990

The semidiscrete Galerkin finite element modeling of compressible viscous flow past an airfoil

The primary project was the numerical simulation, by a finite element/finite difference method, of the viscous flow about an airfoil. The secondary project involved the numerical simulation of the three-dimensional separated and vortex-dominated flow about a hemispherically capped cylinder in the transonic regime. Preliminary calculations were started for the hemisphere-cylinder at 0 and 5 degree angle of attack. The solution of the flow field about airfoils and wings is required to determine the important parameters of lift, moment, and drag. Viscous effects must be accounted for if the drag is to be accurately calculated. At present there are basically two approaches to the numerical simulation of the flow field, the use of fully viscous models and the inviscid/viscous models. The fully viscous models require the solution of an approximation of the Navier-Stokes equations and therefore should simulate most of the physical mechanisms. A fast, accurate, and computationally efficient inviscid flow solver was recently developed by Hartwich. It is thought that Hartwich's program coupled to a fast, accurate, and computationally efficient boundary layer code, will make an excellent tool for airfoil design. The purpose of the primary project was to develop a compressible boundary layer code using the semidiscrete Galerkin finite element method. The numerical scheme employed used the combination of a Dorodnitsyn formulation of the boundary layer equations, with a finite difference/finite element procedure (semidiscrete Galerkin method), in the solution of the compressible two-dimensional boundary layer equations. A laminar compressible boundary layer code was developed and tested for a NACA 0012 airfoil at a Mach number of 0.5, a Reynolds number of 5000, and zero angle of attack. At present the boundary layer program solves up to, but not beyond, separation.

Meade, Andrew J.↗

NASA NTRS1992

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

NASA NTRS1987

Suboptimal shape control for quasi-static distributed-parameter systems

An on-line control approach which will adjust the steady-state shape of a large antenna arbitrarily close to any achievable desired profile is approached. The method makes use of distributed-parameter system theory and allows refocusing using a limited number of control actuators and sensors. The controller gains are calculated by approximating the solution to an infinite-dimensional optimal quasi-static control problem. The controller gain calculation is computationally simpler than that previously proposed. The Galerkin (finite element) approximation method is used for model reduction. It is proved that both gain and state convergence can be achieved by using the proposed approximation scheme.

Balas, M. J.↗

NASA NTRS1977

Three-dimensional modelling in magnetotelluric and magnetic variational sounding

The Galerkin finite-element method is used to obtain approximate solutions for the three-dimensional induction problem. A rectangular conductive prism is considered as an example, and solutions are obtained for linear and circularly polarized incident plane-wave fields. Magnetotelluric tensor impedances and magnetic transfer functions are computed. Polar diagrams of the tensor impedances and magnetic transfer functions along with their amplitude contour maps are presented. The dimensionality parameter, skew, is contoured at the surface of the earth. It is shown that the relative amplitudes and shapes of the additional and principal impedance polar diagrams can be used to determine the dimensionality of geoelectrical structures. Stations with skew values greater than 0.2 are significantly influenced by the three-dimensionality of the geoelectric structure. The amplitudes of the magnetic transfer function and the orientations of its polar diagrams exhibit large anomalies in the vicinity of the intersection of the lateral contacts.

Reddy, I. K.↗

NASA NTRS1979

A split finite element algorithm for the compressible Navier-Stokes equations

An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms.

Baker, A. J.↗

NASA NTRS1988

Evaluation and applicability of transfinite element formulations with particular reference to radiation

The evaluation and applicability of hybrid transfinite element formulations are described for transient nonlinear thermal models with particular reference to account for radiation effects. The formulations are developed from basic concepts and evaluated for both one- and two-dimensional thermal models. The methodology presented is a hybrid approach as it combines finite elements, classical Galerkin schemes, and transform methods. Basic features of the hybrid formulations and the associated solution scheme used are first described. Therein, numerical test models are presented to evaluate the applicability for predicting the nonlinear transient response in one- and two-dimensional thermal models influenced due to radiation effects.

Tamma, Kumar K.↗

NASA NTRS1998

An Astronomical Test of CCD Photometric Precision

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques. we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.

Koch, David↗

NASA NTRS1999

A Posteriori Error Estimation for Discontinuous Galerkin Approximations of Hyperbolic Systems

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.

Larson, Mats G.↗

NASA NTRS2000

A Finite-Element Approach for Modeling Inviscid and Viscous Compressible Flows using Prismatic Grids

The Galerkin finite-element method is used to solve the Euler and Navier-Stokes equations on prismatic meshes. It is shown that the prismatic grid is advantageous for correctly and efficiently capturing the boundary layers in high Reynolds number flows. It can be captured accurately because of the ability to cluster grid points normal to the body. The efficiency derives from the implicit treatment of the normal direction. To treat the normal direction implicitly, a semi-implicit Runge-Kutta time stepping scheme is developed. The semi-implicit algorithm is validated on simple geometries for inviscid and viscous flows and its convergence history is compared to that of the explicit Runge-Kutta scheme. The semi-implicit scheme is shown to be a factor of 3 to 4 faster in terms of CPU time to convergence.

Pandya, S. A.↗

NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗

NASA NTRS1993

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗

NASA NTRS1978

Efficiency trade-offs of steady-state methods using FEM and FDM

The efficiency characteristics of finite element and finite difference approximations for the steady-state solution of the Navier-Stokes equations are examined. The finite element method discussed is a standard Galerkin formulation of the incompressible, steady-state Navier-Stokes equations. The finite difference formulation uses simple centered differences that are O(delta x-squared). Operation counts indicate that a rapidly converging Newton-Raphson-Kantorovitch iteration scheme is generally preferable over a Picard method. A split NOS Picard iterative algorithm for the finite difference method was most efficient.

Gartling, D. K.↗

NASA NTRS1984

Optimal quasi-static shape control for large aerospace antennae

In this paper, an on-line control approach which will adjust the steady-state shape of a large antenna arbitarily close to any achievable desired profile is proposed. The method makes use of distributed parameter system theory and allows refocusing using a limited number of control actuators and sensors. The controller gains are calculated by approximating the solution to an infinite-dimensional optimal quasi-static control problem. A very general convergence result for such quasi-static controllers is proved and applied to the antenna controller to show convergence, using any Galerkin (finite-element) approximation method.

Balas, M. J.↗

NASA NTRS1979

The application of finite element techniques to acoustic transmission in lined ducts with flow

The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.

Astley, R. J.↗