Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element problems”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

An overlapping Domain Decomposition preconditioning method for monolithic solution of shear bands

Metals subjected to high strain rate impact often exhibit a sudden and profound drop in the material’s load bearing capability, a ductile failure phenomenon known as shear banding. Shear bands, characterized as material instabilities, are driven by shear heating and described as narrow regions that have sustained intense plastic deformation and high temperature rise. This coupled thermo-mechanical localization problem can be formulated as a nonlinear system with two balance equations, momentum and energy, and two constitutive laws for elasticity and plasticity. Here in our formulation, mixed finite elements are used to discretize the equations in space and an implicit finite difference scheme is used to advance the system in time, where at every time step, a Newton type method is used to solve the system monolithically. To that end, a block Jacobian matrix is formed analytically using Gâteaux derivatives. The resulting Jacobian is sparse, nonsymmetric and its sparsity pattern and eigenvalue content vary with the different stages of shear bands formation, consequently posing a significant challenge to iterative solvers. To address this issue, an overlapping Domain Decomposition preconditioner that takes into account the physics of shear bands and the domain in which it forms, is proposed. The key idea is to resolve the shear band domain accurately while only solving the remaining domain approximately, with selective updates when necessary. The preconditioner is formed on the basis of an additive Schwartz method that is applied to a Schur complement partition of the system. The proposed preconditioner is implemented in parallel and tested on three different benchmark problems as compared against off-the-shelf solvers. We study the effect of h- and k-refinement that are obtained from Isogeometric discretizations of the system, the overlapping strategy and its parallel performance. Excellent performance is demonstrated in serial and parallel on all benchmark examples.

42 ENGINEERING↗
DOE OSTI2017

Enhancement of flow measurements using fluid-dynamic constraints

Novel experimental modalities acquire spatially resolved velocity measurements for steady state and transient flows which are of interest for engineering and biological applications. One of the drawbacks of such high resolution velocity data is their susceptibility to measurement errors. In this paper, we propose a novel filtering strategy that allows enhancement of the noisy measurements to obtain reconstruction of smooth divergence free velocity and corresponding pressure fields which together approximately comply to a prescribed flow model. The main step in our approach consists of the appropriate use of the velocity measurements in the design of a linearized flow model which can be shown to be well-posed and consistent with the true velocity and pressure fields up to measurement and modeling errors. The reconstruction procedure is then formulated as an optimal control problem for this linearized flow model. The resulting filter has analyzable smoothing and approximation properties. We briefly discuss the discretization of the approach by finite element methods and comment on the efficient solution by iterative methods. The capability of the proposed filter to significantly reduce data noise is demonstrated by numerical tests including the application to experimental data. In addition, we compare with other methods like smoothing and solenoidal filtering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Solving Stochastic Inverse Problems for Property–Structure Linkages Using Data-Consistent Inversion and Machine Learning

Determining process–structure–property linkages is one of the key objectives in material science, and uncertainty quantification plays a critical role in understanding both process–structure and structure–property linkages. In this work, we seek to learn a distribution of microstructure parameters that are consistent in the sense that the forward propagation of this distribution through a crystal plasticity finite element model matches a target distribution on materials properties. This stochastic inversion formulation infers a distribution of acceptable/consistent microstructures, as opposed to a deterministic solution, which expands the range of feasible designs in a probabilistic manner. Furthermore, to solve this stochastic inverse problem, we employ a recently developed uncertainty quantification framework based on push-forward probability measures, which combines techniques from measure theory and Bayes’ rule to define a unique and numerically stable solution. This approach requires making an initial prediction using an initial guess for the distribution on model inputs and solving a stochastic forward problem. To reduce the computational burden in solving both stochastic forward and stochastic inverse problems, we combine this approach with a machine learning Bayesian regression model based on Gaussian processes and demonstrate the proposed methodology on two representative case studies in structure–property linkages.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Classical Field Theory Formulation for the Numerical Solution of Time Harmonic Electromagnetic Fields

Finite element representations of Maxwell's equations pose unusual challenges inherent to the variational representation of the “curl-curl” equation for the fields. We present a variational formulation which solves for the four-potential instead, based on classical field theory. Borrowing from quantum electrodynamics, we modify the Lagrangian by adding an implicit gauge-fixing term. This reformulation explicitly accounts for Gauss’ law through the coupling between $\phi$ and $\rho$ , and enables the use of nodal basis functions instead of edge elements for time-harmonic problems. Here, we demonstrate how this formulation, adhering to the deeper underlying symmetries of the four-dimensional covariant field description, provides a highly general, robust numerical framework.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

PolyMat: an efficient Matlab code for multi-material topology optimization

We introduce a Matlab implementation of topology optimization for compliance minimization on unstructured polygonal finite element meshes that efficiently accommodates many materials and many volume constraints. Leveraging the modular structure of the educational code, PolyTop, we extend it to the multi-material version, PolyMat, with only a few modifications. First, a design variable for each candidate material is defined in each finite element. Next, we couple a Discrete Material Optimization interpolation with the existing penalization and introduce a new parameter such that we can employ continuation and smoothly transition from a convex problem without any penalization to a non-convex problem in which material mixing and intermediate densities are penalized. Mixing that remains due to the density filter operation is eliminated via continuation on the filter radius. To accommodate flexibility in the volume constraint definition, the constraint function is modified to compute multiple volume constraints and the design variable update is modified in accordance with the Zhang-Paulino-Ramos Jr. (ZPR) update scheme, which updates the design variables associated with each constraint independently. The formulation allows for volume constraints controlling any subset of the design variables, i.e., they can be defined globally or locally for any subset of the candidate materials. Borrowing ideas for mesh generation on complex domains from PolyMesher, we conclude which design variables are associated with each local constraint of arbitrary geometry. A number of examples are presented to demonstrate the many material capability, the flexibility of the volume constraint definition, the ease with which we can accommodate passive regions, and how we may use local constraints to break symmetries or achieve graded geometries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Finite element solution of MHD power-law fluid with slip velocity effect and non-uniform heat source/sink

The laminar steady boundary layer flow of an incompressible MHD Power-law fluid past a continuously moving surface is investigated numerically. The study involves the influence of surface slip and non-uniform heat source/sink on the flow and heat transfer. The governing boundary layer flow equations are transformed into non-dimensional, non-linear coupled ordinary differential equations with the help of suitable similarity transformations. The Galerkin finite element method is implemented to crack the resulting system. The impact of different involved physical parameters is exhibited on the dimensionless velocity profile, temperature distributions and rate of heat transfer in graphical and tabular forms for pseudoplastic and dilatant fluids. The local Nusselt number is found to be the decreasing function of slip parameter, temperature and space-dependent heat sink parameter whereas it increases with increasing values of temperature and space dependent heat source parameter. The problem has important application in attaining the sustainable heat transfer rate for the cooling of fluids, especially in heat exchangers used frequently in chemical industry, in order to increase the trustworthiness of a system, as it removes high heat loads from these systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Resolving Rapid Variation in Energy for Particle Transport

Resolving the rapid variation in energy in neutron and thermal radiation transport is needed for the predictive simulation capability in high-energy density physics applications. Energy variation is difficult to resolve due to rapid variations in cross sections and opacities caused by quantized energy levels in the nuclei and electron clouds. In recent work, we have developed a new technique to simultaneously capture slow and rapid variations in the opacities and the solution using homogenization theory, which is similar to multiband (MB) and to the finite-element with discontiguous support (FEDS) method, but does not require closure information. We demonstrated the accuracy and efficiency of the method for a variety of problems. We are researching how to extend the method to problems with multiple materials and the same material but with different temperatures and densities. In this highlight, we briefly describe homogenization theory and some results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order S{sub N} Equations in Near-Void Regions

The common definition of the diffusion coefficient as the inverse of three times the transport cross section is not compatible with voids. Morel introduced a non-local tensor diffusion coefficient that remains finite in voids. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coefficient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction. It was demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors have been applied to modeling the TREAT reactor confirming the findings of Larsen and Trahan. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S{sub N} equations discretized with the discontinuous finite element method (DFEM). This NDA method uses a scalar diffusion coefficient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it suffers from very large and potentially unbounded diffusion coefficients in the low order problem. However, it was noted that the choice of the diffusion coefficient does not influence consistency of the method at convergence and hence the diffusion coefficient is essentially a free parameter. The choice of the diffusion coefficient does, however, affect the convergence behavior of the nonlinear diffusion iterations. Within this work we use Morel's non-local diffusion coefficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant enhancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coefficients, the periodical horizontal interface problem is used. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Iterated Discontinuity Factors Based on Reference Outgoing Currents for Discontinuous Finite Element Diffusion

We present a novel formulation for iteratively computing Discontinuity Factors (DF) designed to reproduce reference partial outgoing currents at the interface between macro-regions using Discontinuous Galerkin Finite Elements to solve the multigroup neutron diffusion equation. Similarly, an expression for the boundary coefficients (BCf) is proposed to match the partial outgoing current at vacuum boundaries. Upon convergence, macro-region averaged reaction and leakage rates are preserved. In addition, we consider a source-free steady-state problem depending on the DFs. The existence of a non-trivial solution is due to the non-linearity of the problem. This form greatly simplifies the implementation within Rattlesnake and efficiently converges by updating the DFs and BCfs at every linear iteration of the Preconditioned Jacobian-Free Newton-Krylov solver. The main challenge then lies in determining an appropriate scaling factor to ensure that the initial values for the DFs and BCfs are not too far from their converged value. While this factor does not change the final solution, it has the potential to significantly affect convergence, for purely numerical reasons. An expression based on the reference fission source integral is proposed and satisfying results are presented for problems and compared to Superhomogénéisation (SPH) results with both reflecting and vacuum boundary conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Iterated Discontinuity Factors Based on Reference Outgoing Currents for Discontinuous Finite Element Diffusion

We present a novel formulation for iteratively computing Discontinuity Factors (DF) designed to reproduce reference partial outgoing currents at the interface between macro-regions using Discontinuous Galerkin Finite Elements to solve the multigroup neutron diffusion equation. Similarly, an expression for the boundary coefficients (BCf) is proposed to match the partial outgoing current at vacuum boundaries. Upon convergence, macro-region averaged reaction and leakage rates are preserved. In addition, we consider a source-free steady-state problem depending on the DFs. The existence of a non-trivial solution is due to the non-linearity of the problem. This form greatly simplifies the implementation within Rattlesnake and efficiently converges by updating the DFs and BCfs at every linear iteration of the Preconditioned Jacobian-Free Newton-Krylov solver. The main challenge then lies in determining an appropriate scaling factor to ensure that the initial values for the DFs and BCfs are not too far from their converged value. While this factor does not change the final solution, it has the potential to significantly affect convergence, for purely numerical reasons. An expression based on the reference fission source integral is proposed and satisfying results are presented for problems and compared to Superhomogénéisation (SPH) results with both reflecting and vacuum boundary conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

FELIX-2.0: New version of the finite element solver for the time dependent generator coordinate method with the Gaussian overlap approximation

The time-dependent generator coordinate method (TDGCM) is a powerful method to study the largeamplitude collective motion of quantum many-body systems such as atomic nuclei. Under the GaussianOverlap Approximation (GOA), the TDGCM leads to a local, time-dependent Schrödinger equation in amulti-dimensional collective space. In this paper, we present the version 2.0 of the code FELIX that solvesthe collective Schrödinger equation in a finite element basis. This new version features: (i) the ability tosolve a generalized TDGCM+GOA equation with a metric term in the collective Hamiltonian, (ii) supportfor new kinds of finite elements and different types of quadrature to compute the discretized Hamiltonianand overlap matrices, (iii) the possibility to leverage the spectral element scheme, (iv) an explicit Krylovapproximation of the time propagator for time integration instead of the implicit Crank–Nicolson methodimplemented in the first version, (v) an entirely redesigned workflow. We benchmark this release on ananalytic problem as well as on realistic two-dimensional calculations of the low-energy fission of240Puand256Fm. Low to moderate numerical precision calculations are most efficiently performed with simplexelements with a degree 2 polynomial basis. Higher precision calculations should instead use the spectralelement method with a degree 4 polynomial basis. We emphasize that in a realistic calculation of fissionmass distributions of240Pu, FELIX-2.0 is about 20 times faster than its previous release (within a numericalprecision of a few percents).

Regnier, D.↗
DOE OSTI2020

Modeling of deformation of battery cells using thick shell element formulation

Here, we describe a new approach for modeling nonlinear deformation and stress distribution of battery cells using a new thick shell finite element formulation with a through-thickness calculation of stresses and strains that satisfy equilibrium conditions. Battery cells are transversely layered materials that contain numerous thin layers in a repeating sequence. The layers are made of materials with significantly different mechanical responses. Explicit discretization of the layers is computationally impractical except for very small domains, while homogenized material models cannot account for stress and strain variations and partition through the thickness. The new formulation allows for calculation of the transverse and interlayer stresses based on the material properties of the individual cell layers. The application to problems where the transverse stresses have strong influence is also possible.

25 ENERGY STORAGE↗
DOE OSTI2017

Modeling 3D PCMI using the Extended Finite Element Method with higher order elements

This report documents the recent development to enable XFEM to work with higher order elements. It also demonstrates the application of higher order (quadratic) elements to both 2D and 3D models of PCMI problems, where discrete fractures in the fuel are represented using XFEM. The modeling results demonstrate the ability of the higher order XFEM to accurately capture the effects of a crack on the response in the vicinity of the intersecting surfaces of cracked fuel and cladding, as well as represent smooth responses in the regions away from the crack.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗