Oscillation criteria for mixed neutral differential equations

Abdulaziz khalid Alsharidi; Ali Muhib; Abdulaziz khalid Alsharidi; Ali Muhib

doi:10.3934/math.2024703

AIMS Mathematics

2024, Volume 9, Issue 6: 14473-14486. doi: 10.3934/math.2024703

Previous Article Next Article

Research article

Oscillation criteria for mixed neutral differential equations

Abdulaziz khalid Alsharidi ^{1
,
,},
Ali Muhib ^{2,3
,
,}

1.
Department of Mathematics and Statistics, College of Science, King Faisal University, Al Ahsa 31982, Saudi Arabia
2.
Department of Mathematics, Faculty of Applied and Educational Sciences, Al-Nadera, Ibb University, Ibb, Yemen
3.
Department of Mathematics, Faculty of Science, Mansoura University, 35516 Mansoura, Egypt

Received: 09 February 2024 Revised: 07 April 2024 Accepted: 11 April 2024 Published: 22 April 2024
MSC : 34C10, 34K11

In this study, we aim to contribute to the increasing interest in functional differential equations by obtaining new theorems for the oscillation of second-order neutral differential equations of mixed type in a non-canonical form. The results obtained here improve and extend those reported in the literature. The applicability of the results is illustrated by several examples.

Keywords:

neutral differential equations of mixed type,
second-order,
non-canonical form

Citation: Abdulaziz khalid Alsharidi, Ali Muhib. Oscillation criteria for mixed neutral differential equations[J]. AIMS Mathematics, 2024, 9(6): 14473-14486. doi: 10.3934/math.2024703

Related Papers:

[1]	Khalil Ur Rehman, Wasfi Shatanawi, Zead Mustafa . Artificial intelligence (AI) based neural networks for a magnetized surface subject to tangent hyperbolic fluid flow with multiple slip boundary conditions. AIMS Mathematics, 2024, 9(2): 4707-4728. doi: 10.3934/math.2024227
[2]	Sami Ul Haq, Saeed Ullah Jan, Syed Inayat Ali Shah, Ilyas Khan, Jagdev Singh . Heat and mass transfer of fractional second grade fluid with slippage and ramped wall temperature using Caputo-Fabrizio fractional derivative approach. AIMS Mathematics, 2020, 5(4): 3056-3088. doi: 10.3934/math.2020198
[3]	Kottakkaran Sooppy Nisar, Muhammad Shoaib, Muhammad Asif Zahoor Raja, Yasmin Tariq, Ayesha Rafiq, Ahmed Morsy . Design of neural networks for second-order velocity slip of nanofluid flow in the presence of activation energy. AIMS Mathematics, 2023, 8(3): 6255-6277. doi: 10.3934/math.2023316
[4]	Muhammad Asif Zahoor Raja, Kottakkaran Sooppy Nisar, Muhammad Shoaib, Ajed Akbar, Hakeem Ullah, Saeed Islam . A predictive neuro-computing approach for micro-polar nanofluid flow along rotating disk in the presence of magnetic field and partial slip. AIMS Mathematics, 2023, 8(5): 12062-12092. doi: 10.3934/math.2023608
[5]	Yousef Jawarneh, Humaira Yasmin, Wajid Ullah Jan, Ajed Akbar, M. Mossa Al-Sawalha . A neural networks technique for analysis of MHD nano-fluid flow over a rotating disk with heat generation/absorption. AIMS Mathematics, 2024, 9(11): 32272-32298. doi: 10.3934/math.20241549
[6]	Nadeem Abbas, Wasfi Shatanawi, Taqi A. M. Shatnawi, Fady Hasan . Theoretical analysis of induced MHD Sutterby fluid flow with variable thermal conductivity and thermal slip over a stretching cylinder. AIMS Mathematics, 2023, 8(5): 10146-10159. doi: 10.3934/math.2023513
[7]	Ali Raza, Umair Khan, Aurang Zaib, Wajaree Weera, Ahmed M. Galal . A comparative study for fractional simulations of Casson nanofluid flow with sinusoidal and slipping boundary conditions via a fractional approach. AIMS Mathematics, 2022, 7(11): 19954-19974. doi: 10.3934/math.20221092
[8]	W. Abbas, Ahmed M. Megahed, Osama M. Morsy, M. A. Ibrahim, Ahmed A. M. Said . Dissipative Williamson fluid flow with double diffusive Cattaneo-Christov model due to a slippery stretching sheet embedded in a porous medium. AIMS Mathematics, 2022, 7(12): 20781-20796. doi: 10.3934/math.20221139
[9]	Abdulaziz H. Alharbi, Ahmed G. Salem . Analytical and numerical investigation of viscous fluid-filled spherical slip cavity in a spherical micropolar droplet. AIMS Mathematics, 2024, 9(6): 15097-15118. doi: 10.3934/math.2024732
[10]	K. Veera Reddy, G. Venkata Ramana Reddy, Ali Akgül, Rabab Jarrar, Hussein Shanak, Jihad Asad . Numerical solution of MHD Casson fluid flow with variable properties across an inclined porous stretching sheet. AIMS Mathematics, 2022, 7(12): 20524-20542. doi: 10.3934/math.20221124

Abstract

Nomenclature

$u\, {\text{and}}\, v$	Fluid velocities corresponding in $x$ and y-directions
$a$	The rate of stretching sheet
$L$	The momentum slip parameter
$V$	The thermal slip factor
$\psi$	Stream function
${r_1}, {r_2}, {r_3}, {r_4}, {r_5}, {r_6}$	Any constant real arbitrary parameter
$\eta$	Similarity variable
$f', \theta, {\phi _1}, {\phi _2}$	Dependent variables
${\alpha _1}$	Momentum slip parameter
${\alpha _2}$	Thermal slip parameter
${\tau _w}$	Shear stress
${q_w}$	Heat flux
${h_{1m}}$ and ${h_{2m}}$	Mass flux
$L{e_1}$ and $L{e_2}$	Lewis numbers
${\phi _1}$	Salt 1 concentration
${\phi _2}$	Salt 2 concentration
${\Lambda _1}$	Thermal buoyancy ratio parameter
${\Lambda _2}$	Solutal buoyancy ratio parameter
$M$	Magnetic field parameter
$T$	Temperature of the fluid
$C$	Concentration of the fluid
${\beta _T}$	Thermal buoyancy
${\beta _{{C_1}}}$	Solutal 1 buoyancy
${\beta _{{C_2}}}$	Solutal 2 buoyancy
${C_{1\infty }}$ and ${C_{2\infty }}$	Solutal 1and 2 at ambient flow
$\nu$	Kinematic viscosity of the fluid
$\Pr$	Prandtl number
$\theta$	Temperature
$g$	Gravity
$\varepsilon$	Small parameter
$\lambda > 0$	Buoyancy assisting flow
$\lambda < 0$	Buoyancy opposing flow
${C_f}$	Local friction factor
$N{u_x}$	Local Nusselt number
$S{h_x}^1$ and $S{h_x}^2$	Local Sherwood numbers
${C_{1w}}$ and ${C_{2w}}$	Concentrations 1 and 2 at the wall
${{Re} _x}$	Reynolds number
$\alpha$	Thermal diffusivity
${D_{{s_1}}}$ and ${D_{{s_2}}}$	Mass 1 and 2 diffusivities

1. Introduction

Lie group transformation or scaling group transformation has elements which are arbitrary close to the identity transformation, which means that the identity transformation is a transformation that changes nothing at all. When two mechanisms contributing to the density diffuse at different rates, convective waves can arise in a steadily stratified liquid. Sophus Lie's ^[1] was developed a differential invariant similarity for the Lie-group. His work was almost forgotten because of the growth of group theory and differential geometry to study differential equations ^[2]. Using Jeffrey and cross fluid models Bhatti et al. ^[3] and Sumer et al. ^[4] analysed importance of momentum with Lie group transformation. Numerical simulation of a thermally enhanced electro-magnetohydrodynamic (EMHD) flow of a heterogeneous micropolar mixture comprising (60%)-ethylene glycol, (40%) water, and copper oxide nanomaterials was studied by Shah et al. ^[5] and MHD heat transfer in thermal slip using semi-infinite domain and Carreau fluid has been discussed by Rehman et al. ^[6]. In recent year, so many researchers made their studies on Lie group transformation ^[7,8,9]. A two-component modeling for non-Newtonian nanofluid slips flow and heat transfer over sheet discussed by Rana et al. ^[10]. They observed that the heat transfer is decreases with increase in Brownian motion and thermophoresis parameters. Amanulla et al. ^[11] studied the numerical analysis of magnetohydrodynamic Williamson nanofluid over an isothermal sphere in the presence of momentum and thermal slip effects. Zhang et al. ^[12] proposed a new multi-view clustering model called Consensus Multi-view Clustering (CMC) model, and it is based on non-negative matrix factorization for predicting the multiple stages of Arithmetic design progression. The proposed CMC model performs multi-view learning to fully capture data features with limited medical images, analyzes similarity relations between different entities, addresses the shortcoming of multi-view fusion that requires manual setting parameters, and further acquires a consensus representation containing shared features and complementary knowledge of multiple view data. Later, Batool et al. ^[13] performed a numerical analysis of heat and mass transfer in micropolar nanofluids flow through a lid-driven cavity by using the finite volume method. They found that a high vortex viscosity parameter produces a weak concentration field and has significant behaviour when the Reynolds number is significant. Maneengam et al. ^[14] analyzed the entropy generation in a two-dimensional Lid-Driven porous container in the presence of obstacles of different shapes and under the influences of buoyancy and Lorentz forces. They concluded that the triangular shape of the baffle is the best in terms of thermal activity and also they observed that increasing the number of lower-wall waves reduces thermal activity. Bhutta et al. ^[15] experimentally investigated the development of novel hybrid two-dimensional and three-dimensional graphene oxide diamond microcomposite polyimide films to ameliorate electrical and thermal conduction. They found that the hybrid fillers with high thermal and low electrical conductivities were responsible for synergistic improvements in the experimental results.Rasool et al. ^[16] performed a numerical investigation of electromagnetohydrodynamic nanofluid flows over a convectively heated Riga pattern positioned horizontally in a Darcy-Forchheimer porous medium. They observed that Darcy-Forchheimer's and Lorentz's forces strengthen the resulting frictional factor at the Riga surface. Rasool et al. ^[17] reported on the Darcy-Forchheimer flow of water conveying multi-walled carbon nanoparticles through a vertical Cleveland z-staggered cavity subject to entropy generation. They found that the velocity of incompressible Darcy-Forchheimer flow at the middle vertical Cleveland Z-staggered cavity declines with a higher Reynolds number.

The phenomena of heat and mass transport associated with entropy production describe the triple diffusion process, which can be found in metallurgy, oceanography, and geophysics. When more than one salt is used in a combined heat and mass transfer, the salts interact and cause coupling effects. In recent years tremendous research activity has been dedicated to developing various fluid models by using the triple diffusive effect: for examples with numerous geometries and multiple body forces, the reader is referred to Arif et al. ^[18], Patil et al. ^[19], Nawaz et al. ^[20] and Manjappa et al. ^[21]. Such techniques used to investigate momentum, thermal and various effects have been discussed by Patil et al. ^[22]; also the roughness of a surface with the nanofluid is studied by Patil et al. ^[23].

The numerical studies incorporated for thermal and momentum slips by Xiong et al. ^[24] and Beg et al. ^[25]. Su et al. ^[26] investigated that increasing flow and convective heat transmission were explored numerically inside circular and parallel-plates micro channels. The slip flow of conductive viscoelastic flow past a permeable radiated shrinking or stretching surface with mass transpiration has been illustrated by Sneha et al. ^[27]. The significance of nanoparticle shape and thermo-hydrodynamic slip constraints on MHD alumina-water nanoliquid flows over a rotating heated disk has been discussed by Sabu et al. ^[28]. Numerous researchers Koriko et al. ^[29], Eleni Seid et al. ^[30], Turkyilmazoglu ^[31], Felipe et al. ^[32], Majeed et al. ^[33] and Zeeshan et al. ^[34] are carried out for different fluid model with various surfaces to analyse the influence of momentum and thermal slips.

No research work has been done using Lie group transformation analysis to analyze the significance of triple diffusive convective flow by considering thermal and momentum slips, as per the authors' knowledge of the literature in most practical uses; however, the components of mass and heat are inextricably linked. The heat and mass diffusion components, on the other hand, are invariably coupled in most real-world applications. Triple diffusion of a fluid has practical applications in many bioengineering and medical sciences. For instance, aqueous suspensions of DNA contain more than two independently diffusing constituents with dissimilar diffusivities and medical drug manufacturing methods etc. However, triple diffusion in the presence of an applied magnetic field has significant applications in metallurgy, solar physics, geophysics, cosmic fluid dynamics, polymer industry, in the motion of the earth's core. This fact prompts us to investigate the combined impact of mass diffusion heat on magnetohydrodynamic free-convective flow in the boundary layer. As far as we can tell, the findings of this work appear to be perfectly consistent with previous research and, due to their simplicity, easily transferable to relevant applications.

2. Mathematical model and formulation

We perform the analysis of the steady triple diffusive hydromagnetic flow of heat and mass transmission of a viscous, electrically conducting fluid. The flow is supposed to be in the x-direction and y-axis perpendicular to it. A uniform magnetic field of strength B is offered perpendicular to the flow direction. It is also considered that all fluid possessions are fixed excluding the impact of the density deviation with temperature and concentration. The surface is retained at a fixed temperature T_w, which is greater than the constant.

For this study, we have considered a triple diffusive free convective MHD incompressible flow with momentum and thermal slip conditions for a fluid model with the help of Lie group transformation analysis. Also, suppose that the flow generation is estimated for a linear widening surface. Additionally, the flow is restricted to the region $y>0$ and equivalent with the plane $y>0$ .

Under the pre-defined assumptions mentioned above, the governing system of continuity, energy, momentum and concentration are expressed as follows:

$\frac{{\partial \bar u}}{{\partial \bar x}} + \frac{{\partial \bar v}}{{\partial \bar y}} = 0 ,$

(1)

$\bar u\frac{{\partial \bar u}}{{\partial \bar x}} + \bar v\frac{{\partial \bar u}}{{\partial \bar y}} = v\frac{{{\partial ^2}\bar u}}{{\partial {{\bar y}^2}}} + \left[ {g{\beta _T}\left( {T - {T_\infty }} \right) + g{\beta _{{C_1}}}\left( {C - {C_1}_\infty } \right) + g{\beta _{{C_2}}}\left( {C - {C_2}_\infty } \right) - \frac{{\sigma {B^2}}}{\rho }\bar u} \right] ,$

(2)

$\bar u\frac{{\partial \bar T}}{{\partial \bar x}} + \bar v\frac{{\partial \bar T}}{{\partial \bar y}} = \alpha \frac{{{\partial ^2}\bar T}}{{\partial {{\bar y}^2}}} ,$

(3)

$\bar u\frac{{\partial {C_1}}}{{\partial \bar x}} + \bar v\frac{{\partial {C_1}}}{{\partial \bar y}} = {D_{{S_1}}}\frac{{{\partial ^2}{C_1}}}{{\partial {{\bar y}^2}}} ,$

(4)

$\bar u\frac{{\partial {C_2}}}{{\partial \bar x}} + \bar v\frac{{\partial {C_2}}}{{\partial \bar y}} = {D_{{S_1}}}\frac{{{\partial ^2}{C_2}}}{{\partial {{\bar y}^2}}} ,$

(5)

Given the following boundary conditions:

$\begin{array}{l} \bar u = a\bar x + L\frac{{\partial u}}{{\partial y}}, \, \bar v = 0, \, T = {T_w} + {D_1}\frac{{\partial T}}{{\partial y}}, {C_1} = {C_{1w}}, {C_2} = {C_{2w}}\, {\text{at}}\, \eta = 0 \hfill \\ \bar u \to {u_\infty }, \, T \to {T_\infty }, \, {C_1} \to C{\, _1}_\infty , \, {C_2} \to C{\, _2}_\infty \, {\text{as}}\, \, \eta \to \infty \hfill \end{array} ,$

(6)

Now, we need to transform the above Eqs (1) to (6) into non-dimensional form. To perform this, the following similarity dimensionless transformations are considered.

$\begin{array}{l} u = \frac{{\bar u}}{{\sqrt {av} }}, v = \frac{{\bar v}}{{\sqrt {av} }}, x = \frac{{\bar x}}{{\sqrt {v/a} }}, y = \frac{{\bar y}}{{\sqrt {v/a} }} \hfill \\ \theta = \frac{{T - {T_\infty }}}{{{T_w} - {T_\infty }}}, {\phi _1} = \frac{{{C_1} - {C_1}_\infty }}{{{C_{1w}} - {C_1}_\infty }}, {\phi _2} = \frac{{{C_2} - {C_2}_\infty }}{{{C_{2w}} - {C_2}_\infty }} \hfill \end{array} ,$

(7)

We substitute the above non-dimensional quantities of (7) in Eqs (1) to (5) along with the boundary conditions of Eq (6) and then complete the simplification results by using the following non-dimensional differential equations :

$\frac{{\partial u}}{{\partial x}} + \frac{{\partial v}}{{\partial y}} = 0 ,$

(8)

$u\frac{{\partial u}}{{\partial x}} + v\frac{{\partial u}}{{\partial y}} = \frac{{{\partial ^2}u}}{{\partial {y^2}}} - \frac{{\sigma {B^2}}}{{\rho a}}u + \lambda \theta + {\Lambda _1}{\phi _1} + {\Lambda _2}{\phi _2} ,$

(9)

$u\frac{{\partial {\phi _1}}}{{\partial x}} + v\frac{{\partial {\phi _1}}}{{\partial y}} = \frac{{{D_{{S_1}}}}}{v}\frac{{{\partial ^2}{\phi _1}}}{{\partial {y^2}}} ,$

(10)

$u \frac{\partial \phi_1}{\partial x}+v \frac{\partial \phi_1}{\partial y} = \frac{D_{S_1}}{\nu } \frac{\partial^2 \phi_1}{\partial y^2},$

(11)

$u\frac{{\partial {\phi _2}}}{{\partial x}} + v\frac{{\partial {\phi _2}}}{{\partial y}} = \frac{{{D_{{S_2}}}}}{v}\frac{{{\partial ^2}{\phi _2}}}{{\partial {y^2}}} ,$

(12)

with the following boundary conditions :

$\begin{array}{l} u = x + L\sqrt {\frac{a}{v}} \frac{{\partial u}}{{\partial y}}, \, v = 0, \, \theta = 1 + {D_1}\sqrt {\frac{a}{v}} \frac{{\partial \theta }}{{\partial y}}, {\phi _1} = 1, {\phi _2} = 1\, {\text{at}}\, \eta = 0 \hfill \\ u \to 0, \theta \to 0, {\phi _1} \to 0, {\phi _2} \to 0\, \, {\text{as}}\, \, \eta \to \infty \hfill \end{array} ,$

(13)

where $\lambda = \frac{{g{\beta _T}\Delta T}}{{a\sqrt {av} }}, {\Lambda _1} = \frac{{g{\beta _{{C_1}}}\Delta {C_1}}}{{a\sqrt {av} }}, {\Lambda _1} = \frac{{g{\beta _{{C_2}}}\Delta {C_2}}}{{a\sqrt {av} }}$ .

3. Scaling transformations

First, the stream function $\psi$ is taken as mentioned below before applying scaling transformations.

$u = \frac{{\partial \psi }}{{\partial y}}\, {\text{and}}\, u = - \frac{{\partial \psi }}{{\partial x}}$

(14)

We introduce the stream functions (14) to (8)–(12) together with the boundary conditions of (13). Obviously, Eq (8) is satisfied by them. And the remaining equations will be converted in terms of the stream function $\psi$ as below period:

$\frac{{\partial \psi }}{{\partial x}}\frac{{{\partial ^2}\psi }}{{\partial x\partial y}} - \frac{{\partial \psi }}{{\partial x}}\frac{{{\partial ^2}\psi }}{{\partial {y^2}}} = \frac{{{\partial ^3}\psi }}{{\partial {y^3}}} - \frac{{\sigma {B^2}}}{{\rho a}}\frac{{\partial \psi }}{{\partial y}} - \lambda \theta + {\Lambda _1}{\phi _1} + {\Lambda _2}{\phi _2} ,$

(15)

$\frac{{\partial \psi }}{{\partial y}}\frac{{\partial \theta }}{{\partial x}} - \frac{{\partial \psi }}{{\partial x}}\frac{{\partial \psi }}{{\partial y}} = \frac{\alpha }{v}\frac{{{\partial ^2}\theta }}{{\partial {y^2}}} ,$

(16)

$\frac{{\partial \psi }}{{\partial y}}\frac{{\partial {\phi _1}}}{{\partial x}} - \frac{{\partial \psi }}{{\partial x}}\frac{{\partial {\phi _1}}}{{\partial y}} = \frac{{{D_{{S_1}}}}}{v}\frac{{{\partial ^2}{\phi _1}}}{{\partial {y^2}}} ,$

(17)

$\frac{{\partial \psi }}{{\partial y}}\frac{{\partial {\phi _2}}}{{\partial x}} - \frac{{\partial \psi }}{{\partial x}}\frac{{\partial {\phi _2}}}{{\partial y}} = \frac{{{D_{{S_2}}}}}{v}\frac{{{\partial ^2}{\phi _2}}}{{\partial {y^2}}} ,$

(18)

with the following boundary conditions:

$\begin{array}{l} \frac{{\partial \psi }}{{\partial y}} = x + L\sqrt {\frac{a}{v}} \frac{{{\partial ^2}\psi }}{{\partial {y^2}}}, \, \frac{{\partial \psi }}{{\partial x}} = 0, \, \theta = 1 + {D_1}\sqrt {\frac{a}{v}\frac{{\partial \theta }}{{\partial y}}} , {\phi _1} = 1, {\phi _2} = 1\, {\text{at}}\, \eta = 0 \hfill \\ \frac{{\partial \psi }}{{\partial y}} \to 0, \, \theta \to 0, {\phi _1} \to 0, {\phi _2} \to 0\, {\text{as}}\, \, \eta \to \infty \hfill \end{array} ,$

(19)

$Also, {e^{\varepsilon \left( {{r_2} - {r_3}} \right)}}\frac{{\partial {\psi ^*}}}{{\partial {y^*}}} = 1, {e^{\varepsilon \left( {{r_1} - {r_3}} \right)}} = 0, \, {e^{\varepsilon {r_4}}}{\theta ^*} = 0$

(20)

Consider $\varepsilon$ as a small parameter of the scaling transformation. Then, the transformation $F$ (a specified set of Lie group analysis) is taken as mentioned below period:

$\begin{array}{l} F:\, {x^*} = x{e^{\varepsilon {r_1}}}, {y^*} = y{e^{\varepsilon {r_2}}}, {\psi ^*} = \psi {e^{\varepsilon {r_3}}} \hfill \\ {\theta ^*} = \theta {e^{\varepsilon {r_4}}}, {\phi _1}^* = {\phi _1}{e^{\varepsilon {r_5}}}, {\phi _2}^* = {\phi _2}{e^{\varepsilon {r_6}}} \hfill \end{array} ,$

(21)

Where, ${r_1}, {r_2}, {r_3}, {r_4}, {r_5}\, {\text{and}}\, {r_6}$ are any real arbitrary constant parameters The point conversion defined through (21) transmuted the co-ordinates $\left({x, y, \psi, {\phi _1}, {\phi _2}} \right)$ into $\left({{x^*}, {y^*}, {\psi ^*}, {\phi _1}^*, {\phi _2}^*} \right)$ . With the help of Lie group analysis mentioned through (15), Eqs (15)–(18) are transformed as

${e^{\varepsilon \left( {{r_1} + {r_2} - {r_3} - {r_5}} \right)}}\left[ {\frac{{\partial {\psi ^*}}}{{\partial {y^*}}}\frac{{{\partial ^2}{\psi ^*}}}{{\partial {x^*}\partial {y^*}}} - \frac{{\partial {\psi ^*}}}{{\partial {x^*}}}\frac{{{\partial ^2}{\psi ^*}}}{{\partial {y^*}^2}}} \right] = {e^{\varepsilon \left( {3{r_2} - {r_3}} \right)}}\frac{{{\partial ^3}{\psi ^*}}}{{\partial {y^*}^3}} - {e^{\varepsilon \left( {{r_2} - {r_3}} \right)}}\frac{{\sigma {B^2}}}{{\rho a}}\frac{{\partial {\psi ^*}}}{{\partial {y^*}}} - \lambda {e^{\varepsilon {r_4}}}\theta$

(22)

${e^{\varepsilon \left( {{r_1} + {r_2} - {r_3} - {r_5}} \right)}}\left[ {\frac{{\partial {\psi ^*}}}{{\partial {y^*}}}\frac{{{\partial ^2}{\theta ^*}}}{{\partial {x^*}\partial {y^*}}} - \frac{{\partial {\psi ^*}}}{{\partial {x^*}}}\frac{{\partial {\theta ^*}}}{{\partial {y^*}}}} \right] = \frac{a}{v}{e^{\varepsilon \left( {2{r_2} - {r_4}} \right)}}\frac{{{\partial ^2}{\theta ^*}}}{{\partial {y^*}^2}}$

(23)

${e^{\varepsilon \left( {{r_1} + {r_2} - {r_3} - {r_5}} \right)}}\left[ {\frac{{\partial {\psi ^*}}}{{\partial {y^*}}}\frac{{\partial {\phi _1}^*}}{{\partial {x^*}}} - \frac{{\partial {\psi ^*}}}{{\partial {x^*}}}\frac{{\partial {\phi _1}^*}}{{\partial {y^*}}}} \right] = \frac{{{D_{{S_1}}}}}{v}{e^{\varepsilon \left( {2{r_2} - {r_5}} \right)}}\frac{{{\partial ^2}{\phi _1}^*}}{{\partial {y^*}^2}}$

(24)

${e^{\varepsilon \left( {{r_1} + {r_2} - {r_3} - {r_6}} \right)}}\left[ {\frac{{\partial {\psi ^*}}}{{\partial {y^*}}}\frac{{\partial {\phi _2}^*}}{{\partial {x^*}}} - \frac{{\partial {\psi ^*}}}{{\partial {x^*}}}\frac{{\partial {\phi _2}^*}}{{\partial {y^*}}}} \right] = \frac{{{D_{{S_2}}}}}{v}{e^{\varepsilon \left( {2{r_2} - {r_6}} \right)}}\frac{{{\partial ^2}{\phi _2}^*}}{{\partial {y^*}^2}} .$

(25)

If the exponent of the above converted equations fulfills the below linear equations, then the set of Eqs (22)–(25) will remain invariant under the folowing transformation

${r_1} = 2{r_2} - 2{r_3} = 3{r_2} - {r_3} = {r_2} - {r_3} = {r_4} = {r_5} = {r_6} ,$

(26)

${r_1} + {r_2} - {r_3} - {r_4} = 2{r_2} - {r_4} ,$

(27)

${r_1} + {r_2} - {r_3} - {r_5} = 2{r_2} - {r_5} ,$

(28)

${r_1} + {r_2} - {r_3} - {r_6} = 2{r_2} - {r_6}$

(29)

By solving the above linear equations (26)–(29) simultaneously, we find

${r_1} = {r_1}, {r_2} = 0, {r_3} = {r_1}, {r_4} = 0, {r_5} = 0, {r_6} = 0 .$

(30)

Replacing the above values of (30) into the scaling transformations given in (21), we get

$F:{x^*} = x{e^{\varepsilon {r_1}}}, {y^*} = y, \, {\psi ^*} = \psi {e^{\varepsilon {r_1}}}, {\theta ^*} = \theta , {\phi _1}^* = {\phi _1}, {\phi _2}^* = {\phi _2} .$

(31)

The Taylor series expansions are

${x^*} - x = x\varepsilon {r_1}, {y^*} - y = 0, \, {\psi ^*} = \psi \varepsilon r, {\theta ^*} - \theta = 0, {\phi _1}^* - {\phi _1} = 0, {\phi _2}^* - {\phi _2} = 0 .$

A simple algebraic expression for the above transformations given by (31) with the help of Taylor series expansion leads to a mono-parametric group of transformations in the form of the characteristic equation given below period

$\frac{{dx}}{{d{r_1}}} = \frac{{dy}}{0} = \frac{{d\psi }}{{\psi {r_1}}} = \frac{{d\theta }}{0} = \frac{{d{\phi _1}}}{0} = \frac{{d{\phi _2}}}{0} .$

(32)

From (32), we can easily obtain new similarity transformations as

$y = \eta , \psi = xf\left( \eta \right), \, \theta = \theta \left( \eta \right), {\phi _1} = {\phi _1}\left( \eta \right), {\phi _2} = {\phi _2}\left( \eta \right) ,$

(33)

where, $\eta$ is the similarity variable and $f, \theta, {\phi _1}, {\phi _2}$ are the dependent variables. Now, we substitute the above quantities specified in (33) into PDEs (16) to (18) then, by applying the boundary conditions mentioned in (19), we attain the set of ordinary differential equations as

$f''' + ff'' - {f'^2} + \lambda \theta + {\Lambda _1}{\phi _1} + {\Lambda _2}{\phi _2} = 0 ,$

(34)

$\theta '' + \Pr \theta ' = 0 ,$

(35)

${\phi _1}^{\prime \prime } + L{e_1}f{\phi _1}^\prime = 0 ,$

(36)

${\phi _2}^{\prime \prime } + L{e_2}f{\phi _2}^\prime = 0 .$

(37)

Together with the boundary situations;

$\begin{array}{l} f = 0, \, f' = 1 + {\alpha _1}f'', \theta = 1 + {\alpha _2}\theta ', {\phi _1} = 1, {\phi _2} = 1\, {\text{at}}\, \eta = 0 \hfill \\ f' \to 0, \theta \to 0, {\phi _1} \to 0, {\phi _2} \to 0\, \, {\text{as}}\, \eta \to \infty \, \hfill \end{array}$

(38)

Here, $\, {\alpha _1} = L\sqrt {\frac{a}{v}}$ is the momentum slip parameter and $\, {\alpha _2} = {D_1}\sqrt {\frac{a}{v}}$ is the thermal slip parameter.

with the following predefined parameters:

$M = \frac{{\sigma {B^2}}}{{\rho a}}, \, \lambda = \frac{{g{\beta _T}\Delta T}}{{a\sqrt {av} }}, {\Lambda _1} = \frac{{g{\beta _{{C_1}}}\Delta {C_1}}}{{a\sqrt {av} }}, {\Lambda _2} = \frac{{g{\beta _{{C_2}}}\Delta {C_2}}}{{a\sqrt {av} }}, \Pr = \frac{v}{\alpha }, L{e_1} = \frac{v}{{{D_{{S_1}}}}}, L{e_2} = \frac{v}{{{D_{{S_2}}}}} .$

The physical measures requiring engineering attention are the local friction factor; local Nusselt number and Sherwood number which are defined as follows:

${C_f} = \frac{{{\tau _w}}}{{\rho {{\left( {ax} \right)}^2}}}, N{u_x} = \frac{{x{q_w}}}{{k\left( {{T_w} - {T_\infty }} \right)}}, N{u_x} = \frac{{x{q_w}}}{{k\left( {{T_w} - {T_\infty }} \right)}}, S{h_x}^1 = \frac{{x{h_m}}}{{{D_{sm}}\left( {{C_{1w}} - {C_{1\infty }}} \right)}}, S{h_x}^2 = \frac{{x{h_m}}}{{{D_{sm}}\left( {{C_{2w}} - {C_{2\infty }}} \right)}} ,$

where, $\tau_{w}, q_{w}, h_{1 m} \& h_{2 m}$ symbolize the shear stress, heat and mass flux change, respectively:

${\tau _w} = - \mu \left( {\frac{{\partial u}}{{\partial y}}} \right), {q_w} = - k\left( {\frac{{\partial T}}{{\partial y}}} \right), {h_{1m}} = - {D_{sm}}\left( {\frac{{\partial {C_1}}}{{\partial y}}} \right), {h_{2m}} = - {D_{sm}}\left( {\frac{{\partial {C_2}}}{{\partial y}}} \right) .$

The dimensionless friction factor, Nusselt number and Sherwood number are defined as

${C_f}{{Re} _x}^{1/2} = f''\left( 0 \right), \, \frac{{N{u_x}}}{{{{{Re} }_x}^{1/2}}} = - \theta '\left( 0 \right), \frac{{S{h_{1x}}}}{{{{{Re} }_x}^{1/2}}} = - {\phi _1}^\prime \left( 0 \right), \frac{{S{h_{2x}}}}{{{{{Re} }_x}^{1/2}}} = - {\phi _2}^\prime \left( 0 \right) ,$

Where, ${{Re} _x} = \sqrt {\frac{a}{{vl}}} x$ .

4. Results and discussions

The collected results show the impact of non-dimensional regulating parameters such as $L{e_1}, L{e_2}, {\Lambda _1}, {\Lambda _2}, {\alpha _1}, {\alpha _2}$ on and $f'\left(\eta \right), \theta \left(\eta \right), {\phi _1}\left(\eta \right), {\phi _2}\left(\eta \right)$ (See –). In the simulations, the pertinent parameter values were $\Pr = 0.72, \, \lambda = - 0.5, L{e_1} = 1, L{e_2} = 1, {\Lambda _1} = 0.5, {\Lambda _2} = 0.5, {\alpha _1} = 0.2, {\alpha _2} = 0.2$ . These variables were kept constant throughout the inquiry, with the exception of the various values given in the associated figures. The results are compared with already divulged results and presents excellent correlation with the Ferdows et al. ^[35] results.

Figure 1. Salt 1 concentration on

$L{e_1}$ .

$\Pr$	Ferdows et al. ^[35]	Present Results
1	0.9547	0.9546
2	1.4715	1.4711
3	1.8691	1.8672

[1]	E. Braverman, B. Karpuz, On oscillation of differential and difference equations with non-monotone delays, Appl. Math. Comput., 218 (2011), 3880–3887. https://doi.org/10.1016/j.amc.2011.09.035 doi: 10.1016/j.amc.2011.09.035
[2]	P. L. Chow, G. Yin, B. Mordukhovich, Topics in stochastic analysis and nonparametric estimation, New York: Springer, 2008. https://doi.org/10.1007/978-0-387-75111-5
[3]	K. Gopalsamy, B. G. Zhang, On a neutral delay logistic equation, Dyn. Stabil. Syst., 2 (1987), 183–195. https://doi.org/10.1080/02681118808806037 doi: 10.1080/02681118808806037
[4]	A. Bellen, N. Guglielmi, A. E. Ruehli, Methods for linear systems of circuit delay differential equations of neutral type, IEEE Trans. Circuits Syst. I Fundam. Theory Appl., 46 (1999), 212–215. https://doi.org/10.1109/81.739268 doi: 10.1109/81.739268
[5]	T. X. Li, S. Frassu, G. Viglialoro, Combining effects ensuring boundedness in an attraction-repulsion chemotaxis model with production and consumption, Z. Angew. Math. Phys., 74 (2023), 109. https://doi.org/10.1007/s00033-023-01976-0 doi: 10.1007/s00033-023-01976-0
[6]	T. X. Li, N. Pintus, G. Viglialoro, Properties of solutions to porous medium problems with different sources and boundary conditions, Z. Angew. Math. Phys., 70 (2019), 1–18. https://doi.org/10.1007/s00033-019-1130-2 doi: 10.1007/s00033-019-1130-2
[7]	Y. Z. Tian, Y. L. Cai, Y. L. Fu, T. X. Li, Oscillation and asymptotic behavior of third-order neutral differential equations with distributed deviating arguments, Adv. Differ. Equ., 2015 (2015), 1–14.
[8]	B. Baculikova, Properties of third-order nonlinear functional differential equations with mixed arguments, Abstr. Appl. Anal., 2011 (2011), 1–15. https://doi.org/10.1155/2011/857860 doi: 10.1155/2011/857860
[9]	H. D. Liu, F. W. Meng, P. C. Liu, Oscillation and asymptotic analysis on a new generalized Emden-Fowler equation, Appl. Math. Comput., 219 (2012), 2739–2748. https://doi.org/10.1016/j.amc.2012.08.106 doi: 10.1016/j.amc.2012.08.106
[10]	H. Ramos, O. Moaaz, A. Muhib, J. Awrejcewicz, More effective results for testing oscillation of non-canonical neutral delay differential equations, Mathematics, 9 (2021), 1–10. https://doi.org/10.3390/math9101114 doi: 10.3390/math9101114
[11]	J. Dzurina, B. Baculikova, Oscillation of third-order quasi-linear advanced differential equations, Differ. Equ. Appl., 4 (2012), 411–421. https://doi.org/10.7153/dea-04-23 doi: 10.7153/dea-04-23
[12]	R. P. Agarwal, M. Bohner, T. X. Li, C. H. Zhang, Oscillation of second-order differential equations with a sublinear neutral term, Carpathian J. Math., 30 (2014), 1–6.
[13]	O. Moaaz, I. Dassios, W. Muhsin, A. Muhib, Oscillation theory for non-linear neutral delay differential equations of third order, Appl. Sci., 10 (2020), 1–16. https://doi.org/10.3390/app10144855 doi: 10.3390/app10144855
[14]	B. Baculikova, J. Dzurina, On the oscillation of odd order advanced differential equations, Bound. Value Probl., 2014 (2014), 1–9. https://doi.org/10.1186/s13661-014-0214-3 doi: 10.1186/s13661-014-0214-3
[15]	B. Almarri, A. H. Ali, A. M. Lopes, O. Bazighifan, Nonlinear differential equations with distributed delay: some new oscillatory solutions, Mathematics, 10 (2022), 1–10. https://doi.org/10.3390/math10060995 doi: 10.3390/math10060995
[16]	O. Bazighifan, O. Moaaz, R. A. El-Nabulsi, A. Muhib, Some new oscillation results for fourth-order neutral differential equations with delay argument, Symmetry, 12 (2020), 1–10. https://doi.org/10.3390/sym12081248 doi: 10.3390/sym12081248
[17]	R. P. Agarwal, M. Bohner, T. X. Li, C. H. Zhang, Even-order half-linear advanced differential equations: improved criteria in oscillatory and asymptotic properties, Appl. Math. Comput., 266 (2015), 481–490. https://doi.org/10.1016/j.amc.2015.05.008 doi: 10.1016/j.amc.2015.05.008
[18]	B. Qaraad, O. Bazighifan, A. H. Ali, A. A. Al-Moneef, A. J. Alqarni, K. Nonlaopon, Oscillation results of third-order differential equations with symmetrical distributed arguments, Symmetry, 14 (2022), 1–14. https://doi.org/10.3390/sym14102038 doi: 10.3390/sym14102038
[19]	C. Cesarano, O. Moaaz, B. Qaraad, N. A. Alshehri, S. K. Elagan, M. Zakarya, New results for oscillation of solutions of odd-order neutral differential equations, Symmetry, 13 (2021), 1–12. https://doi.org/10.3390/sym13061095 doi: 10.3390/sym13061095
[20]	O. Moaaz, R. A. El-Nabulsi, A. Muhib, S. K. Elagan, M. Zakarya, New improved results for oscillation of fourth-order neutral differential equations, Mathematics, 9 (2021), 1–12. https://doi.org/10.3390/math9192388 doi: 10.3390/math9192388
[21]	T. X. Li, Y. V. Rogovchenko, On asymptotic behavior of solutions to higher-order sublinear Emden-Fowler delay differential equations, Appl. Math. Lett., 67 (2017), 53–59. https://doi.org/10.1016/j.aml.2016.11.007 doi: 10.1016/j.aml.2016.11.007
[22]	O. Moaaz, C. Cesarano, A. Muhib, Some new oscillation results for fourth-order neutral differential equations, Eur. J. Pure Appl. Math., 13 (2020), 185–199. https://doi.org/10.29020/nybg.ejpam.v13i2.3654 doi: 10.29020/nybg.ejpam.v13i2.3654
[23]	A. Almutairi, A. H. Ali, O. Bazighifan, L. F. Iambor, Oscillatory properties of fourth-order advanced differential equations, Mathematics, 11 (2023), 1–11. https://doi.org/10.3390/math11061391 doi: 10.3390/math11061391
[24]	M. F. Aktas, A. Tiryaki, A. Zafer, Oscillation criteria for third-order nonlinear functional differential equations, Appl. Math. Lett., 23 (2010), 756–762. https://doi.org/10.1016/j.aml.2010.03.003 doi: 10.1016/j.aml.2010.03.003
[25]	G. E. Chatzarakis, T. X. Li, Oscillation criteria for delay and advanced differential equations with nonmonotone arguments, Complexity, 2018 (2018), 1–18. https://doi.org/10.1155/2018/8237634 doi: 10.1155/2018/8237634
[26]	S. S. Santra, K. M. Khedher, O. Moaaz, A. Muhib, S. W. Yao, Second-order impulsive delay differential systems: necessary and sufficient conditions for oscillatory or asymptotic behavior, Symmetry, 13 (2021), 1–12. https://doi.org/10.3390/sym13040722 doi: 10.3390/sym13040722
[27]	J. Dzurina, I. Jadlovska, A note on oscillation of second-order delay differential equations, Appl. Math. Lett., 69 (2017), 126–132. https://doi.org/10.1016/j.aml.2017.02.003 doi: 10.1016/j.aml.2017.02.003
[28]	J. Dzurina, I. Jadlovska, A sharp oscillation result for second-order half-linear noncanonical delay differential equations, Electron. J. Qual. Theory Differ. Equ., 2020 (2020), 1–14. https://doi.org/10.14232/ejqtde.2020.1.46 doi: 10.14232/ejqtde.2020.1.46
[29]	J. Dzurina, S. R. Grace, I. Jadlovska, T. X. Li, Oscillation criteria for second-order Emden-Fowler delay differential equations with a sublinear neutral term, Math. Nachr., 293 (2020), 910–922. https://doi.org/10.1002/mana.201800196 doi: 10.1002/mana.201800196
[30]	T. X. Li, Y. V. Rogovchenko, Oscillation of second-order neutral differential equations, Math. Nachr., 288 (2015), 1150–1162. https://doi.org/10.1002/mana.201300029 doi: 10.1002/mana.201300029
[31]	T. X. Li, Y. V. Rogovchenko, Oscillation criteria for second-order superlinear Emden-Fowler neutral differential equations, Monatsh. Math., 184 (2017), 489–500. https://doi.org/10.1007/s00605-017-1039-9 doi: 10.1007/s00605-017-1039-9
[32]	G. E. Chatzarakis, I. Jadlovska, Improved oscillation results for second-order half-linear delay differential equations, Hacet. J. Math. Stat., 48 (2019), 170–179. https://doi.org/10.15672/HJMS.2017.522 doi: 10.15672/HJMS.2017.522
[33]	S. R. Grace, J. Dzurina, I. Jadlovska, T. X. Li, An improved approach for studying oscillation of second-order neutral delay differential equations, J. Inequal. Appl., 2018 (2018), 1–13. https://doi.org/10.1186/s13660-018-1767-y doi: 10.1186/s13660-018-1767-y
[34]	I. Jadlovska, New criteria for sharp oscillation of second-order neutral delay differential equations, Mathematics, 9 (2021), 1–23. https://doi.org/10.3390/math9172089 doi: 10.3390/math9172089
[35]	O. Moaaz, M. Anis, D. Baleanu, A. Muhib, More effective criteria for oscillation of second-order differential equations with neutral arguments, Mathematics, 8 (2020), 1–13. https://doi.org/10.3390/math8060986 doi: 10.3390/math8060986
[36]	I. Jadlovska, Oscillation criteria of Kneser-type for second-order half-linear advanced differential equations, Appl. Math. Lett., 106 (2020), 106354. https://doi.org/10.1016/j.aml.2020.106354 doi: 10.1016/j.aml.2020.106354
[37]	G. E. Chatzarakis, J. Dzurina, I. Jadlovska, New oscillation criteria for second-order half-linear advanced differential equations, Appl. Math. Comput., 347 (2019), 404–416. https://doi.org/10.1016/j.amc.2018.10.091 doi: 10.1016/j.amc.2018.10.091
[38]	Y. S. Qi, J. W. Yu, Oscillation of second order nonlinear mixed neutral differential equations with distributed deviating arguments, Bull. Malays. Math. Sci. Soc., 38 (2015), 543–560. https://doi.org/10.1007/s40840-014-0035-7 doi: 10.1007/s40840-014-0035-7
[39]	C. H. Zhang, B. Baculikova, J. Dzurina, T. X. Li, Osillation results for second-order mixed neutral differential equations with distributed deviating arguments, Math. Slovaca, 66 (2016), 615–626. https://doi.org/10.1515/ms-2015-0165 doi: 10.1515/ms-2015-0165
[40]	H. W. Shi, Y. Z. Bai, Oscillatory behavior of a second order nonlinear advanced differential equation with mixed neutral terms, Adv. Differ. Equ., 2019 (2019), 468. https://doi.org/10.1186/s13662-019-2393-9 doi: 10.1186/s13662-019-2393-9
[41]	R. Arul, V. S. Shobha, Oscillation of second order nonlinear neutral differential equations with mixed neutral term, J. Appl. Math. Phys., 3 (2015), 1080–1089. https://doi.org/10.4236/jamp.2015.39134 doi: 10.4236/jamp.2015.39134
[42]	E. Thandapani, S. Padmavathi, P. Pinelas, Oscillation criteria for even-order nonlinear neutral differential equations of mixed type, Bull. Math. Anal. Appl., 6 (2014), 9–22.
[43]	O. Moaaz, A. Muhib, S. S. Santra, An oscillation test for solutions of second-order neutral differential equations of mixed type, Mathematics, 9 (2021), 1634. https://doi.org/10.3390/math9141634 doi: 10.3390/math9141634

AIMS Mathematics

Oscillation criteria for mixed neutral differential equations

Related Papers:

Abstract

1. Introduction

2. Mathematical model and formulation

3. Scaling transformations

4. Results and discussions

5. Concluding remarks

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog