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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

KMT-2019-BLG-1339L: An M Dwarf with a Giant Planet or a Companion near the Planet/Brown Dwarf Boundary

We analyze KMT-2019-BLG-1339, a microlensing event with an obvious but incompletely resolved brief anomaly feature around the peak of the light curve. Although the origin of the anomaly is identified to be a companion to the lens with a low mass ratio q, the interpretation is subject to two different degeneracy types. The first type is the ambiguity in ρ, representing the angular source radius scaled to the angular radius of the Einstein ring, θ {sub E}, and the other is the s ↔ s {sup -1} degeneracy. The former type, “finite-source degeneracy,” causes ambiguities in both s and q, while the latter induces an ambiguity only in s. Here, s denotes the separation (in units of θ {sub E}) in projection between the lens components. We estimate that the lens components have masses (M{sub 1},M{sub 2})-(0.27{sub -0.15}{sup +0.36} M{sub ⊙},11{sub -7}{sup +16} M{sub J}) and -(0.48{sub -0.28}{sup +0.40} M{sub ⊙},1.3{sub -0.7}{sup +1.1} M{sub J}) according to the two solutions subject to the finite-source degeneracy, indicating that the lens comprises an M dwarf and a companion with a mass around the planet/brown dwarf boundary or a Jovian-mass planet. It is possible to lift the finite-source degeneracy by conducting future observations utilizing a high-resolution instrument because the relative lens–source proper motion predicted by the solutions is widely different.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2022

Mechanics of the Ring Tension Test (RTT): A Finite Element-based Investigation

The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted. Two promising test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Hot QCD phase diagram from holographic Einstein–Maxwell–Dilaton models

In this review, we provide an up-to-date account of quantitative holographic descriptions of the strongly coupled quark-gluon plasma (QGP) produced in heavy-ion collisions, based on the class of gauge-gravity Einstein-Maxwell-Dilaton (EMD) models. Holography is employed to tentatively map the QCD phase diagram at finite temperature onto a dual theory of charged, asymptotically AdS black holes in 5D. With a quantitative focus on the hot QCD phase diagram, the EMD models reviewed are adjusted to describe lattice results for the finite-temperature QCD equation of state, with 2+1 flavors and physical quark masses, at zero chemical potential and vanishing electromagnetic fields. The predictive power of EMD models is tested by quantitatively comparing their predictions for the hot QCD equation of state at nonzero baryon density and the corresponding state-of-the-art lattice QCD results. The shear and bulk viscosities predicted by these EMD models are also compared to the corresponding profiles favored by the latest phenomenological multistage models describing different heavy-ion data. We report preliminary results from a Bayesian analysis which provide systematic evidence that lattice results at finite temperature and zero baryon density strongly constrains the free parameters of EMD models. Remarkably, the set of parameters constrained by lattice results at zero chemical potential produces EMD models in quantitative agreement with lattice QCD results also at finite baryon density. We also review results for equilibrium and transport properties from magnetic EMD models, describing the QGP at finite temperatures and magnetic fields. Lastly, we provide a critical assessment of the main limitations and drawbacks of the holographic models reviewed in the present work, and point out some perspectives we believe are of fundamental importance for future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Modeling Isothermal Reduction of Iron Ore Pellet Using Finite Element Analysis Method: Experiments & Validation

Iron ore pellet reduction experiments were performed with pure hydrogen (H2) and mixtures with carbon monoxide (CO) at different ratios. For direct reduction processes that switch dynamically between reformed natural gas and hydrogen as the reductant, it is important to understand the effects of the transition on the oxide reduction kinetics to optimize the residence time of iron ore pellets in a shaft reactor. Hence, the reduction rates were studied by varying experimental parameters such as the temperature (800, 850 & 900 °C), reactant gas flow rate (100, 150 & 200 cm3/min), pellet size and composition of the reactant gas mixture. The rate of reduction was observed to increase with an increase in temperature and reactant gas flow rate, but it decreased with an increase in pellet size. SEM greyscale analysis was performed to analyze the porosity and phase composition of partially reduced pellets. The porosity of the pellets was observed to increase from 0.3 for unreacted pellet to 0.42 for a completely reduced pellet. Energy-dispersive X-ray spectroscopy (EDAX) analysis was performed to identify the phases observed in the SEM images. The fraction of iron phase was observed to increase from the shell region of the pellet to the core region with an increase in the degree of reduction. A 2D-axisymmetric numerical model was developed on COMSOL Multiphysics, and it was validated using the conversion (X) vs. time curves obtained from each experiment. The model was able to accurately predict the total time needed for the complete conversion of a single iron ore pellet for multiple experiments. Effects of changes in the porosity and tortuosity of the pellet on the model were also studied and the rate of reduction was observed to be sensitive to changes in both porosity and tortuosity. The SEM analysis and the model results show that tortuosity is higher for pellets reduced with H2 than for pellets reduced with H2-CO gas mixtures.

08 HYDROGEN↗
DOE OSTI2021

Magnetic properties of CrFeCoNi based high entropy alloys

Monte Carlo simulations are performed on three high entropy alloys: Cr0.25Fe0.25Co0.25Ni0.25, Cr0.2Fe0.2Co0.2Ni0.2Pd0.2, and Cr0.2Mn0.2Fe0.2Co0.2Ni0.2, with exchange interactions extracted from The ab initio Korringa-Kohn-Rostoker method combined with the coherent potential approximation calculations. Using finite size scaling analyses, we estimate the magnetic phase transition temperature for the four component alloy to be 108(2) K, and although the individual critical exponents are different from 3D Heisenberg universality class, the reduced exponent follows Suzuki weak universality. With the additional Palladium component, the transition temperature elevates to about 200 K. In contrast, we find no magnetic order for the five component alloy with Manganese at any finite temperatures.

Yin, Junqi↗
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 OSTI2021

Universal behavior of the bosonic metallic ground state in a two-dimensional superconductor

Anomalous metallic behavior, marked by a saturating finite resistivity much lower than the Drude estimate, has been observed in a wide range of two-dimensional superconductors. Utilizing the electrostatically gated LaAlO 3 /SrTiO 3 interface as a versatile platform for superconductor-metal quantum phase transitions, we probe variations in the gate, magnetic field, and temperature to construct a phase diagram crossing from superconductor, anomalous metal, vortex liquid, to the Drude metal state, combining longitudinal and Hall resistivity measurements. We find that the anomalous metal phases induced by gating and magnetic field, although differing in symmetry, are connected in the phase diagram and exhibit similar magnetic field response approaching zero temperature. Namely, within a finite regime of the anomalous metal state, the longitudinal resistivity linearly depends on the field while the Hall resistivity diminishes, indicating an emergent particle-hole symmetry. The universal behavior highlights the uniqueness of the quantum bosonic metallic state, distinct from bosonic insulators and vortex liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

CFD-DEM and PR-DNS studies of low-temperature densely packed beds

Over the past few decades, granular media is gaining attention as a viable option for heat transfer fluids (HTFs). Several research efforts are studying the use of particle-based heat transfer fluids in a wide variety of applications. With this motivation, the current work focusses on analyzing the different heat transfer mechanisms in low-temperature mono-sized densely packed granular media. To study the heat transfer behavior of granular media at different scales, the current work employs a two-way coupled computational strategy. The motion of particles is solved using the Discrete Element Method (DEM) and the interstitial air is solved using a Finite-Volume (CFD) approach. The Open-Source library CFDEM Coupling® is used in the current study to join the Finite Volume PISO solver of OpenFOAM® and the DEM solver of LIGGGHTS®. Typically, particle-particle contact conduction and particle-air convection are the most popular closure models. But recent research identified a different heat transfer phenomenon in packed beds that cannot be identified by conduction or convection models. While closure models were developed to implement this on a CFD-DEM framework, they did not capture the effect of intra-particulate thermal gradients on this phenomenon. Hence the current work also employs Particle-Resolved Direct Numerical Simulations (PR-DNS) to gain valuable insights allowing for the modification of existing models. A new closure model is then proposed here and is implemented in the CFD-DEM framework. This model provides key insights into the different heat transfer mechanism of packed beds.

42 ENGINEERING↗
DOE OSTI2025

Water Dielectric Function at Finite Wavelength and the Convergence of Its Large Wavelength Limit (Static Dielectric) by Fluctuation Formulas

Here we provide a general analysis of the longitudinal dielectric function ε l (k) via fluctuation formulas at finite wavelength k, including the static dielectric constant ε w = ε l (k → 0), and analyze the different sources of errors, deriving explicit formulas. Simulations with the SPC/E water model show that the convergence of fluctuation formulas to compute the static dielectric constant is slow and requires long simulations. The analysis of ε l (k) allows us to identify the long- and short-wavelength limits and provide a precise determination of the location of the poles as well as the zeros in the complex plane, which determine the short-distance form of the electrostatic (Coulomb) interaction among charged particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2023

Competing magnetic fluctuations and orders in a multiorbital model of doped SrCo 2 ⁢As 2

We revisit the intriguing magnetic behavior of the paradigmatic itinerant frustrated magnet Sr⁢Co 2 ⁢As 2 , which shows strong and competing magnetic fluctuations yet does not develop long-range magnetic order. By calculating the static spin susceptibility Χ⁡(q) within a realistic 16-orbital Hubbard-Hund model, we determine the leading instability to be ferromagnetic (FM). We then explore the effect of doping and calculate the critical Hubbard interaction strength U c that is required for the development of magnetic order. We find that U c decreases under electron doping and with increasing Hund's coupling J, but increases rapidly under hole doping. This suggests that magnetic order could possibly emerge under electron doping but not under hole doping, which agrees with experimental findings. We map out the leading magnetic instability as a function of doping and Hund's coupling and find several antiferromagnetic phases in addition to FM. We also quantify the degree of itinerant frustration in the model and resolve the contributions of different orbitals to the magnetic susceptibility. Lastly, we discuss the dynamic spin susceptibility Χ⁡(q,ω) at finite frequencies, where we recover the anisotropy of the peaks at Q π = (π,0) and (0,π) observed by inelastic neutron scattering that is associated with the phenomenon of itinerant magnetic frustration. By comparing results between theory and experiment, we conclude that the essential experimental features of doped SrCo 2 ⁢As 2 are well captured by an itinerant Hubbard-Hund multiorbital model if one considers a small shift of the chemical potential towards hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Skyrmion lattice manipulation with electric currents and thermal gradients in MnSi

The skyrmion lattice (SkL) in MnSi was studied using small-angle neutron scattering and under the influence of a radial electric current in a Corbino geometry. In response to the applied current, the SkL undergoes an angular reorientation with respect to the MnSi crystal lattice. The reorientation is nonmonotonic with increasing current, with the SkL rotating first in one direction and then the other. The SkL reorientation was studied at different sample locations and found to depend on the local current density as inferred from a finite-element analysis. Here, the nonmonotonic response indicates the presence of two competing effects on the SkL, most likely due to the presence of both radial electric and thermal currents. Such a scenario is supported by micromagnetic simulations, which show how these effects can act constructively or destructively to drive the SkL rotation, depending on the direction of the electric current. In addition, the simulations also suggest how the direction of the skyrmion flow may affect the SkL orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Multi-phase Wireless Power Transfer with High Power Density Inductive Coils for Electric Drone Charging

Wireless charging of unmanned ground vehi-cles and aircraft has been proposed to increase charging reliability and security, allow for autonomous functionality, and either reduce battery size or increase continuous flight time. This paper proposes a three-phase Litz wire primary and a two-phase PCB secondary for high secondary-side power density considering misalignment tolerances, surface and volumetric power density, and coil sizing. Elec-tromagnetic 3D finite element analysis (FEA) simulations are conducted to study variation in mutual inductance and coupling coefficient with different secondary coil sizes and number of turns, horizontal and vertical misalignment between the primary and secondary, and a combination of both for a set of designs with the same rated power. Results indicate that receiver coils of reduced size relative to the primary can deliver rated power for increased power density and higher tolerance to misalignment.

Gastineau, Lucas [University of Kentucky]↗
DOE OSTI2021

Nyx: A Massively Parallel AMR Code for Computational Cosmology

Nyx is a highly parallel, adaptive mesh, finite-volume N-body compressible hydrodynamics solver for cosmological simulations. It has been used to simulate different cosmological scenarios with a recent focus on the intergalactic medium and Lyman alpha forest. Together, Nyx, the compressible astrophysical simulation code, Castro, and the low Mach number code MAESTROeX, make up the AMReX-Astrophysics Suite of open-source, adaptive mesh, performance-portable astrophysical simulation codes. Other examples of cosmological simulation research codes include Enzo, Enzo-P/Cello, RAMSES, ART, FLASH, Cholla, as well as Gadget, Gasoline, Arepo, Gizmo, and SWIFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Advancements Toward ASME Nuclear Code Case for Compact Heat Exchangers

Our research team proposes to advance the state of the ASME section III code (nuclear service) for Compact Heat Exchangers (CHX). This work will improve the technical state of CHXs and lay the foundation necessary for these heat exchangers to be certified for use in nuclear service. During the course of this work, we will advance the understanding of the performance, integrity, and lifetime of the CHXs for use in any industrial application, making their use more attractive and accessible to the industry. We will do this by developing qualification and inspection procedures that utilize Non-destructive evaluation (NDE) and advanced in-service inspection techniques, with insight from the industrial utility leader EPRI. We have enlisted colleagues at MPR Associates (MPR), an elite nuclear code consulting firm, who are experts on the ASME section III code and who, with input from members of the ASME section III committee, will direct the testing and help develop a series of documents that define the rules and regulations for use of the CHX. Colleagues at North Carolina State University (NCSU) and Oregon State University (OSU) will conduct extensive tensile, creep, and fatigue experiments on diffusion bonded samples (manufactured by US-based Vacuum Process engineering) along with modeling using the elastic perfectly plastic assumptions and comprehensive full inelastic finite element analysis (FEA). This work will allow analysis by design and confidence in the strength of different internal structures. To ensure industry acceptance and long term confidence, team members at the University of Wisconsin–Madison (UW), University of Michigan (UM), Georgia Tech (GT), and the University of Idaho (UI) will extensively test prototypic heat exchangers manufactured by US-based manufactures CompRex and Vacuum Process Engineering (VPE), a leader in the development of advanced CHX. This testing will include the use of various working fluids (salt, sodium, helium, and sCO2) to evaluate operational issues as well as structural integrity under the most severe conditions. Post-test analysis of the tested CHXs coupled with pre/in-service/post NDE (ultrasonic and radiography) led by the Electric Power Research Institute (EPRI) will be incorporated into the development of the rules and regulations for their use in nuclear service.

42 ENGINEERING↗
DOE OSTI2023

Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design

Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.

36 MATERIALS SCIENCE↗