Magnetohydrodynamic natural convection heat transfer in a rotating, differentially heated enclosure is studied numerically in this article. The governing equations are in velocity, pressure and temperature formulation and solved using the staggered grid arrangement together with MAC method. The governing parameters considered are the Hartmann number, 0≤Ha≤70, the inclination angle of the magnetic field, 0°≤θ 90°, the Taylor number, 8.9 x 104≤Ta≤1.1 x 106 and the centrifugal force is smaller than the Coriolis force and the both forces were kept below the buoyancy force. It is found that a sufficiently large Lorentz force neutralizes the effect of buoyancy, inertial and Coriolis forces. Horizontal or vertical direction of the magnetic field was most effective in reducing the global heat transfer.

This paper reports an improved constitutive model for the shear creep behavior of offwhite marbles which are selected from slope and underground cavern and contain green schist’s weak structural planes. The shear creep behavior of the samples is characterized using the rheological tests. Based on the experimental measurements on mechanical properties under different normal stress conditions, an improved model is proposed to analyze the experimental results. It is demonstrated from a further discussion that such model can reflect the non-linear creep characteristics of structural planes, and especially, it is suitable for description of the viscoelastic and viscoplastic deformation behavior of structural planes.

This study presents a modification of the central-upwind Kurganov scheme for approximating the solution of the 2D Euler equation. The prototype, extended from a 1D model, reduces substantially less dissipation than expected. The problem arises from over-restriction of some slope limiters, which keep slopes between interfaces of cells to be Total-Variation-Diminishing. This study reports the defect and presents a re-derived optimal formula. Numerical experiments highlight the significance of this formula, especially in long-time, large-scale simulations.

In this paper, numerical sensitivity analysis with respect to the Reynolds number for the flow past obstacle problem is presented. To carry out such analysis, at each time step, we need to solve the incompressible Navier-Stokes equations on irregular domains twice, one for the primary variables; the other is for the sensitivity variables with homogeneous boundary conditions. The Navier-Stokes solver is the augmented immersed interface method for Navier-Stokes equations on irregular domains. One of the most important contribution of this paper is that our analysis can predict the critical Reynolds number at which the vortex shading begins to develop in the wake of the obstacle. Some interesting experiments are shown to illustrate how the critical Reynolds number varies with different geometric settings.

In this paper, new finite difference methods based on the augmented immersed interface method (IIM) are proposed for simulating an inextensible moving interface in an incompressible two-dimensional flow. The mathematical models arise from studying the deformation of red blood cells in mathematical biology. The governing equations are incompressible Stokes or Navier-Stokes equations with an unknown surface tension, which should be determined in such a way that the surface divergence of the velocity is zero along the interface. Thus, the area enclosed by the interface and the total length of the interface should be conserved during the evolution process. Because of the nonlinear and coupling nature of the problem, direct discretization by applying the immersed boundary or immersed interface method yields complex nonlinear systems to be solved. In our new methods, we treat the unknown surface tension as an augmented variable so that the augmented IIM can be applied. Since finding the unknown surface tension is essentially an inverse problem that is sensitive to perturbations, our regularization strategy is to introduce a controlled tangential force along the interface, which leads to a least squares problem. For Stokes equations, the forward solver at one time level involves solving three Poisson equations with an interface. For Navier-Stokes equations, we propose a modified projection method that can enforce the pressure jump condition corresponding directly to the unknown surface tension. Several numerical experiments show good agreement with other results in the literature and reveal some interesting phenomena.