Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Semi‐strong and strong solutions for variable density asymmetric fluids in unbounded domains

We establish the existence of local in time semi‐strong solutions and global in time strong solutions for the system of equations describing flows of viscous and incompressible asymmetric fluids with variable density in general three‐dimensional domains with boundary uniformly of class C³. Under suitable assumptions, uniqueness of local semi‐strong solutions is also proved. Copyright © 2016 John Wiley & Sons, Ltd.     

Non‐homogeneous Navier–Stokes systems with order‐parameter‐dependent stresses

We consider the Navier–Stokes system with variable density and variable viscosity coupled to a transport equation for an order‐parameter c. Moreover, an extra stress depending on c and ?c, which describes surface tension like effects, is included in the Navier–Stokes system. Such a system arises, e.g. for certain models of granular flows and as a diffuse interface model for a two‐phase flow of viscous incompressible fluids. The so‐called density‐dependent Navier–Stokes system is also a special case of our system. We prove short‐time existence of strong solution in L^q‐Sobolev spaces with q>d. We consider the case of a bounded domain and an asymptotically flat layer with a combination of a Dirichlet boundary condition and a free surface boundary condition. The result is based on a maximal regularity result for the linearized system. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of Two Linearized Problems Modeling Viscous Two‐Layer Flows

Two problems that appear in the linearization of certain free boundary value problems of the hydrodynamics of two viscous fluids are studied in the strip‐like domain Π = {x = (x₁, x₂) ∈ ℝ² : x₁ ∈ ℝ¹, (0 < x₂ < h_*) ∨ (h_* < x₂ < 1)}. The first problem arises in the linearization of a two‐layer flow down a geometrically perturbed inclined plane. The second one appears after the linearization of a two‐layer flow in a geometrically perturbed inclined channel with one moving (smooth) wall. For this purpose the unknown flow domain was mapped onto the double strip Π. The arising linear elliptic problems contain additional unknown functions in the boundary conditions. The paper is devoted to the investigation of these boundary problems by studying the asymptotics of the eigenvalues of corresponding operator pencils. It can be proved that the boundary value problems are uniquely solvable in weighted Sobolev spaces with exponential weight. The study of the full (nonlinear) free boundary value problems will be the topic of a forthcoming paper.     

The second‐order upwind finite difference fractional steps method for moving boundary value problem of oil‐water percolation

The research of the three‐dimensional (3D) compressible miscible (oil and water) displacement problem with moving boundary values is of great value to the history of oil‐gas transport and accumulation in basin evolution, as well as to the rational evaluation in prospecting and exploiting oil‐gas resources, and numerical simulation of seawater intrusion. The mathematical model can be described as a 3D‐coupled system of nonlinear partial differential equations with moving boundary values. For a generic case of 3D‐bounded region, a kind of second‐order upwind finite difference fractional steps schemes applicable to parallel arithmetic is put forward. Some techniques, such as the change of variables, calculus of variations, and the theory of a priori estimates, are adopted. Optimal order estimates in l² norm are derived for the errors in approximate solutions. The research is important both theoretically and practically for model analysis in the field, for model numerical method and for software development. Thus, the well‐known problem has been solved.Copyright © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1103–1129, 2014     

On stationary two‐dimensional flows around a fast rotating disk

We study the two‐dimensional stationary Navier–Stokes equations describing flows around a rotating disk. The existence of unique solutions is established for any rotating speed, and qualitative effects of a large rotation are described precisely by exhibiting a boundary layer structure and an axisymmetrization of the flow.     

Classroom Notes

This note gives a simple‐minded approach to the two‐dimensional boundary layer equations. The pressure is eliminated from the equations of motion and the resulting equation is simplified by assuming that certain derivatives in the direction of the boundary are small compared with those at right angles to it. The simplified equation is then integrated to give a single boundary layer equation which, together with the stress rate of strain law and the continuity equation, is sufficient (in theory at least) to predict the flow.

The boundary layer equation as given does not depend on a particular form for the stress rate of strain law and could possibly form the basis for a non‐Newtonian investigation. The viscous boundary layer is given as a special case.     

On uniqueness and time regularity of flows of power‐law like non‐Newtonian fluids

Large class of non‐Newtonian fluids can be characterized by index p, which gives the growth of the constitutively determined part of the Cauchy stress tensor. In this paper, the uniqueness and the time regularity of flows of these fluids in an open bounded three‐dimensional domain is established for subcritical ps, i.e. for p>11/5. Our method works for ‘all’ physically relevant boundary conditions, the Cauchy stress need not be potential and it may depend explicitly on spatial and time variable. As a simple consequence of time regularity, pressure can be introduced as an integrable function even for Dirichlet boundary conditions. Moreover, these results allow us to define a dynamical system corresponding to the problem and to establish the existence of an exponential attractor. Copyright © 2010 John Wiley & Sons, Ltd.     

Well‐posedness in smooth function spaces for moving‐boundary 1‐D compressible euler equations in physical vacuum

The free‐boundary compressible one‐dimensional Euler equations with moving physical vacuum boundary are a system of hyperbolic conservation laws that are both characteristic and degenerate. The physical vacuum singularity (or rate of degeneracy) requires the sound speed $c^2= \gamma \rho^{ \gamma -1}$ to scale as the square root of the distance to the vacuum boundary and has attracted a great deal of attention in recent years. We establish the existence of unique solutions to this system on a short time interval, which are smooth (in Sobolev spaces) all the way to the moving boundary. The proof is founded on a new higher‐order, Hardy‐type inequality in conjunction with an approximation of the Euler equations consisting of a particular degenerate parabolic regularization. Our regular solutions can be viewed as degenerate viscosity solutions. © 2010 Wiley Periodicals, Inc.     

Local Interaction Theory for Laminar Transonic Flows in Slender Channels

Transonic flows through channels so narrow that the classical boundary layer approach fails are considered. As a consequence the properties of the inviscid core and the viscosity dominated boundary layer region can no longer be determined in subsequent steps but have to be calculated simultaneously. The resulting viscous inviscid interaction problem for weakly three dimensional laminar flows is formulated for perfect gases under the requirement that the channel is sufficiently narrow so that the flow outside the viscous wall layers becomes planar in the leading order approximation. Representative solutions for subsonic as well as for supersonic flows disturbed by three dimensional surface mounted obstacles will be presented. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Local well‐posedness and zero‐α limit for the Euler‐α equations

In this paper, we prove the local‐in‐time existence and a blow‐up criterion of solutions in the Besov spaces for the Euler‐α equations of inviscid incompressible fluid flows in . We also establish the convergence rate of the solutions of the Euler‐α equations to the corresponding solutions of the Euler equations as the regularization parameter α approaches 0 in . Copyright © 2016 John Wiley & Sons, Ltd.     

Analytic Study of Magnetohydrodynamic Flow and Boundary Layer Control Over a Wedge

A genuine variational principle developed by Gyarmati, in the field of thermodynamics of irreversible processes unifying the theoretical requirements of technical, environmental and biological sciences is employed to study the effects of uniform suction and injection on MHD flow adjacent to an isothermal wedge with pressure gradient in the presence of a transverse magnetic field. The velocity distribution inside the boundary layer has been considered as a simple polynomial function and the variational principle is formulated. The Euler-Lagrange equation is reduced to a simple polynomial equation in terms of momentum boundary layer thickness. The velocity profiles, displacement thickness and the coefficient of skin friction are calculated for various values of wedge angle parameter m, magnetic parameter ξ and suction/injection parameter H. The present results are compared with known available results and the comparison is found to be satisfactory. The present study establishes high accuracy of results obtained by this variational technique.     

Boundary layers for parabolic perturbations of quasi‐linear hyperbolic problems

In this paper, we study the asymptotic relation between the solutions to the one‐dimensional viscous conservation laws with the Dirichlet boundary condition and the associated inviscid solution. We assume that the viscosity matrix is positive definite, then we prove the existence and the stability of the weak boundary layers by discussing nonlinear well‐posedness of the inviscid flow with certain boundary conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Regularity in Sobolev spaces of steady flows of fluids with shear‐dependent viscosity

The system is considered on a bounded three‐dimensional domain under no‐stick boundary value conditions, where S has p‐structure for some p<2 and D ( u ) is the symmetrized gradient of u . Various regularity results for the velocity u and the pressure π in fractional order Sobolev and Nikolskii spaces are obtained. Copyright © 2006 John Wiley & Sons, Ltd.     

Well-posedness for two-dimensional steady supersonic Euler flows past a Lipschitz wedge

For a supersonic Euler flow past a straight-sided wedge whose vertex angle is less than the extreme angle, there exists a shock-front emanating from the wedge vertex, and the shock-front is usually strong especially when the vertex angle of the wedge is large. In this paper, we establish the L¹ well-posedness for two-dimensional steady supersonic Euler flows past a Lipschitz wedge whose boundary slope function has small total variation, when the total variation of the incoming flow is small. In this case, the Lipschitz wedge perturbs the flow, and the waves reflect after interacting with the strong shock-front and the wedge boundary. We first obtain the existence of solutions in BV when the incoming flow has small total variation by the wave front tracking method and then establish the L¹ stability of the solutions with respect to the incoming flows. In particular, we incorporate the nonlinear waves generated from the wedge boundary to develop a Lyapunov functional between two solutions containing strong shock-fronts, which is equivalent to the L¹ norm, and prove that the functional decreases in the flow direction. Then the L¹ stability is established, so is the uniqueness of the solutions by the wave front tracking method. Finally, the uniqueness of solutions in a broader class, the class of viscosity solutions, is also obtained.     

Stability of global and exponential attractors for a three‐dimensional conserved phase‐field system with memory

We consider a conserved phase‐field system on a tri‐dimensional bounded domain. The heat conduction is characterized by memory effects depending on the past history of the (relative) temperature ?, which is represented through a convolution integral whose relaxation kernel k is a summable and decreasing function. Therefore, the system consists of a linear integrodifferential equation for ?, which is coupled with a viscous Cahn–Hilliard type equation governing the order parameter χ. The latter equation contains a nonmonotone nonlinearity ? and the viscosity effects are taken into account by a term ?αΔ?_tχ, for some α?0. Rescaling the kernel k with a relaxation time ε>0, we formulate a Cauchy–Neumann problem depending on ε and α. Assuming a suitable decay of k, we prove the existence of a family of exponential attractors {?_α,ε} for our problem, whose basin of attraction can be extended to the whole phase–space in the viscous case (i.e. when α>0). Moreover, we prove that the symmetric Hausdorff distance of ?_α,ε from a proper lifting of ?_α,0 tends to 0 in an explicitly controlled way, for any fixed α?0. In addition, the upper semicontinuity of the family of global attractors {??_α,ε} as ε→0 is achieved for any fixed α>0. Copyright © 2009 John Wiley & Sons, Ltd.     

The lower bounds of life span of classical solutions to one‐dimensional initial‐Neumann boundary value problems for general quasilinear wave equations

In this paper, we will study the lower bounds of the life span (the maximal existence time) of solutions to the initial‐boundary value problems with small initial data and zero Neumann boundary data on exterior domain for one‐dimensional general quasilinear wave equations u_tt?u_xx=b(u,Du)u_xx+F(u,Du). Our lower bounds of the life span of solutions in the general case and special case are shorter than that of the initial‐Dirichlet boundary value problem for one‐dimensional general quasilinear wave equations. We clarify that although the lower bounds in this paper are same as that in the case of Robin boundary conditions obtained in the earlier paper, however, the results in this paper are not the trivial generalization of that in the case of Robin boundary conditions because the fundamental Lemmas 2.4, 2.5, 2.6, and 2.7, that is, the priori estimates of solutions to initial‐boundary value problems with Neumann boundary conditions, are established differently, and then the specific estimates in this paper are different from that in the case of Robin boundary conditions. Another motivation for the author to write this paper is to show that the well‐posedness of problem 1.1 is the essential precondition of studying the lower bounds of life span of classical solutions to initial‐boundary value problems for general quasilinear wave equations. The lower bound estimates of life span of classical solutions to initial‐boundary value problems is consistent with the actual physical meaning. Finally, we obtain the sharpness on the lower bound of the life span 1.8 in the general case and 1.10 in the special case. Copyright © 2015 John Wiley & Sons, Ltd.     

Interface boundary value problems of Robin‐transmission type for the Stokes and Brinkman systems on n‐dimensional Lipschitz domains: applications

In this paper, we describe a layer potential analysis in order to show an existence result for an interface boundary value problem of Robin‐transmission type for the Stokes and Brinkman systems on Lipschitz domains in Euclidean setting, when the given boundary data belong to some L^p or Sobolev spaces associated to such domains. Applications related to an exterior three‐dimensional Stokes flow past two concentric porous spheres with stress jump conditions on the fluid‐porous interface are also considered. Copyright © 2013 John Wiley & Sons, Ltd.     

Three‐dimensional internal gravity‐capillary waves in finite depth

We consider three‐dimensional inviscid‐irrotational flow in a two‐layer fluid under the effects of gravity and surface tension, where the upper fluid is bounded above by a rigid lid and the lower fluid is bounded below by a flat bottom. We use a spatial dynamics approach and formulate the steady Euler equations as an infinite‐dimensional Hamiltonian system, where an unbounded spatial direction x is considered as a time‐like coordinate. In addition, we consider wave motions that are periodic in another direction z. By analyzing the dispersion relation, we detect several bifurcation scenarios, two of which we study further: a type of 00(is)(iκ₀) resonance and a Hamiltonian Hopf bifurcation. The bifurcations are investigated by performing a center‐manifold reduction, which yields a finite‐dimensional Hamiltonian system. For this finite‐dimensional system, we establish the existence of periodic and homoclinic orbits, which correspond to, respectively, doubly periodic travelling waves and oblique travelling waves with a dark or bright solitary wave profile in the x direction. The former are obtained using a variational Lyapunov‐Schmidt reduction and the latter by first applying a normal form transformation and then studying the resulting canonical system of equations.     

Time‐weighted blow‐up rates and pointwise profile for single‐point blow‐up solutions in reaction–diffusion equations

This paper deals with asymptotic behavior for blow‐up solutions to time‐weighted reaction–diffusion equations u_t=Δu+e^αtv^p and v_t=Δv+e^βtu^q, subject to homogeneous Dirichlet boundary. The time‐weighted blow‐up rates are defined and obtained by ways of the scaling or auxiliary‐function methods for all α, . Aiding by key inequalities between components of solutions, we give lower pointwise blow‐up profiles for single‐point blow‐up solutions. We also study the solutions of the system with variable exponents instead of constant ones, where blow‐up rates and new blow‐up versus global existence criteria are obtained. Time‐weighted functions influence critical Fujita exponent, critical Fujita coefficient and formulae of blow‐up rates, but they do not limit the order of time‐weighted blow‐up rates and pointwise profile near blow‐up time. Copyright © 2017 John Wiley & Sons, Ltd.