Search NASA⌕ Search

SEARCH · Search NASA

Results for “stable finite element 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS1986

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗
NASA NTRS1985

Structural and aeroelastic analysis of the SR-7L propfan

A structural and aeroelastic analysis of a large scale advanced turboprop rotor blade is presented. This 8-blade rotor is designed to operate at Mach 0.8 at an altitude of 35,000 ft. The blades are highly swept and twisted and of spar/shell construction. Due to the complexity of the blade geometry and its high performance, it is subjected to much higher loads and tends to be much less stable than conventional blades. Four specific analyses were conducted: (1) steady deflection; (2) natural frequencies and mode shapes; (3) steady stresses; and (4) aeroelastic stability. State-of-the-art methods were used to analyze the blades including a large deflection, finite element structural analysis, and an aeroelastic analysis including interblade aerodynamic coupling (cascade effects). The study found the blade to be structurally sound and aeroelastically stable. However, it clearly indicated that advanced turboprop blades are much less robust than conventional blades and must be analyzed and fabricated much more carefully in order to assure that they are structurally sound and aeroelastically stable.

Hirschbein, Murray↗
NASA NTRS1990

Finite-element schemes for extended integrations of atmospheric models

Of the two finite-element models presently used to investigate the effect of the conservation of integral invariants, by means of finite-element discretization schemes of the shallow-water equations, in order to serve as a paradigm of long-term atmospheric-model integrations, the first employs rectangular elements and conserves total energy, while the second uses triangular elements and a high-accuracy two-stage Numerov-Galerkin method. Attention is given to critical times for numerical nonlinear instability, as well as to the determination of the critical degree of dissipation entailed by the achievement of stable long-term integrations. Relative computational efficiency and accuracy comparisons are presented for the two finite-element schemes.

Steppeler, J.↗
NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A stable potential-based time-domain method for wideband elec- tromagnetic analysis

In previous research, the frequency-domain A-ϕ formulation has been validated using the finite element method for electromagnetic simulations of low-frequency and multi- scale problems, demonstrating excellent numerical accuracy, good matrix condition, and high computational efficiency. Time- domain simulations provide significant advantages for modeling wideband problems and are crucial for multiphysics applications. In this paper, the frequency-domain A-ϕ formulation is extended to the time domain. The central difference scheme is employed for temporal discretization to ensure both accuracy and stability. A numerical example is presented to demonstrate the capability of the proposed time-domain method in wideband electromagnetic analysis.

Mekonnen, Minyichil↗
NASA NTRS1986

Analysis of boundary conditions for SSME subsonic internal viscous flow analysis

A study was completed of mathematically proper boundary conditions for unique numerical solution of internal, viscous, subsonic flows in the space shuttle main engine. The study has concentrated on well posed considerations, with emphasis on computational efficiency and numerically stable boundary condition statements. The method of implementing the established boundary conditions is applicable to a wide variety of finite difference and finite element codes, as demonstrated.

Baker, A. J.↗
DOE OSTI2022

A novel post-processing method for progressive failure analysis of brittle composite compression

Finite element analysis of brittle materials in axial compression typically uses element deletion to allow continued global deformation post-element-failure. However, element deletion produces cyclic load-displacement curves that underestimate energy absorption and are not representative of a continuum system. Two key observations support the conclusion that results from an appropriately discretized model can be an adequate representation of a continuum system. Specifically, the frequency of the oscillations in the load-displacement curve is directly dependent upon element length in the loading direction, and the peak amplitudes of oscillations are mesh size independent. A method of post-processing the analysis results, by connecting the peak amplitudes of oscillations, is proposed and applied to a series of continuous carbon fiber composite crush tubes. The load-displacement curve, stable crushing load, and specific energy absorption of the post-processed results compare well to an experimental study of crush tubes with similar layups.

Materials Science↗
NASA NTRS1992

Aeroelastic stability analyses of two counter rotating propfan designs for a cruise missile model

A modal aeroelastic analysis combining structural and aerodynamic models is applied to counterrotating propfans to evaluate their structural integrity for wind-tunnel testing. The aeroelastic analysis code is an extension of the 2D analysis code called the Aeroelastic Stability and Response of Propulsion Systems. Rotational speed and freestream Mach number are the parameters for calculating the stability of the two blade designs with a modal method combining a finite-element structural model with 2D steady and unsteady cascade aerodynamic models. The model demonstrates convergence to the least stable aeroelastic mode, describes the effects of a nonuniform inflow, and permits the modification of geometry and rotation. The analysis shows that the propfan designs are suitable for the wind-tunnel test and confirms that the propfans should be flutter-free under the range of conditions of the testing.

Mahajan, Aparajit J.↗
NASA NTRS2001

Direct Calculation of the Scattering Amplitude Without Partial Wave Analysis

Two new developments in scattering theory are reported. We show, in a practical way, how one can calculate the full scattering amplitude without invoking a partial wave expansion. First, the integral expression for the scattering amplitude f(theta) is simplified by an analytic integration over the azimuthal angle. Second, the full scattering wavefunction which appears in the integral expression for f(theta) is obtained by solving the Schrodinger equation with the finite element method (FEM). As an example, we calculate electron scattering from the Hartree potential. With minimal computational effort, we obtain accurate and stable results for the scattering amplitude.

Shertzer, J.↗
DOE OSTI2021

Variational, stable, and self-consistent coupling of 3D electromagnetics to 1D transmission lines in the time domain

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. Here, we propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗
NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗
DOE OSTI2022

Solid-state synthesis of UV-plasmonic Cr2N nanoparticles

Materials that exhibit plasmonic response in the UV region can be advantageous for many applications, such as biological photodegradation, photocatalysis, disinfection, and bioimaging. Transition metal nitrides have recently emerged as chemically and thermally stable alternatives to metal-based plasmonic materials. However, most free-standing nitride nanostructures explored so far have plasmonic responses in the visible and near-IR regions. Herein, we report the synthesis of UV-plasmonic Cr2N nanoparticles using a solid-state nitridation reaction. The nanoparticles had an average diameter of 9 ± 5 nm and a positively charged surface that yields stable colloidal suspension. The particles were composed of a crystalline nitride core and an amorphous oxide/oxynitride shell whose thickness varied between 1 and 7 nm. Calculations performed using the finite element method predicted the localized surface plasmon resonance (LSPR) for these nanoparticles to be in the UV-C region (100–280 nm). While a distinctive LSPR peak could not be observed using absorbance measurements, low-loss electron energy loss spectroscopy showed the presence of surface plasmons between 80 and 250 nm (or ∼5 to 15 eV) and bulk plasmons centered around 50–62 nm (or ∼20 to 25 eV). Plasmonic coupling was also observed between the nanoparticles, resulting in resonances between 250 and 400 nm (or ∼2.5 to 5 eV).

Chemistry↗
NASA NTRS1982

Finite analytic numerical solution of heat transfer and flow past a square channel cavity

A numerical solution of flow and heat transfer characteristics is obtained by the finite analytic method for a two dimensional laminar channel flow over a two-dimensional square cavity. The finite analytic method utilizes the local analytic solution in a small element of the problem region to form the algebraic equation relating an interior nodal value with its surrounding nodal values. Stable and rapidly converged solutions were obtained for Reynolds numbers ranging to 1000 and Prandtl number to 10. Streamfunction, vorticity and temperature profiles are solved. Local and mean Nusselt number are given. It is found that the separation streamlines between the cavity and channel flow are concave into the cavity at low Reynolds number and convex at high Reynolds number (Re greater than 100) and for square cavity the mean Nusselt number may be approximately correlated with Peclet number as Nu(m) = 0.365 Pe exp 0.2.

Chen, C.-J.↗
NASA NTRS2019

Efficient Entropy Stable Gauss Collocation Methods

The construction of high-order entropy stable collocation schemes on quadrilateral and hexahedral elements has relied on the use of Gauss-Legendre-Lobatto collocation points and their equivalence with summation-by-parts (SBP) finite difference operators. In this work, we show how to efficiently generalize the construction of semidiscrete, entropy stable schemes on tensor product elements to Gauss points and generalized SBP operators. Numerical experiments suggest 8 that the use of Gauss points significantly improves accuracy on curved meshes.

Jesse Chan↗
NASA NTRS2001

The Use of Non-Standard Devices in Finite Element Analysis

A general mathematical description of the response behavior of thin-skin pneumatic envelopes and many other membrane and cable structures produces under-constrained systems that pose severe difficulties to analysis. These systems are mobile, and the general mathematical description exposes the mobility. Yet the response behavior of special under-constrained structures under special loadings can be accurately predicted using a constrained mathematical description. The static response behavior of systems that are infinitesimally mobile, such as a non-slack membrane subtended from a rigid or elastic boundary frame, can be easily analyzed using such general mathematical description as afforded by the non-linear, finite element method using an implicit solution scheme if the incremental uploading is guided through a suitable path. Similarly, if such structures are assembled with structural lack of fit that provides suitable self-stress, then dynamic response behavior can be predicted by the non-linear, finite element method and an implicit solution scheme. An explicit solution scheme is available for evolution problems. Such scheme can be used via the method of dynamic relaxation to obtain the solution to a static problem. In some sense, pneumatic envelopes and many other compliant structures can be said to have destiny under a specified loading system. What that means to the analyst is that what happens on the evolution path of the solution is irrelevant as long as equilibrium is achieved at destiny under full load and that the equilibrium is stable in the vicinity of that load. The purpose of this paper is to alert practitioners to the fact that non-standard procedures in finite element analysis are useful and can be legitimate although they burden their users with the requirement to use special caution. Some interesting findings that are useful to the US Scientific Balloon Program and that could not be obtained without non-standard techniques are presented.

Schur, Willi W.↗