A simulation study of an asymmetric exclusion model with open and periodic boundaries for parallel dynamics

The effect of jumping rate probability on the phase diagram of an asymmetric exclusion model is studied by numerical simulations. Density, current and velocity of particles are calculated for parallel dynamics. In the open boundaries case for one species of particles (particles 1), a passage from first to second order transition occurs by decreasing the jumping rate. In the periodic boundaries case, by introducing another species of particle (particle 2) which plays the role of obstacle for particles 1, the average velocity of particles 1 increases with increasing the jumping rate for small density. While the average velocity of particle 2 decreases for small and intermediate densities. Received: 21 October 1997     

MHD Boundary Layer Flow of a Non-Newtonian Fluid on a Moving Surface with a Power-Law Velocity

A theoretical analysis for MHD boundary layer flow on a moving surface with the power-law velocity is presented. An accurate expression of the skin friction coefficient is derived. The analytical approximate solution is obtained by means of Adomian decomposition methods. The reliability and efficiency of the approximate solutions are verified by numerical ones in the literature.     

Acoustic radiation force on a rigid cylinder near rigid corner boundaries exerted by a Gaussian beam field

Acoustic manipulation is one of the well-known technologies of particle control and a top research in acoustic field. Calculation of acoustic radiation force on a particle nearby boundaries is one of the critical tasks, as it approximates realistic applications. Nevertheless, it is quite difficult to solve the problem by theoretical method when the boundary conditions are intricate. In this study, we present a finite element method numerical model for the acoustic radiation force exerting on a rigid cylindrical particle immersed in fluid near a rigid corner. The effects of the boundaries on acoustic radiation force of a rigid cylinder are analyzed with particular emphasis on the non-dimensional frequency and the distance from the center of cylinder to each boundary. The results reveal that these parameters play important roles in acoustic manipulation for particle-nearby complicated rigid boundaries. This study verifies the feasibility of numerical analysis on the issue of acoustic radiation force calculation close to complex boundaries, which may provide a new idea on analyzing the acoustic particle manipulation in confined space.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

Effect of Rolling Massage on Particle Moving Behaviour in Blood Vessels

The rolling massage manipulation is a classic Chinese massage, which is expected to eliminate many diseases. Here the effect of the rolling massage on the particle moving property in the blood vessels under the rolling massage manipulation is studied by the lattice Boltzmann simulation. The simulation results show that the particle moving behaviour depends on the rolling velocity, the distance between particle position and rolling position. The average values, including particle translational velocity and angular velocity, increase as the rolling velocity increases almost linearly. The result is helpful to understand the mechanism of the massage and develop the rolling techniques.     

Lattice Boltzmann Simulations of Particle-Particle Interaction in Steady Poiseuille Flow

The moving behaviour of two- and three-particles in a pressure-driven flow is studied by the lattice Boltzmann simulation in two dimensions. The time-dependent values, including particles＇ radial positions, translational velocities, angular velocities, and the x-directional distance between the particles are analysed extensively. The effect of flow Reynolds number on particle motion is also investigated numerically. The simulation results show that the leading particle equilibrium position is closer to the channel centre while the trailing particle equilibrium position is closer to the channel wall. If Reynolds number Re is less than 85.30, the larger flow Reynolds number results in the smaller x-directional equilibrium distance, otherwise the x-directional distance increases almost linearly with the increase of time and the particles separate finally. The simulation results are helpful to understand the particle-particle interaction in suspensions with swarms of particles.     

Influence of a high-speed dielectric spherical particle''s movement on its scattering characteristics

The scattered field and differential scattered section (DSS) of a moving spherical particle with a high speed are investigated numerically. The coordinate and vector transformations are used to establish a theoretical basis for studying the laser scattering of a moving particle. The DSS of a moving spherical particle is explained by the electric and magnetic field from Mie scattering theory. Assuming the laser wavelength of 1.06 μm, we compute the ratio of the laser DSS of the moving dielectric spherical particle to that of the static dielectric spherical particle, which changes with radii, speeds and scattering angles of the particle. The numerical results show that the laser DSS of the oving pherical particle is tightly connected with its speed and scattering zenith angle. If a spherical particle moves with high speed, the laser DSS due to movement of the particle could not be neglected. If the speed of the dielectric spherical particle is fluctuating, the Doppler effect and the frequency spectrum expansion play important roles.     

Conformally symmetric vacuum solutions of the gravitational field equations in the brane-world models

A class of exact solutions of the gravitational field equations in the vacuum on the brane are obtained by assuming the existence of a conformal Killing vector field, with non-static and non-central symmetry. In this case, the general solution of the field equations can be obtained in a parametric form in terms of the Bessel functions. The behavior of the basic physical parameters describing the non-local effects generated by the gravitational field of the bulk (dark radiation and dark pressure) is also considered in detail, and the equation of state satisfied at infinity by these quantities is derived. As a physical application of the obtained solutions we consider the behavior of the angular velocity of a test particle moving in a stable circular orbit. The tangential velocity of the particle is a monotonically increasing function of the radial distance and, in the limit of large values of the radial coordinate, tends to a constant value, which is independent on the parameters describing the model. Therefore, a brane geometry admitting a one-parameter group of conformal motions may provide an explanation for the dynamics of the neutral hydrogen clouds at large distances from the galactic center, which is usually explained by postulating the existence of the dark matter.     

二维MLSPH无网格方法

给出MLSPH无网格方法的控制方程、移动最小二乘近似、MLSPH守恒格式等,重点研究MLSPH守恒格式中MLS面积向量的数值求解方法,以提高求解的精度.对激波管问题和平面Noh问题进行模拟,得到较好的结果,确定用三次样条积分公式计算面积向量.     

Moving Finite Element Simulations for Reaction-Diffusion Systems

This work is concerned with the numerical simulations for two reaction-diffusion systems, i.e., the Brusselator model and the Gray-Scott model. The numerical algorithm is based upon a moving finite element method which helps to resolve large solution gradients. High quality meshes are obtained for both the spot replication and the moving wave along boundaries by using proper monitor functions. Unlike [33], this work finds out the importance of the boundary grid redistribution which is particularly important for a class of problems for the Brusselator model. Several ways for verifying the quality of the numerical solutions are also proposed, which may be of important use for comparisons.     

On the vortex formation at the moving front of lightweight granular particles

The formation of vortices at a moving front of lightweight granular particles is investigated experimentally. The particles used in this study are made of polystyrene foam with three different diameters of nearly uniform size. Pairs of vortices are found to emerge at the moving front at regular intervals, thereby forming a wavy pattern. Once the vortices are produced, the flow velocity tends to increase. A simple analysis suggests the existence of a velocity boundary layer at the moving front, whose thickness increases with increasing particle diameter. The frontal radius of each vortex pair is about the size of this boundary layer; when the radius exceeds this size, the front tends to bifurcate into a train of vortices with the size of the boundary layer. The formation of twin vortices leads to a reduction in the air drag force exerted on the system, and thereby the system attains a higher flow velocity, i.e., a higher conversion rate of gravitational potential energy to the kinetic energy of the particle motion. The higher conversion rate of potential energy thus feeds back to the development of the vortex motion, resulting in the twin vortex formation.     

Particle Transport with Asymmetric Unbiased Forces and Entropic Barrier

Transport of a Brownian particle moving along the axis of a three-dimensional （313） symmetric and periodic tube is investigated in the presence of asymmetric unbiased forces. It is found that in the presence of entropic barrier, the asymmetry of the unbiased forces is a way of inducing a net particle current. The particle current is a peaked function of temperature, which indicates that the thermal noise may facilitate the transport even in the presence of entropic barrier. There exists an optimized radius at the bottleneck at which the particle current takes its maximum value. The current can be influenced by the slope of tube walls and there exists an optimized slope for which the particle current takes its maximum value.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier–Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.     

Simulating Free Surface Flows with SPH

The SPH (smoothed particle hydrodynamics) method is extended to deal with free surface incompressible flows. The method is easy to use, and examples will be given of its application to a breaking dam, a bore, the simulation of a wave maker, and the propagation of waves towards a beach. Arbitrary moving boundaries can be included by modelling the boundaries by particles which repel the fluid particles. The method is explicit, and the time steps are therefore much shorter than required by other less flexible methods, but it is robust and easy to program.     

Dynamics in a discontinuous field: The smooth Fermi piston

In this paper we examine a version of the Fermi piston with a discontinuous, but nonimpulsive, periodic driving force. The dynamics of a particle moving in one spatial dimension are studied using a combination of numerical and analytical techniques. The configuration space of the particle is divided into two regions of constant acceleration that are of equal magnitude and opposite direction. The point of discontinuity F(t) dividing the regions changes periodically in time. The method of surface-of-section is used to study the phase space (phi(n), v(n)), where phi(n) is the phase of the driving function and v(n) is the velocity of the particle at the nth encounter between the particle and boundary. We show that it is not possible to stochastically drive up the energy indefinitely except for the cases where F is discontinuous, or dF/dt is not finite everywhere. In addition, we find a new mechanism, other than KAM tori, for segmenting the phase space. As in the KAM picture, the central cause of the new behavior is resonance between the natural period of the particle and the period of the driving force. The boundaries to diffusion for continuous driving functions result from parabolic fixed points that span the entire phase range.     

Scalar Statistics along Inertial Particle Trajectory in Isotropic Turbulence

The statistics of a passive scalar along inertial particle trajectory in homogeneous isotropic turbulence with a mean scalar gradient is investigated by using direct numerical simulation. We are interested in the influence of particle inertia on such statistics, which is crucial for further understanding and development of models in non-isothermal gas-particle flows. The results show that the scalar variance along particle trajectory decreases with the increasing particle inertia firstly; when the particle＇s Stokes number St is less than 1.0, it reaches the minimal value when St is around 1.0, then it increases if St increases further. However, the scalar dissipation rate along the particle trajectory shows completely contrasting behavior in comparison with the scalar variance. The mechanical-to-thermal time scale ratios averaged along particle, （r）p, are approximately two times smaller than that computed in the Eulerian frame r, and stay at nearly 1.77 with a weak dependence on particle inertia. In addition, the correlations between scalar dissipation and flow structure characteristics along particle trajectories, such as strain and vorticity, are also computed, and they reach their maximum and minimum, 0.31 and 0.25, respectively, when St is around 1.0.     

A quasi-elastic regime for vibrated granular gases

Using simple scaling arguments and two-dimensional numerical simulations of a granular gas excited by vibrating one of the container boundaries, we study a double limit of small 1-r and large L, where r is the restitution coefficient and L the size of the container. We show that if the particle density n₀ and (1-r²)(n₀ Ld) where d is the particle diameter, are kept constant and small enough, the granular temperature, i.e. the mean value of the kinetic energy per particle, 〈E 〉/N, tends to a constant whereas the mean dissipated power per particle, 〈D 〉/N, decreases like when N increases, provided that (1-r²)(n₀ Ld)² < 1. The relative fluctuations of E, D and the power injected by the moving boundary, I, have simple properties in that regime. In addition, the granular temperature can be determined from the fluctuations of the power I(t) injected by the moving boundary.