Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical effects”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantum information and quantum simulation of neutrino physics

In extreme astrophysical environments such as core-collapse supernovae and binary neutron star mergers, neutrinos play a major role in driving various dynamical and microphysical phenomena, such as baryonic matter outflows, the synthesis of heavy elements, and the supernova explosion mechanism itself. The interactions of neutrinos with matter in these environments are flavor-specific, which makes it of paramount importance to understand the flavor evolution of neutrinos. Flavor evolution in these environments can be a highly nontrivial problem thanks to a multitude of collective effects in flavor space, arising due to neutrino-neutrino ($v$-$v$) interactions in regions with high neutrino densities. A neutrino ensemble undergoing flavor oscillations under the influence of significant $v$-$v$ interactions is somewhat analogous to a system of coupled spins with long-range interactions among themselves and with an external field (‘long-range’ in momentum-space in the case of neutrinos). As a result, it becomes pertinent to consider whether these interactions can give rise to significant quantum correlations among the interacting neutrinos, and whether these correlations have any consequences for the flavor evolution of the ensemble. In particular, one may seek to utilize concepts and tools from quantum information science and quantum computing to deepen our understanding of these phenomena. In this article, we attempt to summarize recent work in this field. In conclusion, we also present some new results in a three-flavor setting, considering complex initial states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Electron-mediated entanglement of two distant macroscopic ferromagnets within a nonequilibrium spintronic device

Using the nascent concept of quantum spin-transfer torque, we demonstrate that a current pulse can be harnessed to entangle quantum localized spins of two spatially separated ferromagnets (FMs) which are initially unentangled. The envisaged setup is composed of a spin-polarizer (FM p ) and a spin-analyzer (FM a ) FM layers separated by a normal metal (NM) spacer. The injection of a current pulse into the device leads to a time-dependent superposition of many-body states characterized by a high degree of entanglement between the spin degrees of freedom of the two distant FM layers. The nonequilibrium dynamics are due to the transfer of spin angular momentum from itinerant electrons to the localized spins via a quantum spin-torque mechanism that remains active even for collinear but antiparallel arrangements of the FM p and FM a magnetizations (a situation in which the conventional spin torque is absent). We quantify the mixed-state entanglement generated between the FM layers by tracking the time evolution of the full density matrix and analyzing the build-up of the mutual logarithmic negativity over time. The effect of decoherence and dissipation in the FM layers due to coupling to bosonic baths at finite temperature, the use of multielectron current pulses, and the dependence on the number of spins are also considered in an effort to ascertain the robustness of our predictions under realistic conditions. Lastly, we propose a “current-pump–x-ray-probe” scheme, utilizing ultrafast x-ray spectroscopy, that can witness nonequilibrium and transient entanglement of the FM layers by extracting its time-dependent quantum Fisher information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Role of spin-orbit coupling on crystal-field splitting and phase-stability of rare-earth based layered intermetallic

Layered and incommensurate heterostructures have attracted much attention for the occurrence of superconductivity and charge density waves with the possibility of intercalating foreign atoms. However, lanthanide-based heterostructures where (Eu x Sm 1-x )S and TaS 2 are alternatively stacked have been scarcely investigated. In this work, we performed phase stability, bonding behavior and electronic-structure analysis of (Eu x Sm 1-x )TaS 3 using first-principles density-functional theory methods. Our phase stability analysis suggests 50 at.% solubility of Eu in (Eu x Sm 1-x )TaS 3 compared to Eu solubility in all proportions in cubic SmS. The instability of Eu beyond 50 at.% in (Eu x Sm 1-x )TaS 3 was attributed to higher density of Eu-4f states at the Fermi-level. Based on band position calculated from spin-orbit coupling effect, we constructed a qualitative schematic of possible crystal-field analysis. The local change in bond-length and bond-angle around Eu-S in (Eu x Sm 1-x )TaS 3 correlate well with our crystal-field analysis. We believe that quantum mechanical insights provided in this work will be useful to understand other complex heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Defect filtering for thermal expansion induced dislocations in III–V lasers on silicon

Epitaxially integrated III–V semiconductor lasers for silicon photonics have the potential to dramatically transform information networks, but currently, dislocations limit performance and reliability even in defect-tolerant InAs quantum dot (QD)-based lasers. Despite being below the critical thickness, QD layers in these devices contain previously unexplained misfit dislocations, which facilitate non-radiative recombination. We demonstrate here that these misfit dislocations form during post-growth cooldown due to the combined effects of (1) thermal-expansion mismatch between the III–V layers and silicon and (2) mechanical hardening in the active region. By incorporating an additional sub-critical thickness, indium-alloyed “misfit dislocation trapping layer,” we leverage these mechanical hardening effects to our advantage, displacing 95% of misfit dislocations from the QD layer in model structures. Unlike conventional dislocation mitigation strategies, the trapping layer reduces neither the number of threading dislocations nor the number of misfit dislocations. It simply shifts the position of misfit dislocations away from the QD layer, reducing the defects' impact on luminescence. In full lasers, adding a misfit dislocation trapping layer both above and below the QD active region displaces misfit dislocations and substantially improves performance: we measure a twofold reduction in lasing threshold currents and a greater than threefold increase in output power. Overall, our results suggest that devices employing both traditional threading dislocation reduction techniques and optimized misfit dislocation trapping layers may finally lead to fully integrated, commercially viable silicon-based photonic integrated circuits.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic-Field-Driven Electron Dynamics in Graphene

Graphene exhibits unique optoelectronic properties originating from the band structure at the Dirac points. It is an ideal model structure to study the electronic and optical properties under the influence of the applied magnetic field. In graphene, electric field, laser pulse, and voltage can create electron dynamics which is influenced by momentum dispersion. However, computational modeling of momentum-influenced electron dynamics under the applied magnetic field remains challenging. In this work, we perform computational modeling of the photoexcited electron dynamics achieved in graphene under an applied magnetic field. Our results show that magnetic field leads to local deviation from momentum conservation for charge carriers. With the increasing magnetic field, the delocalization of electron probability distribution increases and forms a cyclotron-like trajectory. Our work facilitates understanding of momentum resolved magnetic field effect on non-equilibrium properties of graphene, which is critical for optoelectronic and photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The neutrino gravitational memory from a core collapse supernova: phenomenology and physics potential

General Relativity predicts that the passage of matter or radiation from an asymmetrically-emitting source should cause a permanent change in the local space-time metric. This phenomenon, called the gravitational memory effect, has never been observed, however supernova neutrinos have long been considered a promising avenue for its detection in the future. With the advent of deci-Hertz gravitational wave interferometers, observing the supernova neutrino memory will be possible, with important implications for multimessenger astronomy and for tests of gravity. In this work, we develop a phenomenological (analytical) toy model for the supernova neutrino memory effect, which is overall consistent with the results of numerical simulations. This description is then generalized to several case studies of interest. We find that, for a galactic supernova, the dimensionless strain, h(t), is of order ~ 10 -22 - 10 -21 , and develops over a typical time scale that varies between ~ 0.1 - 10 s, depending on the time-evolution of the anisotropy of the neutrino emission. The characteristic strain, h c (f), has a maximum at a frequency f max ~ Script $\mathcal{O}$(10 -1 ) - Script $\mathcal{O}$(1) Hz. The detailed features of the time- and frequency-structure of the memory strain will inform us of the matter dynamics near the collapsed core, and allow to distinguish between different stellar collapse scenarios. Next generation gravitational wave detectors like DECIGO and BBO will be sensitive to the neutrino memory effect for supernovae at typical galactic distances and beyond; with Ultimate DECIGO exceeding a detectability distance of 10 Mpc

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quantum parton shower with kinematics

Parton showers, which can efficiently incorporate quantum interference effects, have been shown to be run efficiently on quantum computers. However, so far these quantum parton showers have not included the full kinematical information required to reconstruct an event, which in classical parton showers requires the use of a veto algorithm. Here, in this work, we show that adding one extra assumption about the discretization of the evolution variable allows one to construct a quantum veto algorithm, which reproduces the full quantum interference in the event, and allows one to include kinematical effects. We finally show that for certain initial states the quantum interference effects generated in this veto algorithm are classically tractable, such that an efficient classical algorithm can be devised.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

On the Remarkable Superconductivity of FeSe and Its Close Cousins

Emergent electronic phenomena in iron-based superconductors have been at the forefront of condensed matter physics for more than a decade. Much has been learned about the origins and intertwined roles of ordered phases, including nematicity, magnetism, and superconductivity, in this fascinating class of materials. In recent years, focus has been centered on the peculiar and highly unusual properties of FeSe and its close cousins. This family of materials has attracted considerable attention due to the discovery of unexpected superconducting gap structures, a wide range of superconducting critical temperatures, and evidence for nontrivial band topology, including associated spin-helical surface states and vortex-induced Majorana bound states. Here, we review superconductivity in iron chalcogenide superconductors, including bulk FeSe, doped bulk FeSe, FeTe1−xSex, intercalated FeSe materials, and monolayer FeSe and FeTe1−xSex on SrTiO3. We focus on the superconducting properties, including a survey of the relevant experimental studies, and a discussion of the different proposed theoretical pairing scenarios. In the last part of the paper, we review the growing recent evidence for nontrivial topological effects in FeSe-related materials, focusing again on interesting implications for superconductivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Experimental and theoretical studies of the doubly substituted methyl-ethyl Criegee intermediate: Infrared action spectroscopy and unimolecular decay to OH radical products

In this work the infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn ) decay via distinct transition state barriers (16.1 kcal mol –1 and 15.4 kcal mol –1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 –10 7 s –1 . The experimental rates are in good agreement with computed Rice–Ramsperger–Kassel–Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti - and syn -MECI of 473 s –1 and 660 s –1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

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

Extracting doubly excited state lifetimes in helium directly in the time domain with attosecond noncollinear four-wave-mixing spectroscopy

The helium atom, with one nucleus and two electrons, is a prototypical system to study quantum many-body dynamics. Doubly excited states, or quantum states in which both electrons are excited by one photon, showcase electronic-correlation mediated effects. In this paper, the natural lifetimes of the doubly excited 1 P o 2 s n p Rydberg series and the 1 S e 2 p 2 dark state in helium in the 60–65 eV region are measured directly in the time domain with extreme-ultraviolet/near-infrared noncollinear attosecond four-wave-mixing (FWM) spectroscopy. The measured lifetimes agree with lifetimes deduced from spectral linewidths and theoretical predictions, and the roles of specific decay mechanisms are considered. While complex spectral line shapes in the form of Fano resonances are common in absorption spectroscopy of autoionizing states, the background-free and thus homodyned character of noncollinear FWM results exclusively in Lorentzian spectral features in the absence of strong-field effects. The onset of strong-field effects that would affect the extraction of accurate natural lifetimes in helium by FWM is determined to be approximately 0.3 Rabi cycles. This study provides a systematic understanding of the FWM parameters necessary to enable accurate lifetime extractions, which can be utilized in more complex quantum systems in the future. Published by the American Physical Society 2024

Rupprecht, Patrick (ORCID:000000026491793X)↗
DOE OSTI2020

Nuclear-Electronic Orbital Ehrenfest Dynamics

The recently developed real-time nuclear-electronic orbital (RT-NEO) approach provides an elegant framework for treating electrons and selected nuclei, typically pro- tons, quantum mechanically in nonequilibrium dynamical processes. However, the RT-NEO approach neglects the motion of the other nuclei, preventing a complete de- scription of the coupled nuclear-electronic dynamics and spectroscopy. In this work, we describe the dynamical interaction between the other nuclei and the electron-proton subsystem with the mixed quantum-classical Ehrenfest dynamics method. The NEO- Ehrenfest approach propagates the electrons and quantum protons in a time-dependent variational framework, while the remaining nuclei move classically on the correspond- ing average electron-proton vibronic surface. This approach includes the non-Born- Oppenheimer eects between the electrons and quantum protons with RT-NEO and between the classical nuclei and the electron-proton subsystem with Ehrenfest dynam- ics. Spectral features for vibrational modes involving both quantum and classical nuclei are resolved from the time-dependent dipole moments. Moreover, nuclear quantum ef- fects are included directly in nonadiabatic dynamics simulations. The excited state intramolecular proton transfer in the o-hydroxybenzaldehyde molecule is shown to be faster and to exhibit a larger kinetic isotope eect with NEO-Ehrenfest compared to classical Ehrenfest dynamics. This work shows that the NEO-Ehrenfest method is a powerful tool to study dynamical processes with coupled electronic and nuclear degrees of freedom.

Zhao, Luning↗
DOE OSTI2025

Metasurfaces with Multipolar Resonances and Enhanced Light–Matter Interaction

Metasurfaces, composed of engineered nanoantennas, enable unprecedented control over electromagnetic waves by leveraging multipolar resonances to tailor light–matter interactions. This review explores key physical mechanisms that govern their optical properties, including the role of multipolar resonances in shaping metasurface responses, the emergence of bound states in the continuum (BICs) that support high-quality factor modes, and the Purcell effect, which enhances spontaneous emission rates at the nanoscale. These effects collectively underpin the design of advanced photonic devices with tailored spectral, angular, and polarization-dependent properties. This review discusses recent advances in metasurfaces and applications based on them, highlighting research that employs full-wave numerical simulations, analytical and semi-analytic techniques, multipolar decomposition, nanofabrication, and experimental characterization to explore the interplay of multipolar resonances, bound and quasi-bound states, and enhanced light–matter interactions. A particular focus is given to metasurface-enhanced photodetectors, where structured nanoantennas improve light absorption, spectral selectivity, and quantum efficiency. By integrating metasurfaces with conventional photodetector architectures, it is possible to enhance responsivity, engineer photocarrier generation rates, and even enable functionalities such as polarization-sensitive detection. The interplay between multipolar resonances, BICs, and emission control mechanisms provides a unified framework for designing next-generation optoelectronic devices. This review consolidates recent progress in these areas, emphasizing the potential of metasurface-based approaches for high-performance sensing, imaging, and energy-harvesting applications.

Kerker effect↗
DOE OSTI2022

Plasmon-induced excitation energy transfer in silver nanoparticle dimers: A real-time TDDFTB investigation

Here, using real-time quantum dynamics calculations, we perform theoretical investigations of light-induced interactions and electronic excitation transfer in a silver nanoparticle dimer. Real-time time-dependent density functional tight-binding (RT-TDDFTB) calculations provide details of the quantum dynamical processes at an electronic/atomistic level with attosecond resolution. The computational efficiency of RT-TDDFTB allows us to examine electronic dynamics up to picosecond time scales. With time scales varying over six orders of magnitude, we provide insight into interactions between the nanoparticle and laser and between nanoparticles. Our results show that the coupling between nanoparticle monomers is dependent on the separation distance between the nanoparticles in the dimer. As the interparticle distance is varied, the dipole–dipole interactions and electronic excitation transfer mechanisms are markedly different. At large distances (from 50 to 20 Å), the energy transfer from NP1 to NP2 becomes more efficient as the interparticle distance decreases. The total dipole moment of the Ag 14 nanoparticle dimer increases linearly at an interparticle distance of 20 Å and reaches its maximum after 1.2 ps. The electronic excitation transfer is also the most efficient at 20 Å. At short distances, back-transfer effects reduce the ability of the dimer and NP1 to accept energy from the incident electric field. We attribute the distance-dependent features of the nanoparticle dimer to the beating between the laser acting on NP1 and the back transfer from NP2 to NP1.

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

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Nonequilibrium dynamics of two-level systems directly after cryogenic alternating bias

Two-level systems (TLSs) are a dominant decoherence mechanism in superconducting qubits, microwave kinetic inductance detectors, and quantum dots. TLSs are tunneling states commonly found in amorphous materials and interfaces, which interact electromechanically, leading to energy loss. Fluctuations in the frequency of the TLSs are a significant problem for qubit operation and force recalibration every few hours. In this study, we probe the effect of alternating bias at cryogenic temperatures on TLSs in a purpose-built LC oscillator. When an in situ alternating bias is applied, the steady-state spectra of the TLSs disappear. Spectroscopy reveals transient behavior, in which the TLS frequency fluctuates on the order of minutes. Thermal cycling above 10 K reverses these effects, restoring the TLS spectrum to its original state. Importantly, the intrinsic loss tangent of the LC oscillator remains unchanged before and after the application of the alternating bias. We propose that the disappearance of the steady-state spectrum is caused by nonequilibrium energy buildup from strain in the oxide film introduced by the pulsed voltage-bias sequence. Understanding this nonequilibrium energy could shed light on TLS fluctuations and reduce the calibration requirements for qubits.

general physics↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗