A wormhole is an imaginary intuitive concept in general relativity. It acts like a bridge or tunnel to connect two or more asymptotic regions. However, this hypothetical object has become one of the most popular and intensively studied research area in general relativity. The studies so far in this topic can be divided into two classes: static wormholes and dynamic wormholes. Although, there are static wormhole solutions [1–5] since 1973, but, the work by Morris and Thorne [6] has the key role in studying the static wormholes. Usually, the static wormhole space–time is sustained by a single fluid component which requires the violation of the null energy condition (NEC) [7–13]. However, in asymptotically flat space–time, this violation of NEC is a consequence of the topological censorship [14, 15]. Most of the studies in wormholes are related to traversable wormholes which have no horizons, and, as a result, there is two way passage through them. Although, the speed of light is not locally surpasses [16], but due to global space–time topology [7, 16, 17], it is possible to have superluminal travel through these wormholes, and, as a result, there is the idea of time machines [18–20].

Further, it should be noted that, it is possible to construct wormhole space–times with an arbitrarily small violation of the averaged NEC [16]. So, it is speculated that, the wormhole configuration could be realized merely by quantum effects violating the energy conditions.

In general, wormhole geometries are not constructed by solving Einstein field equations, rather one first fixes the form of the space–time metric (i.e., redshift and shape functions) and then matter part is evaluated by computing the field equations. Due to Bianchi identities, the matter part so obtained automatically obey the local conservation equations and violates the NEC [6, 7, 9, 18, 21]. It should be noted that, in modified gravity theories, there are examples of traversable wormhole solutions without any violation of energy conditions, for example in Einstein–Gauss–Bonnet gravity [36, 37] and in higher dimensional Lovelock theories [22–25]. Also, in this context, there are well known non-static Lorentzian wormholes in Einstein gravity where matter component may obey weak energy condition (WEC) but life time may be arbitrarily small, or, large intervals of time [26, 27].

On the other hand, dynamical wormholes (i.e., evolving relativistic wormholes [28–32]) are not as popular as static wormholes and also not well understood. The pioneering work related to dynamical wormholes was done independently by Hochberg and Visser [33] and Hayward [34]. They independently choose quasi local definition of wormhole throat in a dynamical space–time. Essentially, wormhole throat is a trapping horizon [35] of different kind, but, matter in both of them violates the NEC. However, Maeda et al. [36, 37] have shown another class of dynamical wormholes (cosmological wormholes) which are asymptotically Friedmann universe with a big bang singularity at the beginning. This class of wormholes do not need matter which violates NEC rather the dominant energy condition (DEC) is satisfied everywhere. The basic difference between these two class of dynamical wormholes is purely from geometrical aspect. In the former case, wormhole throat is a 2D surface of non-vanishing minimal area on a null hypersurface, while, in the later class of wormholes, due to initial singularity, there is no past null infinity [38]. As a result, the wormhole throat is defined only on a space–like hypersurface. Hence, there is no trapping horizon rather the space–times are trapped everywhere [38]. In recent years, there are works with dynamic wormhole space–time filled with two fluids [39, 40]. Such a matter system is very much relevant in present day cosmology where such two fluid models are widely used to describe the observed accelerated expansion of the universe [41–44].

In the present work, we make an attempt to find dynamical wormhole solutions for two fluid system, where one is a dissipative homogeneous fluid, and the other fluid component is anisotropic and inhomogeneous in nature. We assume that the dissipation arises due to particle creation mechanism in non-equilibrium thermodynamics which for simplicity is assumed to be isentropic in nature. As a result, the dissipative pressure is linearly related to the particle creation rate [45, 46]. The paper is organized as follows: a review of earlier works on wormhole configuration with two non-interacting fluids has been done in Sect. 2. In Sect. 3, particle creation mechanism in the non-equilibrium thermodynamic prescription has been presented. Section 4 shows possible wormhole solutions for different choices of the particle creation parameter. The paper ends with a brief overview in Sect. 5.

The metric ansatz for the dynamic wormhole space–time is given by1ds2=-e2Φ(r,t)dt2+a2(t)dr21-b(r)r-Kr2+r2dΩ22,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm{d}s^2= -\mathrm{e}^{2 \Phi (r, t)} \mathrm{d}t^2+ a^2 (t) \left[ \frac{\mathrm{d}r^2}{1-\frac{b(r)}{r}-K r^2}+r^2 \mathrm{d}\Omega _2 ^2\right] \!,\!\!\nonumber \\ \end{aligned}$$\end{document}where Φ(r,t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi (r, t)$$\end{document} is the redshift function; a(t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a(t)$$\end{document} is the scale factor of the wormhole universe; b(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b(r)$$\end{document} is the usual shape function for the wormhole; dΩ22=dθ2+sin2θdϕ2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}\Omega _2 ^2= \mathrm{d}\theta ^2+ \sin ^2 \theta \mathrm{d}\phi ^2$$\end{document}, and K\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K$$\end{document} takes values 0, ±1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm 1$$\end{document}. In particular, if Φ(r,t)→Φ(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi (r, t)\rightarrow \Phi (r)$$\end{document} and a(t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a(t)$$\end{document}→\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rightarrow $$\end{document}a0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0$$\end{document}, a constant, then the above metric ansatz describes a static wormhole universe, while metric (1) describes a FRW model, if Φ(r,t)=0=b(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi (r, t)= 0= b(r)$$\end{document}.

Suppose the matter distribution of the wormhole universe is described by two non-interacting fluid components (termed as Fluid I and Fluid II) together with a cosmological constant Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}. Fluid I is homogeneous and isotropic, but dissipative in nature having energy–momentum tensor:2TμνI=(ρ+p+Π)uμuν+(p+Π)gμν,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} T_{\mu \nu }^\mathrm{I}= (\rho + p+ \Pi ) u_\mu u_\nu + (p+ \Pi ) g_{\mu \nu }, \end{aligned}$$\end{document}where ρ=ρ(t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho = \rho (t)$$\end{document} is the energy density, p=p(t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p= p (t)$$\end{document} and Π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi $$\end{document} are the isotropic pressure and the pressure due to dissipation respectively and uμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u_\mu $$\end{document} is the four velocity of the fluid. On the other hand, Fluid II is both inhomogeneous and anaisotropic in nature with the energy–momentum tensor:3TμνII=(ρin+ρt)vμvν+ptgμν+(pr-pt)χμχν,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} T_{\mu \nu }^{\mathrm{II}}= (\rho _{\mathrm{in}}+ \rho _t) v_\mu v_\nu + p_t g_{\mu \nu }+ (p_r- p_t) \chi _\mu \chi _\nu , \end{aligned}$$\end{document}where, ρin=ρin(t,r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho _{\mathrm{in}}= \rho _{\mathrm{in}} (t, r)$$\end{document} is the energy density of the inhomogeneous fluid component, the anisotropic pressure is characterized by radial and transversal components by pr=pr(t,r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_r= p_r (t, r)$$\end{document}, and pt=pt(t,r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_t= p_t (t, r)$$\end{document} respectively (pr=pt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_r= p_t$$\end{document} implies that Fluid II is isotropic but inhomogeneous in nature), and vμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_\mu $$\end{document} and χμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi _\mu $$\end{document} are respectively unit time-like and space-like vectors, i.e., vμvμ=-χμχμ=-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_\mu v^\mu = -\chi _\mu \chi ^\mu = -1$$\end{document}, χμvμ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^\mu v_\mu = 0$$\end{document}. Thus, the explicit form of the Einstein’s field equations:4Gμν=-κTμνI+TμνII-Λgμν,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} G_{\mu \nu }= -\kappa \left( T_{\mu \nu }^\mathrm{I}+ T_{\mu \nu }^{\mathrm{II}}\right) - \Lambda g_{\mu \nu }, \end{aligned}$$\end{document}are [39].53e-2Φ(r,t)H2+b′a2r2+3Ka2=κρin+κρ+Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} 3 \mathrm{e}^{-2 \Phi (r, t)} H^2+ \frac{b^\prime }{a^2 r^2}+ \frac{3 K}{a^2}= \kappa \rho _{\mathrm{in}}+ \kappa \rho +\Lambda , \end{aligned}$$\end{document}6-e-2Φ(r,t)2a¨a+H2+Ka2-ba2r3+2e-2Φ(r,t)H∂Φ∂t+2r2a2(r-b)∂Φ∂r=κpr+κ(p+Π)-Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&-\mathrm{e}^{-2 \Phi (r, t)} \left( \frac{2 \ddot{a}}{a}+ H^2\right) + \frac{K}{a^2}- \frac{b}{a^2 r^3}+ 2\mathrm{e}^{-2 \Phi (r, t)} H \frac{\partial \Phi }{\partial t} \nonumber \\&\quad +\frac{2}{r^2 a^2} (r-b) \frac{\partial \Phi }{\partial r}= \kappa p_r+ \kappa (p+ \Pi )- \Lambda , \end{aligned}$$\end{document}7e-Φ(r,t)2a¨a+H2+Ka2+b-rb′2a2r3+2e-Φ(r,t)H∂Φ∂t+2r-b-rb′2a2r2∂Φ∂r+r-ba2r∂Φ∂r2+∂2Φ∂2r=κpt+κ(p+Π)-Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\mathrm{e}^ {-\Phi (r, t)} \left( \frac{2 \ddot{a}}{a} + H^2\right) + \frac{K}{a^2}+ \frac{b-r b^\prime }{2 a^2 r^3}+ 2 \mathrm{e}^{-\Phi (r, t)} H \frac{\partial \Phi }{\partial t}\nonumber \\&\quad \quad + \left( \frac{2r-b-rb^\prime }{2a^2 r^2}\right) \frac{\partial \Phi }{\partial r} +\frac{r-b}{a^2 r} \left[ \left( \frac{\partial \Phi }{\partial r}\right) ^2+ \frac{\partial ^2 \Phi }{\partial ^2 r}\right] \nonumber \\&\quad =\kappa p_t +\kappa (p+ \Pi )- \Lambda , \end{aligned}$$\end{document}and82a˙e-Φ(r,t)r-b(r)r∂Φ(r,t)∂r=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} 2 \dot{a} \mathrm{e}^ {-\Phi (r, t)} \left( \sqrt{\frac{r-b(r)}{r}}\right) \frac{\partial \Phi (r, t)}{\partial r}= 0, \end{aligned}$$\end{document}where κ=8πG\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa = 8 \pi G$$\end{document}, uα=(e-Φ,0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u^\alpha = (\mathrm{e}^{-\Phi }, 0, 0, 0)$$\end{document} is the time-like vector denoting the four velocity of both the fluids, H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H$$\end{document}=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=$$\end{document}a˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{a}$$\end{document}/a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a$$\end{document} is the Hubble parameter, and an ‘overdot’, or, a ‘prime’ denotes the differentiation with respect to the cosmic time ‘t\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t$$\end{document}’, or, the radial co-ordinate r\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r$$\end{document} respectively. From the field Eq. (5), we see that two classes of solutions are possible, namely,(I)a˙=0(static);(II)∂Φ∂r=0(non-static).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} (\mathrm{I})\quad \dot{a}=0\quad (\mathrm{static});\quad ~ (\mathrm{II})\quad \frac{\partial \Phi }{\partial r}= 0\quad ~(\hbox {non-static}). \end{aligned}$$\end{document}So in the present work, we shall consider only the second choice for dynamic wormhole solutions. As a result, (without any loss of generality) we can choose Φ(r,t)=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi (r, t)= 0$$\end{document} (a rescaling of the time co-ordinate). Thus the wormhole metric (1) now simplifies to9ds2=-dt2+a2(t)dr21-b(r)r-Kr2+r2dΩ22.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm{d}s^2= -\mathrm{d}t^2+ a^2(t) \left[ \frac{\mathrm{d}r^2}{1-\frac{b(r)}{r}-K r^2}+r^2 \mathrm{d}\Omega _2 ^2 \right] . \end{aligned}$$\end{document}Now, due to non-interacting nature of the two fluids, both of them satisfy the conservation equations separately as10∂ρ∂t+3H(ρ+p+Π)=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\partial \rho }{\partial t}+ 3H (\rho + p+ \Pi )= 0, \end{aligned}$$\end{document}11∂ρin∂t+H(3ρin+pr+2pt)=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\partial \rho _{\mathrm{in}}}{\partial t}+ H (3 \rho _{\mathrm{in}}+ p_r+ 2 p_t)= 0, \end{aligned}$$\end{document}and12∂pr∂r=2r(pt-pr).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\partial p_r}{\partial r}= \frac{2}{r} (p_t- p_r). \end{aligned}$$\end{document}Here, Eq. (10) is the conservation equation for the homogeneous but dissipative fluid (i.e., Fluid I). Equations (11) and (12) are the conservation equations for the other fluid component (i.e., Fluid II). Note that Eq. (12) is nothing but the relativistic Euler equation.

One may notice that the anisotropic nature (i.e., pt≠pr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_t \ne p_r$$\end{document}) of the inhomogeneous fluid (i.e., Fluid II) is essential, otherwise, the pressure components become homogeneous as seen from Eq. (12), and then the conservation Eq. (11) demands that the density is also homogeneous. So, essentially, we have two non-interacting homogeneous fluid components leading a physically uninteresting situation. Hence Fluid II must be anisotropic (i.e., pt≠pr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_t \ne p_r$$\end{document}), and thus inhomogeneous.

We now write down the simplified form of the field Eqs. (5)–(7) for dynamical wormhole solution as [39, 40]133H2+3Ka2+b′a2r2=κρ+κρin+Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&3 H^2+ \frac{3K}{a^2}+ \frac{b^\prime }{a^2 r^2}= \kappa \rho + \kappa \rho _{\mathrm{in}}+ \Lambda ,\end{aligned}$$\end{document}14-2H˙+3H2+Ka2-ba2r3=κ(p+Π)+κpr-Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&-\left( 2 \dot{H}+ 3 H^2+ \frac{K}{a^2}\right) -\frac{b}{a^2 r^3}= \kappa (p+ \Pi )+ \kappa p_r- \Lambda ,\nonumber \\ \end{aligned}$$\end{document}and15-2H˙+3H2+Ka2+b-rb′2a2r3=κ(p+Π)+κpt-Λ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&-\left( 2 \dot{H}+ 3 H^2+ \frac{K}{a^2}\right) \!+ \!\frac{b-r b^\prime }{2 a^2 r^3}\!=\! \kappa (p+ \Pi )\!+\! \kappa p_t\!-\! \Lambda .\!\!\nonumber \\ \end{aligned}$$\end{document}Now, in order to solve the above non-linear field equations, we assume for simplicity that the radial and the transversal pressure components of Fluid II satisfy barotropic equation of state [40]:16pr(t,r)=ωrρin,andpt(t,r)=ωtρin,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_r (t, r)= \omega _r \rho _{\mathrm{in}},\quad \text {and}\quad ~p_t (t, r)= \omega _t \rho _{\mathrm{in}}, \end{aligned}$$\end{document}where the constants ωr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _r$$\end{document} and ωt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _t$$\end{document} stand for the equation of state parameters. Now, inserting the relations in Eq. (16) to the conservation Eqs. (11) and (12), one immediately gets17ρin(t,r)=ρ0r2(ωt-ωr)ωra3+ωr+2ωt(ρ0=constant of integration).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \rho _{\mathrm{in}} (t, r)= \rho _0 \frac{r^{\frac{2 (\omega _t- \omega _r)}{\omega _r}}}{a^{3+ \omega _r+ 2 \omega _t}} \quad (\rho _0= \text {constant of integration}).\nonumber \\ \end{aligned}$$\end{document}Now, comparing the field Eqs. (14) and (15) and using Eq. (17), the shape function b(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b(r)$$\end{document} can be obtained as18b(r)=K0r3-κρ0ωrr-1ωr,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} b(r)= K_0 r^3-\kappa \rho _0 \omega _r r^{-\frac{1}{\omega _r}}, \end{aligned}$$\end{document}provided, we assume that the equation of state parameters are not independent, rather they are related by the following relation19ωr+2ωt+1=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \omega _r+ 2 \omega _t+ 1= 0. \end{aligned}$$\end{document}In the above solution for the shape function b(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b (r)$$\end{document}, the integration constant K0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_0$$\end{document} behaves as the curvature constant (K\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K$$\end{document}) in the metric shown in Eq. (1) or in Eq. (9). So, without any loss of generality, this integration constant K0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_0$$\end{document} may be absorbed by rescaling the radial co-ordinate ‘r\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r$$\end{document}’ as follows:K+K0=1,when,K+K0>0.K+K0=-1,when,K+K0<0.K+K0=0,for(K=1,K0=-1;K=-1,K0=1;K=0,K0=0).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{array}{ll} K+ K_0= 1, &{}\quad \text {when},\; K+K_0> 0.\\ K+ K_0= -1,&{} \quad \text {when},\; K+K_0< 0.\\ K+ K_0= 0,&{}\quad \text {for } (K= 1, K_0= -1; K= -1,\\ &{}\quad \quad \quad \quad K_0= 1; K=0, K_0= 0). \end{array} \end{aligned}$$\end{document}The gravitational configuration is described by the metric ansatz20ds2=-dt2+a2(t)dr21-Kr2-(rr0)-1+ωrωr+r2dΩ22,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm{d}s^2= -\mathrm{d}t^2+ a^2 (t)\left[ \frac{\mathrm{d}r^2}{1-K r^2-(\frac{r}{r_0})^{-\frac{1+ \omega _r}{\omega _r}}}+ r^2 \mathrm{d}\Omega _2 ^2\right] ,\nonumber \\ \end{aligned}$$\end{document}there are two non-interacting fluid system in which the anisotropic and inhomogeneous matter component has energy density21ρin(t,r)=ρ0r-3-1ωra2=-(rr0)-3-1ωrκa2ωrr02,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \rho _{\mathrm{in}} (t, r)= \frac{\rho _0 r^{-3-\frac{1}{\omega _r}}}{a^2}= -\frac{(\frac{r}{r_0})^{-3-\frac{1}{\omega _r}}}{\kappa a^2 \omega _r r_0 ^2}, \end{aligned}$$\end{document}and the thermodynamic pressures along radial and transverse directions are22pr(t,r)=ωrρin,andpt(t,r)=-12(1+ωr)ρin,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_r (t, r)= \omega _r \rho _{\mathrm{in}},\quad \text {and} \quad p_t (t, r)= -\frac{1}{2} (1+ \omega _r) \rho _{\mathrm{in}}, \end{aligned}$$\end{document}while the homogeneous and isotropic part is described by the Friedmann equations233H2+Ka2=κρ+Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} 3\left( H^2+ \frac{K}{a^2}\right) = \kappa \rho + \Lambda , \end{aligned}$$\end{document}and24-2H˙+3H2+Ka2=κ(p+Π)-Λ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} -\left( 2 \dot{H}+ 3 H^2+ \frac{K}{a^2}\right) = \kappa (p+ \Pi )- \Lambda , \end{aligned}$$\end{document}and the fluid components are related by the conservation Eq. (10). Thus for the present dynamic wormhole universe, the rate of expansion of these evolving wormholes is fully characterized by the homogeneous and isotropic, but dissipative matter component. In the following section we shall take an attempt to find the wormhole solutions when dissipative phenomena is caused by the non-equilibrium thermodynamics due to particle creation mechanism.

As the number of particles is not conserved in non-equilibrium thermodynamic prescription, so the conservation equation for particle number takes the form [47]25n˙+3θn=nΓ,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \dot{n}+ 3 \theta n= n \Gamma , \end{aligned}$$\end{document}where n=NV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n= \frac{N}{V}$$\end{document}, is the particle number density; N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N$$\end{document} is the total number of particles in a co-moving volume V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V$$\end{document}; Nμ=nuμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N^\mu = n u^\mu $$\end{document} is the particle flow vector; θ=u;μμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta = u^\mu _{;\mu }$$\end{document} stands for fluid expansion; Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} represents the particle creation rate, and notationally, n˙=n;μuμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{n}= n_{; \mu } u^\mu $$\end{document}. The sign of Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} indicates creation (Γ>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma > 0$$\end{document}), or, annihilation (Γ<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma < 0$$\end{document}) of particles and Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} represents some dissipative effect to the Cosmic fluid, so that, non-equilibrium thermodynamics comes into picture.

Using Clausius relation, the Gibb’s equation takes the form [47]26Tds=dρn+pd1n,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} T ds= d \left( \frac{\rho }{n}\right) +p d \left( \frac{1}{n}\right) , \end{aligned}$$\end{document}where ‘s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s$$\end{document}’ is the entropy per particle and T\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T$$\end{document} is the fluid temperature. Now, using the conservation relations (10) and (25), the entropy variation can be expressed as [48]27nTs˙=-Πθ-Γρ+p.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} n T \dot{s}= -\Pi \theta - \Gamma \left( \rho + p\right) \!. \end{aligned}$$\end{document}If for simplicity, we assume the thermal process to be ‘adiabatic’ (or, ‘isentropic’, i.e., s˙=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{s}= 0$$\end{document}), then from the above equation we have28Π=-Γθρ+p.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Pi = -\frac{\Gamma }{\theta } \left( \rho + p\right) . \end{aligned}$$\end{document}Hence the dissipative pressure is completely characterized by the particle creation rate for the above isentropic thermodynamical system. Alternatively, the fluid may be considered as perfect fluid with barotropic equation of state: p=(γ-1)ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p= (\gamma - 1) \rho $$\end{document}, and, dissipative phenomena comes into picture through particle creation. Note that, although the entropy per particle is constant but still there is entropy generation due to particle creation, i.e., enlargement of the phase space due to expansion of the universe. So, non-equilibrium configuration is not the conventional one due to the effective bulk pressure, rather a state with equilibrium properties as well (but not the equilibrium era with Γ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma = 0$$\end{document}). Now, we eliminate ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}, p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p$$\end{document} and Π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi $$\end{document} from the Einstein field Eqs. (23) and (24), and the isentropic Eq. (28), and then using barotropic equation of state parameter γ=1+pρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma = 1+ \frac{p}{\rho }$$\end{document}, we obtain [48]29Γ3H=1+23γH˙-Ka2H2+Ka2-Λ3.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\Gamma }{3H}= 1+ \frac{2}{3 \gamma } \left( \frac{\dot{H}-\frac{K}{a^2}}{H^2+ \frac{K}{a^2}-\frac{\Lambda }{3}}\right) . \end{aligned}$$\end{document}So, Eq. (29) helps us to conclude that for adiabatic thermodynamical system, the particle creation rate is related to the evolution of the universe.

In this section, we shall determine evolving wormhole solutions choosing the particle creation rate (Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}) as a function of the Hubble parameter (H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H$$\end{document}). For simplicity, we take the flat model (i.e., K=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K= 0$$\end{document}) of the Universe.

From the field Eqs. (23) and (24), the acceleration equation takes the form30a¨a=-κ6ρ+3p+Π+Λ3.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\ddot{a}}{a}= -\frac{\kappa }{6} \left[ \rho + 3\left( p+ \Pi \right) \right] + \frac{\Lambda }{3}. \end{aligned}$$\end{document}Hence for expansion with constant velocity, we must have31Λ=κ2ρ+3(p+Π),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Lambda = \frac{\kappa }{2} \left[ \rho + 3 (p+ \Pi )\right] \!, \end{aligned}$$\end{document}or, using Eq. (28) for ‘isentropic’ condition, the particle creation rate for uniform velocity is restricted by32Γ=3Hγ(3γ-2)-2Λκρ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Gamma = \frac{3 H}{\gamma } \left[ (3 \gamma - 2)- \frac{2 \Lambda }{\kappa \rho }\right] . \end{aligned}$$\end{document}So, for Λ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda =0$$\end{document} (or, Λ∝H2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \propto H^2$$\end{document}) we have Γ∝ρ12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma \propto \rho ^{\frac{1}{2}}$$\end{document} (or, Γ∝H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma \propto H$$\end{document}). Also, for Λ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda = 0$$\end{document}, positivity of Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} restricts γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} to: γ>23\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma > \frac{2}{3}$$\end{document} or γ<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma < 0$$\end{document}; i.e., the homogeneous and isotropic fluid must be a normal fluid (satisfying strong energy condition) or is in phantom domain (i.e., violating WEC) while there will be particle annihilation if the homogeneous fluid is in the quintessence era.

To evaluate the solutions, we start with the evolution Eq. (24) for K=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K= 0$$\end{document}, Λ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda = 0$$\end{document} and κ=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa = 1$$\end{document}, and using Eq. (28) for Π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi $$\end{document} we have332H˙+3γH2-γHΓ=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} 2 \dot{H}+ 3 \gamma H^2- \gamma H \Gamma = 0, \end{aligned}$$\end{document}which can be integrated to give34H=expγ2∫ΓdtH0+3γ2∫expγ2∫Γdtdt.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H= \frac{\exp \left( {\frac{\gamma }{2} \int \Gamma \mathrm{d}t}\right) }{H_0+ \frac{3 \gamma }{2} \int \exp \left( {\frac{\gamma }{2} \int \Gamma \mathrm{d}t}\right) \mathrm{d}t}. \end{aligned}$$\end{document}Here H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} is the constant of integration. The scale factor evolves as35a=a0H0+3γ2∫expγ2∫Γdtdt23γ;(a0=constant of integration).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} a&= a_0 \left[ H_0+ \frac{3 \gamma }{2} \int \exp \left( {\frac{\gamma }{2} \int \Gamma \mathrm{d}t}\right) \mathrm{d}t\right] ^{\frac{2}{3 \gamma }};\nonumber \\&(a_0= \text {constant of integration}). \end{aligned}$$\end{document}We shall now determine the explicit solutions for the following choices of Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}.

In the present work we deal with a FRW like space–time model which is both inhomogeneous and anisotropic in nature and there is a future singularity at a finite proper time. There are two non-interacting matter components—one is isotropic and homogeneously distributed dissipative fluid, while the other matter component is both inhomogeneous and anisotropic in nature. Here, we have considered dissipation due to particle creation mechanism and for simplicity we restrict ourselves to adiabatic process so that the dissipative pressure is linearly related to the particle creation rate. The solutions presented in the paper describe evolving wormholes which are threaded and sustained by the inhomogeneous and anisotropic fluid component, while the rate of expansion is characterized by the isotropic and homogeneous matter part which in most of the cases is chosen in the phantom domain. Three cosmological models are presented in the paper corresponding to three different choices of the particle creation rates, and in all of them, the wormhole models encounter a big rip singularity in course of its evolution. As all the wormhole solutions are asymptotically flat FRW cosmologies, so all the results on future singularities obtained are also true for flat FRW cosmological models.

Here, the dissipation is chosen as bulk viscosity and the corresponding pressure is related to the particle creation mechanism. As in the literature [49], the bulk viscosity is in the power law form of the Hubble parameter, so, due to the isentropic condition [Eq. (28)], it is reasonable to choose the particle creation rate as some power of the Hubble parameter. For simplicity, we have restricted to: (1) Γ=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma =$$\end{document} constant [50], (2) Γ∝H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma \propto H$$\end{document} [48], and (3) Γ∝1H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma \propto \frac{1}{H}$$\end{document} [51]. The solutions corresponding to the first and third choices, the homogeneous and isotropic fluid is always phantom in nature, but for the second choice, there are two possibilities of which one is similar as the other choices while for the second possibility, there exists evolving wormhole without phantom energy. Hence one can say that phantom energy is not essential for describing evolving wormholes.

Further, it should be noted that in the present work the matter component is chosen such that there is no mixed component of the energy–momentum tensor (i.e., Trt=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{rt}= 0$$\end{document}). As a result, neither there is any radial energy flow, nor there is any accretion onto the wormhole from the Cosmic fluid. So, we may conclude that the present model can not be able to explain the big trip mechanism.