Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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

Fuel compact temperatures are a crucial factor in assessing irradiation performance of the Tristructural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in uncertainty of calculated fuel temperature is the gas gap uncertainty due to the nub-to-shell clearance because of fabrication error. A revision to the thermal model was made to examine the most probable offset position for four different dates of interest. The offset position options varied in magnitude and azimuthal direction of offset of the holder top and bottom position. The best-fit offset of the holder is estimated based on the minimum root mean square error of residuals for all operational thermocouples (TCs). From these results, the following conclusions were made: 1) During earlier cycles (162A – 164A) when numerous TCs were still operational, the best-fit offset distance varied in range [0.002 – 0.0035 in.] for both the top and bottom holder while offset azimuthal direction varied widely, especially for the holder bottom. 2) Holder offset led to slightly lower Capsule 1 average temperature, but wider temperature variation (lower minimum and higher peak fuel temperatures). 3) The maximum offset of 0.006 in. to the northwest direction for both the top and bottom holder yielded a much higher temperature than the original peak fuel temperature, which increased from 1422 to 1557 ? (or 135 ? increase). This could be one of the reasons for massive particle failure near the end of Cycle 165A. 4) Even though, the highest temperature at the TC-1-7 tip reached slightly higher than 1000?, the temperature along the TC wire reached as high as 1335? assuming a holder offset of 0.006 in. to the northwest. Consequently, at this temperature nickel from this TC wire could cause particle failures in Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

AGR-5/6/7 CAPSULE 1 THERMAL MODEL WITH OFFSET GAS GAPS

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tristructural isotropic (TRISO) fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a 3 D finite-element thermal model that was subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty caused by the graphite holder nub-to-capsule clearance. This thermal model was revised to examine the most probable graphite offset position for Capsule 1 for four different days during the irradiation. In the analysis, the offset distance and azimuthal direction at both the top and bottom of the holder were varied. The best-fit offset was estimated based on the minimum root mean square error (RMSE) of the residuals (measured minus calculated thermocouple [TC] temperatures) for the operational TCs in Capsule 1. Near the end of the irradiation, when the graphite had shrunk and the gas gaps had grown, a peak temperature was found with maximum offset along the TC lines, with nickel introducing a potential mechanism for fuel particle failures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

High energy modular bootstrap, global symmetries and defects

We derive Cardy-like formulas for the growth of operators in different sectors of unitary 2 dimensional CFT in the presence of topological defect lines by putting an upper and lower bound on the number of states with scaling dimension in the interval [Δ – δ, Δ + δ] for large Δ at fixed δ. Consequently we prove that given any unitary modular invariant 2D CFT symmetric under finite global symmetry G (acting faithfully), all the irreducible representations of G appear in the spectra of the untwisted sector; the growth of states is Cardy like and proportional to the “square” of the dimension of the irrep. In the Schwarzian limit, the result matches onto that of JT gravity with a bulk gauge theory. If the symmetry is non-anomalous, the result applies to any sector twisted by a group element. For c > 1, the statements are true for Virasoro primaries. Furthermore, the results are applicable to large c CFTs. We also extend our results for the continuous U(1) group.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The estimation of Poisson’s ratio by time-averaging and Cornu’s method for isotropic beams

Prior research suggests that direct static (e.g., uniaxial testing) and dynamic (e.g., ultrasonic wave propagation analysis) measurements differ in the estimation of Poisson’s ratio because anisotropies and heterogeneities in the sample material affect the two types of tests differently. Even assuming isotropic and homogeneous material properties, prior research further suggests that discrepancies between static/dynamic test results will exist because the error of the diagnostic techniques for the measurand are inherently different. Finally, thermodynamic effects are not present in static tests but can significantly affect dynamic test results. Given the potential for all these variables to produce discrepancies, it would be helpful to have the measurement of Poisson’s ratio obtainable from the same theory and experimental measurements by either static or dynamic testing methods. Our finite element calculations show that by combining time-averaged scanning digital holography with Cornu’s method, it is theoretically possible to estimate the effective Poisson’s ratio from the anticlastic contours at the antinode of the first out-of-plane bending mode shape. This is true regardless of frequency and therefore applicable for both static and dynamic measurements. Our results show that the estimate of Poisson’s ratio by Cornu’s method using data from simulations of mode shapes approaches the true value of Poisson’s ratio. Additionally, our research suggests that beam geometry and boundary conditions are fundamental factors limiting the convergence of the estimate of Poisson’s ratio to the true value of Poisson’s ratio regardless of performing a static or dynamic test.

42 ENGINEERING↗
DOE OSTI2023

Flat band separation and robust spin Berry curvature in bilayer kagome metals

Kagome materials have emerged as a setting for emergent electronic phenomena that encompass different aspects of symmetry and topology. It is debated whether the XV 6 Sn 6 kagome family (where X is a rare-earth element), a recently discovered family of bilayer kagome metals, hosts a topologically non-trivial ground state resulting from the opening of spin–orbit coupling gaps. These states would carry a finite spin Berry curvature, and topological surface states. Here we investigate the spin and electronic structure of the XV 6 Sn 6 kagome family. We obtain evidence for a finite spin Berry curvature contribution at the centre of the Brillouin zone, where the nearly flat band detaches from the dispersing Dirac band because of spin–orbit coupling. In addition, the spin Berry curvature is further investigated in the charge density wave regime of ScV 6 Sn 6 and it is found to be robust against the onset of the temperature-driven ordered phase. Utilizing the sensitivity of angle-resolved photoemission spectroscopy to the spin and orbital angular momentum, our work unveils the spin Berry curvature of topological kagome metals and helps to define its spectroscopic fingerprint.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
DOE OSTI2026

Scattering phase shift in quantum mechanics on quantum computers

Here, we investigate the feasibility of extracting infinite volume scattering phase shift on quantum computers in a simple one-dimensional quantum mechanical model, using the formalism established in the work by Guo and Gasparian [Phys. Rev. D 108, 074504 (2023)] that relates the integrated correlation functions for a trapped system to the infinite volume scattering phase shifts through a weighted integral. The system is first discretized in a finite box with periodic boundary conditions, and the formalism in real time is verified by employing a contact interaction potential with exact solutions. Quantum circuits are then designed and constructed to implement the formalism on current quantum computing architectures. To overcome the fast oscillatory behavior of the integrated correlation functions in real-time simulation, different methods of postdata analysis are proposed and discussed. Test results on IBM hardware show that good agreement can be achieved with two qubits, but complete failure ensues with three qubits due to two-qubit gate operation errors and thermal relaxation errors.

Guo, Peng [Dakota State Univ., Madison, SD (United↗
DOE OSTI2023

A Review of the Mechanical Properties of Materials Used in Nb 3 Sn Magnets for Particle Accelerators

Superconducting magnets experience significant thermo-mechanical loads throughout their life cycle. These are introduced by the electro-magnetic forces during powering, but also by the prestress applied in many magnet designs. Further to this, the large thermal excursion that components of different materials experience can generate significant internal forces. The loads are also experienced by the superconducting coils, whose critical current can decrease as a consequence of the applied strain. It is then crucial to predict the overall mechanical behavior and conservatively design a magnet, avoiding failure of the mechanical components and of the superconducting coils. Finite Element Analysis (FEA) is generally used to perform these tasks, but its results rely heavily on the material properties and models used. This is in particular true for the coil composite, which is simplified to allow reasonable model sizes in full magnet models. Here we present the state-of-art knowledge of the mechanical properties of the materials mostly used in superconducting magnet construction. We review elastic and plastic properties at room and cryogenic temperature, thermal contraction, and summarize the state-of-art failure criteria for these materials. Finally, the paper summarizes the present understanding of the mechanical behavior and limits of Nb3Sn coils. For the first time, an orthotropic failure criteria is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗
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

Influence of temperature, doping, and amorphization on the electronic structure and magnetic damping of iron

Hybrid magnonic quantum systems have drawn increased attention in recent years for coherent quantum information processing, but too large magnetic damping is a persistent concern when metallic magnets are used. Their intrinsic damping is largely determined by electron-magnon scattering induced by spin-orbit interactions. In the low scattering limit, damping is dominated by intra-band electronic transitions, which has been theoretically shown to be proportional to the electronic density of states at the Fermi level. In this work, we focus on body-centered-cubic iron as a paradigmatic ferromagnetic material. We comprehensively study its electronic structure using first-principles density functional theory simulations and account for finite lattice temperature, boron (B) doping, and structure amorphization. Our results indicate that temperature induced atomic disorder and amorphous atomic geometries only have a minor influence. Instead, boron doping noticeably decreases the density of states near the Fermi level with an optimal doping level of 6.25%. In addition, we show that this reduction varies significantly for different atomic geometries and report that the highest reduction correlates with a large magnetization of the material. Furthermore, this may suggest materials growth under external magnetic fields as a route to explore in experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Genetic Algorithm-Based Commutation Angle Control for Torque Ripple Mitigation in Switched Reluctance Motor Drives

This work addresses the application of the Genetic Algorithm (GA) technique to optimize the commutation angles of a 2 kW 8/6 switched reluctance machine (SRM). The primary goal is to reduce the well-known drawback of SRMs: the torque ripple. Firstly, the machine was modeled in Matlab /Simulink ® using lookup tables obtained via finite element method (FEM) simulations. Subsequently, the model was used to perform the GA routine aiming to find the optimal phase commutation angles that minimize the torque ripple factor. Notably, the torque performance of the SRM was significantly affected by the commutation angles during the search for the optimal solution. Afterwards, the GA results for four different operation points were verified experimentally through a developed drive platform with digital signal processor-based (DSP) control and an asymmetric bridge converter. As showed by the experiments, the proposed approach was suitable to reduce the torque ripple by more than 50% for one of the evaluated operating points. Furthermore, it was confirmed that the torque ripple mitigation led to acoustic noise improvement.

42 ENGINEERING↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2026

Assessment of Vibration and Shock Loading Environments on Steels Used in Railroad Applications Considering Hydrogen Embrittlement

Studies that examine the combined effects of hydrogen and rail loading environments in the context of hydrogen embrittlement of steels have not been performed. Large-scale hydrogen storage requires either gaseous hydrogen stored in pressure vessels and tubing or liquid hydrogen in cryogenic, insulated tanks. Since the focus of this study is on steels, we examine the influence of hydrogen on stainless steel used for high pressure tubing and cryogenic tanks and Cr-Mo steel liners used in gaseous pressure vessels. Literature studies show that despite hydrogen degradation, these steels can be used effectively and safely if the stresses are managed. The shock and vibration environment on railcars present unique environments where stresses and loading rates can be high. A literature review was conducted to understand the accelerations experienced in rail by examining instrumented railcars of nuclear fuel casks. Hydrogen studies were examined that focused on steel pressure hardware with an emphasis on fatigue, fast loading rates, and low temperatures. A simplified fatigue assessment was performed to examine design cycles for three different hydrogen components: stainless steel tubing, stainless steel cryogenic tanks, and gaseous pressure vessels. The assessments were idealized and used bounding cases of accelerations in both normal rail conditions and coupling events. In normal shock and vibration environments when accelerations were less than 2 g, infinite design life was calculated. In extreme shock or coupling environments, finite design life was calculated which was shown to be dependent on unsupported lengths of pressure components. It was shown that adjusting the unsupported length can reduce the bending stress thereby increasing design life.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling and characterizing the elastodynamic response of octet-truss lattice structures using resonance techniques

Characterization of additively manufactured materials and structures is an ongoing effort for the advanced manufacturing community. This article will investigate an approach using resonant spectral ultrasound (RUS) to measure the effective elastic constants of an Octet Truss lattice and apply the results to continuum bases models representing the lattice regions in different structures. The study is focused on simple lattices structures fabricated from Ti5553 using a laser powder bed fusion process (LBPF). Solid and lattice samples are measured to determine RUS estimates for the elastic properties of the bulk material and the effective elastic properties of the lattice structure. The estimated elastic properties are then incorporated into 3d finite element models using a continuum approximation of the physical AM parts. Comparisons show good agreement between the experimentally measured eigen frequencies of the Ti5553 LBPF parts and the eigen frequencies calculated using the continuum approximations based on the effective RUS elastic properties. The current results suggest there is additional physics and geometrical effects that are accounted for in the RUS continuum approximation of the lattice that are not captured in the full 3d finite element model of the parts utilizing only the base material properties.

36 MATERIALS SCIENCE↗