Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2025

Identification of $^{3}$He–$^{3}$H clusters in the $^{6}$Li+$^{89}$Y experiment using particle-$\gamma$ coincidence measurement

The 6 Li+ 89 Y experiment was performed to explore the reaction mechanism induced by a weakly bound nucleus 6 Li and its cluster configuration. Here, the particle-$\gamma$ coincidence method was used to identify the different reaction channels. The $\gamma$-rays coincident with 3 He/ 3 H indicate that the 3 H/ 3 He stripping reaction plays a significant role in the formation of Zr/Nb isotopes. The obtained results support the existence of a 3 He- 3 H cluster in 6 Li. Direct and sequential transfer reactions are adequately discussed, and the FRESCO code is used to perform precise finite-range cyclic redundancy check calculations. In the microscopic calculation, direct cluster transfer is more predominant than sequential transfer in 3 H transfer. However, the direct cluster transfer is of comparable magnitude to the sequential transfer in the 3 He transfer.

CRC calculations↗
DOE OSTI2021

Fidelity and entanglement entropy for infinite-order phase transitions

Here, we study the fidelity and the entanglement entropy for the ground states of quantum systems that have infinite-order quantum phase transitions. In particular, we consider the quantum O(2) model with a spin-S truncation, where there is an infinite-order Gaussian (IOG) transition for S = 1 and there are Berezinskii-Kosterlitz-Thouless (BKT) transitions for S ≥ 2. We show that the height of the peak in the fidelity susceptibility (χ F ) converges to a finite thermodynamic value as a power law of 1 / L for the IOG transition and as 1 / ln (L) for BKT transitions. The peak position of χ F resides inside the gapped phase for both the IOG transition and BKT transitions. On the other hand, the derivative of the block entanglement entropy with respect to the coupling constant $S^{'}_{vN}$ has a peak height that diverges as ln 2 (L) for S = 1 and ln 3 (L) for S ≥ 2 and can be used to locate both kinds of transitions accurately. We include higher-order corrections for finite-size scalings and obtain the value of the central charge consistent with c = 1 predicted by conformal field theory. The crossing point of χ F between different system sizes is at the IOG point for S = 1 but is inside the gapped phase for S ≥ 2, while those of $S^{'}_{vN}$ are at the phase-transition points for all S truncations. Our work elaborates on how to use the finite-size scaling of χ F or $S^{'}_{vN}$ to detect infinite-order quantum phase transitions and discusses the efficiency and accuracy of the two methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Bipolar Membranes with an Electrospun 3D Junction

Freestanding bipolar membranes (BPMs) with an entirely new and transformative morphology were fabricated and characterized. The key disruptive design element was a water-splitting/water-generation junction layer of finite thickness composed of interwoven and interlocking electrospun anion-exchange polymer (AEP) fibers and cation-exchange polymer (CEP) fibers. Methods were developed to sandwich the 3D junction layer between two pre-formed dense AEP and CEP films., with the closure of all interfiber voids, where the processing steps are easily incorporated into a commercial roll-to-roll membrane manufacturing scheme. A series of membranes were made using different polymers and junction layers to identify the optimum morphology and composition for water-splitting and water-generation applications. For water splitting, the 3D junction BPMs worked remarkably well with operating current densities at/above 1.0 A/cm2 (10X greater than commercial BPMs) at a transmembrane voltage drop of only 1.1 V. In water generation mode, the 3D junction BPMs operated at 0.5 A/cm2, a world-record current density. The membranes were found to exhibit outstanding durability and can be manufactured at scale for low cost. The high operating current densities and stable long-term operation are due to the 3D junction layer design, where there is a high interfacial area for water splitting or water generation reactions and where the interlocking fibers prevent delamination of the outer films. These new BPMs are ideal candidate materials for both existing water-splitting electrodialysis separations and new electrochemical processes, such as large-scale direct air capture of CO 2 , reactors for CO 2 reduction, self-hydrating fuel cells, and redox flow batteries.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Three-body resonances in the φ 4 theory

We study the properties of three-body resonances using a lattice complex scalar φ4 theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2025

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Impact of bulk-edge coupling on observation of anyonic braiding statistics in quantum Hall interferometers

Abstract Quantum Hall interferometers have been used to probe fractional charge and statistics of quasiparticles. We present measurements of a small Fabry–Perot interferometer in which the electrostatic coupling constants which affect interferometer behavior can be determined experimentally. Near the center of the ν = 1/3 state this device exhibits Aharonov–Bohm interference interrupted by a few discrete phase jumps, and Φ 0 oscillations at higher and lower magnetic fields, consistent with theoretical predictions for detection of anyonic statistics. We estimate the electrostatic parameters K I and K I L by two methods: using the ratio of oscillation periods in compressible versus incompressible regions, and from finite-bias conductance measurements. We find that the extracted K I and K I L can account for the deviation of the phase jumps from the theoretical anyonic phase θ a = 2 π /3. At integer states, we find that K I and K I L can account for the Aharonov–Bohm and Coulomb-dominated behavior of different edge states.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A System-Level Model for Estimating Residual Strain and Life of Nuclear Reactor Coolant System Components Under Connected-System-Thermal–Mechanical Boundary Conditions

Background: Environmental-assisted fatigue (EAF) is a major issue for the long-term survival of nuclear power plant fleets in the U.S. and worldwide. Multi-material welded regions (e.g., nozzles) and other high-stress regions of reactor coolant system (RCS) components are prone to EAF-related damage. Objective: The discussed work describes a system-level finite element (FE) model of RCS components of a pressurized water reactor (PWR). This is with the goal of predicting the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Methods: Herein, the FE model was developed considering system-level loading conditions (under connected system thermal-mechanical boundary conditions). Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the FE model results, the strain residuals, strain amplitudes and resulting fatigue lives of RCS components were predicted. Results: The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. Higher residual strain can lead to accelerated cyclic hardening of material and the associated effect of EAF. The simulated component-level strain profile (under realistic multi-axial-multi-physics loading cycle) can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests, which is a focus of future works. Conclusions: Despite similar geometry and material the RCS component can have significantly different strain profiles and resulting fatigue lives. The discussed approach can help to identify and prioritize the RCS components for conducting expensive nondestructive evaluation (NDE) inspections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗
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 Performance-Portable MultiGPU Implementation of 3D Euler Equations using ProtoX and IRIS

Computational scientists often face challenges when developing and optimizing code for high-performance computing (HPC), especially when trying to leverage GPUs. Given the heterogeneity of the nodes that comprise many modern HPC facilities, considerable demand exists for performance portable solutions for the core computational kernels used in many scientific computing libraries. In this work, we demonstrate a fourth-order finite volume method–based implementation of the Euler equations, which are an integral part of computational fluid dynamics. Our performance-portable multiGPU implementation for Euler equations uses ProtoX to generate kernels and IRIS for portability. ProtoX is a domain-specific language that uses a structured-grid partial differential equation library called Proto as its front end and the SPIRAL code generation system as its back end to generate optimized kernels for different architectures. Optimized kernels generated by ProtoX are orchestrated through the IRIS intelligent runtime system to provide portability. Two levels of optimizations within the IRIS runtime— directed acyclic graph fusion and task fusion—are explored to efficiently utilize computing resources in a multiGPU environment. Performance improvement through these optimizations is showcased by comparing the base ProtoX-IRIS implementation on AMD GPUs (Frontier node) and on NVIDIA GPUs (NVIDIA DGX-1).

Mankad, Het↗
DOE OSTI2021

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗
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↗