METHODS AND APPLICATIONS OF ANALYSIS. 2000 International Press Vol. 7, No. 3, pp. 495-510, September 2000 007

In this paper, we present a novel pairwise-force smoothed particle hydrodynamics (PF-SPH) model to allow modeling of various interactions at interfaces in real time. Realistic capture of interactions at interfaces is a challenging problem for SPH-based simulations, especially for scenarios involving multiple interactions at different interfaces. Our PF-SPH model can readily handle multiple kinds of interactions simultaneously in a single simulation; its basis is to use a larger support radius than that used in standard SPH. We adopt a novel anisotropic filtering term to further improve the performance of interaction forces. The proposed model is stable; furthermore, it avoids the particle clustering problem which commonly occurs at the free surface. We show how our model can be used to capture various interactions. We also consider the close connection between droplets and bubbles, and show how to animate bubbles rising in liquid as well as bubbles in air. Our method is versatile, physically plausible and easy-to-implement. Examples are provided to demonstrate the capabilities and effectiveness of our approach.

In this paper, we prove the global existence of solutions with analytic regularity to the 2D magnetohydrodynamic (MHD) boundary layer equations in the mixed Prandtl and Hartmann regime derived by formal multiscale expansion in [D. Gerard-Varet and M. Prestipino, <i>Z. Angew. Math. Phys.</i>, 68 (2017), 76]. The analysis shows that the combined effect of the magnetic diffusivity and transverse magnetic field on the boundary leads to a linear damping on the tangential velocity field near the boundary. And this damping effect yields the global-in-time analytic norm estimate in the tangential space variable on the perturbation of the classical steady Hartmann profile.

In the long-time scale, we consider the fluid dynamical limits for the kinetic equations when the fluctuation is decomposed into even and odd parts with respect to the microscopic velocity with different scalings. It is shown that when the background state is an absolute Maxwellian, the limit fluid dynamical equations are the incompressible Navier-Stokes equations with viscous heating. This is different from the case when the even and odd parts of the fluctuation have the same scaling where the standard incompressible Navier-Stokes equations without viscous heating are obtained. On the other hand, when the background is a local Maxwellian, it is shown that the above even-odd decomposition leads to a non-classical fluid dynamical system without viscous heating which has been used to describe the ghost effect in the kinetic theory. In addition, the above even-odd decomposition is justified rigorously for the Boltzmann equation for the former case when the background is an absolute Maxwellian.

Flow reversals in two-dimensional Rayleigh-Benard convection led by non-Oberbeck-Boussinesq (NOB) effects due to large temperature differences arc studied by direct numerical simulation. Perfect gas is chosen as the working fluid and the Prandtl number is 0.71 for the reference state. 11 NOB effects are included, the flow pattern P-11 with only one dominant roll often becomes unstable by the growth of the cold corner roll, which sometimes results in cession-led flow reversals. By exploiting the vorticity transport equation, it is found that the asymmetries of buoyancy and viscous forces are responsible for the growth of the cold corner roll because both such asymmetries cause an imbalance between the corner rolls and the large-scale circulation (I.SC). The buoyancy force near the cold wall increases and decreases near the hot wall originating from the temperature-dependent isobaric thermal expansion coefficient alpha = 1/T if NOB effects are included; Moreover, the decreased dissipation due to lower viscosity is favourable for the growth of the cold corner roll, while the increased viscosity further suppresses the growth of the hot corner roll. Finally, it is found that the boundary layer near the cold wall plays an important role in the mass transport from LSC to corner rolls subject to mass conservation.