SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗

DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

DOE OSTI2024

Multiphysics Co-Optimization Design and Analysis of Double-Side Cooled Silicon Carbide-Based Power Module: Preprint

With the rapid growth of Electric Vehicles (EVs) and Hybrid Electric Vehicles (HEVs), much more rigorous design targets have been set for automotive power electronics, including high power density, high reliability, and low cost. Novel power module and inverter technologies based on wide bandgap (WEG) semiconductors have been developed to meet these design targets, while providing optimal power semiconductor operating temperature and promising thermomechanical performance. Compared with conventional cooling techniques which are normally applied only on one side of power module, double-side cooling approach is now believed to be the solution to enable high power density and low thermal resistance of WEG semiconductor-based power electronics. In this work, we develop a three-phase power module that is double-sided cooled using dielectric fluid jet impingement. In each phase, four silicon carbide (SiC) power semiconductors are bonded to copper busbars without electrical insulation layers. A finite element analysis (FEA) model is created for thermal and thermomechanical analysis. Based on FEA modeling results, we select particular dimensions for a parametric study to optimize thermal and mechanical performance. Using a multi-objective genetic algorithm (MOGA)-based optimization method, we have minimized the maximum junction temperature and thermal stresses within the power module. The multiphysics co-optimization approach has enabled an efficient design process of power modules with greatly reduced computational cost, as compared to conventional processes that rely on exhaustive numerical simulations and iterations.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Uncertainty Quantification of Fatigue Behavior of Rough AM Surfaces and Microstructures to Enable Hydrogen Gas Turbine Combustion

Modification of fossil-fueled industrial gas turbines to accept no/low carbon fuels (Hydrogen, H2/natural gas blends) is a significant undertaking. Successful deployment of this technology sits at the intersection of three design criteria (1) new functional fuel injectors that can burn these fuels, (2) manufacturability to meet cost and time-to-market targets, and (3) durability in the harsh environment of an operating turbine. Additive manufacturing (AM) provides accelerated product development. However, uncertainty remains around the durability of parts with rough AM surfaces. A fully experimental approach towards quantifying fatigue performance of rough AM microstructures is costly and laborious. Instead, Solar Turbines Incorporated (Solar) proposes the use of a crystal plasticity finite element (CPFE) model to quantify the factors that drive AM surface fatigue behavior. Solar will use the CPFE results, along with targeted experimental data, to train a computationally efficient surrogate model that can be incorporated into existing turbine part lifing methods.

08 HYDROGEN↗

DOE OSTI2020

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Role of bath-induced many-body interactions in the dissipative phases of the Su-Schrieffer-Heeger model

The Su-Schrieffer-Heeger chain is a prototype example of a symmetry-protected topological insulator. Coupling it nonperturbatively to local thermal environments, either through the intercell or the intracell fermion tunneling elements, modifies the topological window. To understand this effect, we employ the recently developed reaction-coordinate polaron transform (RCPT) method, which allows treating system-bath interactions at arbitrary strengths. The effective system Hamiltonian, which is obtained via the RCPT, exposes the impact of the baths on the SSH chain through renormalization of tunneling elements and the generation of many-body interaction terms. By performing exact diagonalization and computing the ensemble geometric phase, a topological invariant, which is applicable even to systems at finite temperature, we distinguish the trivial band insulator (BI) from the topological insulator (TI) phases. Furthermore, through the RCPT mapping, we are able to pinpoint the main mechanism behind the extension of the parameter space for the TI or the BI phases (depending on the coupling scheme, intracell or intercell), which is the bath-induced, dimerized, many-body interaction. In conclusion, we also study the effect of on-site staggered potentials on the SSH phase diagram and discuss extensions of our method to higher dimensions.

1-dimensional systems↗

DOE OSTI2022

Quantifying modeling uncertainty in simplified beam models for building response prediction

The use of simple models for response prediction of building structures is preferred in earthquake engineering for risk evaluations at regional scales, as they make computational studies more feasible. The primary impediment in their gainful use presently is the lack of viable methods for quantifying (and reducing upon) the modeling errors/uncertainties they bear. This study presents a Bayesian calibration method wherein the modeling error is embedded into the parameters of the model. Here, the method is specifically described for coupled shear-flexural beam models here, but it can be applied to any parametric surrogate model. The major benefit the method offers is the ability to consider the modeling uncertainty in the forward prediction of any degree-of-freedom or composite response regardless of the data used in calibration. The method is extensively verified using two synthetic examples. In the first example, the beam model is calibrated to represent a similar beam model but with enforced modeling errors. In the second example, the beam model is used to represent the detailed finite element model of a 52-story building. Both examples show the capability of the proposed solution to provide realistic uncertainty estimation around the mean prediction.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

BISON validation of FeCrAl cladding mechanical failure during simulated reactivity-initiated accident conditions*

Here, a reactivity-initiated accident (RIA) is a postulated design basis accident in light water reactors (LWRs) in which a rapid reactivity insertion induces a fission rate increase and a fuel pellet temperature rise. During RIA, the fuel pellet thermally expands and may cause pellet-cladding mechanical interaction (PCMI). Separate effects PCMI tests were performed on C26M FeCrAl cladding tube samples, introducing biaxial stress via a well understood modified burst test (MBT) system. A high-speed camera in the MBT system captured the projections, covering a 360° view of the cladding deformation and enabling a digital image correlation (DIC) method to quantify the surface strains with high fidelity. Representative hot zero-power RIA mechanical loading conditions were applied to the sample, and the test duration ranged from 20 to 500 ms at an average temperature of 573 K. BISON finite element–based fuel performance code modeling was performed against the high-fidelity DIC data produced from the MBTs. Validation calculations were conducted with 2D models and systematically compared with test data of the burst time, burst pressure, burst hoop strains, and hoop strain rates. Based on the behaviors from the separate effects test, BISON calculations satisfactorily predicted the cladding deformation behaviors. Sensitivity analysis was conducted to identify highly influential mechanical properties responsible for the cladding failure behavior during the MBT experiments. The results highlight the significance of the cladding's mechanical strength in governing cladding strain, followed by the significance of the mechanical interactions between the pellet and the cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Develop a Fast Analysis Solver for Welding Sequence Optimization

During the shipbuilding manufacturing process, materials are exposed to significant stresses, as induced both thermally and mechanically, that alter the intended design and significantly affect the production schedule, labor hours (fitting, welding, rework, etc.), and material structural performance. The type and magnitude of deformation of a given structure depends on many factors such as the material, thickness and quality of components, the process heat input, preheat and inter-pass temperatures, type and size of welds, welding sequence and direction, location, sequence, and degree of fixturing. Numerical simulations using finite element analysis (FEA) have long been used to analyze welding-induced structural distortion. For large assemblies, transient thermal elastic-plastic analysis (TEPA) can take days or weeks to run, and optimization of welding sequence is not feasible. Simplified analysis methods were developed to reduce computational time. However, it is challenging to use these techniques to fully optimize welding sequencing because of their applied simplifications in modeling weld details. A fast analysis solver that could be used by the shipbuilding industry is being developed for optimizing welding sequences by taking full advantage of modern GPU-based HPC hardware and incorporating patented acceleration schemes. The accelerated processing factors are up to 2200 times greater for large, multi-pass welded structures.

Yang, Yu-Ping↗

DOE OSTI2023

Design and Scoping Tests on Alloy 617 Using Notched Specimen Geometry to Validate Methods for Multiaxial Stress Relaxation

In FY 2022, a development effort was initiated at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) to examine the multiaxial stress-relaxation behavior and multiaxial stress-state effect on the creep-fatigue (CF) performance for Alloy 617 at elevated temperatures. This effort supported the development of the design rules in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), Section III, Division 5. In this work, two types of the notch specimen geometries were designed. An inelastic constitutive model was used to investigate the stress triaxiality and elastic follow-up effects on the notch specimens under CF deformation in the finite element (FE) simulations. The numerical study demonstrated that the stress triaxiality and elastic follow-up caused by the notches both play significant roles in the stress-relaxation behavior. In addition, Alloy 617 CF experiments were designed, and the testing on the notch specimens is ongoing. CF test failure data will be generated on the specimens with notches under various conditions in FY 2023, and the results will be used to validate methods for accounting for the multiaxial stress relaxation effect in the design code.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Development and testing of a solar flux heating freeze recovery system for molten salt parabolic troughs

A major obstacle preventing commercialization of molten salt parabolic trough plants is freeze protection and recovery of the solar field. One of the most common nitrate salt compositions under consideration allows for solar field outlet temperatures up to 565°C and begins solidifying above 240°C. Using molten salt directly in the solar field can significantly reduce thermal energy storage costs and the higher temperature compared to oil heat transfer fluids can power a more efficient power cycle, but the very high freezing point is cause for concern. This work builds on a previous modeling effort investigating different methods for melting Solar Salt frozen in the solar field, which suggested the viability of a novel solar heating method to lower costs of the freeze protection system. In this work higher fidelity models are developed to better investigate the melting process with a more detailed 3D geometry which includes the insulated bellows. The thermal-fluid model of the melting process is updated to include temperature dependent density and data for solid phase properties. In addition, a finite element model is developed to resolve thermal stresses at the point with highest thermal gradients. Results with the new model confirm the viability of using controllable solar flux heating to thaw salt frozen in the solar field without damaging the receiver; however, the presence of non-illuminated sections significantly slows down the melting process. Solar heating simulations with pauses off-sun of 120 and 60 s require 18.5 and 11.9 h to melt salt from a night-time temperature of 10°C with 1000 W m2 DNI, which results in maximum thermal stresses of 48.2 and 53.5 MPa, respectively. Adding 150 W m−1 to the solar heating simulation with a 60 s pause off-sun decreases both the melt time and thermal stresses due to the more uniform heating, suggesting a combined heating method may be the best option.

14 SOLAR ENERGY↗

DOE OSTI2020

Dynamic Tensile Characterization of Thin-Sheet Brittle Metallic Materials

Refractory metals are favorable materials in applications where high strength and ductility are needed at elevated temperatures. In some cases, operating temperatures may be near the melting point of the material. However, as temperature drops, refractory metals typically undergo a significant mechanical response change - ductile-to-brittle transition. These materials may be subjected to high strain rate loading at an ambient temperature state, such as an impact or crash. Knowledge of the high rate material properties are essential for design as well as simulation of impact events. The high rate stress-strain behavior of brittle metallic materials at ambient temperature is rarely studied because of experimental challenges, particularly when failure is involved. Failure typically occurs within the non-gage section of the material, which invalidates any collected stress-strain information. Here, a method to determine a specimen geometry which will produce failures in the gage section is presented. Pure tungsten in thin-sheet form was used as a trial material to select a specimen geometry for high rate Kolsky tension bar experiments. A finite element simulation was conducted to derive a strain correction for more accurate results. The room temperature stress-strain behavior of pure tungsten at a strain rate of 24 s –1 is presented. The outcome of this experimental technique can be applied to other brittle materials for dynamic tensile characterization.

36 MATERIALS SCIENCE↗

DOE OSTI2023

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Inverse Design of Plasma Metamaterial Devices for Optical Computing

We apply inverse-design methods to produce two-dimensional plasma metamaterial (PMM) devices. Backpropagated finite-difference frequency-domain (FDFD) simulations are used to design waveguides and demultiplexers operating under both TE and TM modes. Demultiplexing and waveguiding are demonstrated for devices composed of plasma elements with reasonable plasma densities of approximately 7 GHz, allowing for future in situ training and experimental realization of these designs. We also explore the possible applicability of PMMs to nonlinear boolean operations for use in optical computing. Functionally complete logical connectives (or and and) are achieved in the TM mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a densitymatrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. They can enhance individual low-lying transitions, however, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these effects are related to changes in the total nucleon density as the number of neutrons increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Design of High-Power Polyphase PCB Coil Systems for Wireless Power Transfer

Printed circuit board (PCB) coils have been proposed prior for implementation as inductive wireless charging coils to minimize size and cost. Utilization of PCBs can allow for a reduced cost, improved manufacturability, and a wide range of geometric customization options. To circumvent material limitations on insulation and thermal performance, parallel paths can be implemented to divide the current per path accordingly. Within this paper, an unconventional high-power PCB coil is designed employing all possible techniques for wiring with axial and radial parallel paths with equivalent transposition to minimize circulating currents. Design studies are simulated in 3D finite element analysis (FEA) to evaluate imbalance between phases with and without transposition. Two experimental prototype coils were fabricated with measurements for self-inductance and mutual inductance between phases. These measurements were validated to be sufficiently consistent with FEA results. Additionally, a method is proposed for a two-step optimization of coupling coefficient and coil losses.

Lewis, Donovin D.↗

DOE OSTI2020

A reduced-order shifted boundary method for parametrized incompressible Navier–Stokes equations

Here, we investigate a projection-based reduced order model of the steady incompressible Navier–Stokes equations for moderate Reynolds numbers. In particular, we construct an “embedded” reduced basis space, by applying proper orthogonal decomposition to the Shifted Boundary Method, a high-fidelity embedded method recently developed. We focus on the geometrical parametrization through level-set geometries, using a fixed Cartesian background geometry and the associated mesh. This approach avoids both remeshing and the development of a reference domain formulation, as typically done in fitted mesh finite element formulations. Two-dimensional computational examples for one and three parameter dimensions are presented to validate the convergence and the efficacy of the proposed approach.

"embedded" basis↗