SEARCH · Engineering Papers
Results for “finite element methods”
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.
Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity
We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint
AIR Algebraic Multigrid for a Space-Time Hybridizable Discontinuous Galerkin Discretization of Advection(-Diffusion)
This article investigates the efficiency, robustness, and scalability of approximate ideal restriction (AIR) algebraic multigrid as a preconditioner in the all-at-once solution of a space-time hybridizable discontinuous Galerkin discretization of advection-dominated flows. The motivation for this study is that the time-dependent advection-diffusion equation can be seen as a “steady” advection-diffusion problem in $(d+1)$-dimensions and AIR has been shown to be a robust solver for steady advection-dominated problems. Numerical examples demonstrate the effectiveness of AIR as a preconditioner for advection-diffusion problems on fixed and time-dependent domains, using both slab-by-slab and all-at-once space-time discretizations, and in the context of uniform and space-time adaptive mesh refinement. A closer look at the geometric coarsening structure that arises in AIR also explains why AIR can provide robust, scalable, space-time convergence on advective and hyperbolic problems, while most multilevel parallel-in-time schemes struggle with such problems.
Contributions to MoDELib SOFTWARE
The purpose of the current request is to enable LANL employees to contribute computer source code to the existing public repository of the MoDELib software package. This software implements discrete dislocation dynamics (DDD) and finite element (FEM) methods and is currently a vital component of an ongoing DR project at LANL, in collaboration with its original author and maintainer Giacomo Po. Contributions from LANL employees would aim to enhance the reliability, accuracy, and performance of MoDELib simulations using LANL's high performance computing platforms through bug fixes, algorithmic refinements, and parallelization.
HygroThermFEM v1.0
HygroThermFEM is a Finite Element Method-based numerical calculation engine for solving 2-D heat and moisture transfer problems. This numerical engine is used in the THERM software tool, and its primary purpose is for the analysis of building envelopes (e.g., windows, walls, roofs, foundations, etc.). However, the engine can also be used for any heat and moisture transfer problems that require solving fundamental 2-D energy and mass transfer equations. Fluid flow solutions (Navier-Stokes momentum equations) are not included, but the correlations for various convection heat transfer situations are provided, including the translation of complex cavity geometries into those for which correlations are applicable. The calculation engine is written in C++ and includes an API for connecting to third-party tools.
Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation
To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.
EchemFEM: A Firedrake-based Python package for electrochemical transport
The transition from fossil fuels to renewable energy has brought about a rapid increase in the availability of clean electricity. However, electricity generated from sources such as wind and solar are limited to intermittent operation due to daily and seasonal variation. One solution is to utilize electrochemical devices in energy storage and electrochemical manufacturing applications, where they can harness surplus energy and decarbonize chemical industries traditionally reliant on petrochemical feedstocks. Managing the growing prevalence of renewable energy underscores the importance of developing and scaling up these technologies, which can in turn facilitate the achievement of carbon emission reduction commitments of companies and developed economies. Likewise, the electrification of transport creates an increasing need for energy-dense electrochemical energy storage devices such as batteries and supercapacitors. Naturally, simulation tools are required to assist in the design of efficient and industrial-scale electrochemical devices.
A Partitioned Method for Reduced OrderModel-Finite Element Model (ROM-FEM)and ROM-ROM Couplings with SeparateReduced Bases for Interior and InterfaceVariables
Explore the source record for details and available documents.
Large-Scale Welding Process Simulation by GPU Parallelized Computing
The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.
The Biomechanics of Indirect Traumatic Optic Neuropathy Using a Computational Head Model With a Biofidelic Orbit
Indirect traumatic optic neuropathy (ITON) is an injury to the optic nerve due to head trauma and usually results in partial or complete loss of vision. In order to advance a mechanistic understanding of the injury to the optic nerve, we developed a head model with a biofidelic orbit. Head impacts were simulated under controlled conditions of impactor velocity. The locations of impact were varied to include frontal, lateral, and posterior parts of the head. Impact studies were conducted using two types of impactors that differed in their rigidity relative to the skull. The simulated results from both the impactors suggest that forehead impacts are those to which the optic nerve is most vulnerable. The mode and location of optic nerve injury is significantly different between the impacting conditions. Simulated results using a relatively rigid impactor (metal cylinder) suggest optic nerve injury initiates at the location of the intracranial end of the optic canal and spreads to the regions of the optic nerve in the vicinity of the optic canal. In this case, the deformation of the skull at the optic canal, resulting in deformation of the optic nerve, was the primary mode of injury. On the other hand, simulated results using a relatively compliant impactor (soccer ball) suggest that primary mode of injury comes from the brain tugging upon the optic nerve (from where it is affixed to the intracranial end of the optic canal) during coup countercoup motion of the brain. This study represents the first published effort to employ a biofidelic simulation of the full length of the optic nerve in which the orbit is integrated within the whole head.
Failure in Confined Brazilian Tests on Sandstone
Strength of rocks in the confined tension region, where the minimum principal stress is tensile, has only infrequently been measured and is not well understood. Quasi-static confined Brazilian tests under a range of confining stresses (2.76 to 27.58 MPa) where used to determine the strength of sandstone in the confined tension region. The test results indicate that the strength in the confined tension region was a strong function of the intermediate principal stress: increasing the intermediate principal stress significantly increased the strength of the sandstone. The strength data were well fit by the Mogi–Coulomb criterion, which accounts for the intermediate principal stress. Unconfined Brazilian strength data were not well fit to the Mogi–Coulomb criterion derived from the confined Brazilian test data, consistent with a transition from tensile to shear processes dominating failure with increasing confining pressure. Observations of post-failure fracture surfaces reveal more indication of shear processes with increasing confining pressure. Numerical simulations from combined finite-discrete element method are compared to the experimental results and reflect similar conditions for failure compared to the experimental tests in the confined tension region.
Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination
Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.
A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice
The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.
Reduced basis approximations of parameterized dynamical partial differential equations via neural networks
Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. Here, we demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.
Analysis of spent nuclear fuel multipurpose canister dynamics during rail transportation
The primary mode of spent nuclear fuel transportation within the United States will be by railcar. One such system is the Atlas railcar, which is designed to transport 17 different spent nuclear fuel cask systems, including bare fuel systems and canister fuel systems. In the latter configuration, multipurpose canisters containing spent nuclear fuel may be placed within an overpack for storage, or within a cask for transportation. Compared to bare fuel systems, canister fuel systems have additional degrees of freedom for motion during transportation, because clearance between the cask and canister allows for some motion of the canister to occur relative to the cask. This work investigates the effect of canister motion on the shock and vibration imparted to the spent nuclear fuel within. Structural dynamic analyses have been conducted to identify the effects of canister to cask clearance, presence and type of dunnage, and loading direction and frequency. This modeling study calculates anticipated cask motion, canister motion, and spent nuclear fuel structural dynamic response to normal conditions of transportation railcar motion using finite element analysis methods that were developed to model the rail segment of the ENSA/DOE (Equipos Nucleares S.A., U.S. Department of Energy) multimodal transportation test of 2017.
M3 milestone: Advanced contact 2021
This report describes advanced contact approaches and implementations in FY21 for use in emerging reactors with geometries that differ from LWR/EBR-II cylindrical fuel. Contact remains a challenging and crucially important component of our modeling toolkit. While significant improvements have been made through mortar and automatic differentiation, we explored advanced approaches (i.e., dual mortar) that yields better conditioned systems and will improve the robustness and performance of BISON contact problems. This effort is relevant to emerging reactor types, such as micro reactors and the corresponding complex and numerous contact surfaces. In this report, improvements of the mortar-based contact formulation are described. Benchmark problems and BISON assessment cases are included to demonstrate the improvements and advantages that are brought by the recent capability. Remaining issues are discussed and corresponding development plans are outlined.
[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment
In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.
CODE-TO-CODE BENCHMARK STUDY FOR THERMAL STRESS MODELING AND PRELIMINARY ANALYSIS OF THE HIGH-TEMPERATURE SINGLE HEAT PIPE EXPERIMENT
Microreactors are very small nuclear reactors with typical thermal-energy output of up to 20 MWth. The microreactor concept is gaining more and more attention for safe, robust, and reliable supply of electricity for remote locations and to meet industrial process heat needs. Heat-pipe-cooled microreactor is one of the microreactors being investigated at Idaho National Laboratory. In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios for the heat-pipe-cooled microreactor designs. This paper discusses finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory, using two commercial software packages, Abaqus and STAR-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the locations where the largest temperature gradients appeared between heat pipe and electric heater. The temperature fields have good agreements between Abaqus and STAR-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes.