We report on a three-dimensional direct numerical simulation study of flow structure and heat transport in the annular centrifugal Rayleigh-Bdnard convection (ACRBC) system, with cold inner and hot outer cylinders corotating axially, for the Rayleigh number range Ra is an element of [10(6), 10(8)] and radius ratio range eta = R-i/R-o is an element of [0.3, 0.9] (R-i and R-o are the radius of the inner and outer cylinders, respectively). This study focuses on the dependence of flow dynamics, heat transport and asymmetric mean temperature fields on the radius ratio eta. For the inverse Rossby number Ro(-1) = 1, as the Coriolis force balances inertial force, the flow is in the inertial regime. The mechanisms of zonal flow revolving in the prograde direction in this regime are attributed to the asymmetric movements of plumes and the different curvatures of the cylinders. The number of roll pairs is smaller than the circular roll hypothesis as the convection rolls are probably elongated by zonal flow. The physical mechanism of zonal flow is verified by the dependence of the drift frequency of the large-scale circulation (LSC) rolls and the space- and time-averaged azimuthal velocity on eta. The larger eta is, the weaker the zonal flow becomes. We show that the heat transport efficiency increases with eta. It is also found that the bulk temperature deviates from the arithmetic mean temperature and the deviation increases as eta decreases. This effect can be explained by a simple model that accounts for the curvature effects and the radially dependent centrifugal force in ACRBC.

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.

We study the translational and rotational dynamics of neutrally buoyant finite-size spheroids in hydrodynamic turbulence by means of fully resolved numerical simulations. We examine axisymmetric shapes, from oblate to prolate, and the particle volume dependences. We show that the accelerations and rotations experienced by non-spherical inertial-scale particles result from volume filtered fluid forces and torques, similar to spherical particles. However, the particle orientations carry signatures of preferential alignments with the surrounding flow structures, which are reflected in distinct axial and lateral fluctuations for accelerations and rotation rates. The randomization of orientations does not occur even for particles with volume-equivalent diameter size in the inertial range, here up to 60 dissipative units (eta) at Taylor-scale Reynolds number Re-lambda = 120. Additionally, we demonstrate that the role of fluid boundary layers around the particles cannot be neglected in reaching a quantitative understanding of particle statistical dynamics, as they affect the intensities of the angular velocities and the relative importance of tumbling with respect to spinning rotations. This study brings to the fore the importance of inertial-scale flow structures in homogeneous and isotropic turbulence and their impacts on the transport of neutrally buoyant bodies with sizes in the inertial range.

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.

For a fixed geometric configuration, hydrodynamic instabilities and bifurcation processes of laminar isothermal planar opposed jet flows with symmetric and slightly asymmetric inlet boundary conditions are investigated numerically using a high-resolution approach based on spectral element method. In current configuration, when inlet boundary conditions are symmetric, in the range of the Reynolds number considered (Re <= 200), multiple new symmetry-breaking bifurcations are observed and four new flow patterns are identified. Their hydrodynamic characteristics are analyzed, in particular their symmetries. In addition, the case that inlet boundary conditions are slightly asymmetric is investigated. It is found that bifurcation processes are extremely sensitive to this small symmetry-breaking imperfection and much different from those in the symmetric case. Furthermore, model equations are constructed by symmetry consideration to explain the numerical results based on hydrodynamic equations.

In this paper, we are concerned with the motion of electrically conducting fluid governed by the two-dimensional non-isentropic viscous compressible MHD system on the half plane with no-slip condition on the velocity field, perfectly conducting wall condition on the magnetic field and Dirichlet boundary condition on the temperature on the boundary. When the viscosity, heat conductivity and magnetic diffusivity coefficients tend to zero in the same rate, there is a boundary layer which is described by a Prandtl-type system. Under the non-degeneracy condition on the tangential magnetic field instead of monotonicity of velocity, by applying a coordinate transformation in terms of the stream function of magnetic field as motivated by the recent work [27], we obtain the local-in-time well-posedness of the boundary layer system in weighted Sobolev spaces.

We develop an efficient numerical method to compute single slit or double slit diffraction patterns from high frequency wave in inhomogeneous media. We approximate the high frequency asymptotic solution to the Helmholtz equation using the Eulerian Gaussian beam summation proposed in [20, 21]. The emitted rays from a slit are embedded in the phase space using an open segment. The evolution of this open curve is accurately computed using the recently developed Grid Based Particle Method [24] which results in a very efficient computational algorithm. Following the grid based particle method we proposed in [23, 24], we represent the open curve or the open surface by meshless Lagrangian particles sampled according to an underlying fixed Eulerian mesh. The end-points of the open curve are tracked explicitly and consistently with interior particles. To construct the overall wavefield, each of these sampling particles also carry necessary quantities that are obtained by solving advection-reaction equations. Numerical experiments show that the resulting method can model diffraction patterns in inhomogeneous media accurately, even in the occurrence of caustics.

Multiple-fluid interaction is an interesting and common visual phenomenon we often observe. In this paper we present an energybased Lagrangian method that expands the capability of existing multiple-fluid methods to handle various phenomena, including extraction, partial dissolution, etc. Based on our user-adjusted Helmholtz free energy functions, the simulated fluid evolves from high-energy states to low-energy states, allowing flexible capture of various mixing and unmixing processes. We also extend the original Cahn-Hilliard equation to gain abilities of simulating complex fluid-fluid interaction and rich visual phenomena such as motionrelated mixing and position based pattern. Our approach is easy to be integrated with existing state-of-the-art smooth particle hydrodynamic (SPH) solvers and can be further implemented on top of the position based dynamics (PBD) method, improving the stability and incompressibility of the fluid during Lagrangian simulation under large time steps. Performance analysis shows that our method is at least 4 times faster than the state-of-the-art multiple-fluid method. Examples are provided to demonstrate the new capability and effectiveness of our approach.

No abstract uploaded!

A key difficulty in the analysis and numerical approximation of the shallow water equations is the nonconservative product of measures due to the gravitational force acting on a sloped bottom. Solutions may be nonunique, and numerical schemes are not only consistent discretizations of the shallow water equations, but they also determine how to model the physics. Our derivation is based on a subcell reconstruction using infinitesimal singular layers at the cell boundaries, as inspired by S. Noelle, Y. Xing, and C.-W. Shu [J. Comput. Phys., 226 (2007), pp. 29-58]. One key step is to separate the singular measures. Another aspect is the reconstruction of the solution variables in the singular layers. We study three reconstructions. The first leads to the well-known scheme of Audusse et al., [SIAM J. Sci. Comput., 25 (2004), pp. 2050-2065], which introduces the hydrostatic reconstruction. The second is a modification proposed in [T. Morales de Luna, M. J. Castro D´ıaz, and C. Par´es, Appl. Math. Comput., 219 (2013), pp. 9012-9032], which analyzes whether a wave has enough energy to overcome a step. The third is our new scheme, which borrows its structure from the wet-dry front. For a number of cases discussed in recent years, where water runs down a hill, Audusse’s scheme converges slowly or fails. Morales’ scheme gives a visible improvement. Both schemes are clearly outperformed by our new scheme.