Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum fluids”

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

Primordial black hole formation by vacuum bubbles. Part II

The discoveries of LIGO/Virgo black holes in recent years have revitalized the study of primordial black holes. In this work we investigate a mechanism where primordial black holes are formed by vacuum bubbles that randomly nucleate during inflation through quantum tunneling. After inflation, these bubbles typically run into the ambient radiation fluid with a large Lorentz factor. In our previous work, we assumed the bubble fields are strongly coupled to the standard model particles so that the bubble wall is impermeable. Here we complete this picture by considering bubbles interacting with the fluid only through gravity. By studying the scenario in several limits, we found that black holes could form in the either subcritical or supercritical regime. Depending on the model parameters, the resulting mass spectrum of the black holes could be wide or narrow, and may develop two peaks separated by a large mass range. Finally, with different spectra, these black holes may account for the LIGO/Virgo black holes, supermassive black holes, and may play an important role in dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2022

Differential geometric approaches to momentum-based formulations for fluids [Slides]

This SAND report documents CIS Late Start LDRD Project 22-0311, "Differential geometric approaches to momentum-based formulations for fluids". The project primarily developed geometric mechanics formulations for momentum-based descriptions of nonrelativistic fluids, utilizing a differential geometry/exterior calculus treatment of momentum and a space+time splitting. Specifically, the full suite of geometric mechanics formulations (variational/Lagrangian, Lie-Poisson Hamiltonian and Curl-Form Hamiltonian) were developed in terms of exterior calculus using vector-bundle valued differential forms. This was done for a fairly general version of semi-direct product theory sufficient to cover a wide range of both neutral and charged fluid models, including compressible Euler, magnetohydrodynamics and Euler-Maxwell. As a secondary goal, this project also explored the connection between geometric mechanics formulations and the more traditional Godunov form (a hyperbolic system of conservation laws). Unfortunately, this stage did not produce anything particularly interesting, due to unforeseen technical difficulties. There are two publications related to this work currently in preparation, and this work will be presented at SIAM CSE 23, at which the PI is organizing a mini-symposium on geometric mechanics formulations and structure-preserving discretizations for fluids. The logical next step is to utilize the exterior calculus based understanding of momentum coupled with geometric mechanics formulations to develop (novel) structure-preserving discretizations of momentum. This is the main subject of a successful FY23 CIS LDRD "Structure-preserving discretizations for momentum-based formulations of fluids".

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗
DOE OSTI2020

Confinement Controls the Bend Instability of Three-Dimensional Active Liquid Crystals

Spontaneous growth of long-wavelength deformations is a defining feature of active liquid crystals. We investigate the effect of confinement on the instability of 3D active liquid crystals in the isotropic phase composed of extensile microtubule bundles and kinesin molecular motors. When shear aligned, such fluids exhibit finite-wavelength self-amplifying bend deformations. By systematically changing the channel size we elucidate how the instability wavelength and its growth rate depend on the channel dimensions. Experimental findings are qualitatively consistent with a minimal hydrodynamic model, where the fastest growing deformation is set by a balance of active driving and elastic relaxation. Furthermore, our results demonstrate that confinement determines the structure and dynamics of active fluids on all experimentally accessible length scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The αs and θs in Rayleigh–Taylor and Richtmyer–Meshkov instabilities

The alphas (α b , α s ) and thetas (θ b , θ s ) are the essential parameters describing the asymptotic growth behavior of the mixing layer driven by Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities at the interface between the fluids. Extensive simulations and experiments have been conducted in the determination of these parameters. The results depend on numerical algorithms, resolutions, and the accuracy of experimental measurements. In this article, we present a unique relationship between the alphas and thetas by applying a simple buoyancy-drag mixing model. Our first main result is that, for the same fluids, these parameters are not independent because they are fundamentally connected by physics principles and they all are functions of the drag coefficients of the fluids. As a consequence, the alphas and thetas are intrinsically coupled. All are functions of any single one of these. Additionally, our second main result is that the RT bubble growth rate (α b ) is determined by the bubble competition process – a merger model – at the edge of the mixing layer. Together, these two main results determine all the alphas and thetas by strictly theoretical reasoning. Applications to both hydro and laser driven experiments at OMEGA are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
NASA NTRS1991

Extrema Principles Of Dissipation In Fluids

Report discusses application of principle of least action and other variational or extrema principles to dissipation of energy and production of entropy in fluids. Principle of least action applied successfully to dynamics of particles and to quantum mechanics, but not universally accepted that variational principles applicable to thermodynamics and hydrodynamics. Report argues for applicability of some extrema principles to some simple flows.

Horne, W. Clifton↗
DOE OSTI2024

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sneezing and asymptomatic virus transmission

The novel coronavirus disease (COVID-19) spread pattern continues to show that geographical barriers alone cannot contain a virus. Asymptomatic carriers play a critical role in the nature of this virus quickly escalating into a global pandemic. Asymptomatic carriers may transmit the virus unintentionally through sporadic sneezing. A novel Computational Fluid Dynamics (CFD) approach has been proposed with a realistic modeling of a human sneeze achieved by the combination of state-of-the-art experimental and numerical methods. This modeling approach may be suitable for future engineering analyses aimed at reshaping public spaces and common areas, with the main objective to accurately predict the spread of aerosol and droplets that may contain pathogens. This study shows that the biomechanics of a human sneeze, including complex muscle contractions and relaxations, can be accurately modeled by the angular head motion and the dynamic pressure response during sneezing. These have been considered as the human factors and were implemented in the CFD simulation by imposing a momentum source term to the coupled Eulerian–Lagrangian momentum equations. The momentum source was modeled by the measured dynamic pressure response in conjunction with the angular head motion. This approach eliminated the need to create an ad hoc set of inlet boundary conditions. With this proposed technique, it is easier to add multiple fixed and/or moving sources of sneezes in complex computational domains. Additionally, extensive sensitivity analyses based on different environmental conditions were performed, and their impact was described in terms of potential virus spread.

42 ENGINEERING↗
DOE OSTI2024

Emergent normal fluid in the superconducting ground state of overdoped cuprates

The microscopic mechanism for the disappearance of superconductivity in overdoped cuprates is still under heated debate. Here we use scanning tunneling spectroscopy to investigate the evolution of quasiparticle interference phenomenon in Bi 2 Sr 2 CuO 6+δ over a wide range of hole densities. We find that when the system enters the overdoped regime, a peculiar quasiparticle interference wavevector with arc-like pattern starts to emerge even at zero bias, and its intensity grows with increasing doping level. Its energy dispersion is incompatible with the octet model for d-wave superconductivity, but is highly consistent with the scattering interference of gapless normal carriers. The gapless quasiparticles are mainly located near the antinodes and are independent of temperature, consistent with the disorder scattering mechanism. We propose that a branch of normal fluid emerges from the pair-breaking scattering between flat antinodal bands in the quantum ground state, which is the primary cause for the reduction of superfluid density and suppression of superconductivity in overdoped cuprates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The structure of liquid carbon elucidated by in situ X-ray diffraction

Carbon has a central role in biology and organic chemistry, and its solid allotropes provide the basis of much of our modern technology. However, the liquid form of carbon remains nearly uncharted, and the structure of liquid carbon and most of its physical properties are essentially unknown. But liquid carbon is relevant for modelling planetary interiors and the atmospheres of white dwarfs, as an intermediate state for the synthesis of advanced carbon materials, inertial confinement fusion implosions, hypervelocity impact events on carbon materials and our general understanding of structured fluids at extreme conditions. Here we present a precise structure measurement of liquid carbon at pressures of around 1 million atmospheres obtained by in situ X-ray diffraction at an X-ray free-electron laser. Our results show a complex fluid with transient bonding and approximately four nearest neighbours on average, in agreement with quantum molecular dynamics simulations. The obtained data substantiate the understanding of the liquid state of one of the most abundant elements in the universe and can test models of the melting line. The demonstrated experimental abilities open the path to performing similar studies of the structure of liquids composed of light elements at extreme conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Scaling patch analysis of planar turbulent mixing layers

Proper scales for the mean flow and Reynolds shear stress in planar turbulent mixing layers are determined from a scaling patch analysis of the mean continuity and momentum equations. By seeking an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is determined as V ref =(dδ/dx)U ref , where dδ/dx is the growth rate of the mixing layer width and U ref =U h –U l is the difference between the velocity of the high speed stream U h and the velocity of the low speed stream U l . By seeking an admissible scaling for the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is determined as R uv,ref =U avg V ref =[$\frac{1}{2A_u}$$\frac{dδ}{dx}$]$U^{2}_{ref}$ where A u $_{=}^{def}$(U h –U l )/(U h +U l ) is the normalized velocity difference that emerges naturally in the admissible scaling of the mean momentum equation. Self-similar equations for the scaled mean transverse flow V* and Reynolds shear stress $R^{*}_{uv}$=R uv /R uv,ref are derived from the mean continuity and mean momentum equations. Finally, approximate equations for V* and $R^{*}_{uv}$ are developed and found to agree well with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Jet quenching: From theory to simulation

With the explosion of data on jet-based observables in relativistic heavy-ion collisions at the Large Hadron Collider and the Relativistic Heavy-Ion Collider, perturbative Quantum Chromodynamics (pQCD)-based simulations of these processes, often interacting with an expanding viscous fluid dynamical background, have taken center stage. This review is meant to bridge the gap between theory, simulation and phenomenology of jet modification in a dense medium. We will demonstrate how the existence of such end-to-end event generators with semi-realistic or even fully realistic final states allows for the most rigorous comparisons between pQCD-based jet modification theory and experiment. State-of-the-art calculations of several jet-based observables are presented. Extensions of this theory to jets in the small systems of p–A and e–A collisions are discussed.

Physics↗