**1. Introduction**

An air-pocket-incited gas-fluid stream is the premise of smooth motion in numerous compound building gadgets and applications ranging from boilers or evaporators and more than a few stage bubble segment reactors of different structures to enormous-scale vigorous (and sometimes anaerobic) sewage treatment plants. The two-phase transport hypothetical talk and test request are firmly connected. On the other hand, the amalgamation that emerges from this association creates colossal innovative potential for estimations advising and approving unique models. The subsequent innovation develops utility in an expansive range of uses, from cutting-edge atomic hardware and space motors to pharmaceutical assembling, nourishment innovation, vitality, and natural remediation. Sussman et al. [1] considered a level-set methodology for figuring out answers for a incompressible two-stage stream. Their study was about the movement of air bubbles in the water and falling water drops in the air. A single-liquid model for a two-stage stream with variable thickness to the stream-water flow was analyzed by Bankoff [2]. Zuber and Findlay [3] estimated the normal volumetric fixation in two-stage stream frameworks. In their analysis, the outcomes anticipated by the investigation were contrasted with experimental data acquired for different two-phase stream systems, with different fluid gas blends in the adiabatic, vertical stream over a wide weight territory. Picchi and Poesio [4] developed a unified model for both horizontal and slightly inclined fluid

pipes lubricated with two-phase gas/shear-thinning fluid. Sato and Sekoguchi [5] suggested the velocity distribution of liquid in two-phase bubble flow. A more precise analytical procedure was constructed that created the justified foreboding of the liquid velocity dispensation in two-phase bubble flow. Kuwagi et al. [6] investigated the oscillation of bubbly flow through a normally placed cylinder employing a tridimensional system. Picchi and Battiato [7] discussed immiscible two-phase flow in porous media and elaborated the impact of pore-scale flow. Bonzanini et al. [8] simulated 1-D slug and stratified flow in pipes. Sontti and Atta [9] investigated co-flow in microchannels to discuss the viscous effect on Taylor bubble formation. Bhatti et al. [10] broke down the heat and mass exchange of a two-phase stream with an electric twofold layer whose impacts were incited due to the peristaltic impetus within the sight of the transverse attractive field. Haider et al. [11] presented the heat transfer as well as a magnetic field investigation on the peristaltically initiated movement of tiny particles.

Moreover, the mechanism of peristalsis comprises expansion and contraction events that impel an ingredient forward. Examples of some cases of the peristaltic phenomenon are the transport of bile in the bile duct, the transport of urine from the kidney to the gallbladder, the transport of cilia, the vasomotion of small blood vessels, and the mixing of food in the digestive tract, to name a few. The peristaltic phenomenon also has several industrial applications such as in the flow in tube pumps, in the rollers and hoses in heart-lung machines, and in the dialysis process during open-heart surgery. [12,13].

Furthermore, in buoyancy-driven flows, although the difference in inertia is almost negligible, the gravity remains sufficiently strong to make the specific weight significantly different during the flow in multiphase fluids. Tripathi et al. [14] investigated buoyancy effects in the peristaltic flow of nanofluid under the influence of electro-osmosis. Animasaun and Pop [15] numerically explored the effects of buoyancy on the flow driven by catalytic surface reactions. Angirasa et al. [16] reported the buoyancy effects in a fluid saturated with a porous medium. Rashidi et al. [17] studied fluid flow in the presence of buoyancy forces.

In addition to the above, the presence of bubbles has appeared in several applications in a gas-liquid flow. Many theoretical and numerical investigations have been conducting for multiphase bubbly flows in oil, gas, and liquid. It has numerous usages, such as in optical fiber sensing applications [18], sensitive pressure measurement [19], the human bloodstream during decompression sickness, and subcooled flow boiling in macro-channels [20]. Ellhi et al. [21] examined the simulation of bubble through the nozzle of tube. Furthermore, lubricated coatings with bubbles in peristaltic motion have a lot of applications in the biomedical field, and among several of these, in the control of blood pressure. Particle coating with viscous liquids is an essential component in the industry for surface modification purposes in order to induce and improve precise functionalities. Coating with viscous liquids is usually present in very rich industrial trials, which allow the coating of particles under the shear forces exerted in a mixing device. A list of key investigations on peristaltic flows [22–24], multiphase flows [25–29], bubbly flows [30–34], coatings [35–40], elastic medium [41,42], and several other references are available, which provide a more in-depth understanding to the reader.

Due to the immense contribution of two-phase flow structures in many significant fields, this study attempts to trap the structures of gas-liquid bubbly flow inside the elastic walls under the peristaltic mechanism applied over a two-fluid model. Due to the nonlinear model, a powerful and efficient technique called the homotopy perturbation method is used for finding analytic solutions. This method works even without the need of a linearization process of nonlinear differential equations. The parameters affecting the flow prominently have been examined with the help of a graphical illustration.

#### **2. Mathematical Formulation**

We assumed a symmetric channel with flexible walls starting at the origin of a rectangular coordinate system, see Figure 1. Plates were placed parallel to the x-axis, on either side of the origin, separated by a distance of 2h. Only the upper portion of the symmetric channel was taken into

consideration. A continuous wave with long wavelength, λ, and speed, c, traveled on the plates [43] and was defined by

$$z = h(X, t) = a(1 + \eta(X, t))\tag{1}$$

where η(*x*, *t*) = φ*Sin*( <sup>2</sup><sup>π</sup> <sup>λ</sup> (*<sup>X</sup>* <sup>−</sup> *ct*)), <sup>φ</sup> <sup>=</sup> *<sup>b</sup> <sup>a</sup>* having an interval of (0, 1).

**Figure 1.** Geometry of the problem.
