Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Multipulse Phases in k-Mixtures of Bose–Einstein Condensates

For the system

$-\Delta U_i+ U_i=U_i^3-\beta U_i\sum_{j\neq i}U_j^2,\quad i=1,\dots,k,$

(with k ≧ 3), we prove the existence for β large of positive radial solutions on \({\mathbb R^N}\) . We show that as β →  + ∞, the profile of each component U _i separates, in many pulses, from the others. Moreover, we can prescribe the location of such pulses in terms of the oscillations of the changing-sign solutions of the scalar equation  ? ΔW  +  W  =  W³. Within an Hartree–Fock approximation, this provides a theoretical indication of phase separation into many nodal domains for the k-mixtures of Bose–Einstein condensates.

     

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.     

Mechanical Properties of Hard W-C Coating on Steel Substrate Deduced from Nanoindentation and Finite Element Modeling

Mechanical properties of a hard and stiff W-C coating on steel substrate have been investigated using nanoindentation combined with finite element modeling (FEM) and extended FEM (XFEM). The significant pile-up observed around the indents in steel substrate caused an overestimation of hardness and indentation modulus. A simple geometrical model, considering the additional contact surfaces due to pile-up, has been proposed to reduce this overestimation. The presence of W-C coating suppressed the pile-up in the steel substrate and a transition to sink-in behavior occurred. The FEM simulations adequately reproduced the surface topography of the indents in the substrate and coating/substrate systems as well. The maximum principal stresses of the indented W-C/steel coated system were tensile; they were always located in the coating and evolved in 3 stages. Cohesive cracking occurred during loading in the sink-in zone (stage III) when the ultimate tensile strength (σ _max ) of the coating was reached. The obtained hardness (H _c ), indentation modulus (E _c ), yield stress (Y) and strength (σ _max ) of the W-C coating were H _c ? =?20 GPa, E _c ? =?250 GPa, Y?=?9.0 GPa and σ _max ? =?9.35 GPa, respectively. XFEM resulted in fracture energy of the W-C coating of G?=?38.1 J?·?m^-2 and fracture toughness of K _IC ? =?3.5 MPa?·?m^0.5.     

On Properties of the Generalized Wasserstein Distance

The Wasserstein distances W_p (p \({\geqq}\) 1), defined in terms of a solution to the Monge–Kantorovich problem, are known to be a useful tool to investigate transport equations. In particular, the Benamou–Brenier formula characterizes the square of the Wasserstein distance W₂ as the infimum of the kinetic energy, or action functional, of all vector fields transporting one measure to the other. Another important property of the Wasserstein distances is the Kantorovich–Rubinstein duality, stating the equality between the distance W₁(μ, ν) of two probability measures μ, ν and the supremum of the integrals in d(μ ?ν) of Lipschitz continuous functions with Lipschitz constant bounded by one. An intrinsic limitation of Wasserstein distances is the fact that they are defined only between measures having the same mass. To overcome such a limitation, we recently introduced the generalized Wasserstein distances \({W_p^{a,b}}\), defined in terms of both the classical Wasserstein distance W_p and the total variation (or L¹) distance, see (Piccoli and Rossi in Archive for Rational Mechanics and Analysis 211(1):335–358, 2014). Here p plays the same role as for the classic Wasserstein distance, while a and b are weights for the transport and the total variation term. In this paper we prove two important properties of the generalized Wasserstein distances: (1) a generalized Benamou–Brenier formula providing the equality between \({W_2^{a,b}}\) and the supremum of an action functional, which includes a transport term (kinetic energy) and a source term; (2) a duality à la Kantorovich–Rubinstein establishing the equality between \({W_1^{1,1}}\) and the flat metric.     

Measurement of local forcing on a turbulent boundary layer using PIV

An experimental study was carried out to investigate the effect of periodic blowing and suction on a turbulent boundary layer. Particle image velocimetry (PIV) was used to probe the characteristics of the flow. Local forcing was introduced to the boundary layer via a sinusoidally-oscillating jet issuing from a thin spanwise slot. Three forcing frequencies (f⁺=0.44, 0.66 and 0.88) with a fixed forcing amplitude (A⁺=0.6) were employed at Re _θ =690. The effect of three different forcing angles (α=60°, 90° and l20°) was investigated under a fixed forcing frequency (f⁺=0.088). The PIV results showed that the wall-region velocity decreases on imposition of the local forcing. Inspection of the phase-averaged velocity profiles revealed that spanwise large-scale vortices are generated downstream of the slot and persist farther downstream. The highest reduction in skin friction was achieved at the highest forcing frequency (f⁺=0.088) and a forcing angle of α=120°. The spatial fraction of the vortices was examined to analyze the skin friction reduction.     

Simultaneous impacts of MHD and variable wall temperature on transient mixed Casson nanofluid flow in the stagnation point of rotating sphere

A numerical analysis is provided to scrutinize time-dependent magnetohydrodynamics(MHD) free and forced convection of an electrically conducting non-Newtonian Casson nanofluid flow in the forward stagnation point region of an impulsively rotating sphere with variable wall temperature. A single-phase flow of nanofluid model is reflected with a number of experimental formulae for both effective viscosity and thermal conductivity of nanofluid. Exceedingly nonlinear governing partial differential equations(PDEs)subject to their compatible boundary conditions are mutated into a system of nonlinear ordinary differential equations(ODEs). The derived nonlinear system is solved numerically with implementation of an implicit finite difference procedure merging with a technique of quasi-linearization. The controlled parameter impacts are clarified by a parametric study of the entire flow regime. It is depicted that from all the exhibited nanoparticles,Cu possesses the best convection. The surface heat transfer and surface shear stresses in the x-and z-directions are boosted with maximizing the values of nanoparticle solid volume fraction ? and rotation λ. Besides, as both the surface temperature exponent n and the Casson parameter γ upgrade, an enhancement of the Nusselt number is given.     

Entire Solutions to Equivariant Elliptic Systems with Variational Structure

We consider the system Δu ? W _u (u) = 0, where \({u : \mathbb{R}^n \to \mathbb{R}^n}\) , for a class of potentials \({W : \mathbb{R}^n \to \mathbb{R}}\) that possess several global minima and are invariant under a general finite reflection group G. We establish existence of nontrivial G-equivariant entire solutions connecting the global minima of W along certain directions at infinity.     

Regularity and Singularities of Optimal Convex Shapes in the Plane

We focus here on the analysis of the regularity or singularity of solutions Ω ₀ to shape optimization problems among convex planar sets, namely:

$J(\Omega_{0})={\rm min} \{J(\Omega), \Omega \quad {\rm convex},\Omega \in \mathcal{S}_{\rm ad}\},$

where \({\mathcal{S}_{\rm ad}}\) is a set of 2-dimensional admissible shapes and \({J:\mathcal{S}_{\rm ad}\rightarrow\mathbb{R}}\) is a shape functional. Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We prove two types of results:

1. i)
   under a suitable convexity property of the functional J, we prove that Ω ₀ is a W ^2,p -set, \({p\in[1, \infty]}\). This result applies, for instance, with p = ∞ when the shape functional can be written as J(Ω) = R(Ω) + P(Ω), where R(Ω) = F(|Ω|, E _f (Ω), λ₁(Ω)) involves the area |Ω|, the Dirichlet energy E _f (Ω) or the first eigenvalue of the Laplace–Dirichlet operator λ₁(Ω), and P(Ω) is the perimeter of Ω;
    

1. ii)
   under a suitable concavity assumption on the functional J, we prove that Ω ₀ is a polygon. This result applies, for instance, when the functional is now written as J(Ω) = R(Ω) ? P(Ω), with the same notations as above.
    

     

Numerical Simulation of the Negative Magnus Effect of a Two-Dimensional Spinning Circular Cylinder

Based on the finite volume method, the flow past a spinning circular cylinder at a low subcritical Reynolds number (Re =1 × 10 ⁵), high subcritical Reynolds number (Re =1.3 ×10 ⁵), and critical Reynolds number (Re =1.4 ×10 ⁵) were each simulated using the Navier-Stokes equations and the γ-Re _?? transition model coupled with the SST k?ω turbulence model. The system was solved using an implicit algorithm. The freestream turbulence intensity decay was effectively controlled by the source term method proposed by Spalart and Rumsey. The variations in the Magnus force as a function of the spin ratio, α were obtained for the three Reynolds numbers, and the flow mechanism was analyzed. The results indicate that the asymmetric transitions induced by spin affect the asymmetric separations at the top and bottom surfaces of the circular cylinder, which further affects the pressure distributions at the top and bottom surfaces of the circular cylinder and ultimately result in a negative Magnus force, whose direction is opposite to that of the classical Magnus force. This study is the first to use a numerical simulation method to predict a negative Magnus force acting on a spinning circular cylinder. At the low subcritical Reynolds number, the Magnus force remained positive for all spin ratios. At the high subcritical Reynolds number, the sign of the Magnus force changed twice over the range of the spin ratio. At the critical Reynolds number, the sign of the Magnus force changed only once over the range of the spin ratio. For relatively low spin ratios, the Magnus force significantly differed by Reynolds number; however, this variation diminished as the spin ratio increased.     

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

A Simple Approach for Specifying Velocity Inflow Boundary Conditions in Simulations of Turbulent Opposed-Jet Flows

A new methodology is developed to specify inflow boundary conditions for the velocity field at the nozzle exit planes in turbulent counterflow simulations. The turbulent counterflow configuration consists of two coaxial opposed nozzles which emit highly-turbulent streams of varying species compositions depending on the mode considered. The specification of velocity inflow boundary conditions at the nozzle exits in the counterflow configuration is non-trivial because of the unique turbulence field generated by the turbulence generating plates (TGPs) upstream of the nozzle exits. In the method presented here, a single large-eddy simulation (LES) is performed in a large domain that spans the region between the TGPs of the nozzles, and the time series of the velocity fields at the nozzle exit planes are recorded. To provide inflow boundary conditions at the nozzle exit planes for simulations under other conditions (e.g., different stream compositions, bulk velocity, TGP location), transformations are performed on the recorded time series: the mean and r.m.s. (root-mean-square) quantities of velocity, as well as the longitudinal integral length scale on the centerline, at the nozzle exits in simulations are matched to those observed in experiments, thereby matching the turbulent Reynolds number R e _t . The method is assessed by implementing it in coupled large-eddy simulation/probability density function (LES/PDF) simulations on a small cylindrical domain between the nozzle exit planes for three different modes of the counterflow configuration: N ₂ vs. N ₂; N ₂ vs. hot combustion products; and C H ₄/N ₂ vs. O ₂. The inflow method is found to be successful as the first and second moments of velocity from the LES/PDF simulations agree well with the experimental data on the centerline for all three modes. This simple yet effective inflow strategy can be applied to eliminate the computational cost required to simulate the flow field upstream of the nozzle exits. It is also emphasized that, in addition to the predicted time series data, the availability of experimental data close to the nozzle exit planes plays a key role in the success of this method.     

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.     

Assessment of Five SGS Models for Low-Rem MHD Turbulence

Assessment of three regularization-based and two eddy-viscosity-based subgrid-scale (SGS) turbulence models for large eddy simulations (LES) are carried out in the context of magnetohydrodynamic (MHD) decaying homogeneous turbulence (DHT) with a Taylor scale Reynolds number (Re_λ) of 120 and a MHD transition-to-turbulence Taylor-Green vortex (TGV) problems with a Reynolds number of 3000, through direct comparisons to direct numerical simulations (DNS). Simulations are conducted using the low-magnetic Reynolds number approximation (Re_m<<1). LES predictions using the regularization-based Leray- α,LANS- α, and Clark- α SGS models, along with the eddy viscosity-based non-dynamic Smagorinsky and the dynamic Smagorinsky models are compared to in-house DNS for DHT and previous results for TGV. With regard to the regularization models, this work represents their first application to MHD turbulence. Analyses of turbulent kinetic energy decay rates, energy spectra, and vorticity fields made between the varying magnetic field cases demonstrated that the regularization models performed poorly compared to the eddy-viscosity models for all MHD cases, but the comparisons improved with increase in magnitude of magnetic field, due to a decrease in the population of SGS eddies within the flow field.     

Low Reynolds number effects in open-channel turbulent boundary layers

The present study reports measurements of a turbulent boundary layer in an open-channel flow using fiber-optic laser Doppler anemometry. The Reynolds numbers based on momentum thickness and depth of flow are in the range 750≤Re _θ ≤2,400 and 15,300≤Re _h ≤54,200, respectively. It is shown that an accurate estimate of the wall shear stress can be made by fitting a fifth-order polynomial to the near-wall data. The effect of Reynolds number on the mean turbulence intensity and triple correlation is examined using both conventional scaling laws and the recent scaling laws proposed by George and Castillo. The present results show that different scaling laws lead to different conclusions on low Reynolds number effects.     

Effect of Free-Stream Turbulence on the Reduction of Surface Friction in a Turbulent Boundary Layer Behind LEBU-Devices

The effect of increased free-stream turbulence on the reduction of the surface friction coefficient c _f in a turbulent boundary layer behind large-eddy break-up (LEBU) devices is investigated using a gravimetric method. The turbulence level was ε ≈ 1.9–4.9 % and the turbulence scale L _e ≈ 40–110 mm. The boundary layer Reynolds number Re** was varied from 2300 to 7500, with the boundary layer thickness being varied on the range δ = 33–44 mm. It is shown that an increase in the turbulence level ε has almost no impact on the relative reduction of friction behind LEBU-devices, whereas, under similar conditions of elevated free-stream turbulence, for another method, namely, the use of surface riblets, the friction reduction may be more strongly expressed.     

Numerical analysis of the flow pattern and vortex breakdown over a pitching delta wing at supersonic speeds

A supersonic compressible flow over a 60° swept delta wing with a sharp leading edge undergoing pitching oscillations is computationally studied. Numerical simulations are performed by the finite volume method with the use of the k?ω turbulence model for various Mach numbers and angles of attack. Variations of flow patterns in a crossflow plane, hysteresis loops associated with the vortex core location, and vortex breakdown positions during a pitching cycle are investigated. Trends for various Mach numbers, mean angles of attack, pitching amplitudes, and pitching frequencies are illustrated.     

Suspensions of titania nanoparticle networks in nematic liquid crystals: rheology and microstructure

We study the influence of confinement on the rheology and structure of nematic liquid crystals (NLCs). NLCs get confined in networks of titania (TiO₂, primary particle size = 21 nm) nanoparticles in suspensions of TiO₂ and NLC, N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Suspensions with TiO₂ nanoparticle volume fraction (?) of 0.006–0.017, form viscoelastic solids with low elastic modulus (G′) of 10¹ Pa–10² Pa and short relaxation times. Increase in TiO₂ nanoparticle ? leads to a rise in G′ with TiO₂ nanoparticles forming a percolating network at a critical volume fraction (? _c) = 0.023, and G′ of ~10³ Pa. TiO₂/MBBA NLC suspensions at and above ? _c = 0.023 show G′ ~ ω ^x?1 scaling, where ω is the angular frequency and the minimum in loss modulus (G′′) with ω. The effective noise temperature, x decreases and approaches 1 with the increase in the TiO₂ nanoparticle ? from 0.023–0.035, is indicative of an increase in the glassy dynamics. Through the polarized light microscopy and differential scanning calorimetry experiments, we propose that the progressive addition of TiO₂ nanoparticles introduces a quenched random disorder (QRD) in the NLC medium which disturbs the nematic order. This results in metastable TiO₂/MBBA NLC suspensions in which NLC domains get confined in the network of flocs of TiO₂ nanoparticles. We also show that the salient rheological signatures of soft glassy rheology develop only in the presence of NLC MBBA and are absent in the isotropic phase of MBBA.     

Effect of Inflow Turbulence on an Airfoil Flow with Laminar Separation Bubble: An LES Study

The present paper is concerned with numerical investigations on the effect of inflow turbulence on the flow around a SD7003 airfoil. At a Reynolds number Re_c =?60,000, an angle of attack α =?4^° and a low or zero turbulence intensity of the oncoming flow, the flow past the airfoil is known to be dominated by early separation, subsequent transition and reattachment leading to a laminar separation bubble with a distinctive pressure plateau. The objective of the study is to investigate the effect of inflow turbulence on the flow behavior. For this purpose, a numerical methodology relying on a wall-resolved large-eddy simulation, a synthetic turbulence inflow generator and a specific source term concept for introducing the turbulence fluctuations within the computational domain is used. The numerical technique applied allows the variation of the free-stream turbulence intensity (TI) in a wide range. In order to analyze the influence of TI on the arising instantaneous and time-averaged flow field past the airfoil, the present study evaluates the range 0% ≤ TI ≤?11.2%, which covers typical values found in atmospheric boundary layers. In accordance with experimental studies it is shown that the laminar separation bubble first shrinks and finally completely vanishes for increasing inflow turbulence. Consequently, the aerodynamic performance in terms of the lift-to-drag ratio increases. Furthermore, the effect of the time and length scales of the isotropic inflow turbulence on the development of the flow field around the airfoil is analyzed and a perceptible influence is found. Within the range of inflow scales studied decreasing scales augment the receptivity of the boundary layer promoting an earlier transition.