Determination of boundary layers in the plane problem for three-layer strips. Part 2

In the first part of this paper, we considered the exact statement of the plane elasticity problem in displacements for strips made of various materials (problem A, an isotropic material; problem B, an orthotropic material with 2G ₁₂ < √E ₁ E ₂; problem C, an orthotropic material with 2G ₁₂ > √E ₁ E ₂). Further, we stated and solved the boundary layer problem (the problem on a solution decaying away from the boundary) for a sandwich strip of regular structure consisting of isotropic layers (problem AA). In the present paper, we use the solution of the plane problem to consider the problem for sandwich strips of regular structure with isotropic face layers and orthotropic filler (problem AB).     

Revisiting History Effects in Adverse-Pressure-Gradient Turbulent Boundary Layers

The goal of this study is to present a first step towards establishing criteria aimed at assessing whether a particular adverse-pressure-gradient (APG) turbulent boundary layer (TBL) can be considered well-behaved, i.e., whether it is independent of the inflow conditions and is exempt of numerical or experimental artifacts. To this end, we analyzed several high-quality datasets, including in-house numerical databases of APG TBLs developing over flat-plates and the suction side of a wing section, and five studies available in the literature. Due to the impact of the flow history on the particular state of the boundary layer, we developed three criteria of convergence to well-behaved conditions, to be used depending on the particular case under study. (i) In the first criterion, we develop empirical correlations defining the R e _?? -evolution of the skin-friction coefficient and the shape factor in APG TBLs with constant values of the Clauser pressure-gradient parameter β = 1 and 2 (note that β = δ ^?/τ _w dP _e /dx, where δ ^? is the displacement thickness, τ _w the wall-shear stress and dP _e /dx the streamwise pressure gradient). (ii) In the second one, we propose a predictive method to obtain the skin-friction curve corresponding to an APG TBL subjected to any streamwise evolution of β, based only on data from zero-pressure-gradient TBLs. (iii) The third method relies on the diagnostic-plot concept modified with the shape factor, which scales APG TBLs subjected to a wide range of pressure-gradient conditions. These three criteria allow to ensure the correct flow development of a particular TBL, and thus to separate history and pressure-gradient effects in the analysis.     

The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

Improvement for expansion of parabolized stability equation method in boundary layer stability analysis

An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.     

A method of solving the diffusion problem for a vortex layer in a viscoplastic half-plane

A method is proposed to reduce the classical formulation of the problem to a system of two functional equations whose solution can be found numerically. A number of assertions that characterize the behavior of a rigid zone are proved. In particular, the lower estimate h ⁰(t) = 2b√t for the boundary motion is obtained; an explicit expression for b is given as a boundary stress function.     

Limit-analysis of a circular cylinder obeying the Green plasticity criterion and loaded in combined tension and torsion

This short paper presents the limit-analysis of a cylinder with circular basis, made of an ideal-plastic material obeying Green’s yield criterion and subjected to combined tension and torsion. The exact solution of the problem is provided in the form of a statically and plastically admissible stress field and a kinematically admissible velocity field, associated via the normality rule. The overall yield locus, that is the set of pairs [tension force, torsion torque] for which unrestrained plastic flow occurs, is expressed first in parametric form, then explicitly upon elimination of the parameter involved. The explicit expression of this yield locus also entails that of the overall flow rule via the overall normality property. The impact of these results is two-fold. First, they provide a fresh example of a solution to a limit-analysis problem exceptionally combining three generally mutually exclusive features: be non-trivial, exact and explicit. Second, they provide a way of using simple experiments of combined tension and torsion of cylinders to determine the parameter characterizing the influence of the mean stress in Green’s criterion.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

Inverse problem of fracture mechanics for a compound cylinder

We consider the plane problem of fracture mechanics for a compound cylinder. We assume that the sleeve (internal cylinder) is reinforced with negative allowance by the external cylinder and that near the sleeve surface there are N arbitrarily located rectilinear cracks of length 2l _k (k = 1, 2, ..., N). A minimax criterion is used to determine the negative allowance in the junction theoretically minimizing the fracture parameters (the stress intensity coefficients) of the compound cylinder. A simplified method for minimizing the fracture parameters of the compound cylinder is considered separately.     

Inviscid Models Generalizing the Two-dimensional Euler and the Surface Quasi-geostrophic Equations

Any classical solution of the two-dimensional incompressible Euler equation is global in time. However, it remains an outstanding open problem whether classical solutions of the surface quasi-geostrophic (SQG) equation preserve their regularity for all time. This paper studies solutions of a family of active scalar equations in which each component u _j of the velocity field u is determined by the scalar θ through \({u_j =\mathcal{R}\Lambda^{-1}P(\Lambda) \theta}\) , where \({\mathcal{R}}\) is a Riesz transform and Λ = (?Δ)^1/2. The two-dimensional Euler vorticity equation corresponds to the special case P(Λ) = I while the SQG equation corresponds to the case P(Λ) = Λ. We develop tools to bound \({\|\nabla u||_{L^\infty}}\) for a general class of operators P and establish the global regularity for the Loglog-Euler equation for which P(Λ) = (log(I + log(I ? Δ))) ^γ with 0 ≦ γ ≦ 1. In addition, a regularity criterion for the model corresponding to P(Λ) = Λ ^β with 0 ≦ β ≦ 1 is also obtained.     

Bistable Boundary Reactions in Two Dimensions

In a bounded domain \({\Omega \subset \mathbb R^2}\) with smooth boundary we consider the problem

$\Delta u = 0 \quad {\rm{in }}\, \Omega, \qquad \frac{\partial u}{\partial \nu} = \frac1\varepsilon f(u) \quad {\rm{on }}\,\partial\Omega,$

where ν is the unit normal exterior vector, ε > 0 is a small parameter and f is a bistable nonlinearity such as f(u) = sin(π u) or f(u) = (1 ? u ²)u. We construct solutions that develop multiple transitions from ?1 to 1 and vice-versa along a connected component of the boundary ?Ω. We also construct an explicit solution when Ω is a disk and f(u) = sin(π u).

     

Boundary Stabilization of Stokes System in Exterior Domains

The problem of stabilizing a solution to the 2D Stokes system defined in the exterior of the bounded domain with smooth boundary is investigated, i.e. for a given initial velocity field and prescribed positive number k > 0 one has to construct a control function defined on the boundary such that the solution stabilizes to zero at the rate of 1/t ^k .     

Consequences of the Translation Invariance on the Darcy Free Convection Flow Past a Vertical Surface

It is shown that the governing equation for the stream function of the Darcy free convection boundary layer flows past a vertical surface is invariant under arbitrary translations of the transverse coordinate y. The consequences of this basic symmetry property on the solutions corresponding to a prescribed surface temperature distribution T _w (x) are investigated. It is found that starting with a “primary solution” which describes the temperature boundary layer on an impermeable surface, infinitely many “translated solutions” can be generated which form a continuous group, the “translation group” of the given primary solution. The elements of this group describe free convection boundary layer flows from permeable counterparts of the original surface with a transformed temperature distribution \({\tilde {T}_w \left( x \right)}\), when simultaneously a suitable lateral suction/injection of the fluid is applied. It turns out in this way that several exact solutions discovered during the latter few decades are in fact not basically new solutions, but translated counterparts of some formerly reported primary solutions. A few specific examples are discussed in detail.     

New Periodic Solutions for Newtonian n-Body Problems with Dihedral Group Symmetry and Topological Constraints

In this paper, we prove the existence of a family of new non-collision periodic solutions for the classical Newtonian n-body problems. In our assumption, the \({n=2l \geqq 4}\) particles are invariant under the dihedral rotation group D_l in \({\mathbb{R}^3}\) such that, at each instant, the n particles form two twisted l-regular polygons. Our approach is the variational minimizing method and we show that the minimizers are collision-free by level estimates and local deformations.     

Perspectives on the Phenomenology and Modeling of Boundary Layer Transition

The characteristics of the coherent structures in a strongly decelerated large-velocity-defect boundary layer are analysed by direct numerical simulation. The simulated boundary layer starts as a zero-pressure-gradient boundary layer, decelerates under a strong adverse pressure gradient, and separates near the end of the domain, in the form of a very thin separation bubble. The Reynolds number at separation is R e _?? =3912 and the shape factor H=3.43. The three-dimensional spatial correlations of (u, u) and (u, v) are investigated and compared to those of a zero-pressure-gradient boundary layer and another strongly decelerated boundary layer. These velocity pairs lose coherence in the streamwise and spanwise directions as the velocity defect increases. In the outer region, the shape of the correlations suggest that large-scale u structures are less streamwise elongated and more inclined with respect to the wall in large-defect boundary layers. The three-dimensional properties of sweeps and ejections are characterized for the first time in both the zero-pressure-gradient and adverse-pressure-gradient boundary layers, following the method of Lozano-Durán et al. (J. Fluid Mech. 694, 100–130, [2012]). Although longer sweeps and ejections are found in the zero-pressure-gradient boundary layer, with ejections reaching streamwise lengths of 5 boundary layer thicknesses, the sweeps and ejections tend to be bigger in the adverse-pressure-gradient boundary layer. Moreover, small near-wall sweeps and ejections are much less numerous in the large-defect boundary layer. Large sweeps and ejections that reach the wall region (wall-attached) are also less numerous, less streamwise elongated and they occupy less space than in the zero-pressure-gradient boundary layer.     

Pressure-Gradient Turbulent Boundary Layers Developing Around a Wing Section

A direct numerical simulation database of the flow around a NACA4412 wing section at R e _c = 400,000 and 5^° angle of attack (Hosseini et al. Int. J. Heat Fluid Flow 61, 117–128, 2016), obtained with the spectral-element code Nek5000, is analyzed. The Clauser pressure-gradient parameter β ranges from ? 0 and 85 on the suction side, and from 0 to ? 0.25 on the pressure side of the wing. The maximum R e _?? and R e _τ values are around 2,800 and 373 on the suction side, respectively, whereas on the pressure side these values are 818 and 346. Comparisons between the suction side with zero-pressure-gradient turbulent boundary layer data show larger values of the shape factor and a lower skin friction, both connected with the fact that the adverse pressure gradient present on the suction side of the wing increases the wall-normal convection. The adverse-pressure-gradient boundary layer also exhibits a more prominent wake region, the development of an outer peak in the Reynolds-stress tensor components, and increased production and dissipation across the boundary layer. All these effects are connected with the fact that the large-scale motions of the flow become relatively more intense due to the adverse pressure gradient, as apparent from spanwise premultiplied power-spectral density maps. The emergence of an outer spectral peak is observed at β values of around 4 for λ _z ? 0.65δ ₉₉, closer to the wall than the spectral outer peak observed in zero-pressure-gradient turbulent boundary layers at higher R e _?? . The effect of the slight favorable pressure gradient present on the pressure side of the wing is opposite the one of the adverse pressure gradient, leading to less energetic outer-layer structures.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

Effect of tube diameter and capillary number on platelet margination and near-wall dynamics

I investigate the effect of tube diameter D and red blood cell capillary number Ca (i.e. the ratio of viscous to elastic forces) on platelet margination in blood flow at ≈37 % tube haematocrit. The system is modelled as three-dimensional suspension of deformable red blood cells and nearly rigid platelets using a combination of the lattice-Boltzmann, immersed boundary and finite element methods. Results of simulations during the dynamics before the steady state has been reached show that a non-diffusive radial platelet transport facilitates margination. This non-diffusive effect is important near the edge of the cell-free layer, but only for Ca > 0.2, when red blood cells are tank-treading. I also show that platelet trapping in the cell-free layer is reversible for Ca ≤ 0.2. Margination is essentially independent of Ca only for the smallest investigated tube diameter (D = 10 μm). Once platelets have reached the cell-free layer, they tend to slide rather than tumble. The tumbling rate is essentially independent of Ca but increases with D. Strong confinement suppresses tumbling due to the relatively small cell-free layer thickness at ≈ 37 % tube haematocrit.     

Regularity of Solutions to the Navier-Stokes Equations for Compressible Barotropic Flows on a Polygon

The Navier-Stokes system for a steady-state barotropic nonlinear compressible viscous flow, with an inflow boundary condition, is studied on a polygon D. A unique existence for the solution of the system is established. It is shown that the lowest order corner singularity of the nonlinear system is the same as that of the Laplacian in suitable L ^q spaces. Let ω be the interior angle of a vertex P of D. If \(\) and \(\), then the velocity u is split into singular and regular parts near the vertex P. If α < 2 and \(\) or if α > 2 and 2 < q < ∞&;, it is shown that u∈ (H ^{2, q} (D))².