Simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

Fresh cement mortar is a type of workable paste,which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering.In this paper,Papanastasiou’s model for Bingham fluids is solved by using the multiplerelaxation-time lattice Boltzmann model(MRT-LB).Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou’s model provides a good approximation of realistic Bingham plastics for values of m108.For lower values of m,Papanastasiou’s model is valid for fluids between Bingham and Newtonian fluids.The MRT-LB model is validated by two benchmark problems:2D steady Poiseuille flows and lid-driven cavity flows.Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability.We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle.Besides the rich flow structures obtained in this work,the dynamics fluid force on the round particle is calculated.Results show that both the Reynolds number Re and the Bingham number Bn afect the drag coefcients CD,and a drag coefcient with Re and Bn being taken into account is proposed.The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed.Finally,the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields.These results help simulation of fresh concrete flowing in porous media.     

On the lattice Boltzmann method for phonon transport

The lattice Boltzmann method is a discrete representation of the Boltzmann transport equation that has been employed for modeling transport of particles of different nature. In the present work, we describe the lattice Boltzmann methodology and implementation techniques for the phonon transport modeling in crystalline materials. We show that some phonon physical properties, e.g., mean free path and group velocity, should be corrected to their effective values for one- and two-dimensional simulations, if one uses the isotropic approximation. We find that use of the D2Q9 lattice for phonon transport leads to erroneous results in transient ballistic simulations, and the D2Q7 lattice should be employed for two-dimensional simulations. Furthermore, we show that at the ballistic regime, the effect of direction discretization becomes apparent in two dimensions, regardless of the lattice used. Numerical methodology, lattice structure, and implementation of initial and different boundary conditions for the D2Q7 lattice are discussed in detail.     

Non-hydrodynamic modes and general equations of state in lattice Boltzmann equations

The lattice Boltzmann equation is commonly used to simulate fluids with isothermal equations of state in a weakly compressible limit, and intended to approximate solutions of the incompressible Navier–Stokes equations. Due to symmetry requirements there are usually more degrees of freedom in the equilibrium distributions than there are constraints imposed by the need to recover the Navier–Stokes equations in a slowly varying limit. We construct equilibria for general barotropic fluids, where pressure depends only upon density, using the two-dimensional, nine velocity (D2Q9) and one-dimensional, five velocity (D1Q5) lattices, showing that one otherwise arbitrary function in the equilibria must be chosen to suppress instability.     

Exact diagonalization of the S = 1/2 XY ferromagnet on finite bcc lattices to estimate T = 0 properties on the infinite lattice

In this paper finite bcc lattices are defined by a triple of vectors in two different ways - upper triangular lattice form and compact form. In Appendix A are lists of some 260 distinct and useful bcc lattices of 9 to 32 vertices. The energy and magnetization of the S = 1/2 XY ferromagnet have been computed on these bcc lattices in the lowest states for S _z = 0, 1/2, 1 and 3/2. These data are studied statistically to fit the first three terms of the appropriate finite lattice scaling equations. Our estimates of the T = 0 energy and magnetization agree very well with spin wave and series expansion estimates. Received 1st August 2000 and Received in final form 22 December 2000     

Investigation the effect of lattice angle on the band gap width in 3D phononic crystals with rhombohedral(I) lattice

In this paper, the propagation of acoustic waves in the phononic crystal of 3D with rhombohedral(I) lattice is studied theoretically. The crystal composite constituted of nickel spheres embedded in epoxy. The calculations of the band structure and density of states are performed with the plane wave expansion method in the irreducible part of Brillouin zone. In the present work, we have investigated the effect of lattice angle on the band structure and width of the band gap rhombohedral(I) lattice in the irreducible part of the first Brillouin zone and its planes separately. The results show that more than one complete band gape are formed in the four planes of the irreducible part. The most complete band gaps are formed in the (111) plane and the widest complete band gap in (443) with an angle greater than 80 \(^{\circ }\) . So, if the sound passes through the (111) and (443) planes for the lattice angle close to 90 \(^{\circ }\) , the crystal phononic displays the excellent insulation behavior. Moreover, in the other planes, the lattice angle does not affect on the width and the number of band gaps. Also, for the filling fraction 5 %, the widest complete band gap is formed. These results are consistent with the effect of symmetry on the band gap width, because the (111) plane has the most symmetry.     

An implicit high-order spectral difference approach for large eddy simulation

The filtered fluid dynamic equations are discretized in space by a high-order spectral difference (SD) method coupled with large eddy simulation (LES) approach. The subgrid-scale stress tensor is modelled by the wall-adapting local eddy-viscosity model (WALE). We solve the unsteady equations by advancing in time using a second-order backward difference formulae (BDF2) scheme. The nonlinear algebraic system arising from the time discretization is solved with the nonlinear lower–upper symmetric Gauss–Seidel (LU-SGS) algorithm. In order to study the sensitivity of the method, first, the implicit solver is used to compute the two-dimensional (2D) laminar flow around a NACA0012 airfoil at Re = 5 × 10⁵ with zero angle of attack. Afterwards, the accuracy and the reliability of the solver are tested by solving the 2D “turbulent” flow around a square cylinder at Re = 10⁴ and Re =  2.2 × 10⁴. The results show a good agreement with the experimental data and the reference solutions.     

TheO(3) non-linear σ-model on a 2 dimensional random lattice

TheO(3) σ-model on a 2-dimensional random lattice is studied numerically. The comparison of the continuum behaviors of the model on both random and regular lattices is carried out. It is shown that both lattices have the same continuum limit of the model, and the random lattice seems not to have the advantage of the wider scaling window compared to the same sized regular square lattice.     

Implementation of the FCHC lattice gas model on the Connection Machine

The 4-dimensional FCHC lattice gas model has been implemented on a Connection Machine CM-2 with 16K processors. Symmetries are used to reduce the collision table to a size that fits into local memory. This method avoids the degradation of the Reynolds coefficientR _*, but at the price of increased computing time. Bit shuffling between parallel lattices is introduced to reduce the discrepancy between measured viscosities and those predicted from the Boltzmann approximation. Thereby a model with a negative shear viscosity is obtained: a fluid having a uniform initial velocity is unstable and organized nonuniform motions develop. Because of the buildup of very strong correlations between the parallel lattices, the discrepancy with the Boltzmann values decreases only very slowly with the number of parallel lattices.     

Double-component convection due to different boundary conditions in an infinite slot diversely oriented to the gravity

Onset of small-amplitude oscillatory and both small- and finite-amplitude steady double-component convection arising due to component different boundary conditions in an infinite slot is studied for various slot orientations to the gravity. The main focus is on two compensating background gradients of the components. The physical mechanisms underlying steady and oscillatory convection are analyzed from the perspective of a universally consistent understanding of the effects of different boundary conditions. In a horizontal slot with inviscid fluid addressed by Welander [P. Welander, Tellus Ser. A 41 (1989) 66], oscillatory convection sets in with the most unstable wave number and oscillation frequency being zero. Exact expressions for the critical fixed-value background gradient and the respective group velocity at zero wave number are derived from the long-wavelength expansion both for the horizontal slot with independently varying background gradients and for the inclined slot with the compensating gradients. In the horizontal slot with viscous fluid, the dissipation of along-slot perturbation-cell motion reduces efficiency of the oscillatory instability feedback and thus prevents the most unstable wavelength from being infinite. Based on this interpretation, the oscillatory instability of a three-dimensional (3D) nature is predicted for an interval of long two-dimensional (2D) wavelengths in an inclined slot, and such 3D instability is indeed shown to arise. Related general conditions for three-dimensionality of most unstable disturbances are also formulated. As the slot orientation changes from the horizontal by angle θ (?π/2), the oscillatory 2D marginal-stability boundaries in inviscid and viscous fluid are expected to eventually transform into respective steady ones. Oscillatory instability in the vertical slot with viscous fluid, first reported by Tsitverblit [N. Tsitverblit, Phys. Rev. E 62 (2000) R7591], is of a quasi-steady nature. Its (new) mechanism is identified. It is underlain by differential gradient diffusion. As the horizontal slot at θ = π, addressed by Tsitverblit [N. Tsitverblit, Phys. Fluids 9 (1997) 2458], changes its orientation to vertical, the wave number interval of linear steady instability shrinks to the vicinity of the most unstable zero wave number and vanishes. Consistently with the basic nature of finite-amplitude steady convection being the same in the horizontal and vertical slots, the respective convective flows are continuously transformed into each other. The dissimilarity between the nature of finite-amplitude steady convective flows in the horizontal slot with θ = 0, revealed by Tsitverblit [N. Tsitverblit, Phys. Lett. A 329 (2004) 445], and that in the vertical slot is shown to eventually give rise to a region of hysteresis in θ ∈ (0, π/2).     

The bond bending model on triangular lattice and square lattice

The bond bending model is studied using the series expansion method on a triangular lattice and on a square lattice. The elastic splay susceptibility χ_SR and the elastic compressional susceptibility χ_el are calculated up to 11th order for the triangular lattice and up to 14th order for the square lattice. The elastic splay crossover exponent, ζ_SP, is found to be ζ_SP ≈ 1.26 ± 0.05 for the triangular lattice and ζ_SP = 1.30 ± 0.04 for the square lattice which is close to the conductivity exponent, ζ_Re, of the resistor network. From the scaling relation ? _B = dv + ζ_SP, we found that the bulk modulus exponent ? _B = 3.93 ± 0.05 for the triangular lattice and ? _B = 3.97 ± 0.04 for the square lattice which is in good agreement with the result ? _B = 3.96 ± 0.04, obtained by Zabolitzky et al. using a transfer matrix technique on a honeycomb lattice.     

The Landau gauge gluon propagator in lattice QCD

A Monte Carlo calculation of the gluon propagator in the Landau gauge in SU(3) lattice gauge theory is described. The results of calculations at β = 5.6 (200 4³ × 8 lattices), β = 5.8 (400 4³ × 10 lattices and 100 6³ × 12 lattices), and β = 6.0 (100 4³ × 8 lattices) indicate that the gluon propagator resembles a massive particle propagator in which the mass grows with separation. At the largest distances accessible with these lattices, the mass is about 600 MeV.     

The properties of two-dimensional photonic crystals bandgap structure with rhombus lattice

The relative bandgap ω_R of photonic crystals structure with rhombus lattice for both TE and TM modes is studied by the plane wave expansion method. Based on the analysis, a bending waveguide of U type can be easily obtained. The steady-state distribution of the electric field with ω = 0.4765(a/λ) of TM mode is given by FDTD method and the transmission is obtained by FFT of the time domain field. All of the sufficient analysis shows the conceptual framework of the proposed rhombus lattice PCs and its great potential application in photonic integration circuit.     

Structural characterization of thermally evaporated Bi2Te3 thin films

Thermally evaporated Bi₂Te₃ thin films were deposited on glass substrates. X-ray diffraction study confirmed that the growned films are polycrystalline in nature having hexagonal structure. The film exhibits preferential orientation along the [0 1 5] direction for the films of all thickness together with other abundant planes [0 1 1 1] and [1 1 0]. Various structural parameters such as lattice constants, crystallite size, strain, and dislocation density have been calculated and they are found to be thickness dependent. The lattice parameters are found to be a=4.38 Å and c=30.40 Å. The grain size of the films increases with thickness as the dislocation density and the microstrain decreases with thickness. The mean bond energy and the average coordination number of Bi₂Te₃ thin film are found to be 1.72 eV and 2.4, respectively.     

Entropy of a one-dimensional mixed lattice gas

This paper deals with the grand canonical entropy of a lattice gas mixture. The entropy is a function of the multisite densities corresponding to the interaction pattern of the system in question. It is first evaluated for a nearest-neighborinteraction, one-dimensional simple lattice gas to show how the structure of bulk fluid is locally maintained. Generalization requires one set of interrelations among multisite densities presented in closed form for an arbitrary lattice, and one set between Boltzmann factors and multisite densities which is written down for simply connected lattices. Application is made to two-row lattices, which turn out to have local behavior from this viewpoint, as do all single-row or Bethe lattices with complete range-p interactions. Nonlocal examples are also given, and suggestions made for approximation sequences in general lattices.     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.     

A driven three-dimensional electric lattice for polar molecules

Three-dimensional (3D) driven optical lattices have attained great attention for their wide applications in the quest to engineer new and exotic quantum phases. Here we propose a 3D driven electric lattice (3D-DEL) for cold polar molecules as a natural extension. Our 3D electric lattice is composed of a series of thin metal plates in which two-dimensional square hole arrays are distributed. When suitable modulated voltages are applied to these metal plates, a 3D potential well array for polar molecules can be generated and can move smoothly back and forth in the lattice. Thus, it can drive cold polar molecules confined in the 3D electric lattice. Theoretical analyses and trajectory calculations using two types of molecules, ND₃ and PbF, are performed to justify the possibility of our scheme. The 3D-DEL offers a platform for investigating cold molecules in periodic driven potentials, such as quantum computing science, quantum information processing, and some other possible applications amenable to the driven optical lattices.     

Excitation of the orientation transition wave and chaotic dynamics in a lattice of magnetic nanoparticles

The propagation of the front of oscillatory dynamics of magnetic moments caused by the local perturbation of the lattice along a system is studied for planar three-row lattices of magnetic nanoparticles having cubic anisotropy and coupled by the dipole interaction. The propagation of both the transition between two equilibrium planar configurations and chaotic oscillations of magnetic moments in the case of their initial orientation perpendicular to the plane of the lattice is implemented. The possibility of controlling the velocity of propagation of the orientation transition wave and its stop is revealed. The appearance of the moving front of the chaotic dynamics is found at switching on and switching off the local external field and at the action of the alternating field pulse.     

Standing wave instabilities, breather formation and thermalization in a Hamiltonian anharmonic lattice

Modulational instability of travelling plane waves is often considered as the first step in the formation of intrinsically localized modes (discrete breathers) in anharmonic lattices. Here, we consider an alternative mechanism for breather formation, originating in oscillatory instabilities of spatially periodic or quasiperiodic nonlinear standing waves (SWs). These SWs are constructed for Klein-Gordon or Discrete Nonlinear Schr?dinger lattices as exact time periodic and time reversible multibreather solutions from the limit of uncoupled oscillators, and merge into harmonic SWs in the small-amplitude limit. Approaching the linear limit, all SWs with nontrivial wave vectors (0 < Q < π) become unstable through oscillatory instabilities, persisting for arbitrarily small amplitudes in infinite lattices. The dynamics resulting from these instabilities is found to be qualitatively different for wave vectors smaller than or larger than π/2, respectively. In one regime persisting breathers are found, while in the other regime the system thermalizes. Received 6 October 2001 / Received in final form 1st March 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: mjn@ifm.liu.se     

Heat conduction in a three-dimensional momentum-conserving anharmonic lattice

Heat conduction in three-dimensional anharmonic lattices was numerically studied by the Green-Kubo theory. For a given lattice width W, a dimensional crossover is generally observed to occur at a W-dependent threshold of the lattice length. Lattices shorter than W will display a 3D behavior while lattices longer than W will display a 1D behavior. In the 3D regime, the heat current autocorrelation function was found to show a power-law decay as a function of the time lag τ as τ^{β} with β=-1.2. This indicates normal heat conduction. However, the decay exponent deviates significantly from the conventional theoretical value of β=-1.5. A flat power spectrum S(ω) of the global heat current in the low-frequency limit was also observed in the 3D regime. This provides not only an alternative verification of normal heat conduction but also a clear physical insight into its origin.