FEM solution to natural convection flow of a micropolar nanofluid in the presence of a magnetic field

The two-dimensional, laminar, unsteady natural convection flow in a square enclosure filled with aluminum oxide (\(\hbox {Al}_{2} \hbox {O}_{3}\))–water nanofluid under the influence of a magnetic field, is considered numerically. The nanofluid is considered as Newtonian and incompressible, the nanoparticles and water are assumed to be in thermal equilibrium. The mathematical modelling results in a coupled nonlinear system of partial differential equations. The equations are solved using finite element method (FEM) in space, whereas, the implicit backward difference scheme is used in time direction. The results are obtained for Rayleigh (Ra), Hartmann (Ha) numbers, and nanoparticles volume fractions (\(\phi\)), in the ranges of \(10^3 \le Ra \le 10^7\), \(0\le Ha \le 500\) and \(0 \le \phi \le 0.2\), respectively. The streamlines and microrotation contours are observed to show similar behaviors with altering magnitudes. For low Ra values, when \(Ha=0\), symmetric vortices near the walls and a central vortex in opposite direction are observed in vorticity. As Ra increases, the central vortex splits into two due to the circulation in the effect of the buoyant flow. Boundary layer formation is observed when Ha increases for almost all Rayleigh numbers in both streamlines and vorticity. The isotherms have horizontal profiles for high Ra values owing to convective dominance over conduction. As Ha is increased, the convection effect is reduced, and isotherms tend to have vertical profiles. This study presents the first FEM application for solving highly nonlinear PDEs defining micropolar nanofluid flow especially for large values of Rayleigh and Hartmann numbers.     

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.     

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.     

Biorthogonal stretching and shearing of an impermeable surface in a uniformly rotating fluid system

The flow induced by an impermeable flat surface executing orthogonal stretching and orthogonal shearing in a rotating fluid system is investigated. Both the stretching and shearing are linear in the coordinates. An exact similarity reduction of the Navier–Stokes equations gives rise to a pair of nonlinearly-coupled ordinary differential equations governed by three parameters. In this study we set one parameter and analyze the problem which leads to flow for an impermeable surface with shearing and stretching due to velocity u along the x-axis of equal strength a while the shearing and stretching due to velocity v along the y-axis of equal strength b. These solutions depend on two parameters—a Coriolis (rotation) parameter \(\sigma = \Omega /a\) and a stretching/shearing ratio \(\lambda =b/a\). A symmetry in solutions is found for \(\lambda = 1\). The exact solution for \(\sigma = 0\) and the asymptotic behavior of solutions for \(|\sigma | \rightarrow \infty\) are determined and compared with numerical results. Oscillatory solutions are found whose strength increases with increasing values of \(|\sigma |\). It is shown that these solutions tend to the well-known Ekman solution as \(|\sigma | \rightarrow \infty\).     

Two-component modeling for non-Newtonian nanofluid slip flow and heat transfer over sheet: Lie group approach

The present paper deals with the multiple solutions and their stability analysis of non-Newtonian micropolar nanofluid slip flow past a shrinking sheet in the presence of a passively controlled nanoparticle boundary condition. The Lie group transformation is used to find the similarity transformations which transform the governing transport equations to a system of coupled ordinary differential equations with boundary conditions. These coupled set of ordinary differential equation is then solved using the RungeKutta-Fehlberg fourth-fifth order(RKF45) method and the ode15 s solver in MATLAB.For stability analysis, the eigenvalue problem is solved to check the physically realizable solution. The upper branch is found to be stable, whereas the lower branch is unstable. The critical values(turning points) for suction(0 sc s) and the shrinking parameter(χc χ 0) are also shown graphically for both no-slip and multiple-slip conditions. Multiple regression analysis for the stable solution is carried out to investigate the impact of various pertinent parameters on heat transfer rates. The Nusselt number is found to be a decreasing function of the thermophoresis and Brownian motion parameters.     

Error estimate of the modified point-mass trajectory model of an artillery shell

The article deals with the motion of an axially symmetric spinning artillery shell in the gravity field under the action of the system of aerodynamic forces and moments adopted in ballistics. As the starting point, the system of differential equations of motion of the shell is taken, which is obtained from the original “accurate” system by its linearization in the variables describing the angular motion of the symmetry axis and by additional linearization in the angle between the velocity vector of the center of mass and the vertical plane (l-system). This article examines the system of differential equations of the translational motion and axial rotation of the shell which describes its modified point-mass trajectory model as applied to l-system (m-system). By small parameter methods, an estimate is obtained for the difference of the solution of l-system with given initial data and the solution of m-system with the same initial data for the variables of translational motion and axial rotation. This analytical evaluation is built in such a way that it corresponds with certain numerical estimates for components of the translational motion and axial rotation. It is observed that, under accepted assumptions, m-system and l-system determine the translational motion of the shell with the same order of the error as compared to the original “accurate” nonlinear system of equations of motion of the shell. But m-system does not contain rapidly oscillating variables describing the angular motion of the symmetry axis, and so its numerical integration requires tens of times less computational resources than the numerical integration of l-system. Numerical simulation data are represented.     

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.     

Instabilities of a thin viscoelastic liquid film flowing down an inclined plane in the presence of a uniform electromagnetic field

The effect of a uniform electromagnetic field on the stability of a thin layer of an electrically conducting viscoelastic liquid flowing down on a nonconducting inclined plane is studied under the induction-free approximation. Long-wave expansion method is used to obtain the surface evolution equation. The stabilizing role of the magnetic parameter M and the destabilizing role of the viscoelastic parameter Γ as well as the electric parameter E on this flow field are established. A novel result which emerges from our analysis is that the stabilizing effect of M holds no longer true for both viscous and viscoelastic fluids in the presence of electromagnetic field. It is found that when E exceeds a certain critical value depending on Γ, magnetic field exhibits the destabilizing effect on this flow field. Indeed, this critical value decreases with the increase of the viscoelastic parameter Γ since it has a destabilizing effect inherently. Another noteworthy result which arises from the weakly nonlinear stability analysis is that both the subcritical unstable and supercritical stable zones are possible together with the unconditional stable and explosive zones for different values of Γ depending on the wave number k.     

Mixed Convection Boundary-Layer Flow Along a Vertical Cylinder Embedded in a Porous Medium Filled by a Nanofluid

The steady mixed convection boundary-layer flow on a vertical circular cylinder embedded in a porous medium filled by a nanofluid is studied for both cases of a heated and a cooled cylinder. The governing system of partial differential equations is reduced to ordinary differential equations by assuming that the surface temperature of the cylinder and the velocity of the external (inviscid) flow vary linearly with the axial distance x measured from the leading edge. Solutions of the resulting ordinary differential equations for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction ${\phi}$ , the mixed convection or buoyancy parameter ?? and the curvature parameter ??. Results are presented for the specific case of copper nanoparticles. A critical value ?? _c of ?? with ?? _c <?0 is found, with the values of | ?? _c| increasing as the curvature parameter ?? or nanoparticle volume fraction ${\phi}$ is increased. Dual solutions are seen for all values of ?? >??? _c for both aiding, ?? >?0 and opposing, ?? <?0, flows. Asymptotic solutions are also determined for both the free convection limit ${(\lambda \gg 1)}$ and for large curvature parameter ${(\gamma \gg 1)}$ .     

Monte Carlo Assessment of the Impact of Oscillatory and Pulsating Boundary Conditions on the Flow Through Porous Media

The linear stability analysis of vertical throughflow of power law fluid for double-diffusive convection with Soret effect in a porous channel is investigated in this study. The upper and lower boundaries are assumed to be permeable, isothermal and isosolutal. The linear stability of vertical through flow is influenced by the interactions among the non-Newtonian Rayleigh number (Ra), Buoyancy ratio (N), Lewis number (Le), Péclet number (Pe), Soret parameter (Sr) and power law index (n). The results indicate that the Soret parameter has a significant influence on convective instability of power law fluid. It has also been noticed that buoyancy ratio has a dual effect on the instability of fluid flow. Further, it is noticed that the basic temperature and concentration profiles have singularities at \(Pe = 0\) and \(Le = 1\), the convective instability is looked into for the limiting case of \(Pe\rightarrow 0\) and \(Le \rightarrow 1\). For the case of pure thermal convection with no vertical throughflow, the present numerical results coincide with the solution of standard Horton–Rogers–Lapwood problem. The present results for critical Rayleigh number obtained using bvp4c and two-term Galerkin approximation are compared with those available in the literature and are tabulated.     

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.     

The Mixed Convection Boundary Layer on a Vertical Melting Front in a Non-Darcian Porous Medium

The effect of surface melting on the dual solutions that can arise in the problem of the mixed convection boundary-layer flow past a vertical surface embedded in a non-Darcian porous medium is considered. The problem is described by M, melting parameter, \(\lambda \), mixed convection parameter, and \(\gamma \), the flow inertia coefficient, numerical results being obtained in terms of these three parameters. It is seen that the melting phenomenon reduces the heat transfer rate and enhances the boundary-layer separation at the solid–liquid interface. Asymptotic solutions for the forced convection, \(\lambda =0\), and free convection, large \(\lambda \), limits are derived.     

Shear flow of a Newtonian fluid over a quiescent generalized Newtonian fluid

The flow of an upper shear-driven Newtonian fluid above an otherwise still non-Newtonian fluid is considered. The lower fluid is modelled as a generalized Newtonian fluid and set into motion by interfacial shear. By means of similarity transformations, the governing partial differential equations for the two-fluid problem transform exactly into two sets of ordinary differential equations coupled only at the interface. The successful transformation of the two-fluid problem is applied to the particular case when the lower fluid obeys power-law rheology. The resulting three-parameter problem is solved numerically for some different parameter combinations by means of a direct integration approach with the density ratio fixed to unity. We observed that the interfacial velocities decreased with increasing values of the power-law index n in the range from 0.6 to 1.4 whereas the shear-induced motion of the lower fluid penetrates far deeper into a shear-thinning (n < 1) than into a shear-thickening (n > 1) fluid. This phenomenon is ascribed to a corresponding increase of the non-linear viscosity function with lower n-values.     

A numerical study of micropolar flow inside a lid-driven triangular enclosure

A study is carried out to analyze the mixed convection flow and heat transfer inside a lid-driven triangular conduit under the effects of micro-gyration boundary conditions. The micropolar constitutive equation characterizes the fluid inside the cavity. The lower boundary is at a uniform temperature and sliding in its plane with constant velocity u₀, while the inclined walls are cold. Dual cases are considered here, namely the intense concentration (d) and the weak concentration of microelements (\(m = 0.5\)). The governing nonlinear equations are simulated employing the Galerkin finite element method, where the pressure term is handled via the Penalty approach. Using the numerical data, graphical results are produced to illustrate the effects of physical parameters. Specifically, this refers to the effects of the Grashof number (Gr), Prandtl number (Pr), Reynolds number (Re) and vortex viscosity parameter (K) on the streamlines, mid-section velocity profiles, temperature contours, and local and average Nusselt numbers on the cold and heated boundaries of the conduit. Particular emphasis is given on the identification of the set of parameters for which simultaneous symmetry in streamlines and isotherms prevails. The grid independence test is also performed by comparing the average Nusselt numbers (on the hot and cold boundaries of the conduit) for various mesh sizes, and the optimal solution is found. Moreover, the results are also benchmarked with the previously published data.     

Topological Structural Stability of Partial Differential Equations on Projected Spaces

In this paper we study topological structural stability for a family of nonlinear semigroups \(T_h(\cdot )\) on Banach space \(X_h\) depending on the parameter h. Our results shows the robustness of the internal dynamics and characterization of global attractors for projected Banach spaces, generalizing previous results for small perturbations of partial differential equations. We apply the results to an abstract semilinear equation with Dumbbell type domains and to an abstract evolution problem discretized by the finite element method.     

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

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

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

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.     

Pressure-Gradient Turbulent Boundary Layers Developing Around a Wing Section

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

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress,the modified couple stress theory(MCST),and the nonlocal elasticity theories using the differential quadrature method(DQM)is presented.Main advantages of the MCST over the classical theory(CT)are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter.Based on the nonlinear von K′arm′an assumption,the governing equations of equilibrium for the micro-classical plate considering midplane displacements are derived based on the minimum principle of potential energy.Using the DQM,the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained.Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature.A parametric study is conducted to show the effects of the aspect ratio,the side-to-thickness ratio,Eringen’s nonlocal parameter,the material length scale parameter,Young’s modulus of the surface layer,the surface residual stress,the polymer matrix coefficients,and various boundary conditions on the dimensionless uniaxial,biaxial,and shear critical buckling loads.The results indicate that the critical buckling loads are strongly sensitive to Eringen’s nonlocal parameter,the material length scale parameter,and the surface residual stress effects,while the effect of Young’s modulus of the surface layer on the critical buckling load is negligible.Also,considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate.The results show that the critical biaxial buckling load increases with an increase in G12/E2and vice versa for E1/E2.It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude.Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios,it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.