Search NASA⌕ Search

SEARCH · Search NASA

Results for “discontinuous phase transition”

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.

116 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

ULTIMATE Phase I project: Development of Niobium-based alloys for Turbine Applications

Current Ni-based alloys used in turbine blade applications are operating at 1100°C which is approximately 90% of the solidus of Ni-based alloys. Further increases in temperature can be achieved only through the use of alloys with higher solidus temperatures such as refractory alloys which include Mo, Nb, W, Ta alloys. Niobium has a unique combination of physical properties that include a high melting point (2468°C), lowest density (8.57 g/cc), and the lowest modulus amongst refractory alloys, a good thermal conductivity that increases from 45 W/mK at room temperature to 62 W/mK at 1093°C, good ductility with fabricability using traditional techniques, high strength, good creep resistance, and general chemical inactivity. Hence niobium alloys offer many advantages for use at the desired temperature of 1300°C. One of the major deterrents to the use of niobium and its alloys at high temperatures is achieving good oxidation resistance. In addition, achieving balanced properties of room temperature strength, ductility, fracture toughness, high temperature strength, and creep resistance required for 1300°C operation at a density of 9 g/cm 3 is extremely challenging since the addition of W which is a potent strengthener increases the density of the alloy. This project seeks to enable the development of a tri-layered turbine blade that can operate continuously at 1300°C with the core of the turbine blade providing good creep resistance, the second layer providing improved oxidation resistance over the core layer, and the third layer being an environmental barrier providing oxidation resistance. The objective of this Phase 1 project was to use computational modeling and advanced characterization tools to develop two classes of Nb alloys which could be used as the core layer and the transition layer. The first class of alloys developed in this project was a Nb-alloy designed to serve as the turbine blade core with excellent room temperature strength, ductility, and high temperature strength and creep resistance required for 1300°C operation. These alloys were designed to contain sufficient solute solution strengthening elements (W, Mo) but constrained by the density of alloy, along with the addition of elements such as Zr, Hf, Ta, C, and N to achieve a combination of primary and secondary carbide precipitation. A total of 38 “creep-resistant” alloys were designed and cast during the duration of the project. Processing techniques were developed to keep the oxygen contents as low as possible with typical oxygen contents less than 250 ppm. One alloy with a density of < 9.5 g/cc was successful in meeting the Phase 2 intermediate project mechanical property milestone requirements of room temperature ductility greater than 1.0%, 1200°C creep strain of less than 3% at 150 MPa and 100 hours in vacuum, and with solidus temperature greater than 1500ᵒC. The second class of alloys was designed to be a Nb-rich alloy with improved oxidation resistance when compared to the core layer and was designed specifically to be compatible with the core layer and the outer environmental barrier coating. This Nb- alloy will be specifically designed to be microstructurally stable at these temperatures when in contact with the core alloy and provide protection against catastrophic failure of the barrier coating. Two alloys were cast and processed but further development was discontinued to focus on the development of the creep resistant alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

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

Learning protocols for the fast and efficient control of active matter

Exact analytic calculation shows that optimal control protocols for passive molecular systems often involve rapid variations and discontinuities. However, similar analytic baselines are not generally available for active-matter systems, because it is more difficult to treat active systems exactly. Here we use machine learning to derive efficient control protocols for active-matter systems, and find that they are characterized by sharp features similar to those seen in passive systems. We show that it is possible to learn protocols that effect fast and efficient state-to-state transformations in simulation models of active particles by encoding the protocol in the form of a neural network. We use evolutionary methods to identify protocols that take active particles from one steady state to another, as quickly as possible or with as little energy expended as possible. Our results show that protocols identified by a flexible neural-network ansatz, which allows the optimization of multiple control parameters and the emergence of sharp features, are more efficient than protocols derived recently by constrained analytical methods. Our learning scheme is straightforward to use in experiment, suggesting a way of designing protocols for the efficient manipulation of active matter in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRSDate not supplied

Zircon, Baddeleyite, and Reidite Found in Ries Crater Suevite

Introduction: Impact events can generate superheated impact melts and even vapor [1]. In the past decade, evidence of the high temperatures and high pressures of the impact process has been found in the impact melt from terrestrial craters (e.g., Mistastin Lake [2], Meteor Crater [3]) and even on the Moon [4–5]. Often, these studies involve zircon or zirconium-bearing phases. Zircon is a particularly useful mineral, due to its robustness and durability against weathering. It is used principally for chronology, but it has a multitude of geologic applications, including geothermometry and fingerprinting magma sources. Here we focus on its capability of recording impact conditions. Grains of zircon (ZrSiO4) are converted at high temperature and/or pressure during an impact event to reidite (a high-pressure polymorph of ZrSiO4) or to tetragonal-ZrO2 + SiO2 [6]. Reidite and tertragonal-ZrO2 leave identifying relicts in the rocks, markers of the high pressures and/or temperatures these rocks underwent. Ries Crater is a 26-km-diameter peak ring crater in southern Germany that formed approximately 15 Ma [7–9]. Ries is the type locality of the polymict impact breccia known as suevite. Here, we report a microanalytical study of suevite from the Ries impact structure. These analyses are used to inform our understanding of the pressure and temperature conditions involved in creating a polymict impact breccia such as suevite in a peak ring crater. Sample Description: The sample, denoted 16RS08, originates from Otting Quarry, at 48.8777° N, 10.7921° E, approximately 17 km from the center of the Ries Crater in Germany (approximately 4 km outside the crater rim). The material recorded in sample 16RS08 is therefore considered an outer suevite. The outer suevite is a discontinuous layer of polymict impact breccia that occurs outside the central ring of Ries, up to 22 km from the center of the crater [9]. Lithic clasts in the outer suevite consist primarily of crystalline basement rocks (gneiss, granite, amphibolite), with less than 5% of lithic clasts being overlying sedimentary rocks (limestone, sandstone, shale) [9]. The thin section studied contains variably shocked lithic and mineral clasts, impact glass, and interstitial minerals that make up the matrix of the breccia (Fig. 1) [9]. Methods: We used a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona to obtain 15 elemental X-ray maps of 16RS08. We next used the JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC) to obtain electron backscatter diffraction (EBSD) maps and energy dispersive X-ray spectroscopy (EDS) maps of select portions of the section. The EBSD data were collected under beam conditions of 20 kV, and ~9 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We used the elemental X-ray maps to identify the areas of interest in the section, particularly phosphates and Zr-bearing grains. These areas of interest were then targeted for follow up EBSD and EDS analyses. We have identified, through combined EBSD and EDS analysis, the presence of zircon, reidite, and monoclinic-ZrO2 (baddeleyite) in 16RS08. In 16RS08, we have found singular grains of zircon, zircon with a vermicular baddeleyite halo (Fig. 2a), and granular zircon with reidite (Fig. 2b). The different Zr-rich phases and their corresponding textures signifies that this sample underwent a broad spectrum of pressure and temperature conditions during the impact event. For example, zircons surrounded by a vermicular baddeleyite and SiO2 intergrowth (i.e., Fig. 2a) have been shown to preserve evidence of the extremely high temperatures of impact melt, upwards of 2370 °C [2, 6]. Similarly, reidite and granular zircon (i.e., Fig. 2b) have been shown to preserve evidence of high pressure, as the transition to reidite occurs >30 GPa [3, 6]. Future Work: Next, we will process the EBSD and EDS data for these Zr-rich grains, specifically looking for indicators of cubic- or tetragonal-ZrO2 in the baddeleyite remnants and of shock-precursors to the reidite. To further inform our work, we will also obtain BSE images of these grains using a Hitachi S-4800 SEM in KMICF at the University of Arizona, as well as geochemical spot analyses via EPMA. The data collected will be used to constrain the formation conditions of the Ries Crater outer suevite. Acknowledgments: We thank Ken Domanik and Jerry Chang for their support with data collection. This work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. TME thanks A. Cavosie and N. Timms for assistance during field sampling. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Melosh H. J. (1989) Oxf. U. Press. [2] Timms et al. (2017) EPSL 477, 52–58. [3] Cavosie et al. (2016) Geology 44:9, 703–706. [4] White et al. (2020) Nature Astr. 4, 974–978. [5] Crow C. A. et al. (2017) GCA 202, 264–284. [6] Timms et al. (2017) Earth-Sci. Rev. 165, 185–202. [7] Schmieder M. et al. (2018) GCA 220, 146–157. [8] Schwarz W. H. et al. (2020) M&PS 55:2, 312–325. [9] Stöffler et al. (2013) M&PS 43:4, 515–589.

Zircon↗
DOE OSTI2021

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗
DOE OSTI2023

A Graph Dynamical neural network approach for decoding dynamical states in ferroelectrics.

Ferroelectric materials such as BaTiO 3 show tremendous potential for emerging advances in memory devices, particular neuromorphic type devices. High density of memory can be obtained by stabilising polar domain walls at the nanoscale, regions of discontinuity between the well-defined polarization order parameter, but little is known about what controls their structure and dynamics in real nanoscale materials. Indeed, chiral polar domain walls have been observed in heterogeneous ferroelectrics, such as oxygen-deficient BaTiO 3 , but very little is known about how such polar-domains walls interact with defects. Indeed, a critical understanding of how dynamics of domain-walls depend on point-defects is crucial to create engineered ferroelectric memory devices. For this work, we perform large-scale simulations of nansocale domain-wall dynamics in pristine and defective BaTiO 3 using reactive force-field developed by us earlier (Phys. Chem. Chem. Phys., 2019, 21, 18240–18249), and capture their dynamical dependence on point defects using a graph dynamical neural-network approach, which we adapted to interrogate solids with well-defined order-parameters, and implemented using Pytorch based libraries. Our machine learning (ML) approach goes beyond the traditional post-processing methods to capture both spatial and temporal heterogeneities of large-scale molecular dynamics simulations of complex defective ferroelectric oxide materials. We crucially find that isolated oxygen vacancies introduce very localized spatial regions (~1–2 unit-cell in length) that show slow dipole relaxation due to formation of defect-dipoles, and that these defect-dipoles in turn slow the intrinsic dynamics of domain walls. Further, the roughness of domain walls, also influenced by vacancies, introduce dynamic heterogeneity along the domain-wall. As such we find a novel mechanism by which quenched disorder due to defects introduce dynamic heterogeneity thereby influencing response to external fields (particularly time varying fields) in a ferroelectric. Our study also emphasizes the need for creating digital twins of dynamical quantities to achieve autonomous in operando control of nanoscale switching.

42 ENGINEERING↗
NASA NTRS1995

Lead-Lag Control for Helicopter Vibration and Noise Reduction

As a helicopter transitions from hover to forward flight, the main rotor blades experience an asymmetry in flow field around the azimuth, with the blade section tangential velocities increasing on the advancing side and decreasing on the retreating side. To compensate for the reduced dynamic pressure on the retreating side, the blade pitch angles over this part of the rotor disk are increased. Eventually, a high enough forward speed is attained to produce compressibility effects on the advancing side of the rotor disk and stall on the retreating side. The onset of these two phenomena drastically increases the rotor vibratory loads and power requirements, thereby effectively establishing a limit on the maximum achievable forward speed. The alleviation of compressibility and stall (and the associated decrease in vibratory loads and power) would potentially result in an increased maximum forward speed. In the past, several methods have been examined and implemented to reduce the vibratory hub loads. Some of these methods are aimed specifically at alleviating vibration at very high flight speeds and increasing the maximum flight speed, while others focus on vibration reduction within the conventional flight envelope. Among the later are several types passive as well as active schemes. Passive schemes include a variety of vibration absorbers such as mechanical springs, pendulums, and bifilar absorbers. These mechanism are easy to design and maintain, but incur significant weight and drag penalties. Among the popular active control schemes in consideration are Higher Harmonic Control (HHC) and Individual Blade Control (IBC). HHC uses a conventional swash plate to generate a multi-cyclic pitch input to the blade. This requires actuators capable of sufficiently high power and bandwidth, increasing the cost and weight of the aircraft. IBC places actuators in the rotating reference frame, requiring the use of slip rings capable of transferring enough power to the actuators. Both schemes cause an increase in pitch link loads. Trailing Edge Flap (TEF) deployment can also used to generate unsteady aerodynamic forces and moments that counter the original vibratory loads, and thereby reduce rotor vibrations. While the vibrations absorbers, HHC, IBC, and TEF concepts discussed above attempt to reduce the vibratory loads, they do not specifically address the phenomena causing the vibrations at high advance ratios. One passive method that attempts to directly alleviate compressibility and stall, instead of reducing the ensuing vibrations, is the use of advanced tip designs. Taper, sweep, anhedral, and the manipulation of other geometric properties of the blade tips can reduce the severity of stall and compressibility effects , as well as reduce rotor power. A completely different approach to solve these problems is the tiltrotor configuration. As the forward velocity of the aircraft increases, the rotors, in this case, are tilted forward until they are perpendicular to the flow and act as propellers. This eliminates the edgewise flow encountered by conventional rotors and circumvents all the problems associated with flow asymmetry. However, the success involves a tremendous increase in cost and complexity of the aircraft. Another possible approach that has been proposed for the alleviation of vibratory loads at high forward flight speeds involves the use of controlled lead-lag motions to reduce the asymmetry in flow. A correctly phased 1/rev controlled lag motion could be introduced such that it produces a backward velocity on the advancing side and a forward velocity on the retreating side, to delay compressibility effects and stall to a higher advance ratio. Using a large enough lead-lag amplitude, the tip velocities could be reduced to levels encountered in hover. This concept was examined by two groups in the 1950's and early 1960's. In the United States, the Research Labs Division of United Aircraft developed a large lead-lag motion rotor, meant to achieve lag motion amplitudes up to 45 degrees. In order to reduce the required actuation force, the blade hinges were moved to 40% of the blade radius to increase the rotating lag frequency to approximately 1/rev. The blade hinges were redesigned to produce a flap-lag coupling so the large flapwise aerodynamic loads could be exploited to actuate the blades in the lag direction. A wind tunnel test of this rotor concept revealed actuation and blade motion scheduling problems. The project was eventually discontinued due to these problems and high blade stresses. Around the same time, at Boelkow in Germany, a similar lead-lag rotor program was conducted under the leadership of Hans Derschmidt. Here, too, the blade hinges were moved outboard to 34% radius to reduce the actuation loads. The main difference between this and the United Aircraft program was the use of a mechanical actuation scheme with maximum lead-lag motions of 400. This program was also discontinued for unclear reasons. The present study is directed toward conducting a comprehensive analytical examination to evaluate the effectiveness of controlled lead-lag motions in reducing vibratory hub loads and increasing maximum flight speed. Since both previous studies on this subject were purely experimental, only a limited data set and physical understanding of the problem was obtained. With the currently available analytical models and computational resources, the present effort is geared toward developing an in-depth physical understanding of the precise underlying mechanisms by which vibration reduction may be achieved. Additionally, in recognition of the fact that large amplitude lead-lag motions would - (i) be difficult to implement, and (ii) produce very large blade stresses; the present study examines the potential of only moderate-to-small lead-lag motions for reduction of vibratory hub loads. Using such an approach, the emphasis is not on eliminating the periodic variations in tangential velocity at the blade tip, but at best reducing these variations slightly so that compressibility and stall are delayed to slightly higher advance ratios. This study was conducted in two steps. In the first step, a hingeless helicopter rotor was modeled using rigid blades undergoing flap-lag-torsion rotations about spring restrained hinges and bearings. This model was then modified by separating the lead-lag degree of freedom into two components, a free and a prescribed motion. Using this model, a parametric study of the effect of phase and amplitude of a prescribed lead-lag motion on hub vibration was conducted. The data gathered was analyzed to obtain an understanding of the basic physics of the problem and show the capability of this method to reduce vibration and expand the flight envelope. In the second half of the study, the similar analysis was conducted using an elastic blade model to confirm the effects predicted by the simpler model.

Gandhi, Farhan↗