Processing math: 100%
Research article Topical Sections

Investigation of novel turbulator with and without twisted configuration under turbulent forced convection of a CuO/water nanofluid flow inside a parabolic trough solar collector

  • Received: 31 July 2022 Revised: 07 December 2022 Accepted: 13 December 2022 Published: 17 January 2023
  • In this study, we numerically investigated the hydrothermal performance of a parabolic trough solar collector system in which nanofluids are used to transfer thermal energy. The single-phase model has been used to evaluate the respective influences of the spherical shape of nanoparticles with a volume fraction of (φ = 3%), Reynolds number varying between 50,000 ≤ Re ≤ 250,000 and the insertion of a turbulator with and without a twisted configuration on the hydrothermal characteristics created by the turbulent forced convection of a CuO/water nanofluid. The shaped turbulator (+) inserted in the absorber tube had a length turb_L = 2.4 m, a height turb_H = 40 mm and a width turb_t = 2 mm. In the second configuration, the considered turbulator was twisted (N_twist = 5, 10 and 15 twists). The turbulator was positioned at 0.6 m from the inlet of the tube and 1 m from the outlet of the collector. The studied performances included the heat transfer characteristics, pressure drop, friction factor, thermal efficiency, temperature and velocity distribution of the outlet field. The most significant contribution of this study is the proposal of the best parameters to increase the thermal and hydraulic efficiency of parabolic troughs by adding a new turbulator with the considered twists.

    Citation: Omar Ouabouch, Imad Ait Laasri, Mounir Kriraa, Mohamed Lamsaadi. Investigation of novel turbulator with and without twisted configuration under turbulent forced convection of a CuO/water nanofluid flow inside a parabolic trough solar collector[J]. AIMS Materials Science, 2023, 10(1): 112-138. doi: 10.3934/matersci.2023007

    Related Papers:

    [1] Yi Wei . The Riccati-Bernoulli subsidiary ordinary differential equation method to the coupled Higgs field equation. Electronic Research Archive, 2023, 31(11): 6790-6802. doi: 10.3934/era.2023342
    [2] Shao-Xia Qiao, Li-Jun Du . Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, 2021, 29(3): 2269-2291. doi: 10.3934/era.2020116
    [3] Qingming Hao, Wei Chen, Zhigang Pan, Chao Zhu, Yanhua Wang . Steady-state bifurcation and regularity of nonlinear Burgers equation with mean value constraint. Electronic Research Archive, 2025, 33(5): 2972-2988. doi: 10.3934/era.2025130
    [4] Yaning Li, Yuting Yang . The critical exponents for a semilinear fractional pseudo-parabolic equation with nonlinear memory in a bounded domain. Electronic Research Archive, 2023, 31(5): 2555-2567. doi: 10.3934/era.2023129
    [5] Mohammad Kafini, Maher Noor . Delayed wave equation with logarithmic variable-exponent nonlinearity. Electronic Research Archive, 2023, 31(5): 2974-2993. doi: 10.3934/era.2023150
    [6] Abdelhadi Safsaf, Suleman Alfalqi, Ahmed Bchatnia, Abderrahmane Beniani . Blow-up dynamics in nonlinear coupled wave equations with fractional damping and external source. Electronic Research Archive, 2024, 32(10): 5738-5751. doi: 10.3934/era.2024265
    [7] Li Tian, Ziqiang Wang, Junying Cao . A high-order numerical scheme for right Caputo fractional differential equations with uniform accuracy. Electronic Research Archive, 2022, 30(10): 3825-3854. doi: 10.3934/era.2022195
    [8] Tao Zhang, Tingzhi Cheng . A priori estimates of solutions to nonlinear fractional Laplacian equation. Electronic Research Archive, 2023, 31(2): 1119-1133. doi: 10.3934/era.2023056
    [9] Janarthanan Ramadoss, Asma Alharbi, Karthikeyan Rajagopal, Salah Boulaaras . A fractional-order discrete memristor neuron model: Nodal and network dynamics. Electronic Research Archive, 2022, 30(11): 3977-3992. doi: 10.3934/era.2022202
    [10] Adrian Nicolae Branga . Some conditions for the existence and uniqueness of monotonic and positive solutions for nonlinear systems of ordinary differential equations. Electronic Research Archive, 2022, 30(6): 1999-2017. doi: 10.3934/era.2022101
  • In this study, we numerically investigated the hydrothermal performance of a parabolic trough solar collector system in which nanofluids are used to transfer thermal energy. The single-phase model has been used to evaluate the respective influences of the spherical shape of nanoparticles with a volume fraction of (φ = 3%), Reynolds number varying between 50,000 ≤ Re ≤ 250,000 and the insertion of a turbulator with and without a twisted configuration on the hydrothermal characteristics created by the turbulent forced convection of a CuO/water nanofluid. The shaped turbulator (+) inserted in the absorber tube had a length turb_L = 2.4 m, a height turb_H = 40 mm and a width turb_t = 2 mm. In the second configuration, the considered turbulator was twisted (N_twist = 5, 10 and 15 twists). The turbulator was positioned at 0.6 m from the inlet of the tube and 1 m from the outlet of the collector. The studied performances included the heat transfer characteristics, pressure drop, friction factor, thermal efficiency, temperature and velocity distribution of the outlet field. The most significant contribution of this study is the proposal of the best parameters to increase the thermal and hydraulic efficiency of parabolic troughs by adding a new turbulator with the considered twists.



    The introduction of fractional derivatives has significantly expanded the theoretical framework of calculus, enabling the resolution of a broader spectrum of mathematical and applied problems in fields such as physics [1], chemical kinetics systems [2], and control theory [3]. Fractional derivatives play a crucial role in multiple aspects of wave equations. In real-world physical environments, many media exhibit non-uniform properties. The fractional derivative wave equation accounts for these long-term correlations, providing a more precise description of wave behavior and uncovering patterns that traditional integer-order wave equations cannot capture. Therefore, it is of great significance to study the fractional wave equation. With the development of differential equation theory until now, the definition of fractional derivatives is varied. The well-known fractional derivatives include the Caputo derivative [4,5], Caputo-Fabrizio derivative [6], Riemann-Liouville derivative [7], Grünwald-Letnikov derivative [8], Atangana-Baleanu derivative [9], and conformable derivative [10], etc. To date, there have been a large number of methods to analyze and solve fractional nonlinear partial differential equations, such as the F-expansion method [11], generalized finite difference method [12], exp(ϕ(ξ))-expansion and the MAE methods [13], extended (G/G)-expansion method [14], extended trial equation and modified auxiliary equation method [15], the dynamics approach [16], modified generalized Riccati equation mapping method [17], and Hirota method [18].

    The Kadomtsev-Petviashvili (KP) equation [19] was initially introduced by Kadomtsev and Petviashvili in 1970. The equation describes the nonlinear wave phenomena, and can be regarded as a two-dimensional generalization of the Korteweg-de Vries (KdV) equation. The KdV-Sawada-Kotera-Ramani (KdV-SK-R) equation [20], which is a linear combination of the KdV equation and the Sawada-Kotera equation [21,22], can reduce to either the KdV or the Sawada-Kotera equation by eliminating specific terms. In 2022, Guo [23] incorporated a second-order derivative term into the KdV-Sawada-Kotera-Ramani equation, and thereby formulated a new nonlinear partial differential equation known as the (2+1)-dimensional Kadomtsev-Petviashvili-Sawada-Kotera-Ramani (KP-SK-R) equation. This development enables the previous equation to facilitate the description of more complex wave behaviors. Shortly thereafter, Wazwaz et al. [24] extended the (2+1)-dimensional KP-SK-R equation to its (3+1)-dimensional counterpart by introducing additional second-order derivative terms. Recently, Şenol et al. [25] extended the (3+1)-dimensional KP-SK-R equation by introducing the conformable fractional derivative.

    Regarding the KP-SK-R equation, a significant amount of research has been devoted to its analysis. Gu et al. [26] derived soliton, breather, and lump solutions for the (2+1)-dimensional KP-SK-R equation using Hirota's bilinear form and the Duffing chaotic system. Althobaiti et al. [27] investigated the Lax-integrable three-dimensional sixth-order KP-SK-R equation via the Kudryashov method, further extending this equation to an n-dimensional framework further. Ma et al. [28] obtained hybrid solutions for the (3+1)-dimensional KP-SK-R equation through the Hirota bilinear method. Qasim et al. [29] derived exact solutions for the (3+1)-dimensional KP-SK-R equation using the Jacobi elliptic function expansion technique. Şenol et al. [25] formulated the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative by applying the exp(ϕ(ξ))-expansion and Kudryashov R function methods.

    In an effort to further explore the properties of the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative, this paper will continue to study the equation. The polynomial complete discrimination system method is capable of converting non-linear partial differential equations into polynomial equations. Subsequently, it solves these polynomial equations by leveraging the operational rules and discriminant system of polynomials. The solution process features a clear orientation and well-defined steps, rendering it highly amenable to calculation and implementation. Upon examination, it is evident that the (3+1)-dimensional KP-SK-R equation can be transformed into polynomial form via appropriate variable substitution and traveling wave transformation. Therefore, we propose to explore some new traveling wave solutions using the polynomial trial method in conjunction with the polynomial complete discrimination system method [30,31]. Moreover, by introducing external perturbations, we aim to observe the chaotic dynamic behavior inherent in this equation by leveraging the theory of plane dynamic systems [32].

    The following work is divided into four parts: Section 2 gives the mathematical analysis for the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative. Section 3 analyzes chaotic behavior of the equation qualitatively. Section 4 gives some new exact solutions, and part of their graphs using trial method and polynomial complete discrimination system. Section 5 summarizes the whole work.

    Definition 1 ([10]). The αorder conformable derivative of the function f:[0,+)R is defined as

    Dαt(f)(t)=limϵ0f(t+ϵt1α)f(t)ϵ. (2.1)

    If f is α-differentiable for an interval (0,a),a>0, and limt0+Dαt(f)(t) exists, then Dαt(f)(0)=limt0+Dαt(f)(t).

    It has be verified that the fractional-order definition (1) satisfies the following properties:

    Property 1. [10]

    Dαt(tk)=ktkα;

    Dαt(af1+bf2)=aDαt(f1)+bDαt(f2);

    Dαt(f1f2)=f2Dαt(f1)f1Dαt(f2)f22;

    Dαt(f1f2)=f2Dαt(f1)+f1Dαt(f2);

    Dαt(f)=(t1α)df(t)dt

    If f is a constant, then Dαt(f)0.

    Definition ([33]). For a function f with n variables zi,i=1,,n, the α-order conformable derivative of f respect to zi is defined as

    αf(z1,,zn)αzi=limϵ0f(z1,,zi1,zi+ϵz1αi,zi+1,,zn,)f(z1,,zn)ϵ. (2.2)

    According to Definitions 1 and 2, the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative [25] researched in this work is given as follows:

    Dαt(Φx)+κ(3Φ2+Φxx)xx+(15Φ3+15ΦΦxx+Φxxxx)xx+μ1Φyy+μ2Φzz+δ1Φxy+δ2Φyz+δ3Φxz=0, (2.3)

    where x,y,z,t are the wave variables, t>0, κ,μ1,μ2,δ1, δ2, and δ3 are arbitrary constants, and Dαt(Φx) represents the α-order derivative of Φx.

    Given the transformation

    Φ(x,y,z,t)=ϕ(ζ),ζ=ax+by+cz+vtαα, (2.4)

    where v represents the wave velocity, a,b, and c represent the wave width in the x,y, and z directions respectively. Therefore, due to Property 1, we have

    {Dαt(Φx)=avϕ,(3Φ2+Φxx)xx=a2(3ϕ2+a2ϕ),(15Φ3+15ΦΦxx+Φxxxx)xx=a2(15ϕ3+15a2ϕϕ+a4ϕ(4)),μ1Φyy+μ2Φzz=(μ1b2+μ2c2)ϕ,δ1Φxy+δ2Φyz+δ3Φxz=(δ1ab+δ2bc+δ3ac)ϕ. (2.5)

    Subsequently, Eq (2.5) is substituted into Eq (2.3), leading to the ordinary differential form

    κa2(3ϕ2+a2ϕ)+a2(15ϕ3+15a2ϕϕ+a4ϕ(4))+ρϕ=0, (2.6)

    where ρ=av+μ1b2+μ2c2+δ1ab+δ2bc+δ3ac. If Eq (2.6) is integrated twice with respect to ζ and the integration constant is assumed as zero, then we have

    a6ϕ(4)+κa4ϕ+15a4ϕϕ+15a2ϕ3+3κa2ϕ2+ρϕ=0. (2.7)

    By using the polynomial trial method, we take the general trial equation

    ϕ=Nj=0sjϕj, (2.8)

    where sj,j=1,,N, are coefficients to be determined. On the basis of the homogeneous balance principle, it results that N=2. That is to say, ϕ fulfills the trial equation

    ϕ=s2ϕ2+s1ϕ+s0. (2.9)

    Integrating Eq (2.9) with respect to ζ, we have

    (ϕ)2=2s23ϕ3+s1ϕ2+2s0ϕ+d0, (2.10)

    where d0 is the integral constant. From Eqs (2.9) and (2.10), we have

    ϕ(4)=103s22ϕ3+5s1s2ϕ2+(6s2s0+s21)ϕ+s1s0+2s2d0. (2.11)

    New algebraic equations about ϕ can be obtained by substituting both Eqs (2.9) and (2.11) into Eq (2.7) below:

    r3ϕ3+r2ϕ2+r1ϕ+r0=0, (2.12)

    where r3=103a6s22+15a4s2+15a2, r2=5a6s1s2+κa4s2+15a4s1+3κa2, r1=6a6s2s0+a6s21+κa4s1+15a4s0+ρ, and r0=a6s1s0+2a6s2d0+κa4s0. The value of d0,sj,j=1,2,3, can be calculated by solving the equations ri=0,i=1,2,3, which are given as

    s2=3a2,s1 is an arbitrary real number,s0=a6s21+κa4s1+κa4s13a4,d0=(a2s1+κ)(a6s21+κa4s1+κa4s1)18a4, (2.13)

    or

    s2=32a2,s1=κ5a2,s0=4a2κ2+5ρ15a4,d0=4κ(4a2κ2+5ρ)225a4. (2.14)

    It can be seen that s2 is always negative.

    In order to investigate the dynamics of the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative, we analyze Eq (2.9) by employing the theory of plane dynamic systems [32].

    When an autonomous system is subjected to external periodic perturbations, its inherent nonlinear characteristics are activated, leading to a more complex phase space structure and heightened sensitivity to initial conditions. External periodic disturbances are typically modeled using nonlinear periodic functions. Introducing such an excitation function effectively adds an additional degree of freedom to the system, thereby increasing the dimensionality of the phase space. This results in a greater diversity and complexity of possible states and trajectories.

    Let h(ζ)=Q0sin(ωζ) be a periodic external disturbance, where Q0 and ω stands for the amplitude and frequency respectively. When the periodic external disturbance is added to Eq (2.9), this equation is turned into the following form:

    ϕ=s2ϕ2+s1ϕ+s0+Q0sin(ωζ). (3.1)

    Denote g=ϕ and z=ωζ. Then, Eq (3.1) can be transformed into the three-dimensional dynamical system

    {dϕdζ=g,dgdζ=s2ϕ2+s1ϕ+s0+Q0sinz,dzdζ=ω. (3.2)

    By taking the appropriate parameters, some two- and three-dimensional phase portraits of the periodically excited system (3.2) are shown to explain its chaotic behavior. Further, the sensitivity analysis is carried out simultaneously. Several sets of initial values are taken to observe the response of Eq (3.2), see Figure 1.

    Figure 1.  Chaotic behavior of Eq (3.2) with s2=1.1,s1=5.5,s0=1.1, Q0=2.1, ω=1.55, and (ϕ0,g0)=(0.01,0) (red curve) and (0.04,0) (green curve).

    A bifurcation diagram provides a clear visualization of the critical points at which a system abruptly transitions from one state to another as parameters vary. Consequently, it delineates the stable and unstable regions of system parameters, facilitating the observation of various dynamic phenomena and underlying principles. Let s2=14,s1=10,s0=5, ω=2 and the values range of Q0 be from 0 to 5. With the aid of MATLAB software, the bifurcation diagrams in different two directions are given in Figure 2.

    Figure 2.  Bifurcation diagram of Eq (3.2) as Q0 varying with s2=1.1,s1=4.5,s0=1.1, and ω=1.2.

    The Lyapunov exponent quantitatively measures the rate of divergence or convergence of two nearby trajectories in phase space over time. In multidimensional systems, different methods may yield distinct Lyapunov exponents. Positive and negative Lyapunov exponents respectively signify the directions of unstable and stable orbits. The presence of any positive Lyapunov exponent indicates chaotic behavior within the system. In this subsection, the Lyapunov exponents of Eq (3.2) with a set of fixed parameters are calculated by using Wolf algorithm, and its Kaplan-Yorke dimension DKY is obtained. After that, the graph of parameter-dependent Lyapunov exponents are provided when Q0 runs the range of [0,5].

    Let s2=5,s1=15,s0=1, ω=0.9, and Q0=2. From Figure 3, it can be seen that the Lyapunov exponents ξ1,ξ2, and ξ3 converges to 0.032, 0.019, and 0.05, respectively.

    Figure 3.  Dynamics of Lyapunov exponents of Eq (3.2) with s2=5,s1=15,s0=1, and ω=0.9.

    Theorem 1. Under the conditions of s2=5,s1=15,s0=1, ω=0.9, and Q0=2, the Kaplan-Yorke dimension DKY of Eq (3.2) is bigger than 2.

    Proof. The Kaplan-Yorke dimension DKY of Eq (3.2) is computed as follows:

    DKY=2+Σ2i=1ξi|ξ3|3.02. (3.3)

    Therefore, it can be known from Eq (3.3) that the conclusion holds.

    From the above, we can see that both the maximum Lyapunov exponent ξ1>0 and the Kaplan-Yorke dimension DKY>2 verify the chaotic state of Eq (3.2). In Figure 4, it can be seen that as Q0 increases in [0,5], the Lyapunov exponents Le1,Le2, and Le3 present different degrees of chaos, and the maximum Lyapunov exponent Le1 is almost always greater than 0.

    Figure 4.  Lyapunov exponents spectrum of Eq (3.2) depending on Q0 as varying with s2=5,s1=15,s0=1,ω=0.9, and Q0=2.

    In this section, we solve the function Φ for Eq (2.3) by using the third-order polynomial complete discrimination system method [31].

    When applying the substitutions

    u=(23s2)13ϕ, q2=s1(23s2)23, q1=2s0(23s2)13, q0=d0, (4.1)

    Eq (2.10) takes the new form

    (u)2=u3+q2u2+q1u+q0. (4.2)

    The integral form of Eq (4.2) is

    ±(23s2)13(ζζ0)=duu3+q2u2+q1u+q0, (4.3)

    where ζ0 is the integration constant.

    Let F(u)=u3+q2u2+q1u+q0. The complete discrimination system of F(u) is given as follows:

    {Δ=27(2q3227+q0q1q03)24(q1q223)3,H=q1q223. (4.4)

    In the application of the polynomial discriminant system method, the signs of Δ and H can reveal certain characteristics of the roots of the equation. In a cubic equation, when Δ=0, the equation has multiple roots. Specifically, when Δ=0 and H<0, the cubic equation has one real root and a double root. When Δ=0 and H=0, the equation has three equal real roots. When Δ>0 and H<0, the equation has three distinct real roots. When Δ<0, the equation has no real roots. Therefore, the solution of the integral equation (4.3) will be discussed in four cases.

    Theorem 2. Let α1=q2q19q02(q223q1) and α2=q2+q2q19q0q223q1. If Δ=0 and H<0, then the following two conclusions hold:

    (i) When α1>α2, Eq (2.3) has solutions involving a hyperbolic tangent and a hyperbolic cotangent as follows:

    Φ1(x,y,z,t)=(23s2)13{(q23(q2q19q0)2(q223q1))tanh2[12(q23(q2q19q0)2(q223q1))12(23s2)13(ax+by+cz+vtααζ0)]q2+q2q19q0q223q1}; (4.5)
    Φ2(x,y,z,t)=(23s2)13{(q23(q2q19q0)2(q223q1))coth2[12(q23(q2q19q0)2(q223q1))12(23s2)13(ax+by+cz+vtααζ0)]q2+q2q19q0q223q1}, (4.6)

    where α(0,1), a,b,c,v are arbitrary real numbers, s2,s1,s0,d0 are taken as the formulas (2.13) or (2.14), q2,q1,q0 are taken as formulas (4.1) and satisfy that α1>α2.

    (ii) When α1<α2, Eq (2.3) has a tangent solution

    Φ3(x,y,z,t)=(23s2)13{(3(q2q19q0)2(q223q1)q2)tan2[12(3(q2q19q0)2(q223q1)q2)12(23s2)13(ax+by+cz+vtααζ0)]q2+q2q19q0q223q1}, (4.7)

    where the restrictions of α,a,b,c,v,s2,s1,s0,d0 are the same as those of the solution Φ1(x,y,z,t), while q2,q1,q0 are taken as formulas (4.1) and satisfy that α1<α2.

    Proof. Suppose that Δ=0, and H<0. Thereby, we have F(ψ)=(ψα1)2(ψα2). According to the cubic derivation formula, there are α1=q2q19q02(q223q1) and α2=q2+q2q19q0q223q1. If α1>α2, then Eq (4.3) has two solutions by utilizing the third-order polynomial complete discrimination system method as follows:

    u1=(α1α2)tanh2[(α1α2)122(ζζ0)]+α2

    and

    u2=(α1α2)coth2[(α1α2)122(ζζ0)]+α2.

    If α1<α2, then Eq (4.3) has a solution by utilizing the third-order polynomial complete discrimination system method in the following:

    u=(α2α1)tan2[(α1α2)122(ζζ0)]+α2.

    Replacing u,ζ,α1, and α2 in the above three equations by their specific expressions, the solutions of Eq (2.3) are obtained as the forms (4.5)–(4.7), respectively.

    Theorem 3. If Δ=0 and H=0, then Eq (2.3) yields a rational solution given by

    Φ4(x,y,z,t)=4(23s2)23(ax+by+cz+vtααζ0)2q23, (4.8)

    where the restrictions of α,a,b,c,v,s2,s1,s0,d0 are the same as those of the solution Φ1(x,y,z,t), while q2,q1, and q0 are taken as formulas (4.1) and satisfy that α1=q23.

    Proof. Suppose that Δ=0 and H=0. Therefore, F(ψ)=0 has three same real roots, i.e., F(u)=(uα1)3. Because of the cubic derivation formula, there is α1=q23. By using the third-order polynomial complete discrimination system method, Eq (4.3) has a solution as follows:

    u=4(ζζ0)2+α1.

    Substituting u,ζ, and α1 into the above equation with their specific expressions, the solution of Eq (2.3) is obtained as the form (4.8).

    Theorem 4. Let

    {α1=13(q22Acosθ3),α2=13[q2+A(cosθ33sinθ3)],α3=13[q2+A(cosθ3+3sinθ3)], (4.9)

    where θ=arccos2q329q2q1+27q02A32and A=3H=q223q1. If Δ>0 and H<0, then the following two conclusions hold:

    (i) When α1<u=(23s2)13Φ<α3, Eq (2.3) possesses a Jacobian elliptic solution of the form

    Φ5(x,y,z,t)=(23s2)13[13(q22q223q1cosθ3)+(Acosθ333Asinθ3)sn2(12(Acosθ3+33Asinθ3)12(23s2)13(ax+by+cz+vtααζ0),m)], (4.10)

    where sn is the Jacobian elliptic sine function, m2=3Acosθ3Asinθ33Acosθ3+Asinθ3, and the restrictions of α,a,b,c,v,s2,s1,s0,d0 are the same as those of the solution Φ1(x,y,z,t), while q2,q1, and q0 are taken as formulas (4.1) and satisfy that α1<u<α3;

    (ii) When u=(23s2)13Φ>α3, Eq (2.3) possesses a Jacobian elliptic solution of the form

    Φ6(x,y,z,t)=(23s2)13C1+3Asinθ3(C13Asinθ3)sn2[C2(ax+by+cz+vtααζ0),m]3cn2[C2(ax+by+cz+vtααζ0),m], (4.11)

    where cn is the Jacobian elliptic consine function, C1=q2+Acosθ3,C2=12(23s2)13(Acosθ3+33Asinθ3)12, m2=3Acosθ3Asinθ33Acosθ3+Asinθ3, and the restrictions of α,a,b,c,v,s2,s1,s0,d0 are the same as those of the solution Φ1(x,y,z,t), while q2,q1, and q0 are taken as formulas (4.1) and satisfy that u>α3.

    Proof. Suppose that Δ>0 and H<0. Then, F(u)=0 has three real different roots, that is, F(u)=(uα1)(uα2)(uα3). By the cubic derivation formula, we can obtain the specific form of α1,α2, and α3 in Eq (4.9). Furthermore, α1<α2<α3 because of θ(0,π). If α1<u=(23s2)13ϕ<α3, with the aid of the third-order polynomial complete discrimination system method, Eq (4.3) has the following solution:

    u=[α1+(α2α3)sn2((α3α1)122(ζζ0),m)]

    If u=(23s2)13ϕ>α3, then Eq (4.3) has the following solution:

    u=α3α2sn2((α3α1)122(ζζ0),m)cn2((α3α1)122(ζζ0),m).

    Substituting the specific expressions of u,ϕ,ζ,α1,α2, and α3 into the above two equations, two solutions of Eq (2.3) are obtained as the forms (4.10) and (4.11), respectively.

    Theorem 5. Let

    {α1=13(q2+β),β=3γ1+3γ2,γ1=Aq2+32[(q2q19q0)+(q2q19q0)24A(q213q2q0)],γ2=Aq2+32[(q2q19q0)(q2q19q0)24A(q213q2q0)]. (4.12)

    If Δ<0 and u=(23s2)13Φ>α1, then Eq (2.3) has a Jacobian elliptic solution as follows:

    Φ7(x,y,z,t)=(23s2)13[13(q2+β)+229(2q22+q2β+β2)1+cn((29(2q22+q2β+β2))14(23s2)13(ax+by+cz+vtααζ0),m)29(2q22+q2β+β2)], (4.13)

    where cn is the Jacobian elliptic cosine function, m2=12+3(4q2+β)8(2q22+q2β+β2), and the restrictions of α,a,b,c,v,s2,s1,s0,d0 are the same as those of the solution Φ1(x,y,z,t), while q2,q1, and q0 are taken as formulas (4.1) and satisfy that u>α1.

    Proof. Suppose that Δ<0. Then, F(u)=0 has only one real root, that is, F(u)=(uα1)(u2+α2u+α3). By the cubic derivation formula, we have

    α1=13(q2+β),α2=23(q2+12β),α3=19(q2+12β)2+112β2,

    where β is the form in Eq (4.12). If u=(23s2)13ϕ>α1, then Eq (4.3) has the following solution by making use of the third-order polynomial complete discrimination system method

    u=α1+[2(α21+α2α1+α3)121+cn[(α21+α2α1+α3)14(ζζ0),m](α21+α2α1+α3)12].

    Substituting the specific forms of the variables u,ϕ,ζ,α1,α1, and α3 into the above equation and simplifying the expression yields the solution of Eq (2.3) given by (4.13).

    In order to visualize the structure of these traveling wave solutions, the figures of Φ1, Φ5, and Φ6 are shown in Figures 57 under some reasonable parameters by the aid of MATLAB software respectively. To compare the effect of fractional derivatives on the solution, different values of α are taken.

    Figure 5.  The figures of the solution Φ1 with y=3,z=2,a=1.5,b=1.6,c=1.5,v=2,s2=1,q2=3,p1=5,q0=0.5, and ζ0=6.
    Figure 6.  The figures of the solution Φ5 with y=2,z=2,a=1.5,b=2,c=1.5,v=2,s2=1,q2=3,p1=5,q0=0.5, and ζ0=5.
    Figure 7.  The figures of the solution Φ6 with y=2,z=2,a=1.5,b=2,c=1.5,v=2,s2=1,q2=3,p1=5,q0=0.5, and ζ0=5.

    Figure 5 shows a dark soliton of the (3+1)-dimensional KP-SK-R equation. The amplitude of the wave is around 4.3642 to 0.4650, and the width is approximately 9.8997. Throughout the simulation range, as the value of α increases, the wave crest moves to the right while maintaining the amplitude and width unchanged. Figure 6 presents the elliptic periodic solution of the (3+1)-dimensional KP-SK-R equation. This solution is manifested as a waveform that shows periodic changes in space and time, and its shape is similar to the contour of an ellipse. Within one period, the wave completes a full oscillation along the lower x-axis. The period is approximately 3.5338, and the amplitude of the wave varies periodically between 4.6051 and 0.2183. Figure 7 shows the lower part of the graph of a quasi-periodic solution of the (3+1)-dimensional KP-SK-R equation. Some of its wave crests are very high, while some are very low, that is, the amplitude differences of the waves are very large. It can be observed from the figure that the period is approximately 3.5743, but the overall shape of the wave does not precisely repeat within fixed time or spatial intervals as the periodic solution does.

    The (3+1)-dimensional KP-SK-R equation can describe not only the wave phenomenon propagating in weakly nonlinear and weakly dispersive media, but also the shallow water wave phenomenon. Moreover, it finds application in solid-state physics for characterizing the collective vibrational behavior of atoms within the lattice. This study aims to derive a range of new soliton solutions and traveling wave solutions of the (3+1)-dimensional KP-SK-R equation and to probe its behavior under external perturbations. Therefore, four theorems regarding the solution are given. The newly derived exact solutions enhance the understanding of the model structure of the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative. The introduction of external disturbances enables the observation of chaotic behavior within the system, thereby facilitating the exploration of the system's sensitivity and stability characteristics. Two and three-dimensional graphs visually depict the dynamic behavior. Overall, the findings reported herein significantly contribute to the theoretical development of the (3+1)-dimensional KP-SK-R equation with conformable fractional derivative and offer valuable insights into dispersion effects. Additionally, the trial method and polynomial discriminant system utilized in this research hold promise for extension to other fractional wave equations. In future research, building on the work presented in this paper, the stability of the traveling wave solutions of the (3+1)-dimensional KP-SK-R equation can be further investigated, the evolution of solutions under perturbation can be analyzed, and the asymptotic behavior of traveling wave solutions after long-term and long-distance propagation can be explored.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This paper is funded by Opening Fund of Geomathematics Key Laboratory of Sichuan Province (No. scsxdz2025-7).

    The authors declare there is no conflicts of interest.



    [1] Shahbaz M, Raghutla C, Chittedi KR, et al. (2020) The effect of renewable energy consumption on economic growth: Evidence from the renewable energy country attractive index. Energy 207: 118162. https://doi.org/10.1016/j.energy.2020.118162 doi: 10.1016/j.energy.2020.118162
    [2] Bellos E, Tzivanidis C (2019) Alternative designs of parabolic trough solar collectors. Prog Energy Combust Sci 71: 81–117. https://doi.org/10.1016/j.pecs.2018.11.001 doi: 10.1016/j.pecs.2018.11.001
    [3] Ouabouch O, Laasri IA, Kriraa M, et al. (2022) Effects of flow regime and geometric parameters on the performance of a parabolic trough solar collector using nanofluid. Numer Heat Transf Part A App 82: 1–13. https://doi.org/10.1080/10407782.2022.2078632 doi: 10.1080/10407782.2022.2078632
    [4] Tembhare SP, Barai DP, Bhanvase BA(2022) Performance evaluation of nanofluids in solar thermal and solar photovoltaic systems: A comprehensive review. Renew Sustain Energy Rev 153: 111738. https://doi.org/10.1016/j.rser.2021.111738 doi: 10.1016/j.rser.2021.111738
    [5] Qi C, Luo T, Liu M, et al. (2019) Experimental study on the flow and heat transfer characteristics of nanofluids in double-tube heat exchangers based on thermal efficiency assessment. Energy Convers Manag 197: 111877. https://doi.org/10.1016/j.enconman.2019.111877 doi: 10.1016/j.enconman.2019.111877
    [6] Choi SUS, Eastman JA (1995) Enhancing thermal conductivity of fluids with nanoparticles. 1995 International mechanical engineering congress and exhibition, San Francisco, CA (United States.
    [7] Ouabouch O, Kriraa M, Lamsaadi M (2021) Stability, thermophsical properties of nanofluids, and applications in solar collectors: A review. 8: 659–684. https://doi.org/10.3934/matersci.2021040
    [8] Ghasemi SE, Ranjbar AA (2016) Thermal performance analysis of solar parabolic trough collector using nanofluid as working fluid: A CFD modelling study. J Mol Liq 222: 159–166. https://doi.org/10.1016/j.molliq.2016.06.091 doi: 10.1016/j.molliq.2016.06.091
    [9] Ouabouch O, Laasri IA, Kriraa M, et al. (2021) Modelling and comparison of the thermohydraulic performance with an economical evaluation for a parabolic trough solar collector using different nanofluids. Int J Heat Technol 39: 1763–1769. https://doi.org/10.18280/ijht.390609 doi: 10.18280/ijht.390609
    [10] Bellos E, Tzivanidis C, Antonopoulos KA, et al. (2016) Thermal enhancement of solar parabolic trough collectors by using nanofluids and converging-diverging absorber tube. Renew Energy 94: 213–222. https://doi.org/10.1016/j.renene.2016.03.062 doi: 10.1016/j.renene.2016.03.062
    [11] García A, Vicente PG, Viedma A (2005) Experimental study of heat transfer enhancement with wire coil inserts in laminar-transition-turbulent regimes at different Prandtl numbers. Int J Heat Mass Transf 48: 4640–4651. https://doi.org/10.1016/j.ijheatmasstransfer.2005.04.024 doi: 10.1016/j.ijheatmasstransfer.2005.04.024
    [12] Mwesigye A, Bello-Ochende T, Meyer JP (2015) Multi-objective and thermodynamic optimisation of a parabolic trough receiver with perforated plate inserts. Appl Therm Eng 77: 42–56. https://doi.org/10.1016/j.applthermaleng.2014.12.018 doi: 10.1016/j.applthermaleng.2014.12.018
    [13] Khan MS, Yan M, Ali HM, et al. (2020) Comparative performance assessment of different absorber tube geometries for parabolic trough solar collector using nanofluid. J Therm Anal Calorim 142: 2227–2241. https://doi.org/10.1007/s10973-020-09590-2 doi: 10.1007/s10973-020-09590-2
    [14] Liu Y, Chen Q, Hu K, et al. (2016) Flow field optimization for the solar parabolic trough receivers in direct steam generation systems by the variational principle. Int J Heat Mass Transf 102: 1073–1081. https://doi.org/10.1016/j.ijheatmasstransfer.2016.06.083 doi: 10.1016/j.ijheatmasstransfer.2016.06.083
    [15] Saedodin S, Zaboli M, Ajarostaghi SSM (2021) Hydrothermal analysis of heat transfer and thermal performance characteristics in a parabolic trough solar collector with Turbulence-Inducing elements. Sustain Energy Techn 46: 101266. https://doi.org/10.1016/j.seta.2021.101266 doi: 10.1016/j.seta.2021.101266
    [16] Olfian H, Ajarostaghi SSM, Farhadi M, et al. (2021) Melting and solidification processes of phase change material in evacuated tube solar collector with U-shaped spirally corrugated tube. Appl Therm Eng 182: 116149. https://doi.org/10.1016/j.applthermaleng.2020.116149 doi: 10.1016/j.applthermaleng.2020.116149
    [17] Fan F, Qi C, Liu Q, et al. (2020) Effect of twisted turbulator perforated ratio on thermal and hydraulic performance of magnetic nanofluids in a novel thermal exchanger system. Case Stud. Therm Eng 22: 100761. https://doi.org/10.1016/j.csite.2020.100761 doi: 10.1016/j.csite.2020.100761
    [18] Song X, Dong G, Gao F, et al. (2014) A numerical study of parabolic trough receiver with nonuniform heat flux and helical screw-tape inserts. Energy 77: 771–782. https://doi.org/10.1016/j.energy.2014.09.049 doi: 10.1016/j.energy.2014.09.049
    [19] Khanafer K, Vafai K (2011) A critical synthesis of thermophysical characteristics of nanofluids. Int J Heat Mass Transf 54: 4410–4428. https://doi.org/10.1016/j.ijheatmasstransfer.2011.04.048 doi: 10.1016/j.ijheatmasstransfer.2011.04.048
    [20] Xuan Y, Roetzel W (2000) Conceptions for heat transfer correlation of nanofluids. Int J Heat Mass Transf 43: 3701–3707. https://doi.org/10.1016/S0017-9310(99)00369-5 doi: 10.1016/S0017-9310(99)00369-5
    [21] Brinkman HC (1952) The viscosity of concentrated suspensions and solutions. J Chem Phys 20: 571. https://doi.org/10.1063/1.1700493 doi: 10.1063/1.1700493
    [22] Maxwell JC (1891) A Treatise on Electricity and Magnetism, UK: Clarendon Press.
    [23] Turkyilmazoglu M (2017) Condensation of laminar film over curved vertical walls using single and two-phase nanofluid models. Eur J Mech B/Fluids 65: 184–191. https://doi.org/10.1016/j.euromechflu.2017.04.007 doi: 10.1016/j.euromechflu.2017.04.007
    [24] Yuan SW (1967) Foundations of Fluid Mechanics, New York: Prentice-Hall.
    [25] Chakraborty O, Das B, Gupta R, et al. (2020) Heat transfer enhancement analysis of parabolic trough collector with straight and helical absorber tube. Therm Sci Eng Prog 20: 100718. https://doi.org/10.1016/j.tsep.2020.100718 doi: 10.1016/j.tsep.2020.100718
    [26] Petukhov BS (2017) Heat transfer and friction in turbulent pipe flow with variable physical properties. Adv Heat Transf 6: 503–564. https://doi.org/10.1016/S0065-2717(08)70153-9 doi: 10.1016/S0065-2717(08)70153-9
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1870) PDF downloads(84) Cited by(1)

Figures and Tables

Figures(24)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog