Dual solutions of Casson fluid flows over a power-law stretching sheet

A steady boundary layer flow of a non-Newtonian Casson fluid over a power-law stretching sheet is investigated. A self-similar form of the governing equation is obtained, and numerical solutions are found for various values of the governing parameters. The solutions depend on the fluid material parameter. Dual solutions are obtained for some particular range of these parameters. The fluid velocity is found to decrease as the power-law stretching parameter β in the rheological Casson equation increases. At large values of β, the skin friction coefficient and the velocity profile across the boundary layer for the Casson fluid tend to those for the Newtonian fluid.     

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.     

Newtonian,power law,and infinite shear flow characteristics of concentrated slurries using percolation theory concepts

There is a strong interest today in concentrated particulate-filled dispersion and slurries in both polymeric and Newtonian fluids. This paper reviews and extends theoretical approaches using percolation theory concepts to characterize the rheological behavior of fluids filled with particulate solids. First, a previously proposed limiting, zero shear viscosity model based on percolation theory concepts is reviewed. This model has been primarily tested with rigid fillers in a Newtonian carrier and polymeric fluids. Second, all Newtonian fluid-based slurries that have a high concentration of filler become pseudoplastic, shear-thinning slurries at some threshold shear rate. A new theory is reviewed and new data are evaluated that correlate the power law constant, n, to cluster formation of the fillers suspended in the fluids in shear flow. Slurry systems reported here cover a size range from 58 nm to 200 μm. Third, this cluster percolation-based rheological analysis is then extended to a newly proposed model for the calculation of the ratio of infinite shear, η_∞, to the zero shear viscosity, η₀. Using literature data, it is demonstrated that measurements of the viscosity ratio, η_∞/η₀, correlate with the power law through the use of an energy dissipation-based model for Bingham rheological fluids.     

Aerodynamic characteristics of a thin airfoil cascade in an ideal incompressible flow with separation from the leading edges

A new method of solving the problem of separation flow past a cascade of thin airfoils is proposed. The method is based on the model of separated flow past a single airfoil in which vortex wakes shed from the airfoil edges are modeled by vortex layers with time-averaged intensities. Under the assumption of a small deviation of the freestream angle α from the “impactless entry” angle α ₀ the singular integral equation governing the flow kinematics is closed by a dynamic equation. It is shown that the separation effect on the time-averaged aerodynamic characteristics of the cascade is associated with the total pressure loss due to the flow energy expenditure on the vortex wake formation. The aerodynamic characteristics of the cascade calculated with and without account for flow separation differ by the second-order quantity ? = α ? α ₀.     

Some passages to the limit in elasticity theory and “Sapondzhyan’s paradox”

We consider bending of thin plates with polygonal and curvilinear edges and indicate analogies and differences between the boundary conditions and boundary value problems arising in these two cases if the polygon is inscribed in the curvilinear contour and the number k of vertices of the polygon tends to infinity.We believe that the so-called Sapondzhyan paradox that arises when solving the boundary value problems for supported plates with a curvilinear contour and a k-gonal contour inscribed in it as k → ∞ can be called a paradox only by misunderstanding. Sapondzhyan’s paradox was studied in several papers briefly surveyed in the monograph [1]. Apparently, the interpretation of “paradoxes” and the results proposed in the present paper are published for the first time.Sapondzhyan’s paradox can be generalized to the case of bending of the so-called sliding-fixed plates (i.e., the generalized shear force and the rotation angle are zero on the plate contour) with a curvilinear contour and a k-gonal contour inscribed in it as k → ∞.In the case of three-dimensional elasticity problems, we present boundary conditions and boundary value problems similar to those listed above and consider the situations resulting in “paradoxes” similar to those arising in plate bending. We give the corresponding explanations and interpretations.     

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.     

Integral boundary layer relations in the theory of wave flows for capillary liquid films

A generalized method of deriving the model equations is considered for wave flow regimes in falling liquid films. The viscous liquid equations are used on the basis of integral boundary layer relations with weight functions. A family of systems of evolution differential equations is proposed. The integer parameter n of these systems specifies the number of a weight function. The case n = 0 corresponds to the classical IBL (Integral Boundary Layer) model. The case n ≥ 1 corresponds to its modifications called the WIBL (Weighted Integral Boundary Layer) models. The numerical results obtained in the linear and nonlinear approximations for n = 0, 1, 2 are discussed. The numerical solutions to the original hydrodynamic differential equations are compared with experimental data. This comparison leads us to the following conclusions: as a rule, the most accurate solutions are obtained for n = 0 in the case of film flows on vertical and inclined solid surfaces and the accuracy of solutions decreases with increasing n. Hence, the classical IBL model has an advantage over the WIBL models.     

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.     

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.     

Asymptotic Behavior of Type I Blowup Solutions to a Parabolic-Elliptic System of Drift–Diffusion Type

We consider a Cauchy problem for a parabolic-elliptic system of drift–diffusion type. The problem is formally of the form

$ U_t = \nabla \cdot (\nabla U-U \nabla (-\Delta)^{-1}U). $

This system describes a mass-conserving aggregation phenomenon including gravitational collapse and bacterial chemotaxis. Our concern is the asymptotic behavior of blowup solutions when the blowup is type I, in the sense that its blowup rate is the same as the corresponding ordinary differential equation y _t  = y ² (up to a multiple constant). It is shown that all type I blowup is asymptotically (backward) self-similar, provided that the solution is radial, nonnegative when the blowup set is a singleton and the space dimension is greater than or equal to three.

     

Effects of the Reynolds number on a scale-similarity model of Lagrangian velocity correlations in isotropic turbulent flows

A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η≡η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η∝ R_λ~(1.39) obtained from the numerical data.     

Vortex solitons of the (3+1)-dimensional spatially modulated cubic–quintic nonlinear Schrödinger equation with the transverse modulation

Vortex solitons in the spatially modulated cubic–quintic nonlinear media are governed by a (3+1)-dimensional cubic–quintic nonlinear Schrödinger equation with spatially modulated nonlinearity and transverse modulation. Via the variable separation principle with the similarity transformation, we derive two families of vortex soliton solutions in the spatially modulated cubic–quintic nonlinear media. For the disappearing and parabolic transverse modulation, vortex solitons with different configurations are constructed. The similar configurations of vortex solitons exist for the same value of \(l-k\) with the topological charge k and degree number l. Moreover, the number of the inner layer structure of vortex solitons getting rid of the package covering layer is related to \((n-1)/2+1\) with the soliton order number n. For the disappearing transverse modulation, there exist phase azimuthal jumps around their cores of vortex solitons with \(2\pi \) phase change in every jump, and any two jumps one after another realize the change in \(\pi \). For the parabolic transverse modulation, all phases of vortex soliton exist k-jump, and every jump realizes the change in \(2\pi /k\); thus, k-jumps totally realize the azimuthal change in \(2\pi \) around their cores.     

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).

     

Characterization of Viscoelasticity Due to Shear Wave Propagation: a Comparison of Existing Methods Based on Computational Simulation and Experimental Data

Shear wave propagation causes microvibrations within a medium; measuring the wave attenuation coefficient, α, and phase velocity, c _s , the medium shear modulus, μ, and shear viscosity, η, are determined based on a viscoelastic model that includes both c _s and α. The present work compares the performances of nine processing methods, based on cross-correlation and quadrature demodulation, used to extract the motion waveform from a sequence of radio-frequency (RF) echo signals from the medium. Kalman filtering determined the amplitude and the phase of the extracted motion waveform. The comparisons were done with regard to computational simulation and experiments with a gel phantom. Estimates obtained for μ and η of the medium considered different conditions for the vibration amplitude and the signal-to-noise ratio (SNR) of the RF echo signals and the waveform extracted by means of single frequency and shear wave dispersion ultrasound vibration (SDUV) methods. According to the simulated results, the cross-correlation-based processing techniques are more precise and accurate in comparison to quadrature demodulation techniques. The results for c _s , α, μ and η of the phantom and those obtained under the same setup conditions for experimental and computational tests agree with each other. Comparing the estimates based on single frequency and SDUV techniques, they presented similar performances at high SNR of the RF echo signal. On the other hand, the former technique prevailed for low SNR.     

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.     

Stability of plane oscillations of a satellite in a circular orbit

We study motions of a rigid body (a satellite) about the center of mass in a central Newtonian gravitational field in a circular orbit. There is a known particular motion of the satellite in which one of its principal central axes of inertia is perpendicular to the orbit plane and the satellite itself exhibits plane pendulum-like oscillations about this axis. Under the assumption that the satellite principal central moments of inertia A, B, and C satisfy the relation B = A + C corresponding to the case of a thin plate, we perform rigorous nonlinear analysis of the orbital stability of this motion.In the plane of the problem parameters, namely, the oscillation amplitude ε and the inertial parameter, there exist countably many domains of orbital stability of the satellite oscillations in the linear approximation. Nonlinear orbital stability analysis was carried out in thirteen of these domains. Isoenergetic reduction of the system of equations of the perturbed motion is performed at the energy level corresponding to the unperturbed periodic motion. Further, using the algorithm developed in [1], we construct the symplectic mapping generated by the equations of the reduced system, normalize it, and analyze the stability. We consider resonance and nonresonance cases. For small values of the oscillation amplitude, we perform analytic investigations; for arbitrary values of ε, numerical analysis is used.Earlier, numerical analysis of stability of plane pendulum-like motions of a satellite in a circular orbit was performed in several special cases in [1–4].     

Effect of Flow Structures on Turbulence Statistics of Taylor-Couette Flow in the Torque Transition State

Direct numerical simulations of Taylor-Couette flow from Re= 8000 to 25000 have been conducted to investigate changes of turbulence statistics in the transition of the Reynolds number dependency of the mean torque near Re= 10000. The velocity fluctuations are decomposed into the contributions of the Taylor vortex and remaining turbulent fluctuations. Significant Reynolds number dependencies of these components are observed in the radial profiles of the Reynolds stress and the transmission of the mean torque. The contributions of Taylor vortex and turbulent components in the net amount of mean torque are evaluated. The Taylor vortex component is overtaken by the turbulent counterpart around Re= 15000 when they are defined as the azimuthally averaged component and the remnants. The results show that the torque transition can be explained by the competition between the contributions of azimuthally averaged Taylor vortex and the remaining turbulent components.     

An Eulerian–Lagrangian Form for the Euler Equations in Sobolev Spaces

In 2000 Constantin showed that the incompressible Euler equations can be written in an “Eulerian–Lagrangian” form which involves the back-to-labels map (the inverse of the trajectory map for each fixed time). In the same paper a local existence result is proved in certain Hölder spaces \({C^{1,\mu}}\). We review the Eulerian–Lagrangian formulation of the equations and prove that given initial data in H ^s for \({n \geq 2}\) and \({s > \frac{n}{2}+1}\), a unique local-in-time solution exists on the n-torus that is continuous into H ^s and C ¹ into H ^s-1. These solutions automatically have C ¹ trajectories. The proof here is direct and does not appeal to results already known about the classical formulation. Moreover, these solutions are regular enough that the classical and Eulerian–Lagrangian formulations are equivalent, therefore what we present amounts to an alternative approach to some of the standard theory.     

Morse Index and Linear Stability of the Lagrangian Circular Orbit in a Three-Body-Type Problem Via Index Theory

It is well known that the linear stability of the Lagrangian elliptic solutions in the classical planar three-body problem depends on a mass parameter β and on the eccentricity e of the orbit. We consider only the circular case (e = 0) but under the action of a broader family of singular potentials: α-homogeneous potentials, for \(\alpha \in (0, 2)\), and the logarithmic one. It turns out indeed that the Lagrangian circular orbit persists also in this more general setting. We discover a region of linear stability expressed in terms of the homogeneity parameter α and the mass parameter β, then we compute the Morse index of this orbit and of its iterates and we find that the boundary of the stability region is the envelope of a family of curves on which the Morse indices of the iterates jump. In order to conduct our analysis we rely on a Maslov-type index theory devised and developed by Y. Long, X. Hu and S. Sun; a key role is played by an appropriate index theorem and by some precise computations of suitable Maslov-type indices.     

Gradient Bounds for Elliptic Problems Singular at the Boundary

Let Ω be a bounded smooth domain in \({{R}^N, N \geqq 2}\), and let us denote by d(x) the distance function d(x, ?Ω). We study a class of singular Hamilton–Jacobi equations, arising from stochastic control problems, whose simplest model is

$ - \alpha \Delta u+ u + \frac{\nabla u \cdot B (x)}{d (x)}+ c(x) |\nabla u|^2=f (x) \quad {\rm in}\,\Omega, $

where f belongs to \({W^{1,\infty}_{\rm loc} (\Omega)}\) and is (possibly) singular at \({\partial \Omega, c\in W^{1,\infty} (\Omega)}\) (with no sign condition) and the field \({B\in W^{1,\infty} (\Omega)^N}\) has an outward direction and satisfies \({B\cdot \nu\geqq \alpha}\) at ?Ω (ν is the outward normal). Despite the singularity in the equation, we prove gradient bounds up to the boundary and the existence of a (globally) Lipschitz solution. We show that in some cases this is the unique bounded solution. We also discuss the stability of such estimates with respect to α, as α vanishes, obtaining Lipschitz solutions for first order problems with similar features. The main tool is a refined weighted version of the classical Bernstein method to get gradient bounds; the key role is played here by the orthogonal transport component of the Hamiltonian.