Galilean Invariance Research Articles

Galilean relativity for Brownian dynamics and energetics

We study experimentally the impact of inertial reference frame changes on overdamped Brownian motion. The reference frame changes are implemented by inducing, with a laser, laminar convection flows in a column of fluid where Brownian microspheres are dispersed. The convection flow plays the role of the relative velocity between the laboratory and the fluid comoving frames, and enables us to analyse the consequences of Galilean transformations on Brownian diffusion. We verify in particular how the Brownian dynamics remains ‘weakly’ Galilean invariant, in agreement with recent discussions Cairoli et al (2018 Proc. Natl Acad. Sci. USA 115 5714). We also carefully look at the consequences of Galilean relativity on the Brownian energetics. This leads us to derive a Galilean invariant expression of the stochastic thermodynamic first law, consistent with existing theoretical results Speck et al (2008 Phys. Rev. Lett. 100 178302). We finally discuss a potential ambiguity of the Galilean relativistic features of diffusive systems that has obvious practical implications in the context of force measurements in external flows.

The Mystery of the Lorentz Transform: A Reconstruction and Its Implications for Einstein’s Theories of Relativity and cosmology.

The Lorentz Transformation (LT) is an arbitrary and poorly conceived mathematical tool designed to make Maxwell’s electromagnetism conform to Galilean relativity, which formed the basis of classical mechanics and physics. A strange combination of this transform with an axiomatic assumption by Albert Einstein that the velocity of light c is an absolute and universal constant has led to an idealist, geometrical and phenomenological view of the universe, that is at variance with objective reality. This conundrum that has lasted for more than hundred years has led to rampant mysticism and has impaired the development of positive knowledge of the universe. The present reconstruction of LT shows that the gamma term, which fueled mysticism in physics and cosmology is, on the contrary, a natural outcome of the subjective geometrical rendition of the speed of light and the idealist unification of abstract space and time into a 4D “spacetime” manifold; by Minkowski and Einstein. Only a materialist dialectical perspective of space and time can rid physics of all mysticism arising out of LT; from the quantum to the cosmic.

S-Shaped Current-Voltage Characteristics of n+ - i - n - n+ Graphene Field-Effect Transistors due to the Coulomb Drag of Quasiequilibrium Electrons by Ballistic Electrons

We demonstrate that the injection of the ballistic electrons into the two-dimensional electron plasma in lateral n$^+$-i-n-n$^+$ graphene field-effect transistors (G-FET) might lead to a substantial Coulomb drag of the quasi-equilibrium electrons due the violation of the Galilean and Lorentz invariance in the systems with a linear electron dispersion. This effect can result in the S-shaped current-voltage characteristics (IVs). The resulting negative differential conductivity enables the hysteresis effects and current filamentation that can be used for the implementation of voltage switching devices. Due to a strong nonlinearity of the IVs, the G-FETs can be used for an effective frequency multiplication and detection of terahertz radiation.

A COMBINED NEURAL NETWORK AND MULTIPLE MODIFICATION STRATEGY FOR REYNOLDS-AVERAGED NAVIER-STOKES TURBULENCE MODELING

Solving the Reynolds-averaged Navier-Stokes (RANS) equation remains an effective and practical approach in engineering applications, but the uncertainty of Reynolds stress modeling will lead to discrepancies in the prediction accuracy of this approach. With the development of artificial intelligence, the data-driven method of turbulence model combined with machine learning algorithm is more effective than the original RANS model, however, the stability and prediction accuracy of the data-driven method could still be further improved. In the present paper, a fully connected neural network is constructed to predict the eddy viscosity, and this neural network is called as Eddy Viscosity Neural Network (EVNN). Additionally, a tensor-based neural network (TBNN) is also applied to predict the higher-order eddy viscosity relationship between the unclosed quantity and the analytical quantity, and the basis tensors are used to ensure the Galilean invariance. Finally, the closed-loop accuracy of the predicted flow field is realized through multiple modifications. For the method above, the neural network which is combined by EVNN and TBNN, is trained by using the high-fidelity data generated by the large eddy simulation (LES) and the baseline data obtained by the RANS simulation. Compared with the high-fidelity LES results, the results of the modified model exhibit significantly higher accuracy in the posterior velocity field, the mean pressure coefficient, and the mean friction coefficient than the original RANS model. It can be found that the implicit treatment of the linear part of the Reynolds stress can enhance the numerical stability, and the modification of the nonlinear part of the Reynolds stress can better predict the anisotropic characteristics of the flow field. Furthermore, the prediction accuracy is further improved through the multiple modification strategy. Therefore, the combined neural network and multiple modification strategy developed in this paper, have strong potentials in data-driven turbulence modeling and engineering applications in the future.

Enhanced Scattering between Electrons and Exciton-Polaritons in a Microcavity.

The interplay between strong light-matter interactions and charge doping represents an important frontier in the pursuit of exotic many-body physics and optoelectronics. Here, we consider a simplified model of a two-dimensional semiconductor embedded in a microcavity, where the interactions between electrons and holes are strongly screened, allowing us to develop a diagrammatic formalism for this system with an analytic expression for the exciton-polariton propagator. We apply this to the scattering of spin-polarized polaritons and electrons, and show that this is strongly enhanced compared with exciton-electron interactions. As we argue, this counterintuitive result is a consequence of the shift of the collision energy due to the strong light-matter coupling, and hence this is a generic feature that applies also for more realistic electron-hole and electron-electron interactions. We furthermore demonstrate that the lack of Galilean invariance inherent in the light-matter coupled system can lead to a narrow resonancelike feature for polariton-electron interactions close to the polariton inflection point. Our results are potentially important for realizing tunable light-mediated interactions between charged particles.

Study of the high lying electric dipole excitations in Odd-A 153–159Eu isotopes

In this paper, we have studied the electric dipole (E1) transitions in odd-mass 153–159Eu nuclei within Translational+Galilean Invariant Quasiparticle Phonon Nuclear Model (TGI-QPNM). Within the scope, the Giant and Pygmy Dipole Resonances (GDR and PDR) have been investigated in the energy range of 5–20 MeV. The numerical results show that the contribution of PDR and GDR to the total E1 strength is 2% and 98%, respectively. It has been found that the TGI-QPNM predicts the double hump structure and that the agreement between theoretical results and the existing experimental data is perfect, especially near the neutron binding energy Sn and the first hump of the GDR. Also the comparison of the theoretical results and the existing experimental data shows reasonable agreement for the integrated moments of the cross-sections in the GDR region. We can also deduce from the structure analysis of excited states the studied isotopes, that the GDR states can be considered to be more collective than the PDR states.

A simple and general sub-grid drag model for gas–solid fast fluidization

Great efforts have been paid to improve the Energy Minimization Multi-Scale (EMMS) model to better predict heterogeneous gas–solid flow using coarse-grid computational fluid dynamics (CFD) simulation. Nevertheless, the adopted drag correction factor (Hd) in the simulation was generally subject to macroscopic operating conditions (e.g., superficial gas and solid velocities, Ug and Gs), which called for rebuilding the Hd correlation if Ug or Gs varied. In this work we simplified the EMMS sub-grid model by replacing the stability condition with an equivalent cluster voidage correlation, thus making the model possess the Galilean invariance and avoiding the recalculation of Hd in case of varying operating conditions. The pre-constructed drag correction factor can facilitate the coarse-grid CFD simulation of gas–solid fast fluidization within a broad operation range.

Extended Lattice Boltzmann Model.

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.

Free energy fluxes and the Kubo–Martin–Schwinger relation

A general, multi-component Eulerian fluid theory is a set of nonlinear, hyperbolic partial differential equations. However, if the fluid is to be the large-scale description of a short-range many-body system, further constraints arise on the structure of these equations. Here we derive one such constraint, pertaining to the free energy fluxes. The free energy fluxes generate expectation values of currents, akin to the specific free energy generating conserved densities. They fix the equations of state and the Euler-scale hydrodynamics, and are simply related to the entropy currents. Using the Kubo–Martin–Schwinger relations associated to many conserved quantities, in quantum and classical systems, we show that the associated free energy fluxes are perpendicular to the vector of inverse temperatures characterising the state. This implies that all entropy currents can be expressed as averages of local observables. In few-component fluids, it implies that the averages of currents follow from the specific free energy alone, without the use of Galilean or relativistic invariance. In integrable models, in implies that the thermodynamic Bethe ansatz must satisfy a unitarity condition. The relation also guarantees physical consistency of the Euler hydrodynamics in spatially-inhomogeneous, macroscopic external fields, as it implies conservation of entropy, and the local-density approximated Gibbs form of stationarity states. The main result on free energy fluxes is based on general properties such as clustering, and we show that it is mathematically rigorous in quantum spin chains.

Modified Path-conservative HLLEM Scheme for Magnetohydrodynamic Solar Wind Simulations

Abstract The goal of the present work is to solve the magnetohydrodynamic (MHD) system of extended generalized Lagrange multiplier (EGLM) formulation with Galilean invariance (G-EGLM MHD equations) through a modified path-conservative HLLEM finite-volume method. A second-order least-squares reconstruction with Venkatakrishnan limiter is employed for state variables, and a solenoidality-preserving condition is considered for the magnetic field with the purpose of magnetic divergence cleaning. The two-stage Runge–Kutta time-integration method is utilized to advance the MHD governing equations. Compared with the original path-conservative HLLEM method, the modified method in this paper is shock stable and is able to adjust the diffusion according to the smoothness of the physical flow so as to automatically apply more diffusion near strong shocks and less in smooth regions near rarefaction waves and at contact discontinuities. Meanwhile, it can be robustly defined in the low plasma-β region. After several tests of smooth Alfvén wave, strong Lax, odd–even perturbation, and blast-wave problems, the large-scale structures of the solar corona for Carrington Rotation 2185 are numerically modeled in a six-component grid system of spherical coordinates with input from a Carrington rotation synoptic map provided by the Helioseismic and Magnetic Imager. Numerical results show the model’s capability of producing a structured solar wind in agreement with the observations.

A consistent and conservative model and its scheme for N-phase-M-component incompressible flows

In the present work, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

Modeling ternary fluids in contact with elastic membranes.

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modeling approach for the fluid phases, we derive the governing equations for the dynamics of the ternary fluid flow and membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method for the ternary fluid (Eulerian description) and a finite difference representation of the membrane (Lagrangian description). The ternary fluid and membrane solvers are coupled through the immersed boundary method. For validation purposes, we consider the relaxation dynamics of a two-dimensional elastic capsule placed at a fluid-fluid interface. The capsule shapes, resulting from the balance of surface tension and elastic forces, are compared with equilibrium numerical solutions obtained by surface evolver. Furthermore, the Galilean invariance of the proposed model is proven. The proposed approach is versatile, allowing for the simulation of a wide range of geometries. To demonstrate this, we address the problem of a capillary bridge formed between two deformable capsules.

Application of discrete mechanics model to jump conditions in two-phase flows

Discrete mechanics is presented as an alternative to the equations of fluid mechanics, in particular to the Navier-Stokes equation. The derivation of the discrete equation of motion is built from the intuitions of Galileo, the principles of Galilean equivalence and relativity. Other more recent concepts such as the equivalence between mass and energy and the Helmholtz-Hodge decomposition complete the formal framework used to write a fundamental law of motion such as the conservation of accelerations, the intrinsic acceleration of the material medium, and the sum of the accelerations applied to it. The two scalar and vector potentials of the acceleration resulting from the decomposition into two contributions, to curl-free and to divergence-free, represent the energies per unit of mass of compression and shear.The solutions obtained by the incompressible Navier-Stokes equation and the discrete equation of motion are the same, with constant physical properties. This new formulation of the equation of motion makes it possible to significantly modify the treatment of surface discontinuities, thanks to the intrinsic properties established from the outset for a discrete geometrical description directly linked to the decomposition of acceleration. The treatment of the jump conditions of density, viscosity and capillary pressure is explained in order to understand the two-phase flows. The choice of the examples retained, mainly of the exact solutions of the continuous equations, serves to show that the treatment of the conditions of jumps does not affect the precision of the method of resolution.

Метод решёточных уравнений Больцмана: моделирование изотермических низкоскоростных течений

In this work, lattice Boltzmann method on standard lattices was descript as one of the modern method of computation fluid dynamics. The article has main theorems, which prove computational algorithm, different type’s boundary conditions and defect in Galilean invariance. Moreover, the paper has some theoretical background about physical kinetic theory, Hermite polynomials and numeric integration. Here has not any new scientist discoveries, but has explanation of basic lattice Boltzmann theory.

Dependence of Electrodynamics on the Models of the Electromagnetic Medium

This article traces the influence of physical models of the electromagnetic medium on electrodynamics. The connection between the ability of the electrodynamic equations to describe the processes under study and the adequacy of the model of the electromagnetic medium used in their derivation to the known data is shown. It is revealed that the assumption of the absence of the existence of the electromagnetic medium leads to certain methodological difficulties in deriving the electrodynamic equations. The paper highlights the main milestones in the evolution of the physical model of electromagnetic medium. As a rule, the existence of electromagnetic medium in physics has been denied up to the present time due to the lack of its consistent model. A new physical model of the electromagnetic medium is developed and based on it a system of electrodynamic equations is derived that agree with the known data of physical experiments and astronomical observations. On the basis of the proposed physical model, a mathematical description of the method of detecting the electromagnetic medium is obtained. The developed physical model made it possible not only to explain the data of known electromagnetic phenomena and experiments, but also to propose the Galilean relativity principle within the framework of classical mechanics. Possible consequences of detection of the electromagnetic medium are discussed.

Nonlinear energy band structure of spin-orbit coupled Bose-Einstein condensates in optical lattice

<sec>In a recent experiment [Hamner C, et al. <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://doi.org/10.1103/PhysRevLett.114.070401"> 2015 <i>Phys. Rev. Lett.</i> <b>114</b> 070401</ext-link>], spin-orbit coupled Bose-Einstein condensates in a translating optical lattice have been successfully prepared into any Bloch band, and directly proved to be the lack of Galilean invariance in the presence of the spin-orbit coupling. The energy band structure of the system becomes complicated because of the lack of Galilean invariance. At present, the energy band structure of the spin-orbit coupled Bose-Einstein condensates in optical lattice is still an open issue, especially the theoretical evidence for the in-depth understanding of the competition mechanism among the spin-orbit coupling, the Raman coupling, the optical lattice and the atomic interactions of the nonlinear energy band structure has not been clear yet.</sec><sec>In this paper, based on the two-mode approximation and variational analysis, the nonlinear energy band structure and current density of the spin-orbit coupled Bose-Einstein condensates in the one-dimensional optical lattice are investigated. We find that when the spin-orbit coupling, the Raman coupling, the optical lattice, and the atomic interactions satisfy certain conditions, a loop structure in the Brillouin zone edge will emerge. The critical condition for the loop structure emerging in the Brillouin zone edge is obtained in a parameter space. The Raman coupling and the optical lattice suppress the emergence of the loop structure, while the spin-orbit coupling and the atomic interactions promote the emerging of the loop structure and make the energy band structure more complex. Interestingly, the atomic interactions can make the loop structure occur at both the higher-lying bands and the lowest energy band. The energy band structure is closely related to the current density of the system. The spin-orbit coupling causes the current density to be strongly asymmetric and leads the current density distributions of different spin states to be separated from each other in the momentum space near the boundary of the Brillouin zone. The optical lattice strength and the Raman coupling can weaken the asymmetry. The appearance of loop structure breaks the Bloch oscillation and gives rise to the Landau-Zener tunneling. The separation of the current density distributions of different spin states in the momentum space means the emergence of the spin exchange dynamics. Our results are beneficial to the in-depth understanding of the nonlinear dynamics of the spin-orbit coupled Bose-Einstein condensates in optical lattice.</sec>

Dynamical properties of generalized traveling waves of exactly solvable forced Burgers equations with variable coefficients

The initial value problem for a generalized forced Burgers equation with variable coefficients Ut+(μ˙(t)/μ(t))U+UUx=(1/2μ(t))Uxx−a(t)Ux+b(t)(xU)x−ω2(t)x+f(t),x∈R,t>0, is solved using Cole-Hopf linearization and Wei-Norman Lie algebraic approach for finding the evolution operator of the associated linear diffusion type equation. As a result, solution of the initial value problem is obtained in terms of a corresponding linear second-order inhomogeneous ordinary differential equation and a standard Burgers model. Then, using the translation and Galilean invariance of standard Burgers equation, families of generalized nonlinear waves propagating according to a Newtonian type equation of motion are constructed. The influence of the damping, dilatation and forcing terms on the dynamics of shocks, multi-shocks, triangular and N-shaped generalized traveling waves and rational type solutions with moving singularities is investigated. Finally, exactly solvable models with concrete time-variable coefficients are introduced and dynamical properties of certain particular solutions are discussed.

Galilean invariance of shallow cumulus convection large-eddy simulations

In large-eddy simulations (LES) a computational-domain translation velocity can be used to improve performance by allowing longer time-step intervals. The continuous equations are Galilean invariant. However, standard finite-difference-based discretizations are not discretely invariant with the error being proportional to the product of the local translation velocity and the truncation error. Even though such numerical errors are expected to be small, it is shown that in LES of buoyant convection the turbulent large-scale flow organization can modulate and amplify the error. Galilean invariance of global flow statistics is observed in well-resolved direct numerical simulations (DNS). In LES of single-phase convection under an inversion, flow statistics are nearly Galilean invariant and do not depend on the order of accuracy of the finite difference approximation. In contrast, in LES of cloudy convection, flow statistics show strong dependence on the frame of reference and the order of approximation. The error with respect to the frame of reference becomes negligible as the order of accuracy is increased from second to sixth in the present LES. Schemes with low resolving power can produce large dispersion errors in the surface-fixed frame that can be amplified by large-scale flow anisotropies, such as strong updrafts rising in a non-turbulent free troposphere in cumulus-cloud layers. Interestingly, in the present large-eddy simulations, a second-order discretization in the proper Galilean frame can yield comparable accuracy as a high-order scheme in the surface-fixed frame.

Stability in turbulence: The interplay between shocks and vorticity in a superfluid with higher-order dispersion

Higher-order dispersion can lead to intriguing dynamics that are becoming a focus of modern hydrodynamics research. Such systems occur naturally, for example in shallow water waves and nonlinear optics, for which several types of novel dispersive shocks structures have been identified. Here we introduce ultracold atoms as a tunable quantum simulations platform for higher-order systems. Degenerate quantum gases are well controlled model systems for the experimental study of dispersive hydrodynamics in superfluids and have been used to investigate phenomena such as vortices, solitons, dispersive shock waves and quantum turbulence. With the advent of Raman-induced spin-orbit coupling, the dispersion of a dilute gas Bose-Einstein condensate can be modified in a flexible way, allowing for detailed investigations of higher-order dispersion dynamics. Here we present a combined experimental and theoretical study of shock structures generated in such a system. The breaking of Galilean invariance by the spin-orbit coupling allows two different types of shock structures to emerge simultaneously in a single system. Numerical simulations suggest that the behavior of these shock structures is affected by interactions with vortices in a manner reminiscent of emerging viscous hydrodynamics due to an underlying quantum turbulence in the system. This result suggests that spin-orbit coupling can be used as a powerful means to tun the effective viscosity in cold-atom experiments serving as quantum simulators of turbulent hydrodynamics, with applications from condensed matter and optics to quantum simulations of neutron stars.

A decoupled and stabilized lattice Boltzmann method for multiphase flow with large density ratio at high Reynolds and Weber numbers

A decoupled lattice Boltzmann method (LBM) with improved numerical stability is developed from the conventional multiple relaxation time (MRT) model by eliminating the second-order mutual effects of viscosity modes. The dependency between relaxation times and viscosity is decoupled to remove the O(u3) errors in the original model, and thus a stabilization strategy is proposed by adjusting relaxation times without changing the local viscosity under the mass and momentum conservations. This new model still stays in the conventional MRT framework and keeps the advantages of simplicity and efficiency of LBM implementation. The decoupled MRT is demonstrated to be valid for high-speed flow with Galilean invariance. Moreover, the pseudopotential model is employed in this paper to present the improvement of numerical stability in simulating the large-density-ratio multiphase flows at high Reynolds (Re) and Weber (We) numbers. Droplet collision is achieved at density ratio 1153, Re = 34474 and We = 28212. Droplet splashing is achieved at density ratio 6637, Re = 34474 (and 71928) and We = 30538 (and 259). Droplet impacting on dry wall is achieved at density ratio 5720, Re = 11462 (and 22982) and We = 13573 (and 166). All collision regimes are simulated by the decoupled MRT from low Re (We) to high Re (We).