Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Melting”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal energy mitigation and toroidal peaking effects in JET disruptions

Previous investigations on JET suggest half or less of plasma stored thermal energy W th is radiated (f rad,th ≲ 0.5) using either massive gas injection (MGI) or shattered pellet injection (SPI) disruption mitigation. We investigate whether the apparent incomplete f rad,th is explained by radiation peaking near the injection plume. High toroidal peaking throughout the pre-thermal quench is found in argon–deuterium MGI on JET, with typically >3× higher radiation near the injector than toroidally distant. Previously unexplained toroidal bolometry measurements in neon–deuterium SPI are reproduced with similar peaking using the Emis3D radiation analysis code. These observations align with results from Alcator C-Mod and KSTAR. This peaking is not captured by previous JET studies that found poor thermal mitigation. Two sets of neon–deuterium SPI and two sets of argon–deuterium MGI are analyzed using Emis3D. In SPI, f rad,th rises from no-plume estimates of 0.31 and 0.66 to lower bounds of 0.84 and 0.92, respectively, and f rad,th ∼ 1 is possible. In MGI, the toroidal spread of the peaking feature is poorly constrained. f rad,th up to 0.85 and 0.65 are possible using the largest possible spread, increasing from 0.42 and 0.28, although f rad,th ∼ 1 does not appear to be reached. Revised mitigation estimates on JET suggest a lower melt risk to the divertor in mitigated disruptions on ITER and SPARC than previously thought. However, peaking near injectors could increase flash melting risk on nearby plasma facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗
DOE OSTI2026

Chemical bond and phase stability of Ga-doped Sm2Fe17Cx magnet

Sm2Fe17C3 phase (2:17) is metastable and exhibits excellent intrinsic hard magnetic properties. Doping elements such as Ga facilitate the formation of a single-phase 2:17 structure in arc-melted Sm2Fe17Cx alloys, which opens a promising route for fabricating fully dense bulk Sm2Fe17Cx magnets via high-temperature techniques such as melting and sintering. First-principles electronic structure calculation indicates that Ga prefers to partially replace Fe at the 9d and 18h crystallographic sites in Sm2Fe17C3 and Sm2Fe17, respectively. This difference in site preference is attributed to the distinct chemical environments surrounding the Fe atoms in the two compounds. Ga substitution favors the Sm–Ga bonding formation while avoiding Ga–C interactions. Doped Ga atoms result in more negative formation energy in Sm2(Fe, Ga)17C3, indicating improved structural stability. Crystal Orbital Hamilton Population analysis reveals that carbon insertion weakens the bonding of Sm-Fe (18h) and Sm-Fe (18f) in Sm2Fe17C3. Ga doping facilitates electron redistribution across chemical bonds, thereby reinforcing Fe(18h)–Sm and Fe(18f)–Sm interactions and stabilizing the carbon-centered octahedral local structure. This synergistic effect contributes significantly to the observed enhancement in phase stability of Sm2(Fe, Ga)17Cx. These findings suggest that chemical bond engineering through the selective doping of Ga can enhance phase stability and facilitate the synthesis of Sm2Fe17C3, providing a viable strategy for developing advanced magnets.

Liu, Xubo [Critical Materials Innovation Hub, Divi↗
DOE OSTI2020

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core–mantle boundary. In this paper, we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO 3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth’s core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

58 GEOSCIENCES↗
DOE OSTI2021

Metastable–solid phase diagrams derived from polymorphic solidification kinetics

Nonequilibrium processes during solidification can lead to kinetic stabilization of metastable crystal phases. A general framework for predicting the solidification conditions that lead to metastable-phase growth is developed and applied to a model face-centered cubic (fcc) metal that undergoes phase transitions to the body-centered cubic (bcc) as well as the hexagonal close-packed phases at high temperatures and pressures. Large-scale molecular dynamics simulations of ultrarapid freezing show that bcc nucleates and grows well outside of the region of its thermodynamic stability. An extensive study of crystal–liquid equilibria confirms that at any given pressure, there is a multitude of metastable solid phases that can coexist with the liquid phase. We define for every crystal phase, a solid cluster in liquid (SCL) basin, which contains all solid clusters of that phase coexisting with the liquid. Additionally, a rigorous methodology is developed that allows for practical calculations of nucleation rates into arbitrary SCL basins from the undercooled melt. It is demonstrated that at large undercoolings, phase selections made during the nucleation stage can be undone by kinetic instabilities amid the growth stage. On these bases, a solidification–kinetic phase diagram is drawn for the model fcc system that delimits the conditions for macroscopic grains of metastable bcc phase to grow from the melt. We conclude with a study of unconventional interfacial kinetics at special interfaces, which can bring about heterogeneous multiphase crystal growth. A first-order interfacial phase transformation accompanied by a growth-mode transition is examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Direct observation of the collective modes of the charge density wave in the kagome metal CsV 3 Sb 5

A recently discovered group of kagome metals AV 3 Sb 5 (A = K, Rb, Cs) exhibit a variety of intertwined unconventional electronic phases, which emerge from a puzzling charge density wave phase. Understanding of this charge-ordered parent phase is crucial for deciphering the entire phase diagram. However, the mechanism of the charge density wave is still controversial, and its primary source of fluctuations—the collective modes—has not been experimentally observed. Here, we use ultrashort laser pulses to melt the charge order in CsV 3 Sb 5 and record the resulting dynamics using femtosecond angle-resolved photoemission. We resolve the melting time of the charge order and directly observe its amplitude mode, imposing a fundamental limit for the fastest possible lattice rearrangement time. These observations together with ab initio calculations provide clear evidence for a structural rather than electronic mechanism of the charge density wave. Our findings pave the way for a better understanding of the unconventional phases hosted on the kagome lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unveiling the effect of Ni on the formation and structure of Earth’s inner core

Ni is the second most abundant element in the Earth’s core. Yet, its effects on the inner core’s structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe’s melting point at inner core pressures. The melting temperature of Ni is shown to be 700 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe’s crystallization at core pressures. These results suggest that Ni may substantially impact the structure and formation process of the solid inner core.

58 GEOSCIENCES↗
DOE OSTI2025

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗
DOE OSTI2026

Multiscale fatigue crack initiation in hierarchical additively manufactured alloys

Bioinspired hierarchical microstructures offer a route toward engineered fatigue resistance in additively manufactured alloys. However, it remains unclear how discrete structural constituents independently govern damage accumulation, particularly during the critical fatigue initiation regime where short cracks strongly interact with local microstructure. Here, we investigate multiscale fatigue initiation in a dual-phase, nanolamellar AlCoCrFeNi 2.1 high-entropy alloy. By comparing microscale specimens that isolate the nanolamellar structure against macroscale specimens containing the full melt-pool architecture, we identify size-dependent fatigue initiation mechanisms. We find that failure is dictated by nanolamellar interfaces at the microscale, whereas mesoscale melt pool boundaries serve to initiate fatigue at the macroscale. This mechanistic shift is accompanied by a transition from macroscale quasi-brittle failure to microscale plasticity-driven crack extension. Our results provide a physical framework for understanding how structural hierarchy governs the transition from discrete microstructural deformation to continuum fatigue fracture behavior, informing the design of damage-tolerant, additively manufactured alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hot Hydrogen Testing of W-Coated UN Kernels in a Mo30W Matrix

Ceramic uranium mononitride (UN) is being considered as a reactor fuel for nuclear thermal propulsion. To avoid or reduce the dissociation of UN at the high temperatures needed, embedding it in a metallic matrix (cermet) has been proposed. To assess the viability of this concept, hot hydrogen testing of tungsten-coated UN kernels embedded in a Mo-30 wt% W (Mo30W) alloy matrix has been performed at temperatures from 1800°C to 2300°C. Both the isolated kernels and kernels consolidated by spark plasma sintering in the Mo30W matrix were tested. In addition to direct observations and mass loss measurements, the samples were analyzed by X-ray diffraction (XRD) and scanning electron microscopy (SEM)/energy dispersive X-ray spectroscopy (EDS) after each run. The decomposition of UN started at 1800°C despite the coating and matrix, and increased at 2000°C. Uranium seeped through the tungsten grain boundaries of the coating at all temperatures. The consolidated sample expanded irregularly at 2000°C through the formation of voids, and SEM/EDS analysis showed uranium-containing veins in the matrix consisting of U 2 Mo according to the XRD data. The observed pore generation at 2000°C was explained by the formation of water vapor from residual oxides and diffused hydrogen. At 2200°C and above, both the kernels and the consolidated samples melted through the formation of uranium or low–melting point uranium-molybdenum alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

S 8 ring content in liquid sulfur determined from x-ray diffraction

The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .

lambda transition↗
DOE OSTI2020

Degenerate doping in B-Ga 2 O 3 Single Crystals through Hf-doping

n-type conductivity of β-Ga 2 O 3 grown from the melt is typically achieved using Sn and Si. In this paper, we experimentally and computationally investigate Hf doping of β-Ga 2 O 3 single crystals using UV–vis-NIR absorption and Hall effect measurements and hybrid functional calculations. Unintentionally-doped and Hf-doped samples with a nominal concentration of 0.5at% were grown from the melt using vertical gradient freeze and Czochralski method in mixed Ar + O 2 atmosphere. We demonstrate Hf dopants, predicted to incorporate on the octahedral GaII site as a shallow donor, achieve degenerate doping in β-Ga 2 O 3 with a measured electron concentration ~2 × 10 19 cm -3 , mobility 80–65 cm 2 V -1 s -1 , and resistivity down to 5 mΩ cm in our samples. The concentration of Hf was measured to be 1.3 × 10 19 atoms cm -3 using glow discharge mass spectroscopy on doped samples, confirming Hf to be the cause of n-type conductivity (electron concentration ~2 × 10 19 cm -3 ).

36 MATERIALS SCIENCE↗
DOE OSTI2024

The nontrivial effects of annealing on superconducting properties of Nb single crystals

The effect of annealing on the superconducting properties of niobium single crystals was studied using optical, magnetic, and scanning tunneling microscopy (STM) methods. Pieces of the same crystal boule were studied before and after the annealing at 800 ${^\circ}\textrm{C}$, 1400 ${^\circ}\textrm{C}$, and near the melting point of niobium (2477 ${^\circ}\textrm{C}$). The initial samples had a high hydrogen content and low-temperature imaging revealed large hydrides (hundreds of micrometers) appearing below 190 K. The formation of these large precipitates is already completely suppressed by annealing at 800 ${^\circ}\textrm{C}$. However, the overall superconducting properties of the annealed samples did not improve and, in fact, worsened. In particular, the superconducting transition temperature decreased, the upper critical field increased, and the pinning strength increased. In the STM study, the sample was annealed initially at 400 ${^\circ}\textrm{C}$, measured, annealed at 1700 ${^\circ}\textrm{C}$, and measured again. The STM revealed a ‘dirty’ superconducting gap with a significant spatial variation in tunneling conductance after annealing at 400 ${^\circ}\textrm{C}$. The clean gap was recovered after annealing at 1700 ${^\circ}\textrm{C}$. This is likely due to oxygen redistribution near the surface, which is always covered by oxide layers in as-grown crystals. Our results indicate that vacuum annealing at least up to 1400 ${^\circ}\textrm{C}$, while removing a large percentage of hydrogen, introduces additional nanosized defects, likely hydride precipitates, that act as efficient pair-breaking and pinning centers. The dimensionless scattering rate is estimated to have increased from $\Gamma = 0.2$ to about $\Gamma = 0.4$ after annealing at 1400 ${^\circ}\textrm{C}$. These results on single crystals differ drastically from those obtained in polycrystalline bulk niobium (i.e. cut from superconducting radio-frequency cavities), where annealing is known to have a significant positive effect that is attributed to the improvement of the crystalline structure masking the more subtle influence of the hydrides.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Plasma steering to avoid disruptions in ITER and tokamak power plants

Steering tokamak plasmas is commonly viewed as a way to avoid disruptions and runaway electrons. Plasma steering sounds as safe as driving to work but will be shown to more closely resemble driving at high speed through a dense fog on an icy road. The long time required to terminate an ITER discharge compared to time over which dangers can be foreseen is analogous to driving in a dense fog. Here, the difficulty of regaining plasma control if it is lost resembles driving on an icy road. Disruptions and runaways are associated with three issues—a solution to one tends to complicate the solution to the other two: loss of plasma position control, excessive heat deposition, and wall melting due to runaway electrons. All three risks must be addressed for ITER to achieve its mission and essentially eliminated before tokamak power plants can be deployed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Rapid solidified ductile Cu-Al-Mn ribbon and its elastocaloric potential

Abstract Cu-Al-Mn alloys display martensitic transformation over a wide range of temperatures. In addition to low cost, this alloy is known for its low transformation stress with reasonable latent heat favoring elastocaloric applications. However, the ductility of Cu-Al-Mn can be limited owing to ordering and intergranular fracture. Through rapid solidification by melt spinning, we show that Cu-Al-Mn ribbon can be made highly ductile (greater than 8% tensile strain in the as-spun state and 10% tensile strain after heat treatment). The ductility of the melt-spun ribbon is related to the suppression of L2 1 ordering that is characterized through magnetic property measurement. Heat treatment of the ribbon promotes bamboo grain formation, and the latent heat is increased to 6.4 J g −1 . Under tensile conditions, we show that the ribbon exhibited about 4 °C temperature change (4.4 °C on heating and 4.2 °C on cooling from 6.3% strain).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES↗