Vorticity Tensor Research Articles

Micro-Inertia Effects in Material Flow

Abstract The mechanics of understanding a new application of the bracket theory of Non-Equilibrium Thermodynamics that allows for the incorporation of microstructural inertia effects within conformation tensor-based constitutive models of macroscopic material behavior is presented. Introducing inertia effects generally requires the replacement of a first order in time evolution equation for the conformation tensor by a second order one. Through the analysis of a simple damped oscillator we bring forward here the close connection to the structural dissipation brackets present in the two cases, with the weights being inverted as one transitions from the inertialess to the inertial description. Moreover, one may also describe inertial effects in material flow in certain situations through a simple modification of the first order evolution equation for the conformation tensor, which consists of adding a new non-affine term that couples the conformation and the vorticity tensors, as detailed in a recent publication (P. M. Mwasame, N. J. Wagner and A. N. Beris, Phys. Fluids, 30 (2018), 030704). As shown there, when applied to the low particle Reynolds flow of dilute emulsions, this reduced inertial flow model provides predictions consistent with literature-available microscopically based asymptotic results.

A Fully-Mixed Finite Element Method for then-Dimensional Boussinesq Problem with Temperature-Dependent Parameters

AbstractIn this paper, we introduce and analyze a high-order, fully-mixed finite element method for the free convection ofn-dimensional fluids,n∈{2,3}{n\in\{2,3\}}, with temperature-dependent viscosity and thermal conductivity. The mathematical model is given by the coupling of the equations of continuity, momentum (Navier–Stokes) and energy by means of the Boussinesq approximation, as well as mixed thermal boundary conditions and a Dirichlet condition on the velocity. Because of the dependence on the temperature of the fluid properties, several additional variables are defined, thus resulting in an augmented formulation that seeks the rate of strain, pseudostress and vorticity tensors, velocity, temperature gradient and pseudoheat vectors, and temperature of the fluid. Using a fixed-point approach, smallness-of-data assumptions and a slight higher-regularity assumption for the exact solution provide the necessary well-posedness results at both continuous and discrete levels. In addition, and as a result of the augmentation, no discrete inf-sup conditions are needed for the well-posedness of the Galerkin scheme, which provides freedom of choice with respect to the finite element spaces. In particular, we suggest a combination based on Raviart–Thomas, Lagrange and discontinuous elements for which we derive optimal a priori error estimates. Finally, several numerical examples illustrating the performance of the method and confirming the theoretical rates of convergence are reported.

Objective Omega vortex identification method

A new vortex identification method (Liu et al. 2016) was proposed to represent the rotation relative strength and capture and visualize the vortices in our previous study. The basic idea of the Ω method is that a ratio of the vorticity squared over the summation of the vorticity squared and the deformation squared should be used to measure the relative rotation strength. However, the vorticity tensor norm is not objective. Thus, a moving observer will observe different vortex structures in a moving reference frame, which will make people confused with the real vortex structures. In the present study, by the definitions of the net spin tensor and net vorticity vector, an objective omega vortex identification method is presented and the examples are presented to verify the vortex structures will still retain in a moving reference frame.

Effects of rotation and acceleration in the axial current: density operator vs Wigner function

The hydrodynamic coefficients in the axial current are calculated on the basis of the equilibrium quantum statistical density operator in the third order of perturbation theory in thermal vorticity tensor both for the case of massive and massless fermions. The coefficients obtained describe third-order corrections to the Chiral Vortical Effect and include the contribution from local acceleration. We show that the methods of the Wigner function and the statistical density operator lead to the same result for an axial current in describing effects associated only with vorticity when the local acceleration is zero, but differ in describing mixed effects for which both acceleration and vorticity are significant simultaneously.

Rortex based velocity gradient tensor decomposition

Recently, a vector named Rortex was proposed to represent the local fluid rotation [C. Liu et al., “Rortex—A new vortex vector definition and vorticity tensor and vector decompositions,” Phys. Fluids 30, 035103 (2018)]. In this paper, a universal Rortex based velocity gradient tensor decomposition is proposed and the relevant local velocity increment decomposition is provided. Vortex structures in boundary layer transition on a flat plate are analyzed to quantify the local rotational, compression-stretching, and shearing effects. The results demonstrate that vorticity is shearing-dominant, while the rotational part or Rortex in general occupies a small part of vorticity in most areas of this case. In other words, vorticity is a quality representing shearing rather than rotation or vortex in most regions of this case.

Letter: Galilean invariance of Rortex

A new vector named Rortex [C. Liu et al., “Rortex—A new vortex vector definition and vorticity tensor and vector decompositions,” Phys. Fluids 30, 035103 (2018)] was proposed to represent the local fluid rotation in our previous work. However, the Galilean invariance of Rortex is yet to be elaborated. In the present study, we prove that Rortex is invariant under the Galilean transformation and several examples are provided to confirm the conclusion.

Analysis of an augmented fully-mixed formulation for the coupling of the Stokes and heat equations

We introduce and analyse an augmented mixed variational formulation for the coupling of the Stokes and heat equations. More precisely, the underlying model consists of the Stokes equation suggested by the Oldroyd model for viscoelastic flow, coupled with the heat equation through a temperature-dependent viscosity of the fluid and a convective term. The original unknowns are the polymeric part of the extra-stress tensor, the velocity, the pressure, and the temperature of the fluid. In turn, for convenience of the analysis, the strain tensor, the vorticity, and an auxiliary symmetric tensor are introduced as further unknowns. This allows to join the polymeric and solvent viscosities in an adimensional viscosity, and to eliminate the polymeric part of the extra-stress tensor and the pressure from the system, which, together with the solvent part of the extra-stress tensor, are easily recovered later on through suitable postprocessing formulae. In this way, a fully mixed approach is applied, in which the heat flux vector is incorporated as an additional unknown as well. Furthermore, since the convective term in the heat equation forces both the velocity and the temperature to live in a smaller space than usual, we augment the variational formulation by using the constitutive and equilibrium equations, the relation defining the strain and vorticity tensors, and the Dirichlet boundary condition on the temperature. The resulting augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Schauder and Banach theorems, combined with the Lax-Milgram theorem and certain regularity assumptions, are applied to prove the unique solvability of the continuous system. As for the associated Galerkin scheme, whose solvability is established similarly to the continuous case by using the Brouwer fixed-point and Lax–Milgram theorems, we employ Raviart–Thomas approximations of orderkfor the stress tensor and the heat flux vector, continuous piecewise polynomials of order ≤k+ 1 for velocity and temperature, and piecewise polynomials of order ≤ k for the strain tensor and the vorticity. Finally, we derive optimala priorierror estimates and provide several numerical results illustrating the good performance of the scheme and confirming the theoretical rates of convergence.

Vorticity and Stress Tensor

The assumptions of the Navier-Stokes equations are reconsidered. Vorticity is one of the three basic isotropic motions of a fluid element about a point. Taking into consideration the isotropic and symmetric nature of the associated tensor, it was assumed that vorticity could not contribute to the viscous stresses and hence was omitted from subsequent derivations of the Navier-Stokes equations. Deviating from this classical approach, the effect of vorticity is included in the derivation of the stress tensor.The equation hence derived contains additional terms involving vortex viscosity. Analogous to Maxwell stress tensor in electromagnetic fields, another stress tensor is defined in a vorticity field. By treating vortices as physical structures, it is possible to define pressure and the shearing stress that deform the volume element.

A mixed–primal finite element method for the Boussinesq problem with temperature-dependent viscosity

In this paper we focus on the analysis of a mixed finite element method for a class of natural convection problems in two dimensions. More precisely, we consider a system based on the coupling of the steady-state equations of momentum (Navier–Stokes) and thermal energy by means of the Boussinesq approximation (coined the Boussinesq problem), where we also take into account a temperature dependence of the viscosity of the fluid. The construction of this finite element method begins with the introduction of the pseudostress and vorticity tensors, and a mixed formulation for the momentum equations, which is augmented with Galerkin-type terms, in order to deal with the non-linearity of these equations and the convective term in the energy equation, where a primal formulation is considered. The prescribed temperature on the boundary becomes an essential condition, which is weakly imposed, leading us to the definition of the normal heat flux through the boundary as a Lagrange multiplier. We show that this highly coupled problem can be uncoupled and analysed as a fixed-point problem, where Banach and Brouwer theorems will help us to provide sufficient conditions to ensure well-posedness of the problems arising from the continuous and discrete formulations, along with several applications of continuous injections guaranteed by the Rellich–Kondrachov theorem. Finally, we show some numerical results to illustrate the performance of this finite element method, as well as to prove the associated rates of convergence.

Anomalous current from the covariant Wigner function

We consider accelerated and rotating media of weakly interacting fermions in local thermodynamic equilibrium. Kinetic properties of this media are described by covariant Wigner function calculated on the basis of relativistic distribution function of particles with spin. We obtain the formulae for axial current beyond the approximation of smallness of thermal vorticity tensor and chemical potential and calculate its divergence. In the massless limit higher order terms in vorticity and chemical potential compensate each other with only the low-order contributions surviving. It is shown, that axial current gets a topological component along the 4-acceleration vector. The similarity between different approaches to baryon polarisation is established.

Rortex—A new vortex vector definition and vorticity tensor and vector decompositions

A vortex is intuitively recognized as the rotational/swirling motion of the fluids. However, an unambiguous and universally accepted definition for vortex is yet to be achieved in the field of fluid mechanics, which is probably one of the major obstacles causing considerable confusions and misunderstandings in turbulence research. In our previous work, a new vector quantity that is called vortex vector was proposed to accurately describe the local fluid rotation and clearly display vortical structures. In this paper, the definition of the vortex vector, named Rortex here, is revisited from the mathematical perspective. The existence of the possible rotational axis is proved through real Schur decomposition. Based on real Schur decomposition, a fast algorithm for calculating Rortex is also presented. In addition, new vorticity tensor and vector decompositions are introduced: the vorticity tensor is decomposed to a rigidly rotational part and a non-rotationally anti-symmetric part, and the vorticity vector is decomposed to a rigidly rotational vector which is called the Rortex vector and a non-rotational vector which is called the shear vector. Several cases, including the 2D Couette flow, 2D rigid rotational flow, and 3D boundary layer transition on a flat plate, are studied to demonstrate the justification of the definition of Rortex. It can be observed that Rortex identifies both the precise swirling strength and the rotational axis, and thus it can reasonably represent the local fluid rotation and provide a new powerful tool for vortex dynamics and turbulence research.

On the macroscopic modeling of dilute emulsions under flow in the presence of particle inertia

Recently, Mwasame et al. [“On the macroscopic modeling of dilute emulsions under flow,” J. Fluid Mech. 831, 433 (2017)] developed a macroscopic model for the dynamics and rheology of a dilute emulsion with droplet morphology in the limit of negligible particle inertia using the bracket formulation of non-equilibrium thermodynamics of Beris and Edwards [Thermodynamics of Flowing Systems: With Internal Microstructure (Oxford University Press on Demand, 1994)]. Here, we improve upon that work to also account for particle inertia effects. This advance is facilitated by using the bracket formalism in its inertial form that allows for the natural incorporation of particle inertia effects into macroscopic level constitutive equations, while preserving consistency to the previous inertialess approximation in the limit of zero inertia. The parameters in the resultant Particle Inertia Thermodynamically Consistent Ellipsoidal Emulsion (PITCEE) model are selected by utilizing literature-available mesoscopic theory for the rheology at small capillary and particle Reynolds numbers. At steady state, the lowest level particle inertia effects can be described by including an additional non-affine inertial term into the evolution equation for the conformation tensor, thereby generalizing the Gordon-Schowalter time derivative. This additional term couples the conformation and vorticity tensors and is a function of the Ohnesorge number. The rheological and microstructural predictions arising from the PITCEE model are compared against steady-shear simulation results from the literature. They show a change in the signs of the normal stress differences that is accompanied by a change in the orientation of the major axis of the emulsion droplet toward the velocity gradient direction with increasing Reynolds number, capturing the two main signatures of particle inertia reported in simulations.

Numerical simulation for the tip leakage vortex cavitation

The tip leakage vortex (TLV) cavitation is investigated by a commercial Reynolds averaged Navier-Stokes (RANS) solver. A referenced test on a NACA0009 hydrofoil is used to validate the numerical simulation. Considering the local rotation characteristics of the vortical flow, a rotation-curvature corrected Shear-Stress-Transport model (SST-CC model) is applied to simulate the time-averaged turbulent flow. Compared to the original SST model, the SST-CC model improves the prediction of the velocity in TLV on the measured sections in downstream, and the vorticity and pressure features along the TLV trajectory are analysed numerically. In order to increase the prediction accuracy for the TLV cavitation, the empirical condensation coefficient (Fc) in Zwart's cavitation model is calibrated based on the referenced experiment. By introducing a vortex identification parameter (f∗) related to the strain rate tensor and the vorticity tensor, a relationship between the Fc and f∗ is built, and the effects of the rotational motion of the vortex on the cavity are embodied in a modified Zwart's cavitation model. Compared to the conventional Zwart's cavitation model, the modified cavitation model significantly improves the prediction of the TLV cavitation and gets a better agreement with the referenced test on different conditions with various gap widths.

Hydrodynamical model of anisotropic, polarized turbulent superfluids. I: constraints for the fluxes

This work is the first of a series of papers devoted to the study of the influence of the anisotropy and polarization of the tangle of quantized vortex lines in superfluid turbulence. A thermodynamical model of inhomogeneous superfluid turbulence previously formulated is here extended, to take into consideration also these effects. The model chooses as thermodynamic state vector the density, the velocity, the energy density, the heat flux, and a complete vorticity tensor field, including its symmetric traceless part and its antisymmetric part. The relations which constrain the constitutive quantities are deduced from the second principle of thermodynamics using the Liu procedure. The results show that the presence of anisotropy and polarization in the vortex tangle affects in a substantial way the dynamics of the heat flux, and allow us to give a physical interpretation of the vorticity tensor here introduced, and to better describe the internal structure of a turbulent superfluid.

Motion decomposition, frame-indifferent derivatives, and constitutive relations at large displacement gradients from the viewpoint of multilevel modeling

Widespread approaches to generalizing geometrically linear constitutive relations to the case of large displacement gradients have been considered. These approaches are based on the replacement of the material derivatives of stress and strain tensors by frame-indifferent corotational or convective derivatives. The correctness of choosing the indifferent derivatives is analyzed from a more general viewpoint of motion decomposition into rigid and strain-induced motion. It is shown that the use of the Zaremba-Jaumann derivative in constitutive relations corresponds to motion decomposition by the Cauchy-Helmholtz theorem according to which instantaneous rigid rotation of a material particle with small neighborhood is described by the vorticity tensor. The relations derived with the use of the so-called "logarithmic spin" are analyzed. It is noted that the spin tensors entering into these relations are not associated with the material fibers (in particular with the symmetry axes of anisotropic materials) during the entire studied process of deformation. Hence these spins do not describe the rotation of the reference frame (crystallographic one for metals) in which the material property tensor is defined. A new method of motion decomposition is proposed on the basis of a two-level (macro and meso) approach for single and polycrystalline metals. The mesoscopic spin is determined by the rotation rate of the corotational coordinate system associated with the crystallographic direction and crystallographic plane. Mesoscopic constitutive relations are formulated using the proposed spin. The spin of a representative macrovolume is determined by averaging the spins of the crystallites contained in this volume. This spin is used to formulate rate-type elastic constitutive equations. Examples are given to illustrate the stress state determination for loading along closed strain paths and two-segment paths for isotropic and anisotropic (with cubic symmetry, hcp) elastic materials, and an elastoviscoplastic fcc crystallite. The determination is carried out by using the corotational derivatives in the constitutive relations which are obtained by different motion decomposition methods.

Axionic extension of the Einstein-aether theory

We extend the Einstein-aether theory to take into account the interaction between a pseudoscalar field, which describes the axionic dark matter, and a time-like dynamic unit vector field, which characterizes the velocity of the aether motion. The Lagrangian of the Einstein-aether-axion theory includes cross-terms based on the axion field and its gradient four-vector, on the covariant derivative of the aether velocity four-vector, and on the Riemann tensor and its convolutions. We follow the principles of the Effective Field theory, and include into the Lagrangian of interactions all possible terms up to the second order in the covariant derivative. Interpretation of new couplings is given in terms of irreducible parts of the covariant derivative of the aether velocity, namely, the acceleration four-vector, the shear and vorticity tensors, and the expansion scalar. A spatially isotropic and homogeneous cosmological model with dynamic unit vector field and axionic dark matter is considered as an application of the established theory; new exact solutions are discussed, which describe models with a Big Rip, Pseudo Rip and de Sitter-type asymptotic behavior.

A fully-mixed finite element method for the Navier–Stokes/Darcy coupled problem with nonlinear viscosity

AbstractWe propose and analyze an augmented mixed finite element method for the coupling of fluid flow with porous media flow. The flows are governed by a class of nonlinear Navier–Stokes and linear Darcy equations, respectively, and the transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We apply dual-mixed formulations in both domains, and the nonlinearity involved in the Navier–Stokes region is handled by setting the strain and vorticity tensors as auxiliary unknowns. In turn, since the transmission conditions become essential, they are imposed weakly, which yields the introduction of the traces of the porous media pressure and the fluid velocity as the associated Lagrange multipliers. Furthermore, since the convective term in the fluid forces the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin type terms arising from the constitutive and equilibrium equations of the Navier–Stokes equations, and the relation defining the strain and vorticity tensors. The resulting augmented scheme is then written equivalently as a fixed point equation, so that the well-known Schauder and Banach theorems, combined with classical results on bijective monotone operators, are applied to prove the unique solvability of the continuous and discrete systems. In particular, given an integer

A Study of Generalized Quasi Einstein Spacetimes with Applications in General Relativity

The aim of this paper is to investigate some geometric and physical properties of the generalized quasi Einstein spacetime G(Q E)4 under certain conditions. Firstly, we prove the existence of G(Q E)4 by constructing a non trivial example. Then it is proved that the G(Q E)4 spacetime with the conditions \(\mathcal {B}\cdot S=L_{S}Q(g,S)\), where \(\mathcal {B}\) denotes the Ricci tensor or the concircular curvature tensor is an \(N(\frac {a-b}{3})\)-quasi Einstein spacetime and in a G(Q E)4 spacetime with C ⋅ S = 0, where C is the conformal curvature tensor, a − b is an eigenvalue of the Ricci operator. Then, we deal with the Ricci recurrent G(Q E)4 spacetime and prove that in this spacetime, the acceleration vector and the vorticity tensor vanish; but this spacetime has the non-vanishing expansion scalar and the shear tensor. Moreover, it is shown that every Ricci recurrent G(Q E)4 is Weyl compatible, purely electric spacetime and its possible Petrov types are I or D.

Laboratory experiments on counter-propagating collisions of solitary waves. Part 2. Flow field

AbstractIn the companion paper (Chen & Yeh, J. Fluid Mech., vol. 749, 2014, pp. 577–596), collisions of counter-propagating solitary waves were studied experimentally by analysing the measured water-surface variations. Here we study the flow fields associated with the collisions. With the resolved velocity data obtained in the laboratory, the flow fields are analysed in terms of acceleration, vorticity, and velocity-gradient tensors in addition to the velocity field. The data show that flow acceleration becomes maximum slightly before and after the collision peak, not in accord with the linear theory which predicts the maximum acceleration at the collision peak. Visualized velocity-gradient-tensor fields show that fluid parcels are stretched vertically prior to reaching the state of maximum wave amplitude. After the collision peak, fluid parcels are stretched in the horizontal direction. The boundary-layer evolution based on the vorticity generation and diffusion processes are discussed. It is shown that flow separation occurs at the bed during the collision. The collision creates small dispersive trailing waves. The formation of the trailing waves is captured by observing the transition behaviour of the velocity-gradient-tensor field: the direction of stretching of fluid parcels alternates during the generation of the trailing waves.

Einstein-aether theory with a Maxwell field: General formalism

We extend the Einstein-aether theory to include the Maxwell field in a nontrivial manner by taking into account its interaction with the time-like unit vector field characterizing the aether. We also include a generic matter term. We present a model with a Lagrangian that includes cross-terms linear and quadratic in the Maxwell tensor, linear and quadratic in the covariant derivative of the aether velocity four-vector, linear in its second covariant derivative and in the Riemann tensor. We decompose these terms with respect to the irreducible parts of the covariant derivative of the aether velocity, namely, the acceleration four-vector, the shear and vorticity tensors, and the expansion scalar. Furthermore, we discuss the influence of an aether non-uniform motion on the polarization and magnetization of the matter in such an aether environment, as well as on its dielectric and magnetic properties. The total self-consistent system of equations for the electromagnetic and the gravitational fields, and the dynamic equations for the unit vector aether field are obtained. Possible applications of this system are discussed. Based on the principles of effective field theories, we display in an appendix all the terms up to fourth order in derivative operators that can be considered in a Lagrangian that includes the metric, the electromagnetic and the aether fields.