Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamical Study of Adsorbate-Induced Restructuring Kinetics in Bimetallic Catalysts Using the PdAu(111) Model System

The dynamic restructuring of bimetallic catalysts plays a crucial role in their catalytic activity and selectivity. In particular, catalyst pretreatment with species such as carbon monoxide and oxygen has been shown to be an effective strategy for tuning the surface composition and morphology. Mechanistic and kinetic understanding of such restructuring are fundamental to the chemistry and engineering of surface active sites but have remained challenging due to the large structural, chemical, and temporal degrees of freedom. Here, we combine time-resolved temperature-programmed infrared reflection absorption spectroscopy, ab initio thermodynamics, and machine-learning molecular dynamics to uncover previously unidentified timescale and kinetic parameters of in situ restructuring in Pd/Au(111), a highly relevant model system for dilute Pd-in-Au nanoparticle catalysts. The key innovation lies in utilizing CO not only as a chemically sensitive probe of surface Pd, but also as an agent that induces restructuring of the surface. Upon annealing in vacuum, as-deposited Pd islands became encapsulated by Au and partially dissolved into the subsurface, leaving behind isolated Pd monomers on the surface. Subsequent exposure to 0.1 mbar CO enabled Pd monomers to repopulate the surface up to 373 K, above which complete Pd dissolution occurred by 473 K, with apparent activation energies of 0.14 and 0.48 eV, respectively. Furthermore, these restructuring processes occurred over the span of ~1000 s at a given temperature. Such a minute-timescale dynamics not only elucidates the fluxional nature of alloy catalysts but also presents an opportunity to fine-tune the surface at moderate temperature and pressure conditions.

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

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

First Stage Solid Propellant Multi Debris Thermal Analysis

The crew launch vehicle considered for the Constellation (Cx) Program utilizes a first stage solid rocket motor. If an abort is initiated in first stage flight the Crew Module (CM) will separate and be pulled away from the launch vehicle via a Launch Abort System (LAS) in order to safely and quickly carry the crew away from the malfunction launch vehicle. Having aborted the mission, the launch vehicle will likely be destroyed via a Flight Termination System (FTS) in order to prevent it from errantly traversing back over land and posing a risk to the public. The resulting launch vehicle debris field, composed primarily of first stage solid propellant, poses a threat to the CM. The harsh radiative thermal environment induced by surrounding burning propellant debris may lead to CM parachute failure. A methodology, detailed herein, has been developed to address this concern and quantify the risk of first stage propellant debris leading to radiative thermal demise of the CM parachutes. Utilizing basic thermal radiation principles, a software program was developed to calculate parachute temperature as a function of time for a given abort trajectory and debris piece trajectory set. Two test cases, considered worst-case aborts with regard to launch vehicle debris environments, were analyzed using the simulation: an abort declared at Mach 1 and an abort declared at maximum dynamic pressure (Max Q). For both cases, the resulting temperature profiles indicated that thermal limits for the parachutes were not exceeded. However, short duration close encounters by single debris pieces did have a significant effect on parachute temperature, with magnitudes on the order of 10 s of degrees Fahrenheit. Therefore while these two test cases did not indicate exceedance of thermal limits, in order to quantify the risk of parachute failure due to radiative effects from the abort environment, a more thorough probability-based analysis using the methodology demonstrated herein must be performed.

Toleman, Benjamin M.↗
NASA NTRS1989

Vortex creep and the internal temperature of neutron stars - Linear and nonlinear response to a glitch

The dynamics of pinned superfluid in neutron stars is determined by the thermal 'creep' of vortices. Vortex creep can respond to changes in the rotation rate of the neutron star crust and provide the observed types of dynamical relaxation following pulsar glitches. It also gives rise to energy dissipation, which determines the thermal evolution of pulsars once the initial heat content has been radiated away. The different possible regimes of vortex creep are explored, and it is shown that the nature of the dynamical response of the pinned superfluid evolves with a pulsar's age. Younger pulsars display a linear regime, where the response is linear in the initial perturbation and is a simple exponential relaxation as a function of time. A nonliner response, with a characteristic nonlinear dependence on the initial perturbation, is responsible for energy dissipation and becomes the predominant mode of response as the pulsar ages. The transition from the linear to the nonlinear regime depends sensitively on the temperature of the neutron star interior. A preliminary review of existing postglitch observations is given within this general evolutionary framework.

Alpar, M. A.↗
DOE OSTI2022

Comparison of animal and human blood for in vitro dynamic thrombogenicity testing of biomaterials

Abstract Background To determine suitable alternatives to human blood for in vitro dynamic thrombogenicity testing of biomaterials, four different animal blood sources (ovine, bovine, and porcine blood from live donors, and abattoir porcine blood) were compared to fresh human blood. Methods To account for blood coagulability differences between individual donors and species, each blood pool was heparinized to a donor‐specific concentration immediately before testing in a dynamic flow loop system. The target heparin level was established using a static thrombosis pre‐test. For dynamic testing, whole blood was recirculated at room temperature for 1 h at 200 ml/min through a flow loop containing a single test material. Four materials with varying thrombotic potentials were investigated: latex (positive control), polytetrafluoroethylene (PTFE) (negative control), silicone (intermediate thrombotic potential), and high‐density polyethylene (HDPE) (historically thromboresistant). Thrombus weight and surface area coverage on the test materials were quantified, along with platelet count reduction in the blood. Results While donor‐specific heparin levels varied substantially from 0.6 U/ml to 7.0 U/ml among the different blood sources, each source was able to differentiate between the thrombogenic latex and the thromboresistant PTFE and HDPE materials ( p < 0.05). However, only donor ovine and bovine blood were sensitive enough to differentiate an increased response for the intermediate thrombotic silicone material compared to PTFE and HDPE. Conclusions These results demonstrated that multiple animal blood sources (particularly donor ovine and bovine blood) may be suitable alternatives to fresh human blood for dynamic thrombogenicity testing when appropriate control materials and donor‐specific anticoagulation levels are used.

Engineering↗
DOE OSTI2023

Mechanistic origins of methyl-driven Overhauser DNP

The Overhauser effect in the dynamic nuclear polarization (DNP) of non-conducting solids has drawn much attention due to the potential for efficient high-field DNP as well as a general interest in the underlying principles that enable the Overhauser effect in small molecules. We recently reported the observation of 1 H and 2 H Overhauser effects in H 3 C- or D 3 C-functionalized Blatter radical analogs, which we presumed to be caused by methyl rotation. In this work, we look at the mechanism for methyl-driven Overhauser DNP in greater detail, considering methyl librations and tunneling in addition to classical rotation. Here, we predict the temperature dependence of these mechanisms using density functional theory and spin dynamics simulations. Comparisons with results from ultralow-temperature magic angle spinning-DNP experiments revealed that cross-relaxation at temperatures above 60 K originates from both libration and rotation, while librations dominate at lower temperatures. Due to the zero-point vibrational nature of these motions, they are not quenched by very low temperatures, and methyl-driven Overhauser DNP is expected to increase in efficiency down to 0 K, predominantly due to increases in nuclear relaxation times.

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

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an “optimal” value of the charge-transfer energy for the model where T c is largest. Finally, these results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Imaging coherent phonons and precursor dynamics in LaFeAsO with 4D ultrafast electron microscopy

In this work, we used 4D ultrafast electron microscopy (UEM) to directly image femtosecond photoinduced structural dynamics in single-crystal LaFeAsO at initial temperatures of 300 and 100 K, above and below the known structural and magnetic phase-transition temperatures, respectively. With nanometer-picosecond resolution, we resolved an initial (precursor) sigmoidlike response arising from photothermal expansion and lattice reorientation that precedes the onset of propagating coherent acoustic phonons (CAPs). In the specific regions probed, the precursor response at 100 K is shorter than at 300 K ($t$ 0.5;100 K = 11.3ps vs $t$ 0.5;300 K = 17.8ps), and the CAP oscillation frequency is lowered with cooling (ƒ CAP;100 K = 12 GHz vs ƒ CAP;300 K = 21 GHz), correlated to known lattice softening due to the structural phase change. The transient CAP behaviors at 300 K are dispersive, displaying an exponentially decaying phase velocity over the first nanosecond. Further, the CAP symmetry at 300 K matches a first-order antisymmetric shear mode (A 1 ), while at 100 K it is best matched by a mostly nondispersive zero-order symmetric mode (S 0 ). These findings illustrate the sensitivity of UEM imaging to spatially heterogeneous dynamics in the Fe-pnictide materials and more broadly in other quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Real-space observation of ergodicity transitions in artificial spin ice

Abstract Ever since its introduction by Ludwig Boltzmann, the ergodic hypothesis became a cornerstone analytical concept of equilibrium thermodynamics and complex dynamic processes. Examples of its relevance range from modeling decision-making processes in brain science to economic predictions. In condensed matter physics, ergodicity remains a concept largely investigated via theoretical and computational models. Here, we demonstrate the direct real-space observation of ergodicity transitions in a vertex-frustrated artificial spin ice. Using synchrotron-based photoemission electron microscopy we record thermally-driven moment fluctuations as a function of temperature, allowing us to directly observe transitions between ergodicity-breaking dynamics to system freezing, standing in contrast to simple trends observed for the temperature-dependent vertex populations, all while the entropy features arise as a function of temperature. These results highlight how a geometrically frustrated system, with thermodynamics strictly adhering to local ice-rule constraints, runs back-and-forth through periods of ergodicity-breaking dynamics. Ergodicity breaking and the emergence of memory is important for emergent computation, particularly in physical reservoir computing. Our work serves as further evidence of how fundamental laws of thermodynamics can be experimentally explored via real-space imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Single-atom heat engine as a sensitive thermal probe

We propose employing a quantum heat engine as a sensitive probe for thermal baths. In particular, we study a single-atom Otto engine operating in an open thermodynamic cycle. Owing to its cyclic nature, the engine is capable of translating small temperature differences between two baths into a macroscopic oscillation in a flywheel. We present analytical and numerical modeling of the quantum dynamics of the engine and estimate it to be capable of detecting temperature differences as small as 2 μK. This sensitivity can be further improved by utilizing quantum resources such as squeezing of the ion motion. The proposed scheme does not require quantum state initialization and is able to detect small temperature differences in a wide range of base temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2009

Watershed and Hydrodynamic Modeling for Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation and Seagrasses in Mobile Bay

There is a continued need to understand how human activities along the northern Gulf of Mexico coast are impacting the natural ecosystems. The gulf coast is experiencing rapid population growth and associated land cover/land use change. Mobile Bay, AL is a designated pilot region of the Gulf of Mexico Alliance (GOMA) and is the focus area of many current NASA and NOAA studies, for example. This is a critical region, both ecologically and economically to the entire United States because it has the fourth largest freshwater inflow in the continental USA, is a vital nursery habitat for commercially and recreational important fisheries, and houses a working waterfront and port that is expanding. Watershed and hydrodynamic modeling has been performed for Mobile Bay to evaluate the impact of land use change in Mobile and Baldwin counties on the aquatic ecosystem. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use Scenarios in 1948, 1992, 2001, and 2030. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on observed trends. All land use scenarios were developed to a common land classification system developed by merging the 1992 and 2001 National Land Cover Data (NLCD). The LSPC model output provides changes in flow, temperature, sediments and general water quality for 22 discharge points into the Bay. These results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment concentrations on a grid with four vertical profiles throughout the Bay s aquatic ecosystems. The models were calibrated using in-situ data collected at sampling stations in and around Mobile bay. This phase of the project has focused on sediment modeling because of its significant influence on light attenuation which is a critical factor in the health of submerged aquatic vegetation. The impact of land use change on sediment concentrations was evaluated by analyzing the LSPC and EFDC sediment simulations for the four land use scenarios. Such analysis was also performed for storm and non-storm periods. In- situ data of total suspended sediments (TSS) and light attenuation were used to develop a regression model to estimate light attenuation from TSS. This regression model was used to derive marine light attenuation estimates throughout Mobile bay using the EFDC TSS outputs. The changes in sediment concentrations and associated impact on light attenuation in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on seagreasses and Submerged Aquatic Vegetation (SAV) habitat. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of sediments due to land use driven flow changes with the restoration potential of SAVs.

Estes, Maurice G.↗
DOE OSTI2024

A global urban heat island intensity dataset: Generation, comparison, and analysis

The urban heat island (UHI) effect, a phenomenon of local warming over urban areas, is the most well-known impact of urbanization on climate. Globally consistent estimates of the UHI intensity (UHII) are crucial for examining this phenomenon across time and space. However, publicly available UHII datasets are limited and have several constraints: (1) they are for clear-sky surface UHII, not all-sky surface UHII and canopy (air temperature) UHII; (2) the estimation methods often neglect anthropogenic disturbance, introducing uncertainties in the estimated UHII. To address these issues, this study proposes a new dynamic equal-area (DEA) method that can minimize the influence of various confounding factors on UHII estimates through a dynamic cyclic process. Utilizing the DEA method and leveraging various gridded temperature data, we develop a global-scale (>10,000 cities), long-term (over 20 years by month), and multi-faceted (clear-sky surface, all-sky surface, and canopy) UHII dataset. Further, based on these estimates, we provide a comprehensive analysis of the UHII and its trends in global cities. The UHII is found to be greater than zero in >80% of cities, with global annual average magnitudes around 1.0 °C (day) and 0.8 °C (night) for surface UHII, and close to 0.5 °C for canopy UHII. Furthermore, an interannual upward trend in UHII is observed in >60% of cities, with global annual average trends exceeding 0.1 °C/decade (day) and over 0.06 °C/decade (night) for surface UHII, and slightly surpassing 0.03 °C/decade for canopy UHII. Notably, there exists a positive correlation between the magnitude and trend of UHII, suggesting that cities with stronger UHII tend to experience faster growth in UHII. Additionally, discrepancies in UHII are found between different temperature data, stemming not only from distinctions in data types (surface or air temperature) but also from differences in data acquisition times (Terra or Aqua), weather conditions (clear-sky or all-sky), and processing methodologies (with or without gap filling). Overall, our proposed method, dataset, and analysis results have the potential to provide valuable insights for future urban climate studies. The UHII dataset is publicly available at https://doi.org/10.6084/m9.figshare.24821538.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Detached Eddy Simulation of Film Cooling over a GE Flat Plate

The detached eddy simulation of film cooling has been utilized for a proprietary GE plate-pipe configuration. The blowing ratio was 2.02, the velocity ratio was 1.26, and the temperature ratio was 1.61. Results indicate that the mixing processes downstream of the hole are highly anisotropic. DES solution shows its ability to depict the dynamic nature of the flow and capture the asymmetry present in temperature and velocity distributions. Further, comparison between experimental and DES time-averaged effectiveness is satisfactory. Numerical values of span-averaged effectiveness show better prediction of the experimental values at downstream locations than a steady state Glenn HT solution. While the DES method shows obvious promise, there are several issues that need further investigation. Despite an accurate prediction in the hole vicinity, the simulation still falls short in the region x = 10d to 100d. This should be investigated. Also the model used flat plate. Actual turbine blade should be modeled in the future if additional finding is available.

Roy, Subrata↗
DOE OSTI2022

Effects of Temperature on Flow Fouling of Smooth and Nonwetting Surfaces

Mineral fouling (scaling) of heat transfer surfaces is a pervasive problem in heat exchangers, chemical reactors, and other equipment in energy, environmental, and process industrial applications. Many of the applications involve dynamic flow of fluids with the impurity minerals and, furthermore, operate at elevated temperature. Strategies to mitigate fouling under these conditions are of much value in industrial applications. This paper presents a comparative study of temperature-dependent mineral fouling deposition on smooth surfaces and nonwetting superhydrophobic and lubricant-infused surfaces under dynamic flow conditions. The surfaces are represented in a unified manner using the viscosity ratio of the infused material within the porous asperities on a surface to that of the flowing fluid such that the spectrum of surfaces from superhydrophobic to smooth is captured by the range of viscosity ratios from 0 to ∞. Using a forced convection experimental setup, deposition of calcium sulfate on the surfaces is quantified in terms of asymptotic fouling resistance over a range of temperature, Reynolds number, and mineral foulant supersaturation. Through a systematic set of accelerated fouling experiments, an empirical relationship for the asymptotic fouling resistance is developed in terms of Reynolds number, foulant concentration, temperature, and surface type. The empirical model is validated with a comprehensive set of experimental data from this study as well as from the literature. Optimum nonwetting surface designs for minimizing fouling resistance compared to conventional smooth surfaces are developed as a function of temperature. The results of the study offer insight into the temperature-dependent fouling of surfaces under flow conditions and a rational design of fouling-resistant nonwetting surfaces that can be readily translated to practice.

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

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
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↗