Numerical Study of Turbulent Channel Flow over Surfaces with Variable Spanwise Heterogeneities: Topographically-driven Secondary Flows Affect Outer-layer Similarity of Turbulent Length Scales

Wall-bounded turbulent flows over surfaces with spanwise heterogeneous surface roughness – that is, spanwise-adjacent patches of relatively high and low roughness – exhibit mean flow phenomena entirely different to what would otherwise exist in the absence of spanwise heterogeneity. In the outer layer, mean counter-rotating rolls occupy the depth of the flow, and are positioned such that “upwelling” and “downwelling” occurs above the low and high roughness, respectively. It has been comprehensively shown that these secondary flows are Prandtl’s secondary flow of the second kind (Anderson et al., J. Fluid Mech. 768, 316–347 2015). This behaviour indicates that spanwise spacing, s _y , between adjacent patches of high and low roughness is, itself, a problem parameter; in this study, we have systematically assessed how s _y affects turbulence structure in high Reynolds number channel flows via two-point correlations. “High roughness” is imposed with streamwise-aligned pyramid elements with height, h, selected to be ≈ 5% of the channel half height, H. For \(s_{y}/H \gtrsim 1\), we find that the aforementioned domain-scale mean circulations exist and the surface may be regarded as a topography. For s _y /H ? 0.2, turbulence statistics show characteristics very similar to a homogeneous roughness and thus the surface may be regarded as a roughness. For 0.2 ? s _y /H ? 2, the spatial extent of the counter-rotating rolls is controlled by proximity to adjacent rows, and we define such surfaces as being intermediate. We refer to such surfaces as intermediate state.     

Growths of mechanical elasticity and electrical conductance of graphene nanoplatelet/poly(lactic acid) composites under strong electric field: correlation with time evolution of higher order structure of graphene nanoplatelets

For a composite of poly(lactic acid) containing graphene nanoplatelets (GNPs) at a low concentration (0.34 vol%), this study examined growths of mechanical and electrical properties under an alternating current (AC) electric field, focusing on field-induced GNP structures governing those properties. The composite was subjected to the AC field (60 Hz) of various intensities E for various times t _E at 190 °C. A fraction of randomly oriented GNPs was aligned by the field and then connected into columns, as suggested from optical microscopy. This structural evolution led to qualitatively similar growths of low-frequency storage modulus and static electrical conductivity. The key quantity for understanding this growth was a time t _E* for occurrence of short circuit that detected formation of GNP columns conductively bridging the electrodes. The growths of both modulus and conductivity for various E were summarized as functions of a reduced variable, t _E/t _E*, confirming the growths commonly reflected the evolution of the GNP columns. However, the modulus grew fast and leveled off by t _E/t _E* ~ 1, whereas the conductivity kept growing gradually even at t _E/t _E* > 1. This difference was discussed in relation to the matrix chains and leftover GNPs out the column.     

Global Attractor for a Nonlinear Oscillator Coupled to the Klein–Gordon Field

The long-time asymptotics is analyzed for all finite energy solutions to a model\(\mathbf{U}(1)\)-invariant nonlinear Klein–Gordon equation in one dimension, with the nonlinearity concentrated at a single point: each finite energy solution converges as t→ ± ∞ to the set of all “nonlinear eigenfunctions” of the form ψ(x)e^{?iω t}. The global attraction is caused by the nonlinear energy transfer from lower harmonics to the continuous spectrum and subsequent dispersive radiation.We justify this mechanism by the following novel strategy based on inflation of spectrum by the nonlinearity. We show that any omega-limit trajectory has the time spectrum in the spectral gap [ ? m,m] and satisfies the original equation. This equation implies the key spectral inclusion for spectrum of the nonlinear term. Then the application of the Titchmarsh convolution theorem reduces the spectrum of each omega-limit trajectory to a single harmonic \(\omega\in[-m,m]\).The research is inspired by Bohr’s postulate on quantum transitions and Schrödinger’s identification of the quantum stationary states to the nonlinear eigenfunctions of the coupled\(\mathbf{U}(1)\)-invariant Maxwell–Schrödinger and Maxwell–Dirac equations.     

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]}\).     

Does Density Ratio Significantly Affect Turbulent Flame Speed?

In order to experimentally study whether or not the density ratio σ substantially affects flame displacement speed at low and moderate turbulent intensities, two stoichiometric methane/oxygen/nitrogen mixtures characterized by the same laminar flame speed S _L = 0.36 m/s, but substantially different σ were designed using (i) preheating from T _u = 298 to 423 K in order to increase S _L , but to decrease σ, and (ii) dilution with nitrogen in order to further decrease σ and to reduce S _L back to the initial value. As a result, the density ratio was reduced from 7.52 to 4.95. In both reference and preheated/diluted cases, direct images of statistically spherical laminar and turbulent flames that expanded after spark ignition in the center of a large 3D cruciform burner were recorded and processed in order to evaluate the mean flame radius \(\bar {R}_{f}\left (t \right )\) and flame displacement speed \(S_{t}=\sigma ^{-1}{d\bar {R}_{f}} \left / \right . {dt}\) with respect to unburned gas. The use of two counter-rotating fans and perforated plates for near-isotropic turbulence generation allowed us to vary the rms turbulent velocity \(u^{\prime }\) by changing the fan frequency. In this study, \(u^{\prime }\) was varied from 0.14 to 1.39 m/s. For each set of initial conditions (two different mixture compositions, two different temperatures T _u , and six different \(u^{\prime })\), five (respectively, three) statistically equivalent runs were performed in turbulent (respectively, laminar) environment. The obtained experimental data do not show any significant effect of the density ratio on S _t . Moreover, the flame displacement speeds measured at u′/S _L = 0.4 are close to the laminar flame speeds in all investigated cases. These results imply, in particular, a minor effect of the density ratio on flame displacement speed in spark ignition engines and support simulations of the engine combustion using models that (i) do not allow for effects of the density ratio on S _t and (ii) have been validated against experimental data obtained under the room conditions, i.e. at higher σ.     

Decay of Almost Periodic Solutions¶of Conservation Laws

We consider the asymptotic behavior of solutions of systems of inviscid or viscous conservation laws in one or several space variables, which are almost periodic in the space variables in a generalized sense introduced by Stepanoff and Wiener, which extends the original one of H. Bohr. We prove that if u(x,t) is such a solution whose inclusion intervals at time t, with respect to ?>0, satisfy l _epsiv;(t)/t→0 as t→∞, and such that the scaling sequence u ^T (x,t)=u(T x,T t) is pre-compact as t→∞ in L _loc ¹(? ^{d +1} ₊, then u(x,t) decays to its mean value \(\), which is independent of t, as t→∞. The decay considered here is in L ¹ _loc of the variable ξ≡x/t, which implies, as we show, that \(\) as t→∞, where M _x denotes taking the mean value with respect to x. In many cases we show that, if the initial data are almost periodic in the generalized sense, then so also are the solutions. We also show, in these cases, how to reduce the condition on the growth of the inclusion intervals l _?(t) with t, as t→∞, for fixed ? > 0, to a condition on the growth of l _?(0) with ?, as ?→ 0, which amounts to imposing restrictions only on the initial data. We show with a simple example the existence of almost periodic (non-periodic) functions whose inclusion intervals satisfy any prescribed growth condition as ?→ 0. The applications given here include inviscid and viscous scalar conservation laws in several space variables, some inviscid systems in chromatography and isentropic gas dynamics, as well as many viscous 2 × 2 systems such as those of nonlinear elasticity and Eulerian isentropic gas dynamics, with artificial viscosity, among others. In the case of the inviscid scalar equations and chromatography systems, the class of initial data for which decay results are proved includes, in particular, the L ^∞ generalized limit periodic functions. Our procedures can be easily adapted to provide similar results for semilinear and kinetic relaxations of systems of conservation laws.     

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.     

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.     

Two-Scale <Emphasis Type="Italic">Γ</Emphasis>-Convergence of Integral Functionals and its Application to Homogenisation of Nonlinear High-Contrast Periodic Composites

An analytical framework is developed for passing to the homogenisation limit in (not necessarily convex) variational problems for composites whose material properties oscillate with a small period ε and that exhibit high contrast of order \({\varepsilon^{-1}}\) between the constitutive, “stress-strain”, response on different parts of the period cell. The approach of this article is based on the concept of “two-scale Γ-convergence”, which is a kind of “hybrid” of the classical Γ-convergence (De Giorgi and Franzoni in Atti Accad Naz Lincei Rend Cl Sci Fis Mat Natur (8)58:842–850, 1975) and the more recent two-scale convergence (Nguetseng in SIAM J Math Anal 20:608–623, 1989). The present study focuses on a basic high-contrast model, where “soft” inclusions are embedded in a “stiff” matrix. It is shown that the standard Γ-convergence in the L ^p -space fails to yield the correct limit problem as \({\varepsilon \to 0,}\) due to the underlying lack of L ^p -compactness for minimising sequences. Using an appropriate two-scale compactness statement as an alternative starting point, the two-scale Γ-limit of the original family of functionals is determined via a combination of techniques from classical homogenisation, the theory of quasiconvex functions and multiscale analysis. The related result can be thought of as a “non-classical” two-scale extension of the well-known theorem by Müller (Arch Rational Mech Anal 99:189–212, 1987).     

Potentials for tangent tensor fields on spheroids

In three-dimensional Euclidean space let S be a closed simply connected, smooth surface (spheroid). Let \(\hat n\) be the outward unit normal to S, ▽ _S the surface gradient on S, I _S the metric tensor on S, g_ij the four covariant components of I _S (i,j = 1, 2), h _ij the four covariant components of -\(\hat n\)xI _S , and D _i covariant differentiation on S. It is well known that for any tangent vector field u on S there exist scalars ? and ψ on S, unique to within additive constants, such that \(u = \nabla _s \varphi - \hat n \times \nabla _s \psi \); the covariant components of u are \(u_i = D_i \varphi + h_i^j D_j \psi \). This theorem is very useful in the study of vector fields in spherical coordinates. The present paper gives an analogous theorem for real second-order tangent tensor fields F on S: for any such F there exist scalar fields H, L, M, N such that the covariant components of F are

$$F_{ij} = H h{}_{ij} + Lg_{ij} + E_{ij} (M,N),$$     

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.     

Turbulent Scalar Mixing in a Skewed Jet in Crossflow: Experiments and Modeling

Turbulent mixing of an inclined, skewed jet injected into a crossflow is investigated using MRI-based experiments and a high-fidelity LES of the same configuration. The MRI technique provides three-dimensional fields of mean velocity and mean jet concentration. The 30° skew of the jet relative to the crossflow produces a single dominant vortex which introduces spanwise asymmetries to the velocity and concentration fields. The turbulent scalar transport of the skewed jet is investigated in further detail using the LES, which is validated against the experimental measurements. Mixing is found to be highly anisotropic throughout the jet region. Isotropic turbulent diffusivity and viscosity are used to calculate an optimal value of the turbulent Schmidt number, which varies widely over the jet region and lies mostly outside of the typically accepted range 0.7 ≤ S c _t ≤ 0.9. Finally, three common scalar flux models of increasing complexity are evaluated based on their ability to capture the anisotropy and predict the scalar concentration field of the present configuration. The higher order models are shown to better represent the turbulent scalar flux vector, leading to more accurate calculations of the concentration field. While more complex models are better able to capture the turbulent mixing, optimization of model constants is shown to significantly affect the results.     

Families of rational solutions of the <Emphasis Type="Italic">y</Emphasis>-nonlocal Davey–Stewartson II equation

The y-nonlocal Davey–Stewartson II equation is an extension of the usual DS II equation involving a partially parity-time-symmetric potential only with respect to the spatial variable y. By using the Hirota bilinear method, families of n-order rational solutions are obtained, which include lumps in the (x, y)-plane and the (y, t)-plane, growing-and-decaying line waves in the (x, t)-plane, and hybrid solutions of interacting line rogue waves and lumps in the (x, y)-plane.     

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.     

Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

Separation zones in the flow near a small-angle convex corner

On the basis of an asymptotic analysis of the Navier-Stokes system of equations for large Reynolds numbers (Re → ∞), the plane incompressible fluid flow near a surface having a convex corner with a small angle 2θ* is investigated. It is shown that for θ* = O(Re^?1/4), in addition to the known solution that describes a separated flow completely localized in a thin “viscous” sublayer of the interaction region near the corner point, another solution corresponding to a flow with a developed separation zone is possible. For θ ₀ = Re^1/4 θ* = O(1), the longitudinal dimension of this zone varies from finite values up to values of the order of Re^?3/8. The nonuniqueness of the solution is established on a certain range of variation of the parameter θ ₀. The dependence of the drag coefficient on the angle θ* is found.     

Mach waves produced in the supersonic jet mixing layer by shock/vortex interaction

The noise emission of free jets has been extensively investigated for many decades. At subsonic jet velocities, coherent structures of the mixing layer move at subsonic speed and emit sound waves. Free jets blowing at supersonic speeds, however, can emit weak shock waves, called Mach waves. At supersonic speeds, two cases must be distinguished: the structures move either subsonically or supersonically relative to the inside and/or outside speed of sound. In the case of supersonic movement, the Mach waves exist inside as well as outside the jet. At subsonic speeds, no Mach waves appear. Although numerous theories have been established to find the origin of the Mach waves, to the authors’ best knowledge, the mechanism of the Mach wave formation has not yet been clearly explained. Recently another theory of Mach waves in supersonic jets was developed, as described herein, which outlines the causes for the Mach wave production and stability as well as their dynamics. The theory’s principle is that the Mach waves are initiated by vortices which move downstream at three speeds w, \({w}'\) and \({w}''\) inside of the mixing layer. These three types of vortices and Mach waves are described in a comprehensive manner by the theory and are called the “w-, \({w}'\)- and \({w}''\)-vortices” and “w-, \({w}'\)- and \({w}''\)-Mach waves,” respectively.     

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.     

Power-law creep and residual stresses in a carbopol gel

We report on the interplay between creep and residual stresses in a carbopol microgel. When a constant shear stress σ is applied below the yield stress σ _y, the strain is shown to increase as a power law of time, γ(t) = γ ₀ + (t/τ) ^α , with an exponent α = 0.39 ± 0.04 that is strongly reminiscent of Andrade creep in hard solids. For applied shear stresses lower than some typical value σ _c ? 0.2σ _y, the microgel experiences a more complex, anomalous creep behavior, characterized by an initial decrease of the strain, that we attribute to the existence of residual stresses of the order of σ _c that persist after a rest time under a zero shear rate following preshear. The influence of gel concentration on creep and residual stresses are investigated as well as possible aging effects. We discuss our results in light of previous works on colloidal glasses and other soft glassy systems.     

Fine Level Set Structure of Flat Isometric Immersions

A result by Pogorelov asserts that C ¹ isometric immersions u of a bounded domain \({S \subset \mathbb R^2}\) into \({\mathbb {R}^3}\) whose normal takes values in a set of zero area enjoy the following regularity property: the gradient \({f := \nabla u}\) is ‘developable’ in the sense that the nondegenerate level sets of f consist of straight line segments intersecting the boundary of S at both endpoints. Motivated by applications in nonlinear elasticity, we study the level set structure of such f when S is an arbitrary bounded Lipschitz domain. We show that f can be approximated by uniformly bounded maps with a simplified level set structure. We also show that the domain S can be decomposed (up to a controlled remainder) into finitely many subdomains, each of which admits a global line of curvature parametrization.