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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Influence of crack history on the stable tearing behavior of a thin-sheet material with multiple cracks

Fracture tests were conducted on 2.3mm thick, 305mm wide sheets of 2024-T3 aluminum alloy with from one to five collinear cracks. The cracks were introduced (crack history) into the specimens by three methods: saw cutting, fatigue precracking at a low stress range, and fatigue precracking at a high stress range. For the single crack tests, the initial crack history influenced the stress required for the onset of stable crack growth and the first 10mm of crack growth. The effect on failure stress was about 4 percent or less. For the multiple crack tests, the initial crack history was shown to cause differences of more than 20 percent in the link-up stress and 13 percent in failure stress. An elastic-plastic finite element analysis employing the CTOA fracture criterion was used to predict the fracture behavior of the single and multiple crack tests. The numerical predictions were within 7 percent of the observed link-up and failure stress in all the tests.

Dawicke, D. S.↗
NASA NTRS1994

Stable tearing behavior of a thin-sheet material with multiple cracks

Fracture tests were conducted on 2.3mm thick, 305mm wide sheets of 2024-T3 aluminum alloy with 1-5 collinear cracks. The cracks were introduced (crack history) into the specimens by three methods: (1) saw cutting; (2) fatigue precracking at a low stress range; and (3) fatigue precracking at a high stress range. For the single crack tests, the initial crack history influenced the stress required for the onset of stable crack growth and the first 10mm of crack growth. The effect on failure stress was about 4 percent or less. For the multiple crack tests, the initial crack history was shown to cause differences of more than 20 percent in the link-up stress and 13 percent in failure stress. An elastic-plastic finite element analysis employing the Crack Tip Opening Angle (CTOA) fracture criterion was used to predict the fracture behavior of the single and multiple crack tests. The numerical predictions were within 7 percent of the observed link-up and failure stress in all the tests.

Dawicke, D. S.↗
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 is 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 are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been 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.↗
DOE OSTI2020

A Mixed Discontinuous Galerkin Method for Linear Elasticity with Strongly Imposed Symmetry

In this paper, we study a mixed discontinuous Galerkin (MDG) method to solve linear elasticity problem with arbitrary order discontinuous finite element spaces in d -dimension ( \(d=2,3\) ). This method uses polynomials of degree \(k+1\) for the stress and of degree k for the displacement ( \(k\ge 0\) ). The mixed DG scheme is proved to be well-posed under proper norms. Specifically, we prove that, for any \(k \ge 0\) , the \(H(\mathrm{div})\) -like error estimate for the stress and \(L^2\) error estimate for the displacement are optimal. We further establish the optimal \(L^2\) error estimate for the stress provided that the \({\mathcal {P}}_{k+2}-{\mathcal {P}}_{k+1}^{-1}\) Stokes pair is stable and \(k \ge d\) . Finally, we also provide numerical results of MDG showing that the orders of convergence are actually sharp.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2022

In-situ, nanoscale fracture toughness measurements for improved mechanical interfaces

In this project, we demonstrated stable nanoscale fracture in single-crystal silicon using an in-situ wedge-loaded double cantilever beam (DCB) specimen. The fracture toughness KIC was calculated directly from instrumented measurement of force and displacement via finite element analysis with frictional corrections. Measurements on multiple test specimens were used to show K IC = 0.72 ± 0.07 MPa m 1/2 on {111} planes and observe the crack-growth resistance curve in <500 nm increments. The exquisite stability of crack growth, instrumented measurement of material response, and direct visual access to observe nanoscale fracture processes in an ideally brittle material differentiate this approach from prior DCB methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Harriott, G. M.↗
NASA NTRS1996

A Least-Squares Finite Element Method for Electromagnetic Scattering Problems

The least-squares finite element method (LSFEM) is applied to electromagnetic scattering and radar cross section (RCS) calculations. In contrast to most existing numerical approaches, in which divergence-free constraints are omitted, the LSFF-M directly incorporates two divergence equations in the discretization process. The importance of including the divergence equations is demonstrated by showing that otherwise spurious solutions with large divergence occur near the scatterers. The LSFEM is based on unstructured grids and possesses full flexibility in handling complex geometry and local refinement Moreover, the LSFEM does not require any special handling, such as upwinding, staggered grids, artificial dissipation, flux-differencing, etc. Implicit time discretization is used and the scheme is unconditionally stable. By using a matrix-free iterative method, the computational cost and memory requirement for the present scheme is competitive with other approaches. The accuracy of the LSFEM is verified by several benchmark test problems.

Wu, Jie↗
DOE OSTI2025

Feasibility Assessment of Using Electrical Impedance Tomography for Damage Localization in Graphite Microreactor Components

Microreactors have great potential to decrease capital costs and construction timelines, reducing the barriers to implementing advanced nuclear reactor technologies. However, the lower power output of these microreactors introduces economic challenges if the operation and maintenance costs cannot be reduced sufficiently. Many microreactor concepts use graphite materials for in-core moderator and structural purposes, which will require periodic inspection or, ideally, in situ structural health monitoring. This work describes an initial evaluation of the feasibility of leveraging graphite’s semiconducting properties to perform electrical impedance tomography (EIT) for defect localization. First, a study was performed to identify the best methods for bonding electrical contacts to minimize the effects of contact resistance that interfere with impedance measurements of the graphite. After determining the best bonding approach, electrodes were attached to multiple graphite components with varying geometries (e.g., simple block and more representative microreactor hexagonal block). In parallel, finite element analysis approaches were developed and implemented to analyze the impact of defects and contact resistances and to inform an iterative inverse model for reconstructing the conductivity distribution. The results of this study show that defect localization in graphite components is possible using EIT if a sufficient number of electrodes (to improve spatial resolution) can be bonded using a technique with very low contact resistance and those contact resistances remain stable during reactor operation. If the reactor core and vessel design do not allow electrodes to be bonded during reactor operation, then it may be possible to detach/reattach electrodes between reactor operation and outages if a suitable mechanical connection that offers repeatable contact resistance can be identified. Alternatively, future work could focus on identifying contacts that can survive the harsh reactor operating environment so that contacts do not need to be removed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗
NASA NTRS1992

Numerical simulation of conservation laws

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

Chang, Sin-Chung↗
NASA NTRSDate not supplied

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Aeroelastic stability analysis of a high-energy turbine blade

The dynamic analysis for the SSME HPOTP first stage turbine blade is presented wherein the rotor aeroelastic stability is assessed. The method employs normal modes analysis to simulate the coupled blade/fluid system. A three-dimensional finite element model of the blade is used in conjunction with a two-dimensional linearized unsteady aerodynamic theory which accounts for steady aerodynamic loading effects. This unsteady aerodynamic model is applied in stacked axisymmetric strips along the airfoil span. The blade dynamic and aerodynamic behaviors are coupled within modal space by expressing the unsteady aerodynamic forces in the frequency domain. A complex eigenvalue problem is solved to determine the stability of the rotor assuming tuned blades. The present analysis indicates that the HPOTP rotor experiences very low aerodynamic damping in the first four vibrational modes. The edgewise mode was found to be dynamically unstable. This mode of the blade became stable when the effect of mechanical damping was considered.

Smith, Todd E.↗
NASA NTRS1988

Advanced composite materials for precision segmented reflectors

The objective in the NASA Precision Segmented Reflector (PSR) project is to develop new composite material concepts for highly stable and durable reflectors with precision surfaces. The project focuses on alternate material concepts such as the development of new low coefficient of thermal expansion resins as matrices for graphite fiber reinforced composites, quartz fiber reinforced epoxies, and graphite reinforced glass. Low residual stress fabrication methods will be developed. When coupon specimens of these new material concepts have demonstrated the required surface accuracies and resistance to thermal distortion and microcracking, reflector panels will be fabricated and tested in simulated space environments. An important part of the program is the analytical modeling of environmental stability of these new composite materials concepts through constitutive equation development, modeling of microdamage in the composite matrix, and prediction of long term stability (including viscoelasticity). These analyses include both closed form and finite element solutions at the micro and macro levels.

Stein, Bland A.↗