Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element Method”

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 1,153 records · Page 64

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↗

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↗

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↗

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↗

Dynamic energy absorption characteristics of hollow microlattice structures

Hollow microlattice structures are promising candidates for advanced energy absorption and their characteristics under dynamic crushing are explored. The energy absorption can be significantly enhanced by inertial stabilization, shock wave effect and strain rate hardening effect. In this paper we combine theoretical analysis and comprehensive finite element method simulation to decouple the three effects, and then obtain a simple model to predict the overall dynamic effects of hollow microlattice structures. Inertial stabilization originates from the suppression of sudden crushing of the microlattice and its contribution scales with the crushing speed, v. Shock wave effect comes from the discontinuity across the plastic shock wave front during dynamic loading and its contribution scales with e. The strain rate effect increases the effective yield strength upon dynamic deformation and increases the energy absorption density. A mechanism map is established that illustrates the dominance of these three dynamic effects at a range of crushing speeds. Compared with quasi-static loading, the energy absorption capacity a dynamic loading of 250 m/s can be enhanced by an order of magnitude. The study may shed useful insight on designing and optimizing the energy absorption performance of hollow microlattice structures under various dynamic loads. (C) 2014 Elsevier Ltd. All rights reserved.

Liu, YL↗

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↗

Bond-associated deformation gradients for peridynamic correspondence model

Non-ordinary state-based peridynamic correspondence material model is known to have issues with material instability, i.e. the existence of zero-energy modes, due to non-unique mapping between deformation states and force states via the conventional peridynamic deformation gradient. In this paper, an alternative approach in which the deformation gradient hence force state are computed specifically for each individual bond is proposed to eliminate the material instability. Bond-associated deformation gradient is calculated based on deformation states of material points within an individual bond’s proximity, termed here as the bond-associated family, rather than a material point’s whole family. This bond-associated deformation gradient can better represent the force state of each individual bond from the deformation states of material points in its proximity, and hence inherently resolves issues of material instability in the conventional correspondence material model. Parametric study on bond-associated horizon size indicates that the optimal size should be no less than the material point’s horizon size but smaller than two times of that value. Comparisons against reference solutions using finite element method establish the validity and accuracy of the proposed formulation.

42 ENGINEERING↗

Phase-field modeling of the beta to omega phase transformation in Zr–Nb alloys

A three-dimensional elastoplastic phase-field model is developed, using the Finite Element Method (FEM), for modeling the athermal beta to omega phase transformation in Zr–Nb alloys by including plastic deformation and strain hardening of the material. The microstructure evolution during athermal transformation as well as under different stress states, e.g. uni-axial tensile and compressive, bi-axial tensile and compressive, shear and tri-axial loadings, is studied. The effects of plasticity, stress states and the stress loading direction on the microstructure evolution as well as on the mechanical properties are studied. The input data corresponding to a Zr – 8 at.% Nb alloy are acquired from experimental studies as well as by using the CALPHAD method. Our simulations show that the four different omega variants grow as ellipsoidal shaped particles. Our results show that due to stress relaxation, the athermal phase transformation occurs slightly more readily in the presence of plasticity compared to that in its absence. The evolution of omega phase is different under different stress states, which leads to the differences in the mechanical properties of the material. The variant selection mechanism, i.e. formation of different variants under different stress loading directions, is also nicely captured by our model.

36 MATERIALS SCIENCE↗

Coupled strain-induced alpha to omega phase transformation and plastic flow in zirconium under high pressure torsion in a rotational diamond anvil cell

Strain-induced α → ω phase transformation (PT) in the zirconium (Zr) sample under compression and plastic shear in a rotational diamond anvil cell (RDAC) is investigated using the finite element method (FEM). The fields of the volume fraction of the ω phase, all components of the stress tensor, and plastic strain are presented. Before torsion, PT barely occurs. During torsion under a fixed applied force, PT initiates at the center of the sample, where the pressure first reaches the minimum pressure for strain-induced α → ω PT, p d ε , and propagates from the center to the periphery and from the symmetry plane to the contact surface. Salient increase of the shear friction stress and pressure at the center of a sample, so-called pressure self-multiplication effect observed experimentally for some other materials, is predicted here for Zr. It is caused by much higher yield strength of the phase in comparison with the phase. Except at the very center of a sample, the total contact friction stress is equal to the yield strength in shear of the mixture of phases and the plastic sliding occurs there. Due to the reduction in sample thickness and radial material flow during torsion, the ω phase can be observed in the region where pressure is lower than p d ε , which may lead to misinterpretation of the experimental data for p d ε . For the same applied force, torsion drastically promotes PT in comparison with the compression without torsion. However, the PT process in RDAC is far from optimal: (a) due to the pressure self-multiplication effect, the pressure in the transformed region is much higher than that required for PT; (b) the region in which PT occurs is limited by the pressure p d ε and cannot be expanded by increasing a shear under a fixed force; and (c) the significant reduction in thickness during torsion reduces the total mass of the high-pressure phase. These drawbacks can be overcome by placing a sample within a strong gasket with an optimized geometry. It is shown that, due to strong pressure heterogeneity, characterization of α → ω and α → β PTs based on the averaged pressure contains large errors. Furthermore, the obtained results, in addition to providing an improved understanding of the strain-induced PTs, may be beneficial for the optimum design of experiments and the extraction of material parameters, as well as optimization and control of PTs by varying the geometry and loading conditions.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

A transmissive, spectrum-splitting concentrating photovoltaic module for hybrid photovoltaic-solar thermal energy conversion

Here, a spectrum splitting, transmissive concentrating photovoltaic (tCPV) module is proposed and designed for a hybrid photovoltaic-solar thermal(PV/T) system. This module enables the system to more fully utilize the full spectrum of incoming sunlight. By utilizing III-V triple junction solar cells with bandgaps of 2.1eV/1.7eV/1.4eVin the module, ultraviolet (UV) and visible light will be absorbed and converted to electricity, while infrared (IR) light will pass through and be captured by a solar thermal receiver and stored as heat. The stored heat energy may be dispatched as electricity or process heat, as needed.Numerical analysis based on the Finite Elemental Method (FEM), ray tracing, circuit analysis, device modelling, and more is employed to design this tCPV module. According to these simulations, the tCPV module can perform with overall power conversion efficiency exceeding 43.5% for above bandgap (in-band)light under a standard AM1.5D solar spectrum, under an average concentration ratio of 400 suns.Passive and active cooling methods, keeping cells below 110°C, are also investigated and discussed, indicating that a transparent active cooling design could improve the CPV module efficiency by around 1%(absolute), relative to a passive design,by reducing the maximum cell working temperature by ~16°C. Furthermore, cost analysis shows that installation cost of around 1.9/W–2.2/W could be reached for the tCPV based PV/T system,which shows a competitive economic advantage compared to a more conventional PV with battery system.

14 SOLAR ENERGY↗

Multilevel sequential Monte Carlo samplers

Here, we study the approximation of expectations w.r.t. probability distributions associated to the solution of partial differential equations (PDEs); this scenario appears routinely in Bayesian inverse problems. In practice, one often has to solve the associated PDE numerically, using, for instance finite element methods and leading to a discretisation bias, with the step-size level h L . In addition, the expectation cannot be computed analytically and one often resorts to Monte Carlo methods. In the context of this problem, it is known that the introduction of the multilevel Monte Carlo (MLMC) method can reduce the amount of computational effort to estimate expectations, for a given level of error. This is achieved via a telescoping identity associated to a Monte Carlo approximation of a sequence of probability distributions with discretisation levels ${\infty}$ >h 0 >h 1 ...>h L . In many practical problems of interest, one cannot achieve an i.i.d. sampling of the associated sequence of probability distributions. A sequential Monte Carlo (SMC) version of the MLMC method is introduced to deal with this problem. In conclusion, it is shown that under appropriate assumptions, the attractive property of a reduction of the amount of computational effort to estimate expectations, for a given level of error, can be maintained within the SMC context.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

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↗

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗