A note on the hydrodynamics of viscous fluids

The cross-product of the velocity and the vorticity in a viscous incompressible fluid is formulated and its properties investigated. When the cross-product is identically null, either the flow is vortex-free, or the velocity and the vorticity are parallel to each other. The second case yields the following important result for three-dimensional flows: if the velocity and the vorticity are related by a position-independent scalar function, that function must be time-independent as well. (English translation of the seventy-five-years old Czech text — Trkal V.: asopis pro pstování mathematiky a fysiky48 (1919) 302–311.)The paper has been translated (by I. Gregora) from Czech original which appeared in asopis pro pstování mathematiky a fysiky, Vol. 48 (1919), pp. 302–311; see a contextualizing account in the preceding paper in this issue, p. 89. The reprints are available on request from the Editorial Office.I express my warm gratitude to Prof. Dr. F. Závika for his valuable comments.     

基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真

应用基于投影算法的不可压缩光滑粒子动力学(incompressible smoothed particle hydrodynamics, ISPH)法对黏性液滴变形过程进行了数值仿真. 对于张力失稳导致的粒子非物理簇集问题, 采用粒子移位技术加以解决. 为了验证本文ISPH 算法的精度和稳定性, 分别模拟了圆形黏性液滴的拉伸变形过程以及方形液滴的旋转变形过程, 得到了不同时刻液滴内部的压力变化特征, 准确地捕捉了液滴自由面演化过程, 并将数值计算结果与文献中的解析解进行了比较.分析结果表明, 基于投影算法的不可压缩光滑粒子动力学方法结合粒子移位技术, 能够有效地模拟黏性液滴变形过程, 可以得到精确和稳定的结果. 关键词： 不可压缩光滑粒子动力学 黏性液滴 自由面流动 数值仿真     

Second law of thermodynamics in volume diffusion hydrodynamics in multicomponent gas mixtures

We presented the thermodynamic structure of a new continuum flow model for multicomponent gas mixtures. The continuum model is based on a volume diffusion concept involving specific species. It is independent of the observer?s reference frame and enables a straightforward tracking of a selected species within a mixture composed of a large number of constituents. A method to derive the second law and constitutive equations accompanying the model is presented. Using the configuration of a rotating fluid we illustrated an example of non-classical flow physics predicted by new contributions in the entropy and constitutive equations.     

Time-temperature superposition in viscous liquids

Dielectric relaxation measurements on supercooled triphenyl phosphite show that time-temperature superposition (TTS) is obeyed for the primary relaxation process at low temperatures. Measurements on other molecular liquids close to the calorimetric glass transition indicate that TTS is linked to an omega(-1/2) high-frequency decay of the loss, while the loss peak width is nonuniversal.     

A study on the inclusion of body forces in the lattice Boltzmann BGK equation to recover steady-state hydrodynamics

When the lattice Boltzmann (LB) method is used to solve hydrodynamic problems containing a body force term varying in space and/or time, its modelling at the mesoscopic scale must be verified in terms of consistency in order to avoid the appearance of non-hydrodynamic error terms at the macroscopic scale. In the present work it is shown that the modelling of spatially varying steady body force terms in the LB equation must be different from the time-dependent case, when a steady-state flow solution is sought. For that, the Chapman-Enskog analysis is used to derive the LB body force model for the LB BGK equations in a steady-state flow problem. The theoretical findings are supported by numerical tests performed on two different 2D steady-state laminar flows driven by spatially varying body forces with known analytical solutions.     

Longitudinal scaling of elliptic flow in Landau hydrodynamics

This study presents a generalization of the Landau hydrodynamic solution for multiparticle production applied to non-central relativistic heavy-ion collisions. The obtained results show the longitudinal scaling of elliptic flow, v ₂, as a function of rapidity shifted by beam rapidity (y ? y _beam ) for different energies (\(\sqrt {s_{NN} } = 62.4 GeV\) and 200 GeV) and for different systems (Au-Au and Cu-Cu). It is argued that the elliptic flow and its longitudinal scaling are due to the initial transverse energy density distribution and initial longitudinal thickness effect.     

A scaling law of high-order harmonic generation

Using a nonperturbative quantum electrodynamics theory of high-order harmonic generation (HHG), a scaling law of HHG is established. The scaling law states that when the atomic binding energy E_b, the wavelength λ and the intensity I of the laser field change simultaneously to kE_b, λ/k, and k³I, respectively. The characteristics of the HHG spectrum remain unchanged, while the harmonic yield is enhanced k³ times. That HHG obeys the same scaling law with above-threshold ionization is a solid proof of the fact that the two physical processes have similar physical mechanisms. The variation of integrated harmonic yields is also discussed.     

A numerical algorithm for viscous incompressible interfacial flows

We present a new algorithm to numerically simulate two-dimensional viscous incompressible flows with moving interfaces. The motion is updated in time by using the backward difference formula through an iterative procedure. At each iteration, the pseudo-spectral technique is applied in the horizontal direction. The resulting semi-discretized equations constitute a boundary value problem in the vertical coordinate which is solved by decoupling growing and decaying solutions. Numerical tests justify that this method achieves fully second-order accuracy in both the temporal variable and vertical coordinate. As an application of this algorithm, we study the motion of Stokes waves in the presence of viscosity. Our numerical results are consistent with the recently published asymptotic solution for Stokes waves in slightly viscous fluids.     

A semi-empirical scaling law for the β-limit in tokamaks

A limited numerical optimization of the current and pressure profiles in tokamaks of various shapes and aspect ratios has led to a simple scaling law for the maximum volume averaged β which can be reached before the onset of ideal MHD instabilities: β < (2.0?2.5)μ₀I/Ba, where I is the total currrent, B the toroidal field at the centre of the plasma and a the minor plasma radius in the equatorial plane. The limit is set by the appearance of an n = 1 free boundary instability, neglecting the possible stabilizing effect of a shell. The highest β's attained in tokamaks equipped with strong auxiliary heating are shown to be very near the limit given by this scaling law in spite of their very different conditions of operation. Possible implications are discussed.     

A universal scaling law for hadronic reactions

A master scaling law is proposed for arbitrary distributions in arbitrary hadronic processes of which all experimentally established scaling laws (and a host of others, easily deduced as occasion demands) are special cases.     

Integro-differential equation method in the hydrodynamics of an incompressible viscous fluid

An integro-differential equation method is proposed to describe the motion of an incompressible viscous fluid. The method uses an analogy between the hydrodynamic equations for an incompressible viscous fluid and the magnetostatic equations. An analysis is made of the flow of an incompressible viscous fluid round an object as a specific application. The solution automatically satisfies the boundary conditions at the surface of the object and at infinity. Zh. éksp. Teor. Fiz. 112, 1332–1339 (October 1997)     

A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows

We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334–2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case—spectral approximation in space, semi-implicit time-stepping scheme—the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow.     

黏性各向异性磁流体Kelvin-Helmholtz不稳定性:二维数值研究

利用corner transport upwind和constrained transport算法求解非理想磁流体动力学方程组,对匀强平行磁场作用下,黏性各向异性等离子体自由剪切层中的Kelvin-Helmholtz不稳定性进行了数值模拟.从流动结构、涡结构演化、磁场分布、横向磁压力、抗弯磁张力等角度对各向同性和各向异性黏性算例结果进行了讨论,分析了黏性各向异性对Kelvin-Helmholtz不稳定性的影响.结果表明,黏性各向异性比黏性各向同性更利于流动的稳定.其稳定性作用是由于磁感线方向上剪切速率降低导致界面卷起程度和圈数的降低,并使卷起结构中小涡产生增殖、合并,破坏了涡的常规增长,从而导致流动的稳定.黏性各向异性对横向磁压力的影响比对抗弯磁张力更大.     

Universal scaling relations in strongly anisotropic materials

We consider the critical temperature in strongly anisotropic antiferromagnetic materials, with weak coupling between stacked planes, in order to determine the interplane coupling constant from experimentally measured susceptibilities. We present theoretical arguments for a universal relation between interplane coupling and susceptibility shown numerically by Yasuda et al. [Phys. Rev. Lett. 94, 217201 (2005)10.1103/PhysRevLett.94.217201]. We predict a more general scaling function if the system is close to a quantum critical point, a similar relation for other susceptibilities than considered in Yasuda et al., and the validity of these relations for more general phase transitions.     

Kolmogorov scaling of turbulent flow in the vicinity of the wall

How to scale even the simplest of turbulent flows continues to be a cause for considerable controversy. In the present research, a data base compiling results from channel flow direct numerical simulations and turbulent boundary layer experiments is employed to investigate the properties of shear and normal Reynolds stresses very close to the wall. Two types of scaling based on Kolmogorov length and velocity scales are analyzed. It is shown that it is highly likely that large length scales of the order of the channel half-width or the boundary layer thickness play an important role even in the innermost regions of wall-bounded turbulent flows, which hints at the persistence of Reynolds number effects in even high Reynolds number flows.     

A scaling law for molecular orbital X Rays

Absolute molecular-orbital thick-target x-ray yields have been measured for a wide variety of targets using 1.4-MeV/amu⁵⁰Ti,⁵⁸Ni, and⁸⁶Kr and 2.4-MeV/amu⁸⁴Kr beams. A scaling law for calculating two-collision molecular-orbital x-ray yields is derived which fits the shape of the measured yields very well in all collision systems examined, and the intensity within a factor of two for near symmetric ion-atom encounters.