Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 217 records · Page 12

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗

STRUCTURAL MODELING TO SUPPORT POST-YIELD ACCEPTANCE CRITERIA FOR SPENT NUCLEAR FUEL CLADDING

Spent nuclear fuel (SNF) is evaluated for structural failure during storage and transportation scenarios. The U.S. Department of Energy’s Spent Fuel and Waste Science and Technology (SFWST) program has sponsored significant research in quantifying mechanical loads on SNF during storage and transportation scenarios using experimental and modeling methods. The SFWST program has also performed significant research on measuring the mechanical behavior of irradiated SNF as defueled cladding segments and cladding with fuel pellets to measure composite behavior. This paper considers some of the key material data from the Sibling Pin testing and uses structural modeling and analysis methods that have been informed by testing to consider post-yield acceptance criteria for SNF cladding structural analysis. Test data published by Oak Ridge National Laboratory (ORNL) and Pacific Northwest National Laboratory (PNNL) are the foundation for informing the material behavior of the models developed in this study. In particular, four-point bend (4PB) tests of fueled and defueled cladding segments provide significant information about the bending failure mode of SNF. ORNL’s 4PB test data is on fueled cladding segments, so the composite behavior of SNF is demonstrated. This paper describes PNNL’s coincident beam model that was developed to approximate the composite behavior of SNF. This paper also presents PNNL’s structural dynamic finite element models of a cask tip-over scenario, which is predicted to cause the strongest mechanical loads on SNF of all postulated storage and transportation scenarios. SNF bending loads predicted in the cask tip-over scenario and cladding acceptance criteria beyond yield are considered, with justification based on the Sibling Pin test data. ASME Boiler and Pressure Vessel code stress intensity limits are also considered. The ultimate goal of this work is to aid in the justification of structural acceptance criteria for SNF cladding beyond the cladding’s irradiated yield strength for use in structural analysis of all storage and transportation scenarios.

Klymyshyn, Nicholas A.↗

Optimizing the use of pressurized bladders for the assembly of HL-LHC MQXFB magnets

The use of pressurized bladders for stress control of superconducting magnets was firstly proposed at Lawrence Berkeley National Laboratory in the early 2000s. Since then, the so-called ‘bladders and keys’ procedure has become one of the reference techniques for the assembly of high-field accelerator magnets and demonstrators. Exploiting the advantages of this method is today of critical importance for Nb$_{3}$Sn-based accelerator magnets, whose production requires the preservation of tight stress targets in the superconducting coils to limit the effects of the strain sensitivity and brittleness of the conductor. The present manuscript reports on the results of an experimental campaign focused on the optimization of the ‘bladders and keys’ assembly process in the MQXFB quadrupoles. These 7.2 m long magnets shall be among the first Nb$_{3}$Sn cryomagnets to be installed in a particle accelerator as a part of the High Luminosity upgrade of the LHC. One of the main practical implications of the bladders technique, especially important when applied to long magnets like MQXFB, is that to insert the loading keys, the opening of a certain clearance in the support structure is required. The procedure used so far for MQXF magnets involved an overstress in the coils during bladder inflation. The work presented here shows that such an overshoot can be eliminated thanks to additional bladders properly positioned in the structure. This optimized method was validated in a short model magnet and in a full-length mechanical model, becoming the new baseline for the series production at CERN. Furthermore, the results are supported by numerical predictions using finite element models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗

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↗

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

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

14 SOLAR ENERGY↗

Integrating Crack Detection and Pipe Shape Optimization for Enhanced Sewage System Durability

Crack detection in underground reinforced concrete pipes has been essential in determining the state of stormwater infrastructure. Detection models have been implemented for detecting cracks and other defects in pipes using CCTV footage for stormwater drainage systems. In addition, Finite element models have been used to determine optimum shapes and pipe thickness for different boundary conditions such as header pipes in power plants. The concept of shape optimization emerges as a crucial factor in power plant design and operation, with the potential to maximize performance while minimizing the use of materials. Shape optimization not only enhances efficiency but also contributes to reducing the environmental footprint. This paper discusses the integration of both topics by using the cracks detected in underground pipes as boundary conditions for shape optimization of the pipes. A machine learning model has been developed which uses limited data for training and outlines the location of detected cracks. A shape optimization methodology is proposed in which ANSYS modules are used to analyze fluid flow and then optimize the shape of the pipe. The crack detection model developed has been applied to a crack detected in lab setting and machine learning model used has an accuracy of 98% using a random forest algorithm.

20 FOSSIL-FUELED POWER PLANTS↗

Evaluation of Surface Crack Formation in Photovoltaic Backsheets Using Fragmentation and Finite Element Simulations

Backsheet cracking is among the most commonly observed degradation modes of photovoltaic modules in the field. Cracks can reduce the ability of backsheets to fulfill their functions, such as protection of the modules from the environment or electrical insulation. Here, this work presents an evaluation of the degradation and cracking propensity of two backsheets during accelerated IEC TS 62788-7-2 (International Electrotechnical Commission) aging under the A3 condition with a fragmentation test using a coextruded polyamide backsheet (AAA), and a laminated multilayer backsheet with a polyethylene terephthalate core and outer layer and ethylene vinyl acetate inner layer (PPE). Results show the surface embrittlement of the AAA outer layer during exposure. A longer time of exposure causes cracks to form at lower strains during stretching and creates deeper cracks. Accordingly, Young's modulus of the outer layer increases, as measured by cross-sectional nanoindentation. PPE exhibits cracking after exposure as well. While no similar increase of modulus or crack depth can be observed, the outer layer of PPE exhibits more obvious signs of erosion during exposure, including progressive morphological changes and thickness losses. A finite element model is developed to simulate surface crack formation, based on initially zero-thickness decohesion elements. Decohesion criteria define a critical stress, at which these elements grow, and cracks can begin to form. These criteria are obtained via parameter optimization by comparison between simulation and experiment. The model is used to interpret crack formation in both backsheets.

14 SOLAR ENERGY↗

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

Sierra/SolidMechanics 5.0 Goodyear Specific

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 5.0 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 4.56 Goodyear User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.56 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

42 ENGINEERING↗

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

The daily as-run thermal analysis for the Advanced Gas Reactor (AGR)-5/6/7 experiment was documented in ECAR-5633 (Hawkes 2021). A finite element model was created for the entire test train, which consisted of five capsules. The fuel compacts, graphite holders, stainless-steel capsule walls, and all other major components were individually modeled, including all thermocouples (TCs). Daily heat rates for each compact and component in the model were input from daily as-run physics analyses. Daily gas compositions for each capsule were also input. The thermal conductivity of the compacts and graphite holders varied with fast neutron fluence. Gas mixture thermal conductivity was implemented using experimentally attained values from the literature. Fluence and temperature-dependent thermal conductivity was used for the graphite components and the fuel compacts. Radiation heat transfer was implemented. The model was tuned in an attempt to match the TC readings during the first cycle by adjusting the Neolube (graphite lubricant) thickness.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗