Mixed Convection Boundary Layer Flow from a Horizontal Circular Cylinder Embedded in a Porous Medium Filled with a Nanofluid

Steady mixed convection boundary layer flow from an isothermal horizontal circular cylinder embedded in a porous medium filled with a nanofluid has been studied for both cases of a heated and cooled cylinder. The resulting system of nonlinear partial differential equations is solved numerically using an implicit finite-difference scheme. The solutions for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction φ and the mixed convection parameter λ. Three different types of nanoparticles are considered, namely Cu, Al₂O₃ and TiO₂. It is found that for each particular nanoparticle, as the nanoparticle volume fraction φ increases, the magnitude of the skin friction coefficient decreases, and this leads to an increase in the value of the mixed convection parameter λ which first produces no separation. On the other hand, it is also found that of all the three types of nanoparticles considered, for any fixed values of φ and λ, the nanoparticle Cu gives the largest values of the skin friction coefficient followed by TiO₂ and Al₂O₃. Finally, it is worth mentioning that heating the cylinder (λ > 0) delays separation of the boundary layer and if the cylinder is hot enough (large values of λ > 0), then it is suppressed completely. On the other hand, cooling the cylinder (λ < 0) brings the boundary layer separation point nearer to the lower stagnation point and for a sufficiently cold cylinder (large values of λ < 0) there will not be a boundary layer on the cylinder.     

Equivalence between Rank-One Convexity and Polyconvexity for Some Classes of Elastic Materials

Let W(F) = φ(λ ₁ ^s + λ ₂ ^s + λ ₃ ^s ) + ψ(λ ₁ ^r λ ₂ ^r + λ ₁ ^r λ ₃ ^r + λ ₂ ^r λ ₃ ^r ) + f(λ ₁ λ ₂ λ ₃) be a stored energy function. We prove that, for this function, rank-one convexity is equivalent to polyconvexity.under suitable assumptions on φ, ψ and f.     

Existence of the minimal positive solution of some nonlinear elliptic systems when the nonlinearity is the sum of a sublinear and a superlinear term

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｉｓｐａｐｅｒ,ｗｅｃｏｎｓｉｄｅｒｔｈｅｅｌｌｉｐｔｉｃｓｙｓｔｅｍ(1λ)　-Δｕ=ｆ(λ,ｘ,ｕ)-ｖ　　(ｉｎΩ),-Δｖ=δｕ-γｖ(ｉｎΩ),ｕ=ｖ=0(ｏｎΩ),ｗｈｅｒｅΩｉｓａｓｍｏｏｔｈｂｏｕｎｄｅｄｄｏｍａｉｎｉｎＲＮ(Ｎ≥2)ａｎｄλｉｓａｒｅａｌｐａｒａｍｅｔｅｒ.Ｔｈｅｓｏｌｕｔｉｏｎｓ(ｕ,ｖ)ｏｆｔｈｉｓｓｙｓｔｅｍｒｅｐｒｅｓｅｎｔｓｔｅａｄｙｓｔａｔｅｓｏｌｕｔｉｏｎｓｏｆｒｅａｃｔｉｏｎｄｉｆｆｕｓｉｏｎｓｙｓｔｅｍｓｄｅｒｉｖｅｄｆｒｏｍｓｅｖｅｒａｌａｐ…     

A Variational Approach to the Fredholm Alternative for the p-Laplacian near the First Eigenvalue

We are concerned with the existence of a weak solution to the degenerate quasi-linear Dirichlet boundary value problem

On a Crystalline Variational Problem, Part I:¶First Variation and Global L∞ Regularity

Let φ:ℝ ⁿ → [0,+∞[ be a given positively one-homogeneous convex function, and let ?_φ≔{φ≤ 1 }. Pursuing our interest in motion by crystalline mean curvature in three dimensions, we introduce and study the class ?_φ (ℝ ⁿ ) of “smooth” boundaries in the relative geometry induced by the ambient Banach space (ℝ ⁿ , φ). It can be seen that, even when ?_φ is a polytope, ?_φ(ℝ ⁿ ) cannot be reduced to the class of polyhedral boundaries (locally resembling ∂?_φ). Curved portions must be necessarily included and this fact (as well as the nonsmoothness of ∂?_φ) is the source of several technical difficulties related to the geometry of Lipschitz manifolds. Given a boundary δE in the class ?_φ(ℝ ⁿ ), we rigorously compute the first variation of the corresponding anisotropic perimeter, which leads to a variational problem on vector fields defined on δE. It turns out that the minimizers have a uniquely determined (intrinsic) tangential divergence on δE. We define such a divergence to be the φ-mean curvature κ_φ of δE; the function κ_φ is expected to be the initial velocity of δE, whenever δE is considered as the initial datum for the corresponding anisotropic mean curvature flow. We prove that κ_φ is bounded on δE and that its sublevel sets are characterized through a variational inequality.     

Exact analytic solutions for free convection boundary layers on a heated vertical plate with lateral mass flux embedded in a saturated porous medium

 The problem of the self-similar boundary flow of a “Darcy-Boussinesq fluid” on a vertical plate with temperature distribution T _w(x) = T _∞+A·x ^λ and lateral mass flux v _w(x) = a·x ^(λ−1)/2, embedded in a saturated porous medium is revisited. For the parameter values λ = 1,−1/3 and −1/2 exact analytic solutions are written down and the characteristics of the corresponding boundary layers are discussed as functions of the suction/ injection parameter in detail. The results are compared with the numerical findings of previous authors. Received on 8 March 1999     

The mathematical principles of vibration reductor

In engineering and technology, it is often demanded that self-oscillation.be eliminated.so that the equipment or machinery may not be damaged In this paper, a mathematicalmodel for reducing vibration is given by the following equations:(?)_1+(?)((?)_1) +k_1(x_1-x_2) =0, (?)_2+c(?)_1+k_2(x_2-x_1) =0 (*)We have discussed how to choose suitable parameters c_1, k_1,k_2; of equations (*),so as to make the zero solution to be of global stability. Several theorems on the globalstability of the zero solution of equations (*) are also given.     

On a Crystalline Variational Problem, Part II:¶BV Regularity and Structure of Minimizers on Facets

For a nonsmooth positively one-homogeneous convex function φ:ℝ ⁿ → [0,+∞[, it is possible to introduce the class ?_φ (ℝ ⁿ ) of smooth boundaries with respect to φ, to define their φ-mean curvature κ_φ, and to prove that, for E∈?_φ (ℝ ⁿ ), κ_φ∈L ^∞(δE) [9]. Based on these results, we continue the analysis on the structure of δE and on the regularity properties of κ_φ. We prove that a facet F of δE is Lipschitz (up to negligible sets) and that κ_φ has bounded variation on F. Further properties of the jump set of κ_φ are inspected: in particular, in three space dimensions, we relate the sublevel sets of κ_φ on F to the geometry of the Wulff shape ?_φ≔{φ≤ 1 }. Accepted October 11, 2000?Published online 14 February, 2001     

Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Strain energy function of an uncross-linked butadiene rubber estimated from planar extension

The derivatives of the strain energy function u with respect to the invariants of the strain tensor (I₁ and I₂) are estimated for uncross-linked butadiene rubber by using the BKZ constitutive equation. The derivatives at small deformations show anomalous behavior; namely, an upturn for ∂u/∂I₁ and a downturn for ∂u/∂I₂ take place, as is the case of cross-linked rubbers. At large deformations, u is well described by u = A₁(I₁ −3) + A₂(I₂ −3) with numerical constants A₁ and A₂. This behavior is also quite similar to that for cross-linked rubbers. The non-zero positive constant A₂ for the melt suggests that the non-zero value is due to neither the inhomogeneity in network structure nor high extension of constituent polymer chains.     

Axial annular flow of plastic fluids: dead zones and plug-free flow

In axial annular flow, the shear stress decreases from its value τ(κR) at the inner cylinder to 0 at r = λR and increases from then on to τ(R) at the outer cylinder. For plastic fluids with a yield stress τ _c, λ will be such that flow commences when τ(κR) = τ(R) = τ _c. For fluids with position-dependent yield stresses (electro- and magnetorheological fluids are examples), the situation is more complex. While it is possible that yielding and flow occur everywhere, it is also possible that flow occurs only in parts of the fluid-filled space, and a dead zone (region in which the fluid is at rest) close to one of the walls exists. In that case, the fluid will flow no matter how small the applied pressure difference is. If P is large enough, the dead zone ceases to exist and flow without any plug is possible. The fluid flows as if no yield stress exists.

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997     

Twist viscosity of mixtures of low molar mass nematics

Measurements of the twist viscosity, γ₁(DLS) and twist elastic coefficient, K₂₂(DLS) by electric-field-dependent dynamic light scattering (EFDLS) are reported for low molar mass nematics (LMMNs) 4′-heptyl-4-cyanobiphenyl (7CB) and 4′-octyl-4-cyanobiphenyl (8CB), and their binary mixtures at several temperatures in the nematic state. The results are compared with values (γ₁(Rheol)=α₃–α₂) computed from rheological measurements of the Leslie viscosities α₂ and α₃. For the binary mixtures, at each temperature, the measured twist viscosity γ₁(DLS) and corresponding twist elastic constant K₂₂(DLS) show approximately a linear additive dependence on concentration. The calculated twist viscosity, γ₁(Rheol), agrees with γ₁(DLS) for the pure components, but is significantly smaller for the binary mixtures. Our observations appear to be consistent with a recent report of a discrepancy between values of the tumbling parameter λ, determined using a small-strain oscillatory optical technique, vs those measured by a rheological method. These results suggest that, in the rheological measurements at large strains, the rate of director rotation for mixtures may be affected by a flow-induced change in structure, e.g., shear-induced biaxiality. Received: 17 March 2000 Accepted: 17 July 2000     

Steady mixed convection boundary-layer flow over a vertical flat surface in a porous medium filled with water at 4°C: variable surface heat flux

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water at 4°C (maximum density) when the surface heat flux varies as x ^m and the velocity outside the boundary layer varies as x ^(1+2m)/2, where x measures the distance from the leading edge, is discussed. Assisting and opposing flows are considered with numerical solutions of the governing equations being obtained for general values of the flow parameters. For opposing flows, there are dual solutions when the mixed convection parameter λ is greater than some critical value λ _c (dependent on the power-law index m). For assisting flows, solutions are possible for all values of λ. A lower bound on m is found, m > −1 being required for solutions. The nature of the critical point λ _c is considered as well as various limiting forms; the forced convection limit (λ = 0), the free convection limit (λ → ∞) and the limits as m → ∞ and as m → −1.     

Mixed Convection Boundary-Layer Flow in a Porous Medium Filled with Water Close to its Maximum Density

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water close to its maximum density is considered for uniform wall temperature and outer flow. The problem can be reduced to similarity form and the resulting equations are examined in terms of a mixed convection parameter λ and a parameter δ which measures the difference between the ambient temperature and the temperature at the maximum density. Both assisting (λ > 0) and opposing flows (λ < 0) are considered. A value δ₀ is found for which there are dual solutions for a range λ_c < λ < 0 of λ (the value of λ_c dependent on δ) and single solutions for all λ ≥ 0. Another value of δ₁ of δ, with δ₁ > δ₀, is found for which there are dual solutions for a range 0 < λ < λ_c of positive values of λ, with solutions for all λ≤ 0. There is also a range δ₀ <  δ < δ₁ where there are solutions only for a finite range of λ, with critical points at both positive and negative values of λ, thus putting a finite limit on the range of existence of solutions.     

Shear localisation with 2D viscous froth and its relation to the continuum model

Simulations of monodisperse and polydisperse (μ ₂(A) = 0.13±0.002) 2D foam samples undergoing simple shear are performed using the 2D viscous froth (VF) model. These simulations clearly demonstrate shear localisation. The dependence of localisation length on the product λV (shearing velocity V times the wall drag coefficient λ) is examined and is shown to agree qualitatively with published experimental data. A wide range of localisation lengths is found at low λV, an effect which is attributed to the existence of distinct yield and limit stresses. The general continuum model is extended to incorporate such an effect and its parameters are subsequently related to those of the VF model. A Herschel–Bulkley exponent of a = 0.3 is shown to accurately describe the observed behaviour. The localisation length is found to be independent of λV for monodisperse foam samples.     

A proof of the Benjamin-Feir instability

The existence and linear stability problem for the Stokes periodic wavetrain on fluids of finite depth is formulated in terms of the spatial and temporal Hamiltonian structure of the water-wave problem. A proof, within the Hamiltonian framework, of instability of the Stokes periodic wavetrain is presented. A Hamiltonian center-manifold analysis reduces the linear stability problem to an ordinary differential eigenvalue problem on ℝ⁴. A projection of the reduced stability problem onto the tangent space of the 2-manifold of periodic Stokes waves is used to prove the existence of a dispersion relation Λ(λ,σ, I ₁, I ₂)=0 where λ ε ℂ is the stability exponent for the Stokes wave with amplitude I ₁ and mass flux I ₂ and σ is the “sideband’ or spatial exponent. A rigorous analysis of the dispersion relation proves the result, first discovered in the 1960's, that the Stokes gravity wavetrain of sufficiently small amplitude is unstable for F ε (0,F₀) where F ₀ ≈ 0.8 and F is the Froude number.     

Linear viscoelastic behavior of perfluorooctyl sulfonate micelles: effects of counter-ions

Linear viscoelastic behavior was investigated for aqueous solutions of perfluorooctyl sulfonate (C₈F₁₇SO⁻ ₃; abbreviated as FOS) micelles having a mixture of tetraethylammonium (N⁺(C₂H₅)₄; TEA) and lithium (Li⁺) ions as the counter-ions. The solutions had the same FOS concentration (0.1 mol l⁻¹) and various Li⁺ fractions in the counter-ions, φ_Li = 0−0.6, and the FOS micelles in these solutions formed threads which further organized into dendritic networks. At T ≤ 15 °C, the terminal relaxation time τ and the viscosity η, governed by thermal scission of the networks, increased with increasing φ_Li up to 0.55. A further increase of φ_Li resulted in decreases of τ and η and in broadening of the relaxation mode distribution. These rheological changes are discussed in relation to the role of TEA ions in thermal scission: Previous NMR studies revealed that only a fraction of TEA ions were tightly bound to the FOS micellar surfaces and these bound ions stabilized the thread/network structures. The concentration of non-bound TEA ions, C_TEA ^*, decreased and finally vanished on increasing φ_Li up to φ_Li ^* ≅ 0.6, and the concentration of the bound TEA ions significantly decreased on a further increase of φ_Li. The non-bound TEA ions appeared to catalyze the thermal scission of the FOS threads, and the observed increases of τ and η for φ_Li < 0.55 were attributed to the decrease of C_TEA ^*. On the other hand, the decreases of τ and η as well as the broadening of the mode distribution, found for φ_Li > 0.55 (where C_TEA ^* ≅ 0), were related to destabilization of the FOS threads/networks due to a shortage of the bound TEA ions and to the existence of concentrated Li⁺ ions. Viscoelastic data of pure FOSTEA and FOSTEA/FOSLi/TEACl solutions lent support to these arguments for the role of TEA ions in the relaxation of FOSTEA/FOSLi solutions. Received: 12 October 1999/Accepted: 1 November 1999     

Entire Solutions with Merging Fronts to Reaction–Diffusion Equations

We deal with a reaction–diffusion equation u _t = u _xx + f(u) which has two stable constant equilibria, u = 0, 1 and a monotone increasing traveling front solution u = φ(x + ct) (c > 0) connecting those equilibria. Suppose that u = a (0 < a < 1) is an unstable equilibrium and that the equation allows monotone increasing traveling front solutions u = ψ₁(x + c ₁ t) (c ₁ < 0) and ψ₂(x + c ₂ t) (c ₂ > 0) connecting u = 0 with u = a and u = a with u = 1, respectively. We call by an entire solution a classical solution which is defined for all . We prove that there exists an entire solution such that for t≈ − ∞ it behaves as two fronts ψ₁(x + c ₁ t) and ψ₂(x + c ₂ t) on the left and right x-axes, respectively, while it converges to φ(x + ct) as t→∞. In addition, if c > − c ₁, we show the existence of an entire solution which behaves as ψ₁( − x + c ₁ t) in and φ(x + ct) in for t≈ − ∞.     

Initial value problems for a class of nonlinear integro-partial differential equations

１ＰｒｏｂｌｅｍｓａｎｄＭａｉｎＲｅｓｕｌｔｓＩｎｔｈｉｓｐａｐｅｒ,ｗｅｓｔｕｄｙｔｈｅｎｏｎｌｉｎｅａｒｖｉｂｒａｔｉｏｎｓｏｆｉｎｆｉｎｉｔｅｒｏｄｓｗｉｔｈｖｉｓｃｏｅｌａｓｔｉｃｉｔｙ．Ｔｈｅｃｏｎｓｔｉｔｕｔｉｏｎｌａｗｏｆｔｈｅｒｏｄｓ...