Singular perturbations for the exterior three-dimensional slow viscous flow problem

In this paper the rigorous justification of the formal asymptotic expansions constructed by the method of matched inner and outer expansions is established for the three-dimensional steady flow of a viscous, incompressible fluid past an arbitrary obstacle. The justification is based on the series representation of the solution to the Navier-Stokes equations due to Finn, and it involves the reductions of various exterior boundary value problems for the Stokes and Oseen equations to boundary integral equations of the first kind from which existence as well as asymptotic error estimates for the solutions are deduced. In particular, it is shown that the force exerted on the obstacle by the fluid admits the asymptotic representation F = A₀ + A₁Re + O(Re² ln Re⁻¹) as the Reynolds number Re → 0⁺, where the vectors A₀ and A₁ can be obtained from the method of matched inner and outer expansions.     

Simulation of a Nozzle-Jet Configuration Including the Acoustic Near-Field

A numerical simulation setup is presented which allows to study a circular jet flow configuration in which the nozzle is included in the simulation domain. Direct Numerical Simulations (DNS) are performed using up to 10^th order compact finite-difference schemes which are stabilized by applying a mild low-pass filter. A parallelization approach has been implemented which shows good weak and strong scaling behavior. At the inflow the Synthetic Eddy Method is employed to generate turbulent fluctuations in the nozzle boundary layer with prescribed statistics, which are imposed by a sponge (forcing) layer technique. Simulation results for the jet flow field obtained at Reynolds number Re_D = 19000 and a Mach number Ma = 0.9 as well as for the acoustic near-field are found to be in good agreement with recent nozzle-jet simulation results as well as experimental findings. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

DNS and Lie group analysis of turbulent boundarly layer flow

Lie group approach is used to derive new scaling laws for zero-pressure gradient turbulent boundary layer flow. A direct numerical simulation of the flow at Reynolds number Re_θ = 2240 was performed for the verify theoretical results. Navier-Stokes equations were numerically solved using spectral method with up to 160 million grid points. The numerical simulation shows validity of the theoretical results. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Boundary layers for cellular flows at high Péclet numbers

We analyze the behavior of solutions of steady advection‐diffusion problems in bounded domains with prescribed Dirichlet data when the Péclet number Pe ? 1 is large. We show that the solution converges to a constant in each flow cell outside a boundary layer of width O(?^1/2), ? = Pe^?1, around the flow separatrices. We construct an ?‐dependent approximate “water pipe problem” purely inside the boundary layer that provides a good approximation of the solution of the full problem but has ?‐independent computational cost. We also define an asymptotic problem on the graph of streamline separatrices and show that solution of the water pipe problem itself may be approximated by an asymptotic, ?‐independent problem on this graph. Finally, we show that the Dirichlet‐to‐Neumann map of the water pipe problem approximates the Dirichlet‐to‐Neumann map of the separatrix problem with an error independent of the flow outside the boundary layers. © 2004 Wiley Periodicals, Inc.     

Fast Thermoviscous Convection

This paper studies the asymptotic structure of convection in an infinite Prandtl number fluid with strongly temperature-dependent viscosity, in the limit where the dimensionless activation energy 1/ε is large, and the Rayleigh number R, defined (essentially) with the basal viscosity and the prescribed temperature drop, is also large. We find that the Nusselt number N is given by N~CεR^1/5, where C depends on the aspect ratio a. The relative error in this result is O(R^?1/10ε^?1/4, ε^1/2, R^?2/5ε^?2, R^?2/20ε^?1/24), so that we cannot hope to find accurate confirmation of this result at moderate Rayleigh numbers, though it should serve as a useful indicator of the relative importance of R and ε. For the above result to be valid, we require R ? 1/ε⁵ ?1. More important, however, is the asymptotic structure of the flow: there is a cold (hence rigid) lid with sloping base, beneath which a rapid, essentially isoviscous, convection takes place. This convection is driven by plumes at the sides, which generate vorticity due to thermal buoyancy, as in the constant viscosity case (Roberts, 1979). However, the slope of the lid base is sufficient to cause a large shear stress to be generated in the thermal boundary layer which joins the lid to the isoviscous region underneath (though a large velocity is not generated); consequently, the layer does not “see” the shear stress exerted by the interior flow (at leading order), and therefore the thermal boundary layer structure is totally self-determined: it even has a similarity structure (as a consequence). This fact makes it easy to analyse the problem, since the boundary layer uncouples from the rest of the flow. In addition, we find an alternative scaling (in which the lid base is “almost” flat), but it seems that the resulting boundary layer equations have no solution, though this is certainly open to debate: the results quoted above are not for this case. When a free slip boundary condition is applied at the top surface, one finds that there exists a thin “skin” at the top of the lid which is a stress boundary layer. The shear stress changes rapidly to zero, and there exists a huge longitudinal stress (compressive/tensile) in this skin. For earthlike parameters, this stress far exceeds the fracture strength of silicate rocks.     

Free and forced convection flow in a rotating channel bounded below by a permeable bed

The steady flow in a parallel plate channel rotating with an angular velocity Ω and bounded below by a permeable bed is analysed under the effect of buoyancy force. On the porous bed the boundary condition of Beavers and Joseph is applied and an exact solution of the governing equations is found. The solution in dimensionless form contains four parameters: The permeability parameterσ ², the Grashof numberG, the rotation parameterK ² and a dimensionless constantα. The effects of these parameters, specially,σ ², G andK ², on the slip velocities and velocity distributions are studied. For largeK ², there arise thin boundary layers on the walls of the channel.     

Combined forced and free convection heat transfer from a horizontal flat plate

The problem of combined forced and free convection heat transfer in boundary-layer flow past a heated or cooled, horizontal flat plate is studied in the limits of large and small Prandtl number using the technique of matched asymptotic expansions. It is found that the effects of natural convection are confined to small changes in the basic forced convection velocity and temperature fields providedGr/Re ^5/2 Pr ^2/3 1 forPr1, orGr/Re ^5/2 Pr ^1/2 1 forPr1. The governing similarity equations at the first two orders in these small parameters are solved and the results compared with the numerical solutions of Sparrow and Minkowycz (1962).

---

Zusammenfassung Das Problem der kombinierten erzwungenen und freien Konvektion in der Grenzschichtströmung entlang einer geheizten oder gekühlten horizontalen ebenen Platte wurde untersucht im Grenzfall von kleinen und grossen Prandtl-Zahlen, mit Hilfe der angepassten asymptotischen Entwicklungen. Der Effekt der natürlichen Konvektion ergab sich als kleine Korrektur der Geschwindigkeits-und Temperatur-Felder, fallsGr/Re ^5/2 Pr ^2/3 1 fürPr1 oderGr/Re ^5/2 Pr ^1/2 1 fürPr1. Die Ähnlichkeitslösungen der beiden ersten Ordnungen in diesen Parametern wurden bestimmt, und die Resultate wurden mit der numerischen Lösung von Sparrow und Minkowycz (1962) verglichen.

     

DNS of a Turbulent Channel Flow with Streamwise Rotation - Study of the Reverse Effect of the Cross Flow

Modeling of rotating turbulent flows is a major issue in engineering applications. In this work a turbulent channel flow rotating about the streamwise direction is presented. The theory is based on the investigations of [3] employing Lie group analysis. It was found that a cross flow in spanwise direction is induced. A series of direct numerical simulations (DNS) has been conducted for both different rotation rates and different Reynolds numbers to validate the cross flow. In addition some new interesting effects were observed. The averaged profile ū₃ of the cross flow is formed like a ‘S’ that means it exhibits a triple zero-crossing which denotes regions of reverse flow. Alaso a reverse effect is seen which means that for small rotation rates up to Ro=10 the spanwise mean velocity profiles increase and at rotation number Ro=14 this effect appears to reverse. Both effects were observed at two different Reynolds numbers Re = 180 and Re = 270. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent drag reduction at moderate Reynolds numbers via spanwise velocity waves

The question whether large turbulent drag reduction can be achieved at the high values of Re typical of applications is addressed. Answering such question, either by experiments or DNS, is obviously challenging. For DNS, the problem lies in the tremendous increase of the computational cost with Re, that has to be appreciated in view of the need of carrying out an entire parametric study at every Re, owing to the unknown location of the optimal forcing parameters. In this paper we limit ourselves to considering an open-loop technique based on spanwise forcing, the streamwise-traveling waves introduced by [1], and explore via Direct Numerical Simulations (DNS) how the drag reduction varies when the friction Reynolds number is increased from Re_τ = 200 to Re_τ = 2000. To achieve high Re while keeping the computational cost affordable, computational domains of reduced size are employed. We adopted special care to interpret results that are indeed still box-size dependent, as well as strategies to compute the random errors and give the results an error bar. Our results indicate that still R = 0.29 can be obtained at Re_τ = 2000 in the partial region of the parameter space studied. The maximum R is found to decrease as R ˜ Re_τ^−0.22 in the Reynolds range investigated. As most important outcome, we find that the sensitivity of R to Re becomes smaller when far from the low-Re optimum parameters: in this region, we suggest R ˜ Re_τ^−0.08. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Hydromagnetische Stabilität laminarer Grenzschicht an einer längsangeströmten ebenen Platte

Summary The purpose of the paper is to consider the stability for wavelike disturbances in the steady, twodimensional, laminar boundary layer of a magnetic field, which is applied uniformly normal to the flat plate. The results show that the critical Reynolds number (R _c ^* ) increases remarkably with the characteristic parameter (). The increase of the critical Reynolds number depends not only on the shape parameter of the velocity distribution in the boundary layer but also on the peculiarity of the velocity profile. It is also found that the boundary layer holds itself laminar all over the flat plate, when the magnetic parameterN is greater than 1.25×10^–7, then a reduction of the skin-frictin drag might be expeced.     

Direct numerical simulation of transition in MHD duct flow

Transition in the flow of electrically conducting fluid in a square duct with insulating walls is studied by direct numerical simulations. A uniform magnetic field is applied in the transverse direction. Moderate values of the Reynolds (Re = 5000 ) and Hartmann (Ha = 0 … 30 ) numbers are considered that correspond to the classical Hartmann & Lazarus [1] experiments. It is shown that the laminarization begins in the Hartmann layers, whereas the sidewall layers remain turbulent. Complete re-laminarization occurs in the range of R = Re/Ha ≈︁ 220 , which is in agreement with the H. & L. experiments. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of 1‐D incompressible laminar micro‐flow convection behaviour with temperature dependent fluid properties

This paper investigates the qualitative behaviour of single‐phase laminar convection for microchannels and conventionallysized channels formed between two parallel plates, captured by a numerical simulation on water flow. The convection parameters are obtained by separate numerical calculations on a series of parallel plates at constant temperatures. The pairs of parallel plates are maintained at progressively greater temperatures, to simulate the condition of increasing fluid temperature in a channel. The governing one‐dimensional (1‐D) momentum and energy equations are formulated to incorporate the dependence on temperature of both fluid viscosity (μ) and thermal conductivity (k). The qualitative behaviour of Nusselt number (Nu) decreasing with increasing Reynolds number (Re), exhibited by reported experimental data in literature, is simulated. Results show that it is practically dif_cult to observe this behaviour in the conventionally‐sized channels, but the effect easily surfaces in microchannels for practical lengths of flow and allowable high heat flux (q^″_W). (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Longitudinal and transversal two-point correlations in a turbulent flow

Based on measurements with hot-wire anemometers in a fully developed turbulent flow in the wake of a cylinder we investigate different aspects of longitudinal and transversal two-point correlations. In particular we use Taylor's hypothesis as well as simultaneous measurements with two X-hot-wire probes, and investigate how these correlations change with the relative orientation of their separation vector with respect to the mean flow direction. For the wake at a distance of eighty cylinder diameters at a Reynolds number Re = 28 000 (Taylor-based Reynolds number Re_λ = 290), we find evidence of strong anisotropy. We analyze the data using increment statistics, in particular increment distributions and structure functions. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic description of incipient separation bubble bursting

The appearance of short laminar separation bubbles in high Reynolds number (Re) wall bounded flows due to appropriate adverse pressure gradient conditions is usually associated with minor effects on global flow properties (e.g. lift force). However, localized reverse flow regions are known to react very sensitively to perturbations and in further consequence may trigger the laminar-turbulent transition process or even cause global separation. The present investigation of marginally separated boundary layer flows is based on an asymptotic approach Re → ∞. Special emphasis is placed on solutions of the corresponding model equations which blow up within finite time indicating the ejection of a vortical structure and the emergence of shorter spatio-temporal scales reminiscent of the early transition scenario (‘ bubble bursting’ ). Within the framework of marginal separation theory, an alternative adjoint operator method is used to formulate evolution equations governing the viscous-inviscid interaction process in leading and higher order correction required for the study of later stages of the flow development. Their blow up structure specifies the initial condition of and the match to the subsequent triple deck stage. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Boundary value problems with eigenvalue depending boundary conditions

We investigate some classes of eigenvalue dependent boundary value problems of the form where A ? A⁺ is a symmetric operator or relation in a Krein space K, τ is a matrix function and Γ₀, Γ₁ are abstract boundary mappings. It is assumed that A admits a self‐adjoint extension in K which locally has the same spectral properties as a definitizable relation, and that τ is a matrix function which locally can be represented with the resolvent of a self‐adjoint definitizable relation. The strict part of τ is realized as the Weyl function of a symmetric operator T in a Krein space H, a self‐adjoint extension à of A × T in K × H with the property that the compressed resolvent P_K (à – λ)^–1|_K k yields the unique solution of the boundary value problem is constructed, and the local spectral properties of this so‐called linearization à are studied. The general results are applied to indefinite Sturm–Liouville operators with eigenvalue dependent boundary conditions (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new approach to the nonlinear stability of viscous flow in a coplanar magnetic field

We present a new Lyapunov function for laminar flow, in the x‐direction, between two parallel planes in the presence of a coplanar magnetic field for three‐dimensional perturbations with stress‐free boundary planes that provides conditional nonlinear stability for all Reynolds numbers(R_e) and magnetic Reynolds numbers(R_m) below π²/2M. Compared with previous results on the nonlinear stability of this problem, the radius of stability ball and the energy decay rate obtained in this paper are independent of the magnetic field. Copyright © 2016 John Wiley & Sons, Ltd.     

On fractional K‐factors of random graphs

Let K be a graph on r vertices and let G = (V,E) be another graph on ∣V ∣ = n vertices. Denote the set of all copies of K in G by 𝒦. A non‐negative real‐valued function f : 𝒦→ ℝ₊ is called a fractional K‐factor if ∑ _{K:v∈K∈𝒦}f(K) ≤ 1 for every v ∈ V and ∑ _K∈𝒦f(K) = n/r. For a non‐empty graph K let d(K) = e(K)/v(K) and d⁽¹⁾(K) = e(K)/(v(K) ‐ 1). We say that K is strictly K₁‐balanced if for every proper subgraph K^′ ⊊ K, d⁽¹⁾(K^′) < d⁽¹⁾(K). We say that K is imbalanced if it has a subgraph K^′ such that d(K^′) > d(K). Considering a random graph process on n vertices, we show that if K is strictly K₁‐balanced, then with probability tending to 1 as n →∞, at the first moment τ₀ when every vertex is covered by a copy of K, the graph has a fractional K‐factor. This result is the best possible. As a consequence, if K is K₁‐balanced, we derive the threshold probability function for a random graph to have a fractional K‐factor. On the other hand, we show that if K is an imbalanced graph, then for asymptotically almost every graph process there is a gap between τ₀ and the appearance of a fractional K‐factor. We also introduce and apply a criteria for perfect fractional matchings in hypergraphs in terms of expansion properties. © 2006 Wiley Periodicals, Inc. Random Struct. Alg., 2007