Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary finite element”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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 OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Strain localization criteria for viscoplastic geomaterials

This report presents a viscoplastic localization criterion to detect quasi-instantaneous (i.e., load-induced) and delayed (creep-induced) strain localization in rate-dependent solids. The study is based on the theory of controllability and a viscoplastic description of the mechanical response. Analytical precursors of unstable states are defined through systems of ordinary differential equations (OEDs). The use of the proposed criteria is illustrated at the material point level through a set of strain localization analyses simulating active strain localization of a porous rock. In addition, full-field finite element simulations of compression tests conducted under various pressures are reported to demonstrate the role of local unstable viscoplasticity in the spontaneous propagation of deformation bands under stationary boundary conditions. The study shows that the viscoplastic localization criterion maintains a negative sign as long as the behavior is unstable, that is, the rate of deformation is accelerating. The sign switch coincides with the transition to decelerating deformation. The analyses revealed that pulses of overstress always emerge in correspondence with the growth of unstable behavior, and the peak matches the transition to stable behavior. The local responses recovered from full-field analyses were consistent with those observed in analyses at material point level and the predictions of the presented theory.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗
DOE OSTI2024

Reduction of Methane Leaks through Corrosion Mitigation Pre-treatments for Pipelines with Field Applied Coatings

Corrosion of buried, coated steel pipelines transporting natural gas is a significant source of methane emissions, from pipeline venting required for maintenance and repairs and from pipeline leaks and incidents. Corrosion of steel under field applied coatings is an important safety concern for the pipeline industry. This project investigates the application of a field applied alloy over girth welds to mitigate external corrosion of buried coated steel pipelines. Various metallic coating options were considered, which were required to meet several criteria: (1) it must resist corrosion under open-circuit or mild cathodic protection conditions, (2) it must protect the substrate steel, and (3) it must not negatively affect the adhesion of the field coating. Finite element models and lab testing were performed of alloy coating compositions to identify promising alloy types underneath disbonded coatings. Polarization curves of coating alloys were generated to provide the boundary conditions for the COMSOL model to compute potential and current distributions around coated areas. Sacrificial and corrosion-resistant metal alloy coatings were evaluated and optimized using corrosion modeling and laboratory electrochemical testing, where aluminum alloy 5356 (5% Mg) and steel alloy B9 (9% Cr) were selected. Corrosion test coupons were designed and fabricated using thermal spray aluminum 5356 and welded B9 steel overlays on API 5L grade X42 line pipe steel. The corrosion test coupons, with simulated pipe coating damage, were tested in a laboratory soil box and a field pipeline site in Texas for 3-months. Corrosion test coupons were then tested for 6-months at field pipeline sites in Texas and Tennessee to quantify corrosion rates and performance of the aluminum and steel alloys under polyethylene tape and 2-part epoxy coatings, various coating holidays, and with and without cathodic protection.

03 NATURAL GAS↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗
DOE OSTI2022

Crystallographic effects on transgranular chloride-induced stress corrosion crack propagation of arc welded austenitic stainless steel

Abstract The effect of crystallography on transgranular chloride-induced stress corrosion cracking (TGCISCC) of arc welded 304L austenitic stainless steel is studied on >300 grains along crack paths. Schmid and Taylor factor mismatches across grain boundaries (GBs) reveal that cracks propagate either from a hard to soft grain, which can be explained merely by mechanical arguments, or soft to hard grain. In the latter case, finite element analysis reveals that TGCISCC will arrest at GBs without sufficient mechanical stress, favorable crystallographic orientations, or crack tip corrosion. GB type does not play a significant role in determining TGCISCC cracking behavior nor susceptibility. TGCISCC crack behaviors at GBs are discussed in the context of the competition between mechanical, crystallographic, and corrosion factors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

SC Wall-to-RC Basemat Over-Strength Connection: Behavior and Design

This paper presents results from experimental and analytical investigations conducted to evaluate the lateral load behavior and capacity of steel-plate composite (SC) wall-to-reinforced concrete (RC) basemat connections. Two SC wall-to-reinforced concrete basemat connection specimens were tested. These SC wall specimens had a height-to-length ratio of 0.6 and did not include boundary elements. The experimental results include the lateral force-displacement (V-Δ) responses of the specimens and observations of local damage such as steel plate local buckling and concrete crushing. 3D finite element models were developed and benchmarked using the experimental results. The benchmarked models were used to conduct analytical parametric studies, expand the database, and gain additional insights into the behavior of SC wall-to-RC basemat connections. The parameters included in analytical investigations were the wall aspect ratio (h/lw), reinforcement ratio (ρ), and wall thickness (T).

Kurt, Efe G.↗
DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part II: Entropy Stability

In the hyperbolic research community, there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability. In the first part of the series, the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly. By applying this technique, the authors demonstrate that a pure continuous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way. In this work, we extend this investigation to the nonlinear case and focus on entropy conservation. Here, by switching to entropy variables, we provide an estimation of the boundary operators also for nonlinear problems, that guarantee conservation. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

University of Missouri Research Reactor (MURR) Design Demonstration Element End Fitting Structural Rigidity Analysis

The primary objective of this work is to assess the extent to which the stiffness (measured by means of maximum displacement) of the end fittings in the DDE contributes to the stiffness of the entire element, and how it compares to the equivalent stiffness of the end fittings in the LEU element. Structural analysis of both the LEU element and the DDE were performed using COMSOL 5.3a finite element software. Supporting combs are used on the leading and trailing edges of fuel plates for both the LEU element and the DDE. Therefore, simulations with and without combs are performed as two bounding boundary conditions on the leading edge of the fuel plates. Three types of loads are analyzed in this work: the hydraulic load due to the channel flow disparity-induced pressure differential, the thermal load due to the thermal expansion of the fuel plates, and a point load equal in magnitude to the LEU element’s weight.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Three-dimensional nucleation and growth of deformation twins in magnesium

At two-thirds the weight of aluminum, magnesium alloys have the potential to reduce the fuel consumption of transportation vehicles. These advancements depend on our ability to optimize the desirable versus undesirable effects of deformation twins, which are three-dimensional (3D) microstructural domains that form under mechanical stresses. Previously only characterized through surface or thin-film measurements, we here present 3D in situ characterization of deformation twinning inside an embedded grain over mesoscopic fields of view using dark-field x-ray microscopy supported by crystal plasticity finite element analysis. The results revealed the role of triple junctions on twin nucleation and the sequence and irregularity of twin growth and showed that twin-grain junctions, twin-twin junctions, and twin boundaries were the sites of localized dislocation accumulation.

Lee, Sangwon [Univ. of Michigan, Ann Arbor, MI (Un↗
DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗