Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Experimental and numerical investigation of inclined air/SF_6 interface instability under shock wave

The shock tube experiments of inclined air/SF6 interface instability under the shock wave with the Mach numbers 1.23 and 1.41 are conducted. The numerical simulation is done with the parallel algorithm and the multi-viscous-fluid and turbulence （MVFT） code of the large-eddy simulation （LES）. The developing process of the interface accelerated by the shock wave is reproduced by the simulations. The complex wave structures, e.g., the propagation, refraction, and reflection of the shock wave, are clearly revealed in the flows. The simulated evolving images of the interface are consistent with the experimental ones. The simulated width of the turbulent mixing zone （TMZ） and the displacements of the bubble and the spike also agree well with the experimental data. Also, the reliability and effectiveness of the MVFT in simulating the problem of interface instability are validated. The more energies are injected into the TMZ when the shock wave has a larger Mach number. Therefore, the perturbed interface develops faster.

Large eddy simulation of boundary layer flow under cnoidal waves

Water waves in coastal areas are generally non- linear, exhibiting asymmetric velocity profiles with different amplitudes of crest and trough. The behaviors of the bound- ary layer under asymmetric waves are of great significance for sediment transport in natural circumstances. While pre- vious studies have mainly focused on linear or symmetric waves, asymmetric wave-induced flows remain unclear, par- ticularly in the flow regime with high Reynolds numbers. Taking cnoidal wave as a typical example of asymmetric waves, we propose to use an infinite immersed plate oscillat- ing cnoidally in its own plane in quiescent water to simulate asymmetric wave boundary layer. A large eddy simulation approach with Smagorinsky subgrid model is adopted to investigate the flow characteristics of the boundary layer. It is verified that the model well reproduces experimental and theoretical results. Then a series of numerical experiments are carried out to study the boundary layer beneath cnoidal waves from laminar to fully developed turbulent regimes at high Reynolds numbers, larger than ever studied before. Results of velocity profile, wall shear stress, friction coeffi- cient, phase lead between velocity and wall shear stress, and the boundary layer thickness are obtained. The dependencies of these boundary layer properties on the asymmetric degree and Reynolds number are discussed in detail.

Reynolds number dependence of particle transport in a 90°bend with electrostatic effects

In this work,particle transport in a fully developed turbulent 90°bend flow at the"electrostatic equi-librium"state is simulated using large eddy simulation coupled with Lagrangian particle tracking technique.The flow Reynolds numbers(based on bulk velocity)considered is from 34000 to 58000.Three particle size 5,10 and 50μm are considered and their corresponding St number are from 2.2 to 547.Simulation results of the bend flow agree well with experimental data.The electrostatic field formed in the bend is symmetric in the spanwise direction but asymmetric in the vertical direction and radial direction,which is independent of Reynolds number.The minimum electrostatic field strength occurs at r/ra=0.25 near the inner wall of the bend.Particles transported in a bend gradually accumulate near the wall due to turbophoresis,such trend is improved by electrostatics.In addition,under the effect of electrostatics,the plume pattern of particle distribution disappeared.Particle concentration at the inner wall of the bend is higher than that at the outer wall,which depends on the combined effect of elec-trostatics and Dean vortices in the bend.