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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Scale–Resolving Simulations of a Fundamental Trailing–Edge Cooling Slot Using a Discontinuous–Galerkin Spectral–Element Method

The accurate prediction of turbulent mixing in high-pressure turbines that incorporate various airfoil surface-cooling strategies is becoming increasing critical to the design of modern gas turbine engines where the quest for improved efficiency is driving compressor overall pressure ratios and turbine inlet temperatures to much higher levels than ever before. In the present paper, a recently developed computational capability for accurate and efficient scale–resolving simulations of turbomachinery is extended to study the turbulent mixing mechanism of a simplified abstraction of an airfoil trailing-edge cooling slot - a plane wall jet with finite lip thickness discharging into an ambient flow. The computational capability is based on an entropy– stable, discontinuous–Galerkin approach that extends to arbitrarily high orders of spatial and temporal accuracy. The numerical results show that the present simulations capture the trends observed in the experiments. Discrepancies between the simulations and experiments are believed to be due to differences in the inflow profiles and tunnel side–wall effects. The thick lip configuration leads to a thicker wake and higher unsteadiness in the wall jet compared to the thin lip. A detailed comparison of the turbulent flowfields is presented to highlight differences arising due to lip thickness variations.

Scale Resolving↗
DOE OSTI2022

Pit Stability Predictions of Additively Manufactured SS316 Surfaces Using Finite Element Analysis

Stainless steels are susceptible to localized forms of corrosion attack, such as pitting. The size and lifetime of a nucleated pit can vary, depending on a critical potential or current density criterion, which determines if the pit repassivates or continues growing. This work uses finite element method (FEM) modeling to compare the critical pit radii predicted by thermodynamic and kinetic repassivation criteria. Experimental electrochemical boundary conditions are used to capture the active pit kinetics. Geometric and environmental parameters, such as the pit shape and size (analogous to additively manufactured lack-of-fusion pores), solution concentration, and water layer thickness were considered to assess their impact on the pit repassivation criterion. The critical pit radius (the transition point from stable growth to repassivation) predicted for a hemispherical pit was larger when using the repassivation potential (E rp ) criteria, as opposed to the current density criteria (pit stability product). Including both the pit stability product and E rp into its calculations, the analytical maximum pit model predicted a critical radius two times more conservative than the FEA approach, under the conditions studied herein. The complex pits representing lack-of-fusion pores were shown to have minimal impact on the critical radius in atmospheric conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Reduced modeling of flexible structures for decentralized control

Based upon the modified finite element-transfer matrix method, this paper presents a technique for reduced modeling of flexible structures for decentralized control. The modeling decisions are carried out at (finite-) element level, and are dictated by control objectives. A simply supported beam with two sets of actuators and sensors (linear force actuator and linear position and velocity sensors) is considered for illustration. In this case, it is conjectured that the decentrally controlled closed loop system is guaranteed to be at least marginally stable.

Yousuff, A.↗
NASA NTRS1993

High-order cyclo-difference techniques: An alternative to finite differences

The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy.

Carpenter, Mark H.↗
NASA NTRS1990

Adaptive quadrilateral and triangular finite-element scheme for compressible flows

The development of an adaptive mesh refinement procedure for analyzing high-speed compressible flows using the finite-element method is described. This new adaptation procedure, which uses both quadrilateral and triangular elements, was implemented with two explicit finite-element algorithms - the two-step Taylor-Galerkin and the multistep Galerkin-Runge-Kutta schemes. A von Neumann stability analysis and a rotating 'cosine hill'problem demonstrate the instability of the Taylor-Galerkin scheme when coupled with the adaptation procedure. For the same adaptive refinement scheme, the Galerkin-Runge-Kutta procedure yields stable solutions within its explicit stability limit. The utility of this new adaptation procedure for the prediction of compressible flow features is illustrated for inviscid problems involving strong shock interactions at hypersonic speeds.

Ramakrishnan, R.↗
NASA NTRS1985

Optimal design of a thermally stable composite optical bench

The Lidar Atmospheric Sensing Experiment will be performed aboard an ER-2 aircraft; the lidar system used will be mounted on a lightweight, thermally stable graphite/epoxy optical bench whose design is presently subjected to analytical study and experimental validation. Attention is given to analytical methods for the selection of such expected laminate properties as the thermal expansion coefficient, the apparent in-plane moduli, and ultimate strength. For a symmetric laminate in which one of the lamina angles remains variable, an optimal lamina angle is selected to produce a design laminate with a near-zero coefficient of thermal expansion. Finite elements are used to model the structural concept of the design, with a view to the optical bench's thermal structural response as well as the determination of the degree of success in meeting the experiment's alignment tolerances.

Gray, C. E., Jr.↗
NASA NTRS1992

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallelize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet -- and perhaps will never -- be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
NASA NTRS1973

Free vibration analysis of spinning structural systems.

This article presents an efficient digital computer procedure, along with the complete listing of the associated computer program, which may be conveniently utilized for the accurate solution of a wide range of practical eigenvalue problems. Important applications of the present work are envisaged in the natural frequency analysis of spinning structures discretized by the finite element technique, and in the determination of transfer functions associated with the dynamic blocks of control systems of spacecraft utilizing gas jets or reaction wheels for attitude control, as well as of spin-stabilized and dual-spin-stabilized satellites. The validity of the Sturm sequence property is first established for the related matrix formulation involving Hermitian and real symmetric, positive-definite matrices, both being usually of highly banded configuration. A numerically stable algorithm based on the Sturm sequence method is then developed which fully exploits the banded form of the associated matrices.

Gupta, K. K.↗
NASA NTRS1992

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

Aeroelastic stability analyses were performed to insure structural integrity of two counterrotating propfan blade designs for a NAVY/Air Force/NASA cruise missile model wind tunnel test. This analysis predicted if the propfan designs would be flutter free at the operating conditions of the wind tunnel test. Calculated stability results are presented for the two blade designs with rotational speed and Mach number as the parameters. A aeroelastic analysis code ASTROP2 (Aeroelastic Stability and Response of Propulsion Systems - 2 Dimensional Analysis), developed at LeRC, was used in this project. The aeroelastic analysis is a modal method and uses the combination of a finite element structural model and two dimensional steady and unsteady cascade aerodynamic models. This code was developed to analyze single rotation propfans but was modified and applied to counterrotating propfans for the present work. Modifications were made to transform the geometry and rotation of the aft rotor to the same reference frame as the forward rotor, to input a non-uniform inflow into the rotor being analyzed, and to automatically converge to the least stable aeroelastic mode.

Mahajan, Aparajit J.↗
NASA NTRS1993

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics, part 2

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗
NASA NTRS2004

Fracture Analysis of Semi-Elliptical Surface Cracks in Ductile Materials

Accurate life assessment of structural components may require advanced life prediction criteria and methodologies. Structural components often exhibit several different types of defects, among the most prevalent being surface cracks. A semi-elliptical surface crack subjected to monotonic loading will exhibit stable crack growth until the crack has reached a critical size, at which the crack loses stability and fracture ensues (Newman, 2000). The shape and geometry of the flaw are among the most influential factors. When considering simpler crack configurations, such as a through-the-thickness crack, a three-dimensional (3D) geometry may be modeled under the approximation of two-dimensional (2D) plane stress or plane strain. The more complex surface crack is typically modeled numerically with the Finite Element Method (FEM). A semi-elliptical surface crack is illustrated in Figure 1-1.

Daniewicz, S. R.↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
NASA NTRS1982

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗
NASA NTRS1977

Stability of numerical integration techniques for transient rotor dynamics

A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.

Kascak, A. F.↗
NASA NTRS1979

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗
NASA NTRSDate not supplied

Solar Cruiser TRAC Boom Development

The TRAC boom studied is a 30 meter, high-aspect ratio and highly nonlinear structural element planned for the NASA Solar Cruiser spacecraft where it would serve as the solar sail’s skeletal system designed to deploy the packaged reflective sail membrane as well as support operational sail tension loads for flight. Key design drivers for this boom are understanding and characterizing load-deformation behavior and quantifying buckling performance–under flight conditions (zero gravity). Knowing that the full-scale solar sail system cannot be feasibly tested for flight conditions under Earth’s gravity, a building block approach is methodically followed to verify finite element analysis (FEA) predictions and to evolve the finite element model (FEM) accordingly, beginning with the lowest level complexity to match global elastic behavior to the highest level complexity to bound theoretical responses under conservative expected flight conditions. Although the testing campaign for the TRAC boom has not developed yet to produce a full-scale test data point that loads the boom to failure, the fidelity of simulations have evolved with success and warrant further interest and work on this high-aspect ratio, but highly stable, deployable structural boom.

Lee Nguyen↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗