Forced two layer beta-plane quasigeostrophic flow. Part I: Long-time existence and uniqueness of weak solutions

We consider a model of quasigeostrophic turbulence that has proven useful in theoretical studies of large scale heat transport and coherent structure formation in planetary atmospheres and oceans. The model consists of a coupled pair of hyperbolic PDEs with a forcing which represents domain-scale thermal energy source. Although the use to which the model is typically put involves gathering information from very long numerical integrations, little of a rigorous nature is known about long-time properties of solutions to the equations. In this first paper we define a notion of weak solution, and show using Galerkin methods the long-time existence and uniqueness of such solutions.     

Calculation of critical properties for the anisotropic two-layer Ising model on the Kagome lattice: Cellular automata approach

The critical point (K_c) of the two-layer Ising model on the Kagome lattice has been calculated with a high precision, using the probabilistic cellular automata with the Glauber algorithm. The critical point is calculated for different values of the inter- and intra-layer couplings (K₁≠K₂≠K₃≠K_z), where K₁, K₂ and K₃ are the nearest-neighbor interactions within each layer in the 1, 2 and 3 directions, respectively, and K_z is the intralayer coupling. A general ansatz equation for the critical point is given as a function of the inter- and intra-layer interactions, ξ=K₃/K₁,σ=K₂/K₁ and ω=K_z/K₁ for the one- and two-layer Ising models on the Kagome lattice.     

The temperature-dependent phase diagrams of the spin-3/2 Ising model on a FM/AFM two-layer lattice with a crystal field

The temperature-dependent phase diagrams of the spin-3/2 Ising model on a two-layer Bethe lattice with ferromagnetic (FM)/antiferromagnetic (AFM) intra-layer and either FM or AFM type inter-layer interactions are investigated under a constant magnetic field (H) and in the presence of a crystal field (D) by using exact recursion equations in a pairwise approach for coordination numbers q=3,4 and 6, in detail. In the light of the ground-state (GS) phase diagrams, the temperature-dependent phase diagrams of the model are obtained by studying the thermal variations of the order parameters, response functions and free energy. Then, they are illustrated on the (kT/J₁,J₃/J₁) and (kT/J₁,J₂/J₁) planes for the given system parameters. It is observed that the system exhibits first- and second-order phase transitions for all q values, and hence, in some cases, tricritical points. The existence of critical-end points and that of isolated points are also observed. The re-entrant behavior owes its presence to the two Néel temperatures, T_N, that are present for all q.     

Complex turbulent compressible flow computation using a two-layer approach

A two-layer approach is proposed to compute complex flows including separations. For high- and low-Reynolds-number regions we use a two-equation k-? model and a one-equation k-L model respectively. A robust algorithm is proposed for the treatment of the convective part of the turbulence equations. Several complex configurations including separations are computed.     

Numerical simulation and experimental verification of silicone oil flow over magnetic fluid under applied magnetic field

Two-layer flow of magnetic fluid and non-magnetic silicone oil was simulated numerically. The continuity equation, momentum equations, kinematic equation, and magnetic potential equation were solved in two-dimensional Cartesian coordinate. PLIC （piecewise linear integration calculation） VOF （volume of fluid） scheme was employed to track the free interface. Surface tension was treated via a continuous surface force （CSF） model that ensures robustness and accuracy. The influences of applied magnetic field, inlet velocity profile, initial surface disturbance of interface and surface tension were analyzed. The computed interface shapes at different conditions were compared with experimental observation.     

Suppression or enhancement of interfacial instability in two-layer plane Couette flow of FENE-P fluids past a deformable solid layer

The linear stability of two-layer plane Couette flow of FENE-P fluids past a deformable solid layer is analyzed in order to examine the effect of solid deformability on the interfacial instability due to elasticity and viscosity stratification at the two-fluid interface. The solid layer is modeled using both linear viscoelastic and neo-Hookean constitutive equations. The limiting case of two-layer flow of upper-convected Maxwell (UCM) fluids is used as a starting point, and results for the FENE-P case are obtained by numerically continuing the UCM results for the interfacial mode to finite values of the chain extensibility parameter. For the case of two-layer plane Couette flow past a rigid solid surface, our results show that the finite extensibility of the polymer chain significantly alters the neutral stability boundaries of the interfacial instability. In particular, the two-layer Couette flow of FENE-P fluids is found to be unstable in a larger range of nondimensional parameters when compared to two-layer flow of UCM fluids. The presence of the deformable solid layer is shown to completely suppress the interfacial instability in most of the parameter regimes where the interfacial mode is unstable, while it could have a completely destabilizing effect in other parameter regimes even when the interfacial mode is stable in rigid channels. When compared with two-layer UCM flow, the two-layer FENE-P case is found in general to require solid layers with relatively lower shear modulii in order to suppress the interfacial instability. The results from the linear elastic solid model are compared with those obtained using the (more rigorous) neo-Hookean model for the solid, and good agreement is found between the two models for neutral stability curves pertaining to the two-fluid interfacial mode. The present study thus provides an important extension of the earlier analysis of two-layer UCM flow [V. Shankar, Stability of two-layer viscoelastic plane Couette flow past a deformable solid layer: implications of fluid viscosity stratification, J. Non-Newtonian Fluid Mech. 125 (2005) 143–158] to more accurate constitutive models for the fluid and solid layers, and reaffirms the central conclusion of instability suppression in two-layer flows of viscoelastic fluids by soft elastomeric coatings in more realistic settings.     

Roll waves structure in two-layer Hele–Shaw flows

The mathematical model of inhomogeneous fluid motion in a Hele–Shaw cell is proposed. Based on this model the equations for describing two-layer flows and development of roll waves at the interface are derived. Conditions of roll waves existence are formulated in terms of Whitham criterion. Numerical calculations of the interface position are provided. It is shown that small perturbations of the interface in the inlet section of the channel lead to the roll waves for certain parameters of the flow. Two-parametric class of exact solutions corresponding to the roll waves regime is obtained. Diagrams of critical depths of roll waves development are constructed.     

Two-layer film flow over a rotating disk

Unsteady two-layer liquid film flow on a horizontal rotating disk is analyzed using asymptotic method for small values of Reynolds number. This analysis of non-linear evolution equation elucidates how a two-layer film of uniform thickness thins when the disk is set in uniform rotation. It is observed that the final film thickness attains an asymptotic value at large time. It is also established that viscous force dominates over centrifugal force and upper layer thins faster than lower layer at large time.     

Water wave interaction with a sphere in a two-layer fluid flowing through a channel of finite depth

Using linear water wave theory, we consider a three-dimensional problem involving the interaction of waves with a sphere in a fluid consisting of two layers with the upper layer and lower layer bounded above and below, respectively, by rigid horizontal walls, which are approximations of the free surface and the bottom surface; these walls can be assumed to constitute a channel. The effects of surface tension at the surface of separation is neglected. For such a situation time-harmonic waves propagate with one wave number only, unlike the case when one of the layers is of infinite depth with the waves propagating with two wave numbers. Method of multipole expansions is used to find the particular solutions for the problems of wave radiation and scattering by a submerged sphere placed in either of the upper or lower layer. The added-mass and damping coefficients for heave and sway motions are derived and plotted against various values of the wave number. Similarly the exciting forces due to heave and sway motions are evaluated and presented graphically. The features of the results find good agreement with previously available results from the point of view of physical interpretation.     

Prediction of turbulent wall shear flows directly from wall

The fully elliptic Reynolds-averaged Navier–Stokes equations have been used together with Lam and Bremhorst's low-Reynolds-number model, Chen and Patel's two-layer model and a two-point wall function method incorporated into the standard k-? model to predict channel flows and a backward-facig step flow. These flows enable the evaluation of the performance of different near-wall treatments in flows involving streamwise and normal pressure gradients, flows with separation and flows with non-equilibrium turbulence characteristics. Direct numerical simulation (DNS) of a channel flow with Re =3200 further provides the detailed budgets of each modelling term of the k and ?-transport equations. Comparison of model results with DNS data to evaluate the performance of each modelling term is also made in the present study. It is concluded that the low-Reynolds-number model has wider applicability and performs better than the two-layer model and wall function approaches. Comparison with DNS data further shows that large discrepancies exist between the DNS budgets and the modelled production and destruction terms of the ? equation. However, for simple channel flow the discrepancies are similar in magnitude but opposite in sign, so they are cancelled by each other. This may explain why, even when employing such an inaccurately modelled ?-equation, one can still predict satisfactorily some simple turbulent flows.