Gas In Thermal Equilibrium Research Articles

Numerical investigations on flow and heat transfer characteristics of a high-enthalpy double-cone in thermal and chemical non-equilibrium for hypersonic propulsion

There are strong interests in hypersonic propulsion communities in predicting the flow and heat transfer characteristics over a hypersonic vehicle. However, this is quite challenging due to the emergence of thermal and chemical out-of-equilibrium effects as prominent in a high-temperature flow. In this work, we conduct 2D numerical investigations on the hypersonic double cone flow in a surrounding of the stagnation enthalpy of 15.23 MJ/kg by solving the laminar Navier–Stokes equation systems with the vibrational non-equilibrium processes and seven species air reactions. We then conduct the sensitivities analyses on the flow and heat transfer features, as different transport models and different thermochemical models are applied. It is shown that the aerodynamic heat is more sensitive to the transport models than the flow, and the transport coefficients computed by the Gupta-Yos model are more suitable for predicting the double cone flow at the high enthalpy. On the other hand, the flow field of the high enthalpy double cone is more sensitive to the thermochemical models than the heat transfer on the wall excluding the peak position. Behind the detached shock, chemical reactions are found to enhance thermal or vibrational out-of-equilibrium effects. However, vibrational non-equilibrium effects are revealed to weaken the degree of molecule dissociation reactions. Two flow separations occur in the flow field computed by using the thermal non-equilibrium gas models, but the flows of thermal equilibrium gas models separate only once. The maximum values of wall parameters without chemical reactions are found to deviate greatly from the experimental measurements. The thermochemical out-of-equilibrium gas model can predict the separation region size closer to the experiment than the two reported numerical results. The wall static pressure and heat flux are also shown to well match the experimental results. The present work should shed lights on better understanding thermochemical non-equilibrium effects of complex flow phenomenon like the shock wave/boundary layer interaction in hypersonic propulsion with high enthalpy surroundings.

Thermal rarefied gas flow simulations with moving boundaries based on discrete unified gas kinetic scheme and immersed boundary method

Rarefied gas flows in complex geometries with heat transfer and moving boundaries have received remarkable attention due to the prevailing of micro/nano science and engineering applications. This study presents an immersed boundary (IB)-based discrete unified gas kinetic scheme (DUGKS) for the simulation of complex moving boundary thermal rarefied flows. The fluid and temperature fields are solved according to the DUGKS with Shakhov collision model (SDUGKS), which is capable of modeling non-unit Prandtl number problems in all Knudsen regimes. The IB approach, featured by its superior performance to incorporate complicated boundary constraints on Cartesian grids, is applied here to handling the interaction between the micro gas and curved boundaries. Specifically, for the IB adopted in this study, we use the least square method for the reconstruction of the distribution functions near the boundary. The accuracy of the IB-SDUGKS is explored comprehensively in six benchmark problems, including the micro lid-driven cavity flow, thermal equilibrium gas in a circular domain, thermal transpiration in an enclosure, moving shuttle in a rarefied gas, movement of a piston with pressure differences and particle free motion in a lid-driven cavity. The present results agree well with the corresponding analytical and direct simulation Monte Carlo (DSMC) solutions, illustrating the versatility of the present IB-SDUGKS complex thermal rarefied flows.

Complex Langevin approach to interacting Bose gases

Quantitative numerical analyses of interacting dilute Bose-Einstein condensates are most often based on semiclassical approximations. Since the complex-valued field-theoretic action of the Bose gas does not offer itself to standard Monte Carlo techniques, simulations beyond their scope almost exclusively rely on quantum-mechanical techniques, which are inherently limited to small particle numbers. Here we explore an alternative approach based on a Langevin-type sampling in an extended state space, known as complex Langevin (CL) algorithm. While the use of the CL technique has a long-standing history in high-energy physics, in particular in the simulation of QCD at finite baryon density, applications to ultracold atoms are still in their infancy. Here we examine the applicability of the CL approach for a one- and two-component, three-dimensional non-relativistic Bose gas in thermal equilibrium, below and above the Bose-Einstein phase transition. By comparison with analytic descriptions at the Gaussian level, including Bogoliubov and Hartree-Fock theory, we find that the method allows computing reliably and efficiently observables in the regime of experimentally accessible parameters. Close to the transition, quantum corrections lead to a shift of the critical temperature which we reproduce within the numerical range known in the literature. With this work, we aim to provide a first test of CL as a potential out-of-the-box tool for the simulation of experimentally realistic situations, including trapping geometries and multicomponent/-species models.

Interacting Loop Ensembles and Bose Gases

We study interacting Bose gases in thermal equilibrium on a lattice. We establish convergence of the grand-canonical Gibbs states of such gases to their mean-field (classical field) and large-mass (classical particle) limits. The former is a classical field theory for a complex scalar field with quartic self-interaction. The latter is a classical theory of point particles with two-body interactions. Our analysis is based on representations in terms of ensembles of interacting random loops, the Ginibre loop ensemble for Bose gases and the Symanzik loop ensemble for classical scalar field theories. For small enough interactions, our results also hold in infinite volume. The results of this paper were previously sketched in Fröhlich et al. (J Stat Phys 180(1–6):810–831, 2020).

Constraint andli, work andqi: Deacon and Zhu Xi

AbstractIn this article, I bring together the famous American semiotician Terrence Deacon and the most famous proponent of Neo-Confucianism, Zhu Xi 朱熹 (1130–1200). I focus on two of Deacon’s central concepts, constraints and work. All systems are constrained in some way, i.e. they behave in certain ways, but not in all ways. And “work” means that in their existence they make some difference: a gas in thermal equilibrium does not produce work. I bring these notions together with Zhu Xi’s concepts ofli理 andqi氣.Lican be understood as internal articulation of a system or an event, its “veins.” Andqiis the power of existence of a system or an event, its “energy.” In this light, I discuss the topics of the priority ofli, the coagulation ofqi, the normativity implied in theli, and self-cultivation. A connection to Deacon can give us new tools to make sense of those ancient topics of Chinese philosophy, and a connection to Zhu Xi can give semiotics in general and Deacon’s theory in particular an extension to certain fields that have been underdeveloped in Western thought, for instance self-cultivation.

Zero absolute vorticity state in thermal equilibrium as a hydrodynamic analog of the quantum harmonic oscillator ground state

Here, we point out on an intriguing mapping between the quantum harmonic oscillator ground state and the zero absolute vorticity plane Couette flow of a Boltzmann-like density distributed ideal gas in thermal equilibrium. The mapping is obtained when the gas is embedded in a rotating frame whose rotation rate is equal to half of the frequency of the quantum harmonic oscillator.

A Path-Integral Analysis of Interacting Bose Gases and Loop Gases

We review some recent results on interacting Bose gases in thermal equilibrium. In particular, we study the convergence of the grand-canonical equilibrium states of such gases to their mean-field limits, which are given by the Gibbs measures of classical field theories with quartic Hartree-type self-interaction, and to the Gibbs states of classical gases of point particles. We discuss various open problems and conjectures concerning, e.g., Bose-Einstein condensation, polymers and $\vert \boldsymbol{\phi} \vert^{4}$-theory.

Tuning the distance to equipartition by controlling the collision rate in a driven granular gas experiment.

In a granular gas experiment of magnetized particles confined in a thin layer, the rate of dissipative collisions is tuned by adjusting the amplitude of an external magnetic field. The velocity statistics are analyzed using the dynamic and static structure factors of transverse velocity modes. Using the fluctuating hydrodynamics theory, we measure the deviation from kinetic energy equipartition in this out-of-equilibrium system as a function of the dissipative collision rate. When the collision rate is decreased, the distance to equipartition becomes smaller, meaning that the dynamical properties of this granular gas approach by analogy those of a molecular gas in thermal equilibrium.

Direct measurement of polariton-polariton interaction strength in the Thomas-Fermi regime of exciton-polariton condensation

Bosonic condensates of exciton polaritons (light-matter quasiparticles in a semiconductor) provide a solid-state platform for studies of non-equilibrium quantum systems with a spontaneous macroscopic coherence. These driven, dissipative condensates typically coexist and interact with an incoherent reservoir, which undermines measurements of key parameters of the condensate. Here, we overcome this limitation by creating a high-density exciton-polariton condensate in an optically-induced "box" trap. In this so-called Thomas-Fermi regime, the condensate is fully separated from the reservoir and its behaviour is dominated by interparticle interactions. We use this regime to directly measure the polariton-polariton interaction strength, and reduce the existing uncertainty in its value from four orders of magnitude to within three times the theoretical prediction. The Thomas-Fermi regime has previously been demonstrated only in ultracold atomic gases in thermal equilibrium. In a non-equilibrium exciton-polariton system, this regime offers a novel opportunity to study interaction-driven effects unmasked by an incoherent reservoir.

The Near-infrared CO Absorption Band as a Probe to the Innermost Part of an AGN-obscuring Material

Abstract We performed a systematic analysis of the 4.67 μm CO ro-vibrational absorption band toward nearby active galactic nuclei (AGNs) and analyzed the absorption profiles of 10 nearby galaxies collected from the AKARI and Spitzer spectroscopic observations that show the CO absorption feature by fitting a plane-parallel local thermal equilibrium gas model. We found that CO gas is warm (200–500 K) and has a large column density ( ). The heating of the gas is not explicable by either UV heating or shock heating because these processes cannot represent the large column densities of the warm gas. Instead, X-ray photons from the nuclei, which can produce large columns of warm gas with up to , are the most convincing power source. The hydrogen column density estimated from the CO band is smaller than that inferred from X-ray observations. These results indicate that the region probed by the near-infrared CO absorption is in the vicinity of the nuclei and is located outside the X-ray emitting region. Furthermore, the covering factors of nearly unity required by the observed deep absorption profiles suggest that the probed region is close to the continuum source, which can be designated as the inner rim of the obscuring material around the AGN.

On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes

Current Monte Carlo codes use one of three models: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S(α,β) model, depending on the neutron energy and the specific Monte Carlo code. This paper addresses the consequences of using the free gas scattering model, which assumes that the neutron interacts with atoms in thermal motion in a monatomic gas in thermal equilibrium at material temperature T. Most importantly, the free gas model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei, where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that the exact resonance scattering model is temperature dependent, and neglecting the resonances in the lower epithermal range can underpredict resonance absorption due to the upscattering phenomenon mentioned above, leading to an overprediction of keff by several hundred pcm. Existing methods to address this issue involve changing the neutron weights or implementing an extra rejection scheme in the free gas sampling scheme, and these all involve performing the collision analysis in the center-of-mass (CM) frame, followed by a conversion back to the laboratory frame to continue the random walk of the neutron.The goal of this paper was to develop a sampling methodology that (1) accounted for the energy-dependent scattering cross sections in the collision analysis and (2) was performed in the laboratory frame, avoiding the conversion to the CM frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials (second and fourth order) to approximate the scattering cross section in Blackshaw’s equations for the moments of the differential scattering probability distribution functions. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on the fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using methods developed in this paper showed very close comparison to results using the reference Doppler-broadened rejection correction scheme.

Quantum theory of the dissipative Josephson parametric amplifier

SummaryRecent research in superconducting quantum circuits operating close to the quantum limit results in the need of a quantum mechanical treatment of losses. Of special interest is the dynamic behaviour of an open quantum system. As an example, a negative‐resistance Josephson parametric amplifier is treated. The DC bias voltage is chosen such that a strong interaction between the Josephson junction and the two resonant circuits, the signal and the idler circuit, is achieved. Power exchange occurs between the two considered resonator modes and also between the resonator modes and the DC power supply. Losses in the resonators are modeled by the quantum Langevin method, which describes the losses by coupling the resonators to a heat bath representing a photon gas in thermal equilibrium. The derived dynamic behaviour does not provide signal energy saturation, like classically expected for parametric amplifiers. Introducing a phenomenological multi‐photon coupling approach, saturation of the amplified signal is ensured. The time evolution of the signal and noise energy is calculated and numerically evaluated for a specific example in cases of both, the quantum Langevin method with and without the phenomenological multi‐photon coupling approach. Copyright © 2017 John Wiley & Sons, Ltd.

Effective action and electromagnetic response of topological superconductors and Majorana-mass Weyl fermions

Motivated by an apparent paradox in [X L. Qi, E.Witten, S-C. Zhang, Phys. Rev. B 87 134519 (2013)] we use the method of gauged Wess-Zumino-Witten functionals to construct an effective action for a Weyl fermion whose Majorana mass arises from coupling to a charged condensate. We obtain expressions for the current induced by an external gauge field and observe that the topological part of the current is only one-third of that that might have been expected from the gauge anomaly. The anomaly is not changed by the induced mass gap however. The topological current is supplemented by a conventional supercurrent that supplies the remaining two-thirds of the anomaly once the equation of motion for the Goldstone mode is satisfied. We apply our formula for the current to resolve the apparent paradox, and also to the chiral magnetic effect (CME) where it predicts a reduction of the CME current to one third of its value for a free Weyl gas in thermal equilibrium. We attribute this reduction to a partial cancelation of the CME by a chiral vortical effect (CVE) current arising from the persistent rotation of the fluid induced by the external magnetic field.

Relativistic Rotating Boltzmann Gas Using the Tetrad Formalism

Abstract We consider an application of the tetrad formalism introduced by Cardall et al. [Phys. Rev. D 88 (2013) 023011] to the problem of a rigidly rotating relativistic gas in thermal equilibrium and discuss the possible applications of this formalism to rel- ativistic lattice Boltzmann simulations. We present in detail the transformation to the comoving frame, the choice of tetrad, as well as the explicit calculation and analysis of the components of the equilibrium particle ow four-vector and of the equilibrium stress-energy tensor.

Search for a Lorentz invariant velocity distribution of a relativistic gas

We examine the problem of the relativistic velocity distribution in a 1-dim relativistic gas in thermal equilibrium. We use numerical simulations of the relativistic molecular dynamics for a gas with two components, light and heavy particles. However in order to obtain the numerical data our treatment distinguishes two approaches in the construction of the histograms for the same relativistic molecular dynamic simulations. The first, largely considered in the literature, consists in constructing histograms with constant bins in the velocity variable and the second consists in constructing histograms with constant bins in the rapidity variable which yields Lorentz invariant histograms, contrary to the first approach. For histograms with constant bins in the velocity variable the numerical data are fitted accurately by the Jüttner distribution which is also not Lorentz invariant. On the other hand, the numerical data obtained from histograms constructed with constant bins in the rapidity variable, which are Lorentz invariant, are accurately fitted by a Lorentz invariant distribution whose derivation is discussed in this paper. The histograms thus constructed are not fitted by the Jütter distribution (as they should not). Our derivation is based on the special theory of relativity, the central limit theorem and the Lobachevsky structure of the velocity space of the theory, where the rapidity variable plays a crucial role. For v2/c2≪1 and 1/β≡kBT/m0c2≪1 the distribution tends to the Maxwell–Boltzmann distribution.

Assessment of Heat Flux Prediction Capabilities of Residual Distribution Method: Application to Atmospheric Entry Problems

AbstractIn the present contribution we evaluate the heat flux prediction capabilities of second-order accurate Residual Distribution (RD) methods in the context of atmospheric (re-)entry problems around blunt bodies. Our departing point is the computation of subsonic air flows (with air modeled either as an inert ideal gas or as chemically reacting and possibly out of thermal equilibrium gas mixture) around probe-like geometries, as those typically employed into high enthalpy wind tunnels. We confirm the agreement between the solutions obtained with the RD method and the solutions computed with other Finite Volume (FV) based codes.However, a straightforward application of the same numerical technique to hypersonic cases involving strong shocks exhibits severe deficiencies even on a geometry as simple as a 2D cylinder. In an attempt to mitigate this problem, we derive new variants of RD schemes. A comparison of these alternative strategies against established ones allows us to derive a diagnose for the shortcomings observed in the traditional RD schemes.

Maxwell speed distribution and analogue Hawking–Unruh temperature in an ontological model of a Harmonic oscillator ground state

Within an ontological (hidden variable) model of quantum fluctuation, one can discuss the actual properties of a system regardless (independent) of measurement. Here we apply an ontological model proposed earlier to investigate a Harmonic oscillator in the quantum mechanical ground state. We first show that the actual speed of the oscillator fluctuates randomly following the Maxwell–Boltzmann distribution. On the other hand, the actual energy obeys a broad Gamma distribution with an average 3ħω/2, where ω is the classical angular frequency, so that one may conclude that the outcome of a single energy measurement reveals the average of the actual energy. The distribution of actual speed (energy) thus formally resembles the distribution of speed (energy) of an ideal gas in thermal equilibrium of temperature Tg=ħω/2. We shall then argue that Tg can be written in a form analogous to the Hawking temperature for a Schwarzschild black hole in which the average distance of the oscillator from the origin plays the analogous role of the radius of the black hole event horizon. It can also be written in a form analogous to the Unruh temperature experienced by a body moving with a uniform acceleration. In the analogy, the oscillator suffers an effective acceleration which balances the attractive force of the trapping Harmonic potential, thus keeps its average position away from the origin.

Is the mean free path the mean of a distribution?

We bring attention to the fact that Maxwell's mean free path for a dilute hard-sphere gas in thermal equilibrium, (2σn)−1, which is ordinarily obtained by multiplying the average speed by the average time between collisions, is also the statistical mean of the distribution of free path lengths in such a gas.

Synchrotron Radiation Photoelectron Spectroscopy Study on Oxide Evolution during Oxidation of a Si(111)-7 × 7 Surface at 300 K: Comparison of Thermal Equilibrium Gas and Supersonic Molecular Beams for Oxygen Adsorption

The translational energy of gas molecules is a fundamental physical quantity which dominates adsorption processes in gas–solid reactions, and for this reason molecular beam experiments have been widely conducted. Although gas in thermal equilibrium (thermal gas) plays an important role in many beneficial chemical reactions supporting the most advanced technologies, understanding the adsorption processes of thermal gas in terms of translational energy remains an unsolved issue. Here, for thermal-O2 and supersonic O2 molecular beams, we present a comparative study of oxide evolution at early stages of oxidation of a Si(111)-7 × 7 surface at room temperature. Real-time observations of oxides at the surface were conducted by employing high energy resolution photoelectron spectroscopy using high brilliance synchrotron radiation. It was found that ins structure which has one oxygen atom at the backbond of a Si adatom is the common first product for oxidation by thermal-O2 and supersonic O2 molecular beams. Similarities between thermal-O2 and supersonic O2 molecular beams were found for oxide evolution and illustrate the notion that translational energy of thermal O2 can be ascribed to the average molecular kinetic energy defined based on the most probable speed of the Maxwell–Boltzmann distribution. This experimental result suggests that translational energy is the unifying reaction parameter that rationalizes adsorption mechanisms for both thermal gas and supersonic molecular beams.

Extension of the quantum-kinetic model to lunar and Mars return physics

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high-mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. A recently introduced molecular-level chemistry model, the quantum-kinetic, or Q-K, model that predicts reaction rates for gases in thermal equilibrium and non-equilibrium using only kinetic theory and fundamental molecular properties, is extended in the current work to include electronic energy level transitions and reactions involving charged particles. Like the Q-K procedures for neutral species chemical reactions, these new models are phenomenological procedures that aim to reproduce the reaction/transition rates but do not necessarily capture the exact physics. These engineering models are necessarily efficient due to the requirement to compute billions of simulated collisions in direct simulation Monte Carlo (DSMC) simulations. The new models are shown to generally agree within the spread of reported transition and reaction rates from the literature for near equilibrium conditions.