Engineering Papers⌕ Search

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.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Constitutive modeling of δ-phase zircaloy hydride based on strain rate dependent nanoindentation and nano-scale impact dataset

In this work, a constitutive model for zircaloy-4 embedded with δ-phase hydride rim/blister structures (δ-phase zircaloy hydride) is presented. Nanoindentation and nano-scale impact experiments were performed on the δ-hydride blisters as a function of temperature in order to obtain strain rate dependent location specific constitutive response. A power law viscoplastic constitutive model was chosen to represent the strain rate dependent behavior. The strength coefficient and the strain hardening exponent in the power law relation were found to decrease with increase in temperature indicating softening behavior at higher temperatures. The strain rate exponent remained constant with increase in temperature indicating that the rate dependent ductility does not change significantly at higher temperatures. A finite element method based model was created to analyze the failure behavior and the ductile to brittle transition in zircaloy hydride. The model predicted that fracture strength of δ-phase zircaloy hydride decreases linearly with increase in temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Weight function procedure for reduced order fracture analysis of arbitrary flaws in cylindrical pressure vessels

Fracture mechanics calculations using the finite element method can be computationally expensive, which makes them challenging to use in engineering evaluations of flaws in pressure vessels. The weight function (WF) technique is a reduced order fracture modeling method that is widely used to greatly reduce these computational costs. Its computational efficiency is essential for use in probabilistic fracture mechanics evaluations of embrittled nuclear reactor pressure vessels (RPVs), due to the large number of sampled flaws that must be evaluated. Although the WF technique is general, it is typically used only for axis-aligned flaws. Recent discoveries of off-axis flaws in operating nuclear reactors have necessitated detailed simulations of such flaws and the interactions between neighboring flaws. This paper presents a generalized WF approach applicable for analyzing arbitrary flaw geometries in thick-walled cylindrical vessels, including surface-breaking and subsurface flaws, which can either axis-aligned or off-axis, and can account for interactions with other flaws. Herein, this approach is demonstrated on representative simulations of multiple flaw geometries in a RPV subjected to transient loading conditions. In all cases, the WF technique gives good comparison with benchmark results from direct simulation with greatly reduced computational effort.

42 ENGINEERING↗
DOE OSTI2022

A unified finite strain gradient-enhanced micropolar continuum approach for modeling quasi-brittle failure of cohesive-frictional materials

Here, in this work, a novel framework for modeling quasi-brittle crack propagation and shear band dominated failure of cohesive-frictional materials like concrete, mortar, rock, tough ceramics, energetic materials, but also granular materials like sands or powders in terms of a unified continuum approach is proposed. It is based on a combination of the gradient-enhanced continuum with gradients of internal variables for representing quasi-brittle cracking, and the micropolar continuum, accounting for the deformation of the microstructure. For developing the gradient-enhanced micropolar framework, the set of balance equations and the kinematic relations are derived, and the constitutive relations are established in a general manner. The framework is formulated in a geometrically exact setting, based on the thermodynamically sound theory of hyperelasto-plasticity, and the numerical implementation by means of the finite element method is discussed. For assessing the approach, realizations of this new approach in terms of constitutive models for particular materials are developed. They are applied to numerical benchmark examples, investigating various loading conditions, and the obtained results are validated by means of a comparison with experiments from the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating the bond strength and fracture mechanisms of cold-sprayed zinc coating on AZ91 magnesium substrate via a combined experimental and computational approach

Magnesium (Mg) alloys are ideal candidates for automotive applications due to their high strength to weight ratio, castability, recyclability etc., however, they lack corrosion and oxidation resistance. Solid-state deposition techniques, such as cold spray, have been demonstrated to enhance their corrosion resistance as it relies on the severe plastic deformation of powder particles upon impact with the substrate to form a metallurgical bond with the substrate and within the coating. At cold sprayed interfaces, a heterogeneous microstructure is formed that includes some porosity, oxides and intermetallics which can significantly affect coating performance. Thus, establishing a direct correlation between the interface microstructure and its properties can aid in designing optimal cold spray parameters. In this study, we investigated the microstructure and mechanical properties of a zinc (Zn) coating deposited on a high pressure die cast (HPDC) AZ91 Mg substrate via high resolution scanning transmission electron microscopy, in situ micro-tensile testing, and finite element method (FEM) modeling. Micro-tensile pillars fabricated using the plasma focused ion beam (PFIB) successfully isolates the coating-substrate interface within the gauge length. The average bond strength of Zn-Mg interface was determined to be ∼140 MPa with failure occurring partially at the interface and mostly into the coatings. A detailed microstructural characterization revealed evidence of a strong metallurgical bonding at the Zn-Mg interface and formation of the C14 MgZn 2 laves phase interlayer resulting in a mixed mode of fracture during the micro-tensile experiments. FEM modeling reveals the stress distribution along the interfaces and suggests that a MgZn 2 layer thickness between 200–400 nm is optimum to increase the bond strength and minimize the triaxiality. Such a site-specific interfacial analysis with correlative computational modeling provides crucial insight into the overall performance of cold spray interfaces.

Bond strength↗
DOE OSTI2021

Micro-mechanical investigation of the thermo-mechanical properties of micro-architectured tungsten coatings

The thermo-mechanical response of micro-architectured tungsten coatings is characterized in the temperature range of 293 to 673 K using both in situ micro-compression experiments inside a scanning electron microscope (SEM) as well as image-based crystal plasticity finite element method (CPFEM) simulations. The experiments were conducted on micropillar-like specimens that were focus ion beam milled into the coatings, while the simulations were conducted on columnar-grained micropillar simulation cells constructed based on the statistics of the coating’s microstructure. The experimental results show that the stress–strain response and deformation mode exhibit a strong temperature-dependence. At room temperature, catastrophic failure is observed shortly after yield and is manifested in the form of intergranular fracture and buckling of individual columnar grains. With increasing temperature, this catastrophic failure is gradually suppressed and the material exhibits a steadier strain hardening response at 693 K. The CPFEM simulations are also shown to be in good agreement with the experimental results, and these simulations indicate that the material response is strongly influenced by the local crystallographic anisotropy and microstructure inhomogeneity. Furthermore, the simulations capture the underlying mechanisms that control the temperature-dependent transition in deformation mode. The current results highlight that the micro-architectured microstructure offers a great combination of excellent mechanical strength and structural integrity at elevated temperatures, which is of importance for high temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
DOE OSTI2021

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Mesoscale modeling of microstructure-dependent thermal conductivity in U-Zr fuels

In uranium-zirconium (U-Zr) based metallic fuels, different phases can form at different compositions and temperatures. Typically, lamellar δ-UZr 2 and α-U phases are the dominant microstructures in U-rich U-Zr alloys at temperatures below 880 K. In this work, a finite element method based mesoscale modeling technique is used to calculate the effective thermal conductivities of such heterogeneous microstructures, using the thermal conductivities of two individual phases and their interphase thermal resistance (Kapitza resistance) as input parameters. The Kapitza resistance between δ-UZr 2 and α-U is determined at different temperatures, which shows an approximately T 3 dependence in the temperature range between 300 and 800 K. In addition, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr 2 phase. Further, an analytical model is therefore developed to quantify the effects of both temperature and δ-UZr 2 aspect ratio on the Kapitza resistance. Using this newly developed Kapitza resistance model, the effective thermal conductivities of a number of δ-UZr 2 + α-U heterogeneous microstructures in U-Zr alloys, including non-lamellar microstructures, can be estimated accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement

Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗
DOE OSTI2021

Impact-induced twinning in a magnesium alloy under different stress conditions

We investigate deformation twinning in a highly textured magnesium alloy plate of Mg-3.1%Al-0.9%Zn-0.4%Mn under edge-on impact with in situ, synchrotron-based, ultrafast X-ray diffraction measurements, and corresponding stress states are simulated with the finite element method. Deformation twinning and its anisotropy under the triaxial stress condition are explained by a statistical analysis of resolved shear stress. The critical resolved shear stress criterion is applicable under complicated stress conditions induced by high strain rate impact loading. Three typical simple stress conditions are further explored as verification and application cases: uniaxial-stress, uniaxial-strain, and plane-stress. Extension twinning in the magnesium alloy is prone to occur for impact loading applied perpendicular to the crystallographic c-axis, regardless of the exact stress conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of scan strategy on porosity and microstructure of Ti-6Al-4V fabricated by electron beam powder bed fusion

Ti-6Al-4V is of high value in biomedical, aerospace, and industrial sectors and is ideal for additive manufacturing (AM) because of its poor machinability. In electron beam powder bed fusion of Ti-6Al-4V, a raster scan is commonly used and often results in porosity that is detrimental to the fatigue strength. In this paper, we demonstrate that spot melting can significantly minimize porosity but produces a coarser microstructure compared to raster melting. Process parameters can be optimized to achieve minimal porosity with tailored microstructures. A finite element method was used to rationalize these observations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2024

Convergent laser beam shapes: Unveiling the dynamics of Laser-induced elastic waves in composite materials

Overcoming the low signal-to-noise ratio (SNR) in laser ultrasonic testing of composite materials remains a significant challenge. Current efforts focus on enhancing SNR by inserting more energy into the material through temporal and/or spatial modulation of the laser beam. However, potential SNR improvements through wave convergence and wave energy manipulation have been overlooked. This paper addresses this gap by demonstrating the convergence of different wave types to a designated point and by showing the feasibility of directing absorbed laser energy into a specific wave type through spatial modulation of the laser beam. To achieve this, mathematical expressions for the convergent laser beams are derived. Various laser beam profiles are then introduced to the thermoelastic equations and solved using the finite element method. The sample under investigation is a transversely isotropic unidirectional carbon fiber reinforced plastic, characterized by anisotropic thermal expansion coefficients and thermal conductivities. Results reveal pronounced convergence of the intended wave type at the center due to laser beam shaping. This study showcases the ability to direct absorbed laser energy toward a specific wave type through spatial modulation of the laser beam and highlights the role of material anisotropy in energy focusing.

composite materials↗
DOE OSTI2024

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗