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 523 records · Page 29

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2022

Reduced-order modeling of neutron transport separated in energy by Minimax Proper Generalized Decomposition

In this article, we demonstrate a Petrov-Galerkin Proper Generalized Decomposition (PGD) known as Minimax PGD for modeling neutron transport separated in energy. To compare the Minimax with the classical Galerkin PGD, we assess both on a model problem of UO{sub 2} or Mixed Oxide (MOX) fuel pins in an infinite lattice with 3 industry-standard energy meshes. We find the Minimax PGD achieves a superior decomposition to Galerkin PGD, both with and without update of the energy modes. This suggests Minimax PGD may be more computationally efficient, provided this reduction in modes (to achieve a given accuracy) outweighs the cost of solving the necessary adjoint problems. In either case, we note that PGD offers an a priori Reduced-Order Model (ROM) which may be dramatically cheaper to solve than the full-order model, especially in problems with fine to ultrafine energy meshes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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

DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1983

Hybrid perturbation/Bubnov-Galerkin technique for nonlinear thermal analysis

A two step hybrid analysis technique to predict the nonlinear steady state temperature distribution in structures and solids is presented. The technique is based on the regular perturbation expansion and the classical Bubnov-Galerkin approximation. The functions are obtained by using the regular perturbation method. These functions are selected as coordinate functions and the classical Bubnov-Galerkin technique is used to compute their amplitudes. The potential of the proposed hybrid technique for the solution of nonlinear thermal problems is discussed. The effectiveness of this technique is demonstrated by the effects of conduction, convection, and radiation modes of heat transfer. It is indicated that the hybrid technique overcomes the two major drawbacks of the classical techniques: (1) the requirement of using a small parameter in the regular perturbation method; and (2) the arbitrariness in the choice of the coordinate functions in the Bubnov-Galerkin technique. The proposed technique extends the range of applicability of the regular perturbation method and enhances the effectiveness of the Bubnov-Galerkin technique.

Noor, A. K.↗

NASA NTRS1987

An adaptive finite element method for high speed flows

The solution of the equations of compressible high speed flow, on unstructured triangular grids in 2D and tetrahedral grids in 3D, is considered. Solution methods based upon both Taylor-Galerkin and Runge-Kutta time-stepping techniques are presented and the incorporation of the ideas of flux corrected transport (FCT) is discussed. These methods are combined with an adaptive mesh regeneration procedure and are employed in the solution of several examples, consisting of Euler flows in both 2D and 3D and Navier-Stokes flows in 2D.

Peraire, J.↗

DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗

NASA NTRS1987

A hybrid transfinite element approach for nonlinear transient thermal analysis

A new computational approach for transient nonlinear thermal analysis of structures is proposed. It is a hybrid approach which combines the modeling versatility of contemporary finite elements in conjunction with transform methods and classical Bubnov-Galerkin schemes. The present study is limited to nonlinearities due to temperature-dependent thermophysical properties. Numerical test cases attest to the basic capabilities and therein validate the transfinite element approach by means of comparisons with conventional finite element schemes and/or available solutions.

Tamma, Kumar K.↗

NASA NTRS1988

Application of finite element and remeshing technique to shock interference on a cylindrical leading edge

The problem of planar oblique shock impingement on a cylindrical leading edge in hypersonic flow is modeled using a Galerkin-Runge Kutta finite element method. The method utilizes a four stage Runge-Kutta time stepping scheme to solve the compressible Euler equations. Freestream Mach numbers of 6.5, 8.0 and 16.0 are studied. The computed surface pressure distributions consistently agree well with available experimental data. The peak pressure amplification ranges from 5.45 at M = 6.5 to approximately 17.0 at M = 16.0. Stagnation point heat transfer rate amplifications are calculated from the inviscid solution using the method of Fay and Riddell. The value and wall location of the peak pressure and heat transfer rate amplifications are extremely sensitive to the location of the impinging shock/bow shock intersection point.

Stewart, James R.↗

NASA NTRS1992

Convection in superposed fluid and porous layers

Thermal convection due to heating from below in a porous layer underlying a fluid layer has been analyzed using the Navier-Stokes equations for the fluid layers and the extended Darcy equation (including Brinkman and Forchheimer terms) for the porous layer. The flow is assumed to be two-dimensional and periodic in the horizontal direction. The numerical scheme used is a combined Galerkin and finite-difference method, and appropriate boundary conditions are applied at the interface. Results have been obtained for depth ratios of 0, 0.1, 0.2, 0.5, and 1.0, where this ratio is defined as the ratio of the thickness of the fluid layer to that of the porous layer. For the depth ratio of 0.1, the convection is dominated by the porous layer, similar to the situation at onset, even though the Rayleigh number for the fluid layer is well into the supercritical regime.

Chen, Falin↗