Womersley Number Research Articles

Oscillating viscous flow past a streamwise linear array of circular cylinders.

This paper addresses the viscous flow developing about an array of equally spaced identical circular cylinders aligned with an incompressible fluid stream whose velocity oscillates periodically in time. The focus of the analysis is on harmonically oscillating flows with stroke lengths that are comparable to or smaller than the cylinder radius, such that the flow remains two-dimensional, time-periodic and symmetric with respect to the centreline. Specific consideration is given to the limit of asymptotically small stroke lengths, in which the flow is harmonic at leading order, with the first-order corrections exhibiting a steady-streaming component, which is computed here along with the accompanying Stokes drift. As in the familiar case of oscillating flow over a single cylinder, for small stroke lengths, the associated time-averaged Lagrangian velocity field, given by the sum of the steady-streaming and Stokes-drift components, displays recirculating vortices, which are quantified for different values of the two relevant controlling parameters, namely, the Womersley number and the ratio of the inter-cylinder distance to the cylinder radius. Comparisons with results of direct numerical simulations indicate that the description of the Lagrangian mean flow for infinitesimally small values of the stroke length remains reasonably accurate even when the stroke length is comparable to the cylinder radius. The numerical integrations are also used to quantify the streamwise flow rate induced by the presence of the cylinder array in cases where the periodic surrounding motion is driven by an anharmonic pressure gradient, a problem of interest in connection with the oscillating flow of cerebrospinal fluid around the nerve roots located along the spinal canal.

Fluid–structure interaction and flow sensing of primary cilia in oscillating fluid flows

This study systematically investigates the interaction between an oscillating flow and primary cilia using numerical simulations. The primary cilia are modeled as elastic filaments with rotatable basal ends to mimic real ciliary deflections. How some governing parameters [i.e., the peak Reynolds number (Repeak), Womersley number (Wo), cilium length (L*), and streamwise spacing interval (Ld*)] regulate fluid–cilia interaction is explored. Our results indicate that within a certain range, both the span of deflection (SD) and the maximal curvature increase with the increase in Repeak, L*, and Ld*, while they decrease as the Wo increases. Compared with other parameters, Ld* affects ciliary deflection less significantly and its impact becomes nearly negligible when the cilia are separated over twice their length. Three typical stretch states are captured. For primary cilia with a short or medium length, an increase in the SD is accompanied by a greater propagation distance of the location of the maximal tensile stress (LMTS). However, this is not the case for long cilia that protrude into 1/3 of the lumen, as the arising third stretch state may greatly suppress the LMTS propagation. Our study further confirms the role of primary cilia in decreasing the wall shear stress (WSS) and altering its oscillating feature. The WSS decrease is more significant for cilia undergoing a larger SD and/or when Ld* is reduced. For a constant Ld*, a larger SD corresponds to a more uneven oscillatory shear index distribution, and the affected (i.e., less oscillatory) region appears to greatly depend on Ld*.

Effect of a bend on vortex formation and evolution in a three-dimensional stenosed geometry during pulsatile flow

Stenosis at arterial bends alters hemodynamics and instigates abnormal disease progression. This configuration is addressed numerically by exploring pulsatile flow (Reynolds number Re = 300–1200; Womersley number Wo = 7.62–15.24) in arteries encountering bend angles of θ = 20°–60°. Individual influences of stenosis and bend on flow dynamics are investigated. Validations against particle image velocimetry experiments for Re = 800 and Wo = 7.62 are carried out in straight and 60° bend stenosed models. For Re = 300–800, the shear layer along the stenosis rolls up into a primary vortex, that is, constrained by the outer wall forming a secondary vortex. At Re = 1200, shear layers undergo instabilities along the post-stenotic region and develop new vortices that promote disturbances and induce asymmetries over the cross-plane flow structures. These features are not present in a straight stenosed tube, showing that the bend is responsible for flow distortion. During the pulsatile cycle, increasing bend angles intensify the size and strength of vortices, while these are suppressed at higher frequencies. A higher bend of 60° experiences large time-averaged wall shear stress and oscillatory loads. In time, wall loading spatially circumscribes the post-stenotic region followed by wall loading during cycle deceleration. These features are consistent with the skewing of a three-dimensional ring structure formed in a stenosed tube that evolves into disintegrated structures in the post-stenotic region. Overall, simulations reveal that strongly bent stenosed arteries experience aggravated oscillatory loading. In the biomedical context, such arterial geometries will require special attention.

Computational Modeling of Thermodynamical Pulsatile Flow with Uncertain Pressure in Hydrocephalus

The watery cerebrospinal fluid that flows in the subarachnoid space (SAS) surrounds the entire central nervous system via symmetrical thermo-solute flow. The significance of this study was to present a flexible simulation based on theoretical vivo inputs onto a mathematical framework to describe the interaction of cerebrospinal fluid circulation restricted to a pathological disorder. The pathophysiology disorder hydrocephalus is caused by an enormous excess of asymmetric fluid flow in the ventricular region. This fluid imposition increases the void space of its boundary wall (Pia mater). As a result, the dumping effect affects an inertial force in brain tissues. A mathematical model was developed to impose the thermal dynamics of hydrocephalus, in which solute transport constitutes the excess watery CSF fluid caused by hydrocephalus, in order to demonstrate perspective changes in ventricular spaces. This paper investigated brain porous spaces in order to strengthen the acceleration and thermal requirements in the CNS mechanism. To characterize neurological activities, a unique mathematical model that includes hydrodynamics and nutrient transport diffusivity was used. We present the analytical results based on physical experiments that use the novel Laplace method to determine the nutrients transported through permeable pia (brain) parenchyma with suitable pulsatile boundary conditions. This causes high CSF pressure and brain damage due to heat flux over the SAS boundary wall. As a result of the increased Schmidt number, the analysis of the hydrocephalus problem revealed an increase in permeability and drop in solute transport. A high-velocity profile caused a rise in thermal buoyancy (Grashof number). When the CSF velocity reached an extreme level, it indicated a higher Womersley number. Additionally, the present study compared a number of clinical studies for CSF amplitude and pressure. We validated the results by providing a decent justification with the clinical studies by appropriate field references.

Multi-scale analysis of concentration distribution in unsteady Couette–Poiseuille flows through a porous channel

A multiple-scale perturbation analysis is presented to analyse the two-dimensional concentration distribution of passive contaminant released in an incompressible viscous fluid flowing between two parallel plates filled with a porous medium. The flow is driven by the combined effect of the upper plate oscillation in its own plane moving with a constant velocity, and the periodic pressure gradient. Mei’s homogenization technique is used to find the concentration distribution up to third order, complemented with the dispersion coefficients for four different situations, namely, steady, pulsatile, oscillatory and the combined effect of all these. We observe that when the flow is under the combined effect of wall oscillation and pressure pulsation, then the respective frequency (Womersley number) and amplitude parameters oppose each other while influencing the dispersion coefficient. Our analysis reveals that for a fixed amplitude of oscillation and pulsation, the frequency of pressure pulsation has a stronger effect on the dispersion coefficient compared with the wall oscillation. On the other hand, when the Womersley number is kept fixed, amplitude of the wall oscillation dominates the pressure pulsation. This behaviour is more prominent for higher values of the Darcy number. The transverse concentration distribution and its dependency on porous medium parameters are also discussed in detail.

Visco-Elastic Fluid Model in an Inclined Porous Stenosed Artery with Slip Effect and Body Acceleration

The present paper analyzes an unsteady magnetohydrodynamic blood flow model of an visco-elastic fluid through an inclined porous stenosed artery with body acceleration and slip effect. Navier-Stokes equations have been used to describe the blood flow model. The governing equation of blood flow is solved by an analytic method by considering blood as an incompressible, visco-elastic fluid, and suspension of RBC’s in plasma. Axial velocity, blood acceleration, flow rate, and shear stress are derived numerically by using the finite Laplace and Hankel transformation and their inverse. The effect of parameters such as the visco-elasticity parameter, Womersley number, Hartmann number, inclination angle, parameter of slip, and body acceleration frequency is analyzed. Axial velocity reduces as the Hartmann number and visco-elasticity parameter enhance and it enhances with the enhancement of the slip parameter and inclination angle. The study is beneficial for finding the effect of slip parameter, porosity factor and Hartmann number when a human body is exposed to MRI and CT scan.

Effects of diffusive Reynolds number on electro-osmotic pulsating nanofluid flow

We examine pulsating electro-osmotic nanofluid flow phenomena in a microchannel with porous walls. The combined effect of the injected nanofluid velocity and ion diffusion coefficients on the electrical potential formation is considered. A novel boundary condition is introduced so as to examine the effects of electro-osmosis and frictional forces on thermal profiles and nanoparticle volume fractions of nanofluids. Being motivated by the experimental works of Kong et al. [Phys. Chem. Chem. Phys. 19, 7678 (2017).], this paper aims to extend the study of ion diffusivity in terms of diffusive Reynolds number on nanofluid temperature in the pulsating pressure gradient setting. The semi-analytic differential transform method is used to solve the physical equations, represented as coupled ordinary differential equations, with a special emphasis on the convergence of solutions, which is presented in terms of tables and graphs. The study shows that the nanofluid velocity, temperature, and mass concentration are strongly influenced by the ion diffusion coefficient and the frequency of pulsating pressure gradient. The diffusive Reynolds number significantly influences the electric potential distribution. The velocity and temperature show an increasing trend in terms of diminishing sensitivity parameter. However, nanoparticle concentration increases with an enhancement of the sensitivity parameter. Finally, velocity and temperature increase with a diminution of the Womersley number.

Analysis of pulsatile shear-thinning flows in rectangular channels

We present an analysis of pulsatile shear-thinning flows in two-dimensional channels. Characteristic viscosity is determined based on steady-state analysis of non-Newtonian flows and is used to non-dimensionalize the flow system by introducing the non-Newtonian Womersley number. Numerical analyses on various Carreau fluids reveal the existence of master curves related to the amplitude and phase lag of the flows, where the shape of the master curve is determined by the degree of shear-thinning. Furthermore, the pulsatile flow dynamics can be predicted accurately using pre-computed master curves, thus predicting shear-thinning pulsatile flow dynamics without explicit transient computations.

Womersley's solution for the measurement of volume flow rates in transient laminar flow tubes

The characterization of transient flows within the Reynolds number range of Re = 10–100 and the Womersley number range α = 0.8–10 is required for the ongoing development of green monopropellant thrusters. However, at the ml/min scale flow rates of interest, these measurements are outside the capabilities of current commercial flow meters. It is proposed here that transient flows under the required dynamic conditions can be characterized via Womersley's solution for transient flow in a rigid tube. This solution method requires only the measured transient pressure gradient within a controlled laminar flow section and can be accomplished using existing commercial pressure measurement hardware. Experiments were performed where flow similarity was maintained with the ultimate thruster characterization application, but the radius of the flow passage was increased so that flows could be simultaneously characterized by both the proposed solution method and a commercial ultrasonic flow meter. It was shown that across the range of interest, applying Womersley's solution to a measured pressure gradient was an effective method of transient flow characterization. Additionally, it was shown that non-periodic flows can be characterized except for the initial flow startup transient with solution convergence times approximating the analytical solution to starting flow in a pipe. While these results were expected due to an experimental design matching the assumptions required for Womersley's analytic solution, this work demonstrates that this method is practically feasible as novel instrumentation enabling previously unobtainable measurements.

Electroosmotic effect on two immiscible (conducting–non-conducting) fluids flowing in the porous channel under magnetic field

In this article, we investigate the motion of two-phase immiscible fluids in the channel with the porous medium in the presence of a magnetic field. Depending on the different levels of permeability, the channel is classified into two regions. The lower region (region-I) consists of non-conducting modified Casson fluid and the upper region (region-II) consist of conducting Newtonian fluid. The external electric field is applied in the region-II and the electroosmotic effect is produced at the wall and the interface region. Due to the interfacial effect, the non-conducting fluid was driven by conducting fluid. To compute the flow dynamic and electric potential nature, the momentum equation and Poisson–Boltzmann equation have been solved analytically in the state of periodic vibrations. The impact of physiological parameters such as Debye–Hückel parameter, electric field, magnetic field, Darcy numbers, viscosity ratio, Casson parameter and Womersley number has been investigated and the various behaviours were characterised. Based on numerical computations performed on the aforementioned parameters, we found that the Casson parameter, Darcy numbers and electric double layer (EDL) effects are predominant features. In the velocity profiles, we identified the flow of the fluid was enhanced by the augmented values of Darcy numbers, electric field (wall zeta potential), Debye–Hückel parameter and Casson parameter. The retardation flow occurred for the higher values of the magnetic field, viscosity ratio and electric field (interface zeta potential). Based on the conditional cases, the Womersley number possesses inconsistent behaviour. The depletion nature was found in the percentage flow rate of the Casson parameter. By increasing the value of the electric field and Casson parameter, the magnitude of wall shear stress to Darcy number increases and decreases for higher values of the magnetic field, viscosity ratio and pulsatile Reynolds number. It is first noticed that for any value of M1 and M2, Darcy number (Da2) plays a dual role in the upper wall shear stress. The present study leads to the pumping effect of gas pipelines or designing the microfluidic devices which are manipulated in biosensors/biochips.

Performance analysis of thermoacoustic plasma MHD generation

We develop a quasi-one-dimensional theoretical model describing the oscillatory plasma MHD conversion in a channel with parallel electrodes. Based on the model, a thermoacoustic plasma MHD generation system with a looped configuration is proposed and the characteristics of the operation of the system are investigated. Our results show that, for efficient oscillatory MHD generation, the channel dimension of the MHD duct should be much larger than the penetration depth (a proper Womersley number is above 100), otherwise the thermoacoustic effect will play a major role and consume considerable energy. Further, the electromagnetic induction parameter should be larger than 100S⋅T2/m, so that a majority of the kinetic energy of the oscillatory plasma can be harnessed for energy conversion. Moreover, a load ratio between 0.6 to 0.8 is suggested for the presented system, with which an efficiency above 20% can be achieved when the electromagnetic induction parameter is 500S⋅T2/m. In the optimal case, an efficiency of 24% can be obtained, with the corresponding output electricity of 1644 W. These results shed light on the feasibility of the closed-cycle MHD generation without any moving part, which endows the system with advantages of high reliability and long lifetime.

The Influence of Womersley Number on Non-Newtonian Effects: Transient Computational Study of Blood Rheology

Abstract A properly validated computational fluid dynamics methodology is a valuable predictive tool capable of aiding in the development of new methods for treating cardiovascular disease (CVD). Although blood is a shear-thinning non-Newtonian fluid, a key assumption that remains highly contested is whether non-Newtonian blood can be approximated as a Newtonian fluid. Previously, a preliminary link was established between the effects of non-Newtonian viscosity and the Womersley number, α, which could lend an explanation to the varied conclusions from previous comparison studies. Building upon this foundation, computational fluid dynamics was utilized to perform an in-depth investigation into the link between blood rheology and α for multiple geometries. For the first time in the open literature, the present research sheds definitive light on the source of the diverse results from previous studies. It demonstrates how α can affect the severity of non-Newtonian effects when compared to Newtonian viscosity, while otherwise maintaining the same boundary conditions. These results show that an increase in α reduces the peak global importance factor, a measure of the difference between Newtonian and non-Newtonian models, by upwards of 90%. Additionally, this results in a decrease in the relative difference for disturbed flow factors, parameters linked to the initiation and progression of CVD, from upwards of 34% down to approximately 5%. This study proves that there is a significant relationship between α and blood rheology, with higher α shifting the apparent viscosity of non-Newtonian models further toward the constant Newtonian viscosity.

Acoustic viscous boundary layer beyond laminar regime: Particle Image Velocimetry measurements

The current study reports an experimental investigation of flow regimes for acoustically oscillating flows in the near-wall region of a resonator. Particle Image Velocimetry (PIV) is used to measure axial and transverse velocity components inside an empty square resonator over the acoustic cycle for a wide range of Reynolds numbers and at a high Womersley number. The measurements are compared with LDV measurements in a circular duct. A good agreement between the PIV and LDV measurements is observed. The measured transverse distribution of the axial velocity agrees well with the Stokes boundary layer expectations up to the Reynolds number (based on acoustic velocity and viscous boundary layer depth) of 240. At higher Reynolds numbers, the measured velocity distribution deviates from the theoretical expectations, suggesting a change of regime at a Reynolds number between 240 and 272. To assert the critical Reynolds number, the velocity fluctuations and Reynolds stress are estimated. At Reynolds number of 272, an increase in the value of the Reynolds stress is observed at the edge of the viscous penetration depth which is regarded as an indication of changing the flow regime. To better qualify the high amplitude acoustic flow, the spectrum of the turbulent kinetic energy and evolution of the integral length scale over the acoustic cycle are investigated. The decay of turbulent kinetic energy with the wavenumber complies with the universal slope (−5/3). The integral length scale increases during the deceleration stage and decreases at the beginning of the acceleration stage.

On Stokes' second problem solutions in cylindrical and Cartesian domains

It is well known that drag created by turbulent flow over a surface can be reduced by oscillating the surface in the direction transverse to the mean flow. Efforts to understand the mechanism by which this occurs often apply the solution for laminar flow in the infinite half-space over a planar, oscillating wall (Stokes' second problem) through the viscous and buffer layer of the streamwise turbulent flow. This approach is used for flows having planar surfaces, such as channel flow, and flows over curved surfaces, such as the interior of round pipes. However, surface curvature introduces an additional effect that can be significant, especially when the viscous region is not small compared to the pipe radius. The exact solutions for flow over transversely oscillating walls in a laminar pipe and planar channel flow are compared to the solution of Stokes' second problem to determine the effects of wall curvature and/or finite domain size. It is shown that a single non-dimensional parameter, the Womersley number, can be used to scale these effects and that both effects become small at a Womersley number of greater than about 6.51, which is the Womersley number based on the thickness of the Stokes' layer of the classical solution.

An experimental study on heat transfer of pulsating air flow in metal foam subjected to constant heat flux

In this study, the effect of pulsating air flow on heat transfer was investigated experimentally, while constant heat flux was applied with an electric heater to the side surface of the cylindrical duct with metal foam inside. Reynolds numbers in the range of 5932–8592 were studied at three different fan speeds. Pulsating flow experiments were carried out in the range of 12.0–27.9 Womersley numbers. It has been observed that the surface and air temperatures are not affected by the flow frequencies in this range and are independent of time. With the data obtained, Nusselt number values were calculated for the fully developed region, and the Nusselt number ratio was defined in order to compare the steady flow with the pulsating flow. Nusselt correlation was established for the fully developed region depending on the Reynolds number and the Womersley number in the studied range. It has been determined that the effect of frequency on the heat transfer in the case of pulsating flow is weak, but the heat transfer tends to decrease as the frequency increases.

Linear instabilities of pulsatile plane channel flow between compliant walls

The linear dynamics of perturbations developing in an infinite channel with compliant walls is investigated for pulsatile flow conditions. Two-dimensional modal perturbations are considered for Womersley-type pulsating base flows and the wall motion is only allowed in the normal direction. It is found that the flow dynamics is mainly governed by four control parameters: the Reynolds number $Re$ , the reduced velocity $V_R$ , the Womersley number $Wo$ and the amplitude of the base-flow modulation $\tilde Q$ . Linear stability analyses are carried out within the framework of Floquet theory, implementing an efficient approach for removing spurious eigenmodes. The characteristics of flow-based (Tollmien–Schlichting) and wall-based (both travelling-wave flutter and divergence) modes are investigated over a large control-parameter space. It is shown that travelling-wave flutter (TWF) modes are predominantly influenced by the reduced velocity and that the Reynolds number has only a marginal effect. The critical reduced velocity (corresponding to onset of linear instability) is demonstrated to depend both on the Womersley number and modulation amplitude for a given set of wall parameters. Similarly to the steady flow case, the Tollmien–Schlichting (TS) mode is also found to be only weakly affected by the flexibility of the wall in pulsatile flow conditions. Finally, the classification given by Benjamin (J. Fluid Mech., vol. 16, 1963, pp. 436–450) is found to be too restrictive in the case of pulsatile base flows. In particular, a new type of mode is identified that shares characteristics of two distinct Floquet eigenmodes: TS and TWF modes. Due to coupling of the different Floquet harmonics, a phenomenon specific to time-periodic base flows, this two-wave mode exhibits a beating over the intracyclic dynamics.

A computational conjugate heat transfer study of a rectangular minichannel undergoing sinusoidal flow pulsations

The development of current and next generation high performance electronic devices has led to smaller components in more densely packed spaces. The increasing power levels have resulted in ever-increasing heat flux densities which necessitates the evolution of new liquid-based heat exchange technologies. Pulsating flow in single-phase cooling systems is viewed as a potential solution to the problems involving high heat flux densities. A review of published literature indicates a lack of time-resolved and space-resolved links between the hydrodynamic pulsating characteristics and associated heat transfer perturbations. The scope of this work involves the development of a validated three-dimensional conjugate heat transfer computational model to investigate hydrodynamically and thermally fully developed pulsating flows in a heated rectangular minichannel. Simulations were performed for a sinusoidal waveform with a fixed pulsation amplitude for varying pulsation frequencies in the range of 0.02 Hz to 25 Hz, corresponding to Womersley numbers in the range of 0.5 ≤Wo≤ 18.33. Low pulsation frequencies exhibited the well known parabolic profile for the fluctuating hydrodynamic and thermal parameters, i.e., velocity, wall shear, and wall temperature. As a result, the axial pressure gradient, velocity, and wall shear stress profiles were in phase and similar results were obtained for the oscillating wall and bulk fluid temperatures. For the inertia dominated high frequency flows, an increase in axial pressure gradient leads to a phase lag of π/2 when compared with the velocity and wall shear profiles. The shorter time period pulsations exhibit unique attributes in the form of flow reversal effects at local near wall regions. High near wall thermal gradients were observed as a result of stronger viscous effects due to a narrowing thermal boundary layer; consequently the transverse diffusion of heat was ineffective. A phase lag and a subsequent drop in the peak magnitudes existed between the oscillating bulk and wall temperatures for high frequency flows. Fluctuations in near wall heat flux profiles showed a dependency on the imposed pulsation frequencies. For the chosen pulsation profile and frequencies, the overall time averaged thermal performance indicates that pulsating flow performs worse than steady flow for a flow rate amplitude of 1. The highest thermal performance was achieved for Wo=5.1 while maintaining a low friction factor.

Onset of nonlinearity in oscillatory flow through a hexagonal sphere pack

We simulated laminar flow through a hexagonal sphere pack driven by a sinusoidal volume force using direct numerical simulation. We vary two independent parameters, the Hagen and Womersley numbers, representing the amplitude and frequency of the forcing. First, we determine for which regions in the parameter space nonlinear effects have to be considered. We judge the presence of nonlinear effects from the departure of the superficial velocity and kinetic energy from a linear behaviour as well as from the presence of higher harmonics in the discrete Fourier transform of the velocity field. We discuss the asymptotic behaviour of the onset of nonlinearity in the limits of low and high Womersley number, and we delineate approximately the parameter region that can be described using the linear theory. Second, we document the changes of instantaneous velocity fields with Hagen and Womersley numbers. We show that the onset of nonlinearity is accompanied by a loss of fore–aft symmetry of the flow, and subsequently, we employ the deviation from this symmetry to quantify the strength of nonlinear effects in the instantaneous velocity fields. Based on this analysis, we demonstrate that for higher Womersley numbers, the strongest nonlinear effects occur during the deceleration of the superficial velocity; consequently, the development of the nonlinearity is not in phase with the superficial velocity. Finally, we describe the leading-order nonlinear effects in the frequency domain and the interaction among the nonlinear Fourier modes that leads to a temporal variation in the strength of nonlinear effects.

Pulsatile Flow Through a Curved Pipe – An Analytical Study

In this paper, the unsteady flow of Newtonian fluid through a curved pipe due to a pulsatile pressure gradient has been considered. The flow is three dimensional and the partial differential equations governing the flow are highly coupled and non-linear. Approximate analytical solutions of the governing partial differential equations have been obtained without neglecting any term containing the curvature ratio. Perturbation series in terms of curvature ratio has been used for obtaining the solutions. It is interesting to note that the solutions are valid for large as well as small values of Womersley number. The effect of different parameters such as, Reynolds number, Womersley number and curvature ratio on the flow through curved pipe is discussed in this paper. It is found that the axial velocity is qualitatively periodic in nature, as expected. The Reynolds number and curvature ratio are found to shift the axial velocity towards the outer boundary of the curved pipe.

Mass Transfer Effects on the Mucus Fluid with Pulsatile Flow Influence of the Electromagnetic Field

The influence of pulsatile flow on the oscillatory motion of an incompressible conducting boundary layer mucus fluid flowing through porous media in a channel with elastic walls is investigated. The oscillatory flow is treated as a cyclical time-dependent flux. The Laplace transform method using the Womersley number is used to solve non-linear equations controlling the motion through porous media under the influence of an electromagnetic field. The theoretical pulsatile flow of two liquid phase concurrent fluid streams, one kinematic and the other viscoelastic, is investigated in this study. To extend the model for various physiological fluids, we postulate that the viscoelastic fluid has several distinct periods. We also apply our analytical findings to mucus and airflow in the airways, identifying the wavelength that increases dynamic mucus permeability. The microorganism’s thickness, velocity, energy, molecular diffusion, skin friction, Nusselt number, Sherwood number, and Hartmann number are evaluated. Discussion is also supplied in various sections to investigate the mucosal flow process.