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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Stability of Panels of a Thermoelastic Shell Heated at the Edge by a Mobile Point Source

The thermostability of thin-walled structures is considered by means of a thermoelasticity problem in a three-dimensional quasi-static formulation. Finite-element models of panels in the form of a circular cylinder and a plate reinforced by ribs are constructed. Numerical solutions are obtained for the stability loss of thermoelastic thin-walled reinforced structures with intense edge heating by a mobile heat source.

42 ENGINEERING↗
DOE OSTI2017

Heat Source Characterization in a TREAT Fuel Particle Using Coupled Neutronics-BCMC Calculations

This work presents a multi-physics, multi-scale approach to modeling the Transient Test Reactor (TREAT) currently prepared for restart at the Idaho National Laboratory. TREAT fuel is made up of microscopic fuel grains (r ˜ 20µm) dispersed in a graphite matrix. The novelty of this work is in coupling a binary collision Monte-Carlo (BCMC) model to the Finite Element based code Moose for solving a microsopic heat-conduction problem whose driving source is provided by the BCMC model tracking fission fragment energy deposition. This microscopic model is driven by a transient, engineering scale neutronics model coupled to an adiabatic heating model. The macroscopic model provides local power densities and neutron energy spectra to the microscpic model. Currently, no feedback from the microscopic to the macroscopic model is considered. TREAT transient 15 is used to exemplify the capabilities of the multi-physics, multi-scale model, and it is found that the average fuel grain temperature differs from the average graphite temperature by 80 K despite the low-power transient. The large temperature difference has strong implications on the Doppler feedback a potential LEU TREAT core would see, and it underpins the need for multi-physics, multi-scale modeling of a TREAT LEU core.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Time Scale Effects on the Temperature Coefficient of Reactivity

For a nuclear reactor, the temperature coefficient of reactivity (TCR) is of significance for both power control and safety analysis, but it is not always a constant. For example, the TCR may be dependent upon temperature due to Doppler broadening, and may be also dependent upon time scale due to heat transfer and thermal expansion. A fast burst reactor is a favorable nuclear facility to simulate the time scale effects on the TCR. Fast burst reactors are generally unmoderated research reactors of enriched uranium metal or uranium alloys. The neutrons are fast and the Doppler broadening effects can be neglected, and the TCR is mainly contributed by fuel thermal expansion. Fast burst reactors can operate on both steady state power mode and burst power mode. In steady state power mode, the systems are assumed to be in thermal equilibrium. In burst mode, the full width at half maximum (FWHM) of the power is from tens of microseconds to a few milliseconds, and the heat transfer is generally neglected. If FWHM is short enough, thermal expansion will have a time delay due to inertial effect. This paper is aimed to study the exact time scale effects on the TCR due to thermal expansion. The TCR contributed by the Doppler broadening effects was neglected, because it is relatively small (The k{sub eff} as a function of temperature was also calculated using the Monte Carlo particle transport code JMCT by our group, and no definite difference due to the temperature change was obtained near room temperature), and there is no method to evaluate the time scale dependence of the Doppler broadening effects yet. In the time scale from 1 x 10{sup -6} s to 1 x 10{sup 5} s, the TCR of a Godiva-like spherical fast burst reactor was calculated with perturbation theory. A Godiva-like spherical fast burst reactor was studied, the calculation of the time scale effects on TCR involved thermal conduction, thermoelasticity and reactivity variation analysis. Both the thermal and the thermoelastic problems were solved using ANSYS software. ANSYS is software implementation of finite element simulations, which is widely used in engineering field. The reactivity variation due to thermal expansion was calculated with perturbation theory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
DOE OSTI2023

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A modified moment-fitted integration scheme for X-FEM applications with history-dependent material data

Here, we present a strategy for the numerical integration of partial elements with the eXtended finite element method (X-FEM). The new strategy is specifically designed for problems with propagating cracks through a bulk material that exhibits inelasticity. Following a standard approach with the X-FEM, as the crack propagates into new regions of the domain, elements are split into several copies that contain pieces of the domain on either side of the crack. We examine quadrature rules that have sufficient accuracy to calculate stiffness matrices regardless of the orientation of the crack with respect to the element. This permits the number of integration points within elements to remain constant as a crack propagates, and for state data to be easily transferred between successive discretizations. In order to maintain weights that are strictly positive, we propose an approach that blends moment-fitted weights with volume-fraction based weights. To demonstrate the efficacy of this simple approach, we present results from numerical tests and examples with both elastic and plastic material response.

42 ENGINEERING↗
DOE OSTI2021

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2021

Mechanics of vitrimer particle compression and fusion under heat press

The compression and fusion of vitrimer particles is a fundamental problem underlying the recycling of vitrimers. Incomplete particle compression can lead to voids in the fused vitrimer, thereby impacting its mechanical property. In this work, we present a two-dimensional finite element model to capture the interplay of multiple complex mechanisms including random particle packing, large deformation, inter-particle contact, thermally activated bulk stress relaxation, and interface healing. Specifically, we focus on understanding two key components of this problem: i) evolution of porosity during compression and fusion of randomly packed vitrimer particles, and ii) effective tensile modulus and strength of fused vitrimer. Using the finite element model, we first show that using smaller particles can reduce the porosity in the fused vitrimer, which only slightly increases the effective modulus, but the increase in tensile strength is more pronounced. In addition, under the same processing conditions, particles with mixed sizes can achieve better densification and higher tensile strength after fusion than uniformly sized particles with the same average radius. Regarding the effects of processing conditions, we discuss how the porosity and effective modulus of the fused vitrimer depend on the processing time, temperature and pressure, and find a general trend consistent with experimental observations, i.e. longer processing time, higher temperature or higher pressure can lead to lower porosity and higher modulus of the fused vitrimer. Finally, the theoretical insights towards the compression and fusion process as well as the simulation framework can be useful in optimizing the powder-based heat press process of vitrimer and its composites.

42 ENGINEERING↗
DOE OSTI2017

Stencil computations for PDE-based applications with examples from DUNE and hypre

Here, stencils are commonly used to implement efficient on–the–fly computations of linear operators arising from partial differential equations. At the same time the term “stencil” is not fully defined and can be interpreted differently depending on the application domain and the background of the software developers. Common features in stencil codes are the preservation of the structure given by the discretization of the partial differential equation and the benefit of minimal data storage. We discuss stencil concepts of different complexity, show how they are used in modern software packages like hypre and DUNE, and discuss recent efforts to extend the software to enable stencil computations of more complex problems and methods such as inf–sup–stable Stokes discretizations and mixed finite element discretizations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Generalized Multiscale Finite Element Methods with energy minimizing oversampling

In this work, we propose a general concept for constructing multiscale basis functions within Generalized Multiscale Finite Element Method, which uses oversampling and stable decomposition. The oversampling refers to using larger regions in constructing multiscale basis functions and stable decomposition allows estimating the local errors. The analysis of multiscale methods involves decomposing the error by coarse regions, where each error contribution is estimated. In this estimate, we often use oversampling techniques to achieve a fast convergence. We demonstrate our concepts in the mixed, the Interior Penalty Discontinuous Galerkin, and Hybridized Discontinuous Galerkin discretizations. One of the important features of the proposed basis functions is that they can be used in online Generalized Multiscale Finite Element Method, where one constructs multiscale basis functions using residuals. In these problems, it is important to achieve a fast convergence, which can be guaranteed if we have a stable decomposition. In our numerical results, we present examples for both offline and online multiscale basis functions. Our numerical results show that one can achieve a fast convergence when using online basis functions. Moreover, we observe that coupling using Hybridized Discontinuous Galerkin provides a better accuracy compared with Interior Penalty Discontinuous Galerkin, which is due to using multiscale glueing functions.

42 ENGINEERING↗
DOE OSTI2018

Unique Determination of 3-Connected Plane Domains by Relative Conformal Moduli of Pairs of Boundary Components

We study the boundary values of conformal mappings of plane finitely connected domains. We present a complete exposition of a theorem giving an exhaustive description of the boundary values of such mappings in terms of conformal moduli (extremal lengths) of pairs of boundary components of a domain of connectedness at most 3. Bibliography: 12 titles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A Finite-Element Approach for the Partitioning of Carbon in Q&P Steel

This paper reports a Galerkin finite-element analysis of carbon partitioning from martensite into austenite during the quenching and partitioning (Q&P) processes in steels. In contrast to classical or sophisticated diffusion field models, an alternative nonlinear governing equation based on chemical potential and composition is considered. The model is applied to simulate the carbon partitioning of modified 22MnB5 alloys assuming an immobile austenite–martensite phase boundary and, in turn, is compared with experimental measurements of the volume fraction and carbon content of retained austenite. The simulations show outstanding results on the influence of Si and C content on the partitioning parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Mesoscale Numerical Characterization of Kevlar and Carbon–Kevlar Hybrid Plain-Woven Fabric Compression Behavior

A mesoscale finite element model (FEM) of the plain-woven fabric representative elementary volume (REV/unit cell) is developed. The developed FEM is used to predict the nonlinear compaction behavior of monolithic and intraply hybrid fabric made of high-strength carbon and Kevlar fibers. For geometrical modeling of the fabric REV, TexGen, an open-source geometrical modeling software (developed at Nottingham University), is utilized, and for mechanical characterization, FE simulation tool Abaqus{sup ®} is used. Geometrical nonlinearity of the fabric, including periodic boundary conditions and yarn interactions, is incorporated in the REV FE model. Material nonlinearity caused by transversely isotropic behavior of fibrous yarns is implemented through the user material subroutine form in FE simulation along with local orientation assigned to each yarn in the REV. Material properties of fibers are extracted from the literature and supplier data. The compression behavior of monolithic Kevlar and hybrid carbon–Kevlar (carbon warp and Kevlar weft yarns, C–K) fabric is predicted using the proposed FE model. FEM-predicted results were found to have good agreement with the experimentally obtained ones. Results depict that C–K hybrid fabrics have improved compression strength. Hybridization effect is governed by yarns mutual interaction, yarn area at crossovers, and transverse in-plane properties of the interyarn hybridized fabric. The proposed FE model can be used to predict the mechanical behavior of woven fabrics with different weaving patterns, unit cell geometry, material properties, and loading conditions.

36 MATERIALS SCIENCE↗