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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Sinc-Galerkin estimation of diffusivity in parabolic problems

A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise.

Smith, Ralph C.↗

NASA NTRS1992

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗

NASA NTRS1993

Double-diffusive fingering convection in a porous medium

The characteristics of nonlinear two-dimensional horizontally periodic double-diffusive fingering convection in a saturated porous medium is investigated, using the Darcy equation including Brinkman and Forchheimer terms for the momentum equation. To solve the equations and the corresponding initial and boundary conditions, a Galerkin method is applied in the horizontal direction, and a hybrid finite difference method is used in the vertical direction. The developed computer code was used to compute the thermal convection case, and the results were found to be in good agreement with existing results.

Chen, Falin↗

NASA NTRS1993

Simulation of tail buffet using delta wing-vertical tail configuration

Computational simulation of the vertical tail buffet problem is accomplished using a delta wing-vertical tail configuration. Flow conditions are selected such that the wing primary-vortex cores experience vortex breakdown and the resulting flow interacts with the vertical tail. This multidisciplinary problem is solved successively using three sets of equations for the fluid flow, aeroelastic deflections and grid displacements. For the fluid dynamics part, the unsteady, compressible, full Navier-Stokes equations are solved accurately in time using an implicit, upwind, flux-difference splitting, finite-volume scheme. For the aeroelastic part, the aeroelastic equation for bending vibrations is solved accurately in time using the Galerkin method and the four-stage Runge-Kutta scheme. The grid for the fluid dynamics computations is updated every few time steps using a third set of interpolation equations. The computational application includes a delta wing of aspect ratio 1 and a rectangular vertical tail of aspect ratio 2, which is placed at 0.5 root-chord length downstream of the wing trailing edge. The wing angle of attack is 35 deg and the flow Mach number and Reynolds number are 0.4 and 10,000, respectively.

Kandil, Osama A.↗

NASA NTRS1950

Critical stress of ring-stiffened cylinders in torsion

A chart in terms of nondimensional parameters is presented for the theoretical critical stress in torsion of simply supported cylinders stiffened by identical equally spaced rings of zero torsional stiffness. The results are obtained by solving the equation of equilibrium by means of the Galerkin method. Comparison of the theoretical results with experimental results indicates that ring-stiffened cylinders buckle, on the average, at a buckling stress about 15 percent below the theoretical buckling stress. (author)

Stein, Manuel↗

NASA NTRS1993

Wavelet multiresolution analyses adapted for the fast solution of boundary value ordinary differential equations

We present ideas on how to use wavelets in the solution of boundary value ordinary differential equations. Rather than using classical wavelets, we adapt their construction so that they become (bi)orthogonal with respect to the inner product defined by the operator. The stiffness matrix in a Galerkin method then becomes diagonal and can thus be trivially inverted. We show how one can construct an O(N) algorithm for various constant and variable coefficient operators.

Jawerth, Bjoern↗

NASA NTRS1995

Scattering Of Nonplanar Acoustic Waves

Report presents theoretical study of scattering of nonplanar acoustic waves by rigid bodies. Study performed as part of effort to develop means of predicting scattering, from aircraft fuselages, of noise made by rotating blades. Basic approach was to model acoustic scattering by use of boundary integral equation to solve equation by the Galerkin method.

Gillman, Judith M.↗

NASA NTRS1995

An hp-adaptivity and error estimation for hyperbolic conservation laws

This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Bey, Kim S.↗

NASA NTRS1998

Detection of Non-Symmetrical Damage in Smart Plate-Like Structures

A two-dimensional model for in-plane vibrations of a cantilever plate with a non-symmetrical damage is used in the context of defect identification in materials with piezoelectric ceramic patches bonded to their surface. These patches can act both as actuators and sensors in a self-analyzing fashion, which is a characteristic of smart materials. A Galerkin method is used to approximate the dynamic response of these structures. The natural frequency shifts due to the damage are estimated numerically and compared to experimental data obtained from tests on cantilever aluminum plate-like structures damaged at different locations with defects of different depths. The damage location and extent are determined by an enhanced least square identification method. Efficacy of the frequency shift based algorithms is demonstrated using experimental data.

Blanks, H. T.↗

NASA NTRS1998

Effect of Apex Flap Deflection on Vertical Tail Buffeting

A computational study of the effect of vortex breakdown location on vertical tail buffeting is conducted. The position of the breakdown is modified by employing an apex flap deflected by an experimentally determined optimal angle. The delayed breakdown flow and buffeting response is then compared to the nominal undeflected case. This multidisciplinary problem is solved sequentially for the fluid flow, the elastic tail deformations and the grid displacements. The fluid flow is simulated by time accurately solving the unsteady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, upwind, flux-difference splitting finite volume scheme. The elastic vibrations of the tails are modeled by uncoupled bending and torsion beam equations. These equations are solved accurately in time using the Galerkin method and a five-stage Runge-Kutta-Verner scheme. The grid for the fluid dynamics calculations is continuously deformed using interpolation functions to disperse the displacements smoothly throughout the computational domain. An angle-of-attack of 35 deg.is chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on tile vertical tails. The dimensions and material properties of the vertical tails are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. Results are presented for a baseline uncontrolled buffeting case and a delayed breakdown case in which the apex flap has been deflected 15 deg. The flap was found to be very effective in delaying the breakdown, increasing the location from 50%c to 94%c, which resulted in a 6% increase in lift coefficient and pitching moment. However, the integrated buffet loads and tip responses were roughly equivalent for the two cases.

Massey, Steven J.↗

NASA NTRS2000

Discontinuous dual-primal mixed finite elements for elliptic problems

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Bottasso, Carlo L.↗

NASA NTRS1997

Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology

The compact form of the discontinuous Galerkin method allows for a detailed local analysis of the method in the neighborhood of the shock for a non-linear model problem. Insight gained from the analysis leads to new flux formulas that are stable and that preserve the compactness of the method. Although developed for a model equation, the flux formulas are applicable to systems such as the Euler equations. This article presents the analysis for methods with a degree up to 5. The analysis is accompanied by supporting numerical experiments using Burgers' equation and the Euler equations.

Atkins, H. L.↗

NASA NTRS2006

Low-Velocity Impact Response of Sandwich Beams with Functionally Graded Core

The problem of low-speed impact of a one-dimensional sandwich panel by a rigid cylindrical projectile is considered. The core of the sandwich panel is functionally graded such that the density, and hence its stiffness, vary through the thickness. The problem is a combination of static contact problem and dynamic response of the sandwich panel obtained via a simple nonlinear spring-mass model (quasi-static approximation). The variation of core Young's modulus is represented by a polynomial in the thickness coordinate, but the Poisson's ratio is kept constant. The two-dimensional elasticity equations for the plane sandwich structure are solved using a combination of Fourier series and Galerkin method. The contact problem is solved using the assumed contact stress distribution method. For the impact problem we used a simple dynamic model based on quasi-static behavior of the panel - the sandwich beam was modeled as a combination of two springs, a linear spring to account for the global deflection and a nonlinear spring to represent the local indentation effects. Results indicate that the contact stiffness of thc beam with graded core Increases causing the contact stresses and other stress components in the vicinity of contact to increase. However, the values of maximum strains corresponding to the maximum impact load arc reduced considerably due to grading of thc core properties. For a better comparison, the thickness of the functionally graded cores was chosen such that the flexural stiffness was equal to that of a beam with homogeneous core. The results indicate that functionally graded cores can be used effectively to mitigate or completely prevent impact damage in sandwich composites.

Apetre, N. A.↗

NASA NTRSDate not supplied

Aeroelastic Trim Drag Optimization of Mach 0.8 TransonicTruss-Braced Wing Aircraft with Variable Camber ContinuousTrailing Edge Flap

This paper presents an aeroelastic trim drag optimization study of the Mach 0.8 Tran-sonic Truss-Braced Wing (TTBW) aircraft with the Variable Camber Continuous TrailingEdge Flap (VCCTEF). An aero-structural analysis solver VSPAERO with transonic small dis-turbance, integral boundary-layer, and wing-strut interference corrections coupled to modeshapes computed by NASTRAN using the Galerkin method is developed to provide a rapidaircraft aeroelastic performance evaluation. Aeroelastic trim drag optimization studies areconducted for a VCCTEF configuration with 6-spanwise sections. Three different flight condi-tions corresponding to Mach 0.8 are selected for the aeroelastic trim drag optimization at thedesign and off-design cruise lift coefficients. The preliminary optimization results show thatthe TTBW aircraft with the optimized VCCTEF deflection achieves a drag reduction of about 9.2 counts, 9.6 counts, and 12.3 counts corresponding to the lift coefficients 0.661, 0.695, and0.729, respectively. When accounting for the actuator weight penalty, the corresponding dragreductions are 1.77%, 1.82%, 2.41%. A high-fidelity CFD solver FUN3D is used to validatethe aeroelastic trim drag optimization.

Aeroelastic trim drag↗

NASA NTRSDate not supplied

Aeroelastic Trim Drag Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft using High-Lift Devices and Control Surfaces

This paper presents an aeroelastic trim drag optimization study of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using the High-lift devices and control surfaces. An aero-structural analysis solver VSPAERO with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections coupled to mode shapes computed by NASTRAN using the Galerkin method is developed to provide a rapid aircraft aeroelastic performance evaluation. Three different flight conditions corresponding to Mach 0.8 are selected for the aeroelastic trim drag optimization at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized deflection of the high-lift devices and control surfaces achieves a drag reduction of about 8.2 counts, 6.5 counts, and 9.7 counts corresponding to the lift coefficients 0.661, 0.695, and 0.729, respectively. A high-fidelity CFD solver FUN3D is used to verify the aeroelastic trim drag optimization.

TTBW↗

DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗

DOE OSTI2025

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗

NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗