Search

The exploration addresses the effect of variable viscosity and thermal conductivity on the peristaltic mechanism of Bingham fluid. A two-dimensional non-uniform porous channel is considered for the fluid flow, which is assumed to be inclined. The impact of heat, slip conditions, wall properties, homogeneous and heterogeneous reactions are examined. The resulting nonlinear differential equations are solved by employing the perturbation method. The solutions acquired are analyzed and sketched through graphs that show that the variable viscosity renders a critical role in regulating the velocity of the fluid in the channel's central part. The stream function has been analyzed to observe the trapping phenomenon. Further, the obtained results find its application in understanding the flow of blood in micro arteries.

The rheological behaviors of highly concentrated fine particle suspensions (clay-silt-water mixtures) and coarse particle suspensions (coarse particles within a fine particle suspension) were investigated in this study. Experimental results demonstrated that the Bingham Fluid Model with two rheological parameters, Bingham yield stress and viscosity, well characterized the rheological behavior of fine particle suspensions at shear rates between 4 and 20 s−1. The inclusion of coarse particles within a fine particle suspension induced an enhancement to the rheological parameters. The rheological parameters of a coarse particle suspension not only depend on its total particle fraction but also on its relative fine/coarse particle fractions. Empirical equations of these two parameters were proposed, quantitatively related to both fine and coarse particle fractions. Results indicated that the Bingham yield stress and viscosity are much more (an order larger) sensitive to changes in fine particle fraction than to changes in coarse particle fraction.

An artery narrowing referred to as atherosclerosis or stenosis causes a reduction in the diameter of the artery. When blood flow through an artery consists of stenosis, the issue of solute dispersion is more challenging to solve. A mathematical model is developed to examine the unsteady solute dispersion in an overlapping stenosed artery portraying blood as Bingham fluid model. The governing of the momentum equation and the constitutive equation is solved analytically. The generalized dispersion model is imposed to solve the convective-diffusion equation and to describe the entire dispersion process. The dispersion function at steady-state decreases at the center of an artery as the stenosis height increase. A reverse behavior is shown at an unsteady-state. As the plug core radius, time and stenosis height increase, the dispersion function decreases at the center of an artery. There is a high amount of red blood cells at the center of the artery but no influences near the wall. Hence, this model is useful in transporting the drug or nutrients to the targeted stenosed region in the treatment of diseases and in understanding many physiological processes.

A crucial part in a two-dimensional discrete fracture network (DFN) model of two-phase flow in a naturally fractured rock is the updating algorithm for the fluid interface positions inside the fractures. In this study, a single fracture situated in an impermeable rock matrix is considered. Displacements of a Newtonian fluid by another Newtonian fluid, by a power-law fluid, or by a Bingham fluid are analysed. It is shown that a significant improvement in accuracy and computational efficiency in all three types of displacement can be achieved by updating the interface position with predictor–corrector scheme instead of forward Euler scheme currently used in such simulators. In order to achieve the same accuracy in the breakthrough time as with the forward Euler method, predictor–corrector can be run with a timestep that is 3 to 50 times larger. The gain in accuracy (or, equivalently, increase of the adjustable timestep) achieved with predictor–corrector is shown to outweigh the extra computational cost of this scheme. It is also shown that further increase in accuracy, e.g. by using the Runge–Kutta 4 scheme, does not pay off. In order to have control over the error in the interface positions, a viable strategy is to choose the timestep equal to a given fraction of the shortest single-fracture breakthrough time in the network. The value of this fraction can be chosen based on the required accuracy in the breakthrough time. In the case of Newtonian/Newtonian displacement, the optimal ratio of timestep to the single-fracture breakthrough time is found to depend only on the mobility ratio and the prescribed accuracy in the numerically computed breakthrough time. In case of a single interface in the fracture, numerical experimentation yields simple equations for recommended timestep in Newtonian/Newtonian displacement based on the fluid mobility ratio and the maximum allowed error in the breakthrough time.

In this paper aims, we found the fluid concentration after calculating the velocity and temperature of the fluid with a variable viscosity that depends on the fluid moving through an inclined porous channel. We examined the influences of certain parameters that are active on fluid velocity by analyzing the graphs obtained after we reached the momentum equation solution, and used the MATHEMATICA program for plot the velocity and temperature of the fluid for two types of flow (Poiseuille and Couette).

Recently, Zou et al. (2020a) published a theoretical analysis on the radial flow of a Bingham fluid, where they argued that the classical analysis by Dai and Bird (1981) violates the mass conservation. The present discussion aims to clarify this conflict between those two studies. It is noted that Zou et al. (2020a) presumed the gap-wise mass flux is negligible in the mass conservation equation, while Dai and Bird (1981) did not require so in their model, and this is found to be the origin of the conflict. In fact, Dai and Bird (1981)’s model is shown to not violate the mass conservation. Therefore, those two models should be viewed as separate models derived from different perspectives. Details of the major difference between the two models are discussed.

In this study, we performed direct simulations of particle suspensions in a simple shear flow with a viscoplastic model (Bingham fluid) using a finite-element/fictitious-domain method, as a model system for cuttings transport with a drilling mud in the oil and gas exploration. Bulk rheology and particle interaction were investigated via several example problems including single-particle, two-interacting particle, and many-particle problems for various solid fractions in both Bingham and Newtonian fluids. We report that, at low shear rate, huge reductions occur in both the relative bulk viscosity and angular velocity of particle suspension in a Bingham fluid compared to that in a Newtonian fluid. This indicates that, although particle incorporation increases bulk viscosity in a Bingham fluid, the amount of increase appears much smaller than that with the Newtonian medium, and particle motions were restricted to some extend at a low shear rate with a yield stress fluid. Moreover, the suppression in viscosity increase with particles disappears at a high shear rate, just like Newtonian particle suspension. This phenomenon has been interpreted by the suppressed particle angular rotation in a Bingham fluid. In addition, flow fields with low shear regions, where the viscosity of a Bingham fluid is high and the fluid mobility decreased significantly, on the left and right sides of a particle also confirm this interpretation. The two-interacting particle problem reveals delayed particle interaction with increased time period in particle trajectories (drafting, kissing, tumbling and separation) in a Bingham fluid. Similar behaviors with the suppression in viscosity increase can be also observed in the many-particle problem.

Avalanche is a natural disaster in the snow-covered mountainous area in winter, which may cause great disasters to human life and property. It is also a danger for skiers and climbers. This paper presents a new physically based algorithm to simulate the dynamic avalanches under position-based dynamics framework. To realistically simulate avalanches’ dynamic characteristics, we introduce the Bingham plastic model from geodynamics to model snow flow motion in avalanches. The interaction between snow flow in the avalanche and the surrounding objects is simulated by a level set-based two-way fluid–solid coupling model. We also propose static and kinetic friction mixed model to determine the accumulated transition of the avalanche. To create an avalanche scene with more realistic details, we employ an aerodynamics-based snow drag force model to generate snow fog effect. Finally, by choosing different criterion shear rate and friction parameters, different kinds of wet and dry avalanche scenes are realistically rendered. Compared with the real photographs of avalanches, our simulated results are quite satisfactory.

In nature, many fluid-like materials exhibit a yield stress below which they behave like a solid. The Bingham model aims to describe such materials. This paper draws some mathematical considerations on the flow of a Bingham fluid in a vertical channel. The situation due to the presence of an external magnetic field and natural convection is analyzed: the external magnetic field, which is orthogonal to the walls of the channel, generates the Lorentz forces that influence the motion through the Hartmann number. The behavior of the velocity, the induced magnetic field and the thickness of the plug regions are discussed and presented graphically. We find that the velocity is a decreasing function of the Bingham and Hartmann numbers. In particular, the presence of the external magnetic field increases the thickness of the plug region. The modulus of the induced magnetic field is not monotone when the Hartmann number changes, but it is a decreasing function of the Bingham number.

The paper examined a theoretical investigation of the stimulus of mass and heat transfer on the channel's peristaltic utilization of the MHD Bingham liquid. The research study is motivated to explore blood circulation in the little vessels by taking the slip, variable thermal conductivity, and thickness of the wall features into account. The leading constitutive equations are established based on low Reynolds number and approximations for long wavelengths. The solution for the resulting nonlinear energy and momentum equations is obtained using a semi-analytic method, while the exact solution for the concentration field is obtained. Using the MATLAB software, the influences of different constraints on the interest of physiological quantities are represented graphically. One of the considerable outcomes of the current model exposes that the existence of variable fluid properties boosts the rate as well as temperature level areas. The rheological and flow properties of various biological fluids can be derived from this model as a particular case. Moreover, the formation of stuck bolus diminishes for larger values of the magnetic and velocity slip constraints.

Purpose The purpose of this study is to solve convective diffusion equation analytically by considering appropriate boundary conditions and using the Taylor-Aris method to determine the solute concentration, the effective and relative axial diffusivities. Design/methodology/approach >An analysis has been conducted on how body acceleration affects the dispersion of a solute in blood flow, which is known as a Bingham fluid, within an artery. To solve the system of differential equations analytically while validating the target boundary conditions, the blood velocity is obtained. Findings The blood velocity is impacted by the presence of body acceleration, as well as the yield stress associated with Casson fluid and as such, the process of dispersing the solute is distracted. It graphically illustrates how the blood velocity and the process of solute dispersion are affected by various factors, including the amplitude and lead angle of body acceleration, the yield stress, the gradient of pressure and the Peclet number. Originality/value It is witnessed that the blood velocity, the solute concentration and also the effective and relative axial diffusivities experience a drop when either of the amplitude, lead angle or the yield stress rises.

Accurate pressure drop estimation is important for drill string and bit nozzles design and optimized fluid circulations as well as identifying the drilling problems such as bit nozzle(s) washout or plugging. In this study, the Bingham Plastic model has been modified by applying a coefficient to its turbulent pressure loss calculations. This coefficient encompasses the effects of the drill pipe tool joints and other effects in estimation of pressure losses. The range of the coefficient was determined in field applications for different hole sizes and mud types. The results showed that applying a correction coefficient of 1.08–1.12 to turbulent pressure loss equations (depending on borehole size and mud type) improves the pressure loss estimation. By applying this coefficient, the estimated pressure losses are increased to compensate the under-estimation of the Bingham Plastic model. This is considered a significant contribution to accurate calculation of borehole hydraulics and in-time detection and identification of borehole problems and reduction of invisible lost time. The findings also showed that this enhanced effect is independent of the mud type. The use of this coefficient removes the necessity of using rather complex mud rheological models such as the Herschel–Bulkley model.

Present analysis is conducted for Bodewadt flow along a permeable plane surface immersed in a Bingham fluid, which represents a class of non-Newtonian fluids exhibiting yield stress. The resulting heat transfer phenomenon is analyzed in the existence of viscous dissipation. After implementing boundary layer approximations, von-Karman substitutions are employed resulting in a coupled non-linear system. For the numerical solution, a reliable numerical routine bvp4c based on collocation method, is adopted. Flow model contains a parameter Bn , measuring the relative importance of fluid yield stress. Boundary layer is markedly risen as a consequence of fluid yield stress. Viscous dissipation also remarkably influences the resulting heat transfer process. Subtle fluid dynamics quantities including normalized skin friction, disk cooling rate and volumetric flow rate are scrutinized under the influence of fluid yield stress. The results predict a deceleration in entrained flow whenever fluid exhibits yield stress.

We have developed and tested a numerical model for turbulence resolving simulations of dense mud–water mixtures in oscillatory bottom boundary layers, based on a low Stokes number formulation of the two-phase equations. The resulting non-Boussinesq equation for the fluid momentum is coupled to a transport equation for the mud volumetric concentration, giving rise to a volume-averaged fluid velocity that is non-solenoidal, and the model was implemented as a new compressible flow solver. An oscillating pressure gradient force was implemented in the correction step of the standard semi-implicit method for pressure linked equations (SIMPLE), for consistency with the treatment of other volume forces (e.g., gravity). The flow solver was further coupled to a new library for Bingham plastic materials, in order to model the rheological properties of dense mud mixtures using empirically determined concentration-dependent yield stress and viscosity. We present three direct numerical simulation tests to validate the new MudMixtureFoam solver against previous numerical solutions and experimental data. The first considered steady flow of Bingham plastic fluid with uniform concentration around a sphere, with Bingham numbers ranging from 1 to 100 and Reynolds numbers ranging from 0.1 to 100. The second considered the development of turbulence in oscillatory bottom boundary layer flow, and showed the formation of an intermittently turbulent layer with peak velocity perturbations exceeding 10 percent of the freestream flow velocity and occurring at a distance from the bottom comparable to the Stokes boundary layer thickness. The third considered the effects of density stratification due to resuspended sediment on turbulence in oscillatory bottom boundary layer flow with a bulk Richardson number of 1×10−4 and a Stokes–Reynolds number of 1000, and showed the formation of a lutocline between 20 and 40 Stokes boundary layer depths. In all cases, the new solver produced excellent agreement with the previous results.

Quick-setting grout has time-varying characteristics, and the space porosity of the injected medium changes after being grouted. The above two points need to be considered when loose and broken coal rock masses are grouted with rapid-setting grout. In other words, the coupling effect of the slurry-coal rock mass causes the slurry to exhibit the temporal and spatial characteristics of slurry viscosity (TSCSV). Based on the Bingham fluid constitutive model, a capillary group model with a linear increase in diameter was created, and a one-dimensional penetration grout diffusion model was established considering the TSCSV. The one-dimensional visualization splitting-permeation grout diffusion simulation test system was designed to obtain the space porosity of the broken coal rock mass under different final grouting pressures to obtain the law of the changes in grouting pressure and diffusion height with time under the condition of a fixed grouting rate. Furthermore, the mechanism of splitting-permeation grouting in loose and broken coal rock masses is analysed in depth considering the TSCSV. The results showed that the space porosity of coal rock masses increases linearly. In addition, the calculated value of the slurry diffusion distance was 0.93 – 1.1 times the experimental value considering the TSCSV. The calculated value was in good agreement with the experimental value. The calculated value of the final pressure of orifice grouting was approximately 1.5 times the test value without considering the unevenness of the space porosity distribution (USPD). The grouting pressure was overestimated without considering the USPD. Therefore, the TSCSV of the quick-setting slurry should be fully considered in the grouting design.

In this work, the momentum and heat transfer characteristics of a time-dependent flow of Bingham plastic fluid over a heated isothermal sphere have been investigated numerically over wide ranges of conditions: Reynolds number, 5 ≤ Re ≤ 120, Bingham number, 0.1 ≤ Bn ≤ 100, Prandtl number, 0.7 ≤ Pr ≤ 100, frequency, π/4 ≤ ω* ≤ π and amplitude, 0 ≤ A ≤ 0.8. The influence of the fluid yield stress and inertia due to the imposed flow pulsations on the flow and thermal fields have been examined in detail. Detailed structure of the flow and temperature fields are analysed in terms of the instantaneous streamlines, isothermal contours, yielded/unyielded zones, surface pressure profiles, time-average drag coefficient and surface- and time-average Nusselt number. The influence of the frequency and amplitude of pulsations on the size of yielded (fluid-like) and unyielded (solid-like) zones is considered in order to understand convective heat transport. The size of yielded zones is seen to be in phase with the imposed pulsating velocity. However, the yield stress effects suppress the influence of flow pulsations. The temporal evolution of the drag coefficient and Nusselt number lag the imposed pulsating flow by different degrees thereby indicating the inherently different evolution of the momentum and thermal boundary layers. Broadly, the pulsating flow conditions may lead to up to 20% augmentation in heat transfer provided there is a moderate degree of advection and/or when the fluid yield stress effects are weak, i.e., small Bingham numbers. Thus, the maximum benefits of pulsating flow accrue in Newtonian fluids only. Finally, the present values of the time-average Nusselt number have been consolidated in the form of a predictive expression.

Grouting is a commonly used technique in rock engineering to enhance the joint strength and improve the stability of surrounding rock. Grout penetration characteristic is controlled by grouting parameters and has a significant role on practice. In the study, a numerical manifold method (NMM) for grout penetration process simulation in fractured rock mass is firstly proposed. The fluid flow behaviour of the grout is assumed to be a Bingham fluid and control equations are established using discrete fracture network model. The global discretization equation, element sub-matrixes and NMM simulation algorithm for grouting are presented. Then, numerical tests for grouting process in a single fracture and a regular fracture network are conducted firstly to verify the proposed NMM grouting model by comparing with analytical solutions and experimental results. Furthermore, the effects of mesh size, fracture and slurry parameters on the grouting performance are systematically investigated using a random fracture network example. The numerical results indicate that the grouted zone and propagation depth decrease as the mesh size of numerical model and yield strength increases, while it increases as initial fracture aperture and grouting pressure increases.

The superimposition of mechanical oscillations can lead to a remarkable modification in the friction factor (Cf) with a fixed flowrate for viscous non-Newtonian fluids in a pipe but cannot affect the flow of Newtonian types. Therefore, less driving pressure is required for fluids with a shear-thinning nature, while more for shear-thickening types. A validated computational fluid dynamics (CFD) model is employed to compare the efficacy of three modes of oscillation: longitudinal, transverse, and rotational within the laminar flow region. Rheologies considered here are power-law fluids and Bingham plastic types. At a given Reynolds number (Re), Cf is found to be a function of rheological parameters (i.e., flow behavior index, n and Bingham number, Bn) and the vibration intensity, λ. Longitudinal oscillation is the most effective, with transverse oscillation also producing comparable results. Rotational oscillation causes the least effect on the viscosity and therefore, the friction factor is even nearly unchanged while λ ≤ 20 for power-law fluids, and λ ≤ 5 for Bingham plastic rheologies. Expressions for Cf are proposed, giving a quantitative estimation on the effect of mechanical vibration on the driving pressure.

This paper discusses the well posedness of an initial value problem describing the motion of a Bingham fluid in a basin with a degenerate bottom topography. A physical interpretation of such motion is discussed. The system governing such motion is obtained from the Shallow Water-Bingham models in the regime where the Froude number degenerates, i.e taking the limit of such equations as the Froude number tends to zero. Since we are considering equations with degenerate coefficients, then we shall work with weighted Sobolev spaces in order to establish the existence of a weak solution. In order to overcome the difficulty of the discontinuity in Bingham's constitutive law, we follow a similar approach to that introduced in [G. DUVAUT and J.-L. LIONS, Springer-Verlag, 1976]. We study also the behavior of this solution when the yield limit vanishes. Finally, a numerical scheme for the system in 1D is furnished.

Accurate determination of hydraulic parameters such as pressure losses, equivalent circulation density (ECD), etc. plays profound roles in drilling, cementing and other well operations. Hydraulics characterization requires that all factors are considered as the neglect of any could become potential sources of errors that would be detrimental to the overall well operation. Drilling Hydraulics has been extensively treated in the literature. However, these works almost entirely rely on the assumption that the drill string lies perfectly at the center of the annulus—the so-called “concentric annulus”. In reality, concentricity is almost never achieved even when centralizers are used. This is because of high well inclination angles and different string geometries. Thus, eccentricity exists in practical oil and gas wells especially horizontal and extended reach wells (ERWs) and must be accounted for. The prevalence of drillstring (DS) eccentricity in the annulus calls for a re-evaluation of existing hydraulic models. This study evaluates the effect of drilling fluid rheology types and DS eccentricity on the entire drilling hydraulics. Three non-Newtonian fluid models were analyzed, viz: Herschel Bulkley, power law and Bingham plastic models. From the results, it was observed that while power law and Bingham plastic models gave the upper and lower hydraulic values, Herschel Bulkley fluid model gave annular pressure loss (APL) and ECD values that fall between the upper and lower values and provide a better fit to the hydraulic data than power law and Bingham plastic fluids. Furthermore, analysis of annular eccentricity reveals that APLs and ECD decrease with an increase in DS eccentricity. Pressure loss reduction of more than 50% was predicted for the fully eccentric case for Herschel Bulkley fluids. Thus, DS eccentricity must be fully considered during well planning and hydraulics designs.

In this paper, experimental experiments were conducted using advanced computerized rheometer. The lean solid slurry (10 to 40 percent by weight) requires an immense amount of water and requires high-energy pumping power. The rheological properties are highly dependent on shear strength, solid density and particle size distribution. The plastic viscosity and yield stress enhanced with a solid concentration. At higher concentrations of solids, the experimental results fitted well with the Bingham plastic model. The results reveal that the plastic viscosity and yield stress of the slurry samples are increased due to finer particles. Bingham plastic viscosity depends primarily on solid concentration and particle size.

The paper is devoted to the study of convective motions of viscoplastic fluids in closed two-dimensional rectangular domains with different aspect ratios under lateral heating. The problem was solved numerically using the ANSYS Fluent software. The Herschel-Bulkley model was chosen for describing the rheological behavior of the fluid. Under certain rheological parameters, this model was transformed into a Newtonian fluid model, whose behavior was also simulated as a limiting case. Based on the calculated results, the dependences of the maximum value of the stream function in the cavity on the Rayleigh number are plotted. It was found that at small Rayleigh numbers the intensity of motion is close to zero. At a certain threshold value of the Rayleigh number there is a sharp increase in the intensity of motion, and a further increase in the Rayleigh number causes an almost linear growth of the maximum value of the stream function. For each of the considered ratios of cavity sides, the "threshold" values of the Rayleigh number, at which a sharp increase in the intensity of motion of liquid is observed, were determined. The obtained values of Rayleigh numbers turned out to be close to the "threshold" values of the Rayleigh number for the Bingham fluid found in the works of other authors earlier. The fields of the stream function and the square root of the second invariant of the viscous stress tensor were obtained for different values of the Rayleigh number and different aspect ratios. The scenarios of rearrangement of quasisolid motion zones with increasing Rayleigh number were compared with the available results for the Bingham fluid.

Modeling of fluids with complex rheology in the lattice Boltzmann method (LBM) is typically realized through the introduction of an effective viscosity. For fluids with a yield stress behavior, such as so-called Bingham fluids, the effective viscosity has a singularity for low shear rates and may become negative. This is typically avoided by regularization such as Papanastasiou’s method. Here we argue that the effective viscosity model can be re-interpreted as a generalized equilibrium in which no violation of the stability constraint is observed. We implement a Bingham fluid model in a three-dimensional cumulant lattice Boltzmann framework and compare the direct analytic effective viscosity/generalized equilibrium method to the iterative approach first introduced by Vikhansky which avoids the singularity in viscosity that can arise in the analytic method. We find that both methods obtain similar results at coarse resolutions. However, at higher resolutions the accuracy of the regularized method levels off while the accuracy of the direct method continuously improves. We find that the accuracy of the proposed direct method is not limited by the singularity in viscosity indicating that a regularization is not strictly necessary.

This framework presents heat transfer analysis for swirling flow of viscoplastic fluid bounded by a permeable rotating disk. Problem formulation is made through constitutive relations of Bingham fluid model. Viscous dissipation effects are pre-served in the mathematical model. Entropy production analysis is made which is yet to be explored for the von-Karman flow of non-Newtonian fluids. Having found the similarity equations, these have been dealt numerically for broad parameter values. The solutions are remarkably influenced by wall suction parameter and Bingham number which measures the fluid yield stress. Akin to earlier numerical results, thermal boundary-layer suppresses upon increasing wall suction velocity. Thermal penetration depth is much enhanced when fluid yield stress becomes large. Higher heat transfer rate can be accomplished by employing higher suction velocity at the disk. However, deterioration in heat transfer is anticipated as fluid yield stress enlarges. Current numerical results are in perfect line with those of an existing article in limiting sense.

We study the onset of convection within a rectangular porous cavity which is saturated with a Bingham fluid and subjected to a uniform internal heat generation. When such a cavity is saturated by a Newtonian fluid, convection takes place at all nonzero values of the Darcy–Rayleigh number, Ra. In such cases convection takes the form of two contrarotating cells with flow down the cold sidewalls when Ra first increases from zero. However, when the enclosure is saturated by a Bingham fluid, then we find that the cavity remains stagnant until the Darcy–Rayleigh number is sufficiently large that buoyancy overcomes the yield threshold. Numerical solutions are obtained using a second-order accurate finite-difference methodology where convergence is accelerated using line relaxation. The presence of the yield surfaces, which mark the boundaries of stagnant regions, is modelled by means of a regularization of the yield threshold. It is found that the critical value of Ra above which convection arises depends roughly linearly on the value of Rb, which may be described as a convective porous Bingham number. When the cavity has a sufficiently large aspect ratio, the fluid admits more than one stable steady-state solution.

We present rheological and mechanical behaviors of a partially cured photopolymer. When an ultraviolet (UV)-curable resin is exposed to UV light in atmospheric conditions, a partially cured layer is formed on the top of the resin owing to inhibitory effects of oxygen. Interestingly, such a partially cured resin behaves like a Bingham plastic with a yield stress, being a rigid solid at low shear stress and a viscous liquid at high stress. Unlike typical Bingham plastic materials, however, deformation rate saturation is observed with an increase in applied stress, which is attributed to the gradient in the degree of photopolymerization of the resin (termed “gradient Bingham plastic”). This gradient Bingham plastic can be utilized for the robust fabrication of diverse 3D, multiscale structures.

Understanding flow of non-Newtonian fluids in porous media is critical to successful operation of several important processes, including polymer flooding, filtration, and food processing. In some cases, such as when a non-Newtonian fluid can be represented by the power-law model, simulation of its flow in a porous medium is a straightforward extension of that of Newtonian fluids. In other cases, such as flow of Bingham fluids, there is a minimum external threshold pressure below which there would be no macroscopic flow in the porous medium. Computing the threshold pressure is a difficult problem, however. We present a new algorithm for determining the threshold pressure for flow of a Bingham fluid through a porous medium, modeled by a pore-network (PN) model. The algorithm, the ant colony optimization (ACO), is described in detail and together with the PN model is used to determine the minimum pressure for flow of Bingham fluids in a heterogeneous porous medium, the Mt. Simon sandstone, whose PN and morphological properties were extracted from the sandstone’s image. To assess the accuracy and computational efficiency of the ACO algorithm, we also carry out the same computations with two previous methods, namely invasion percolation with memory (IPM) and the path of minimum pressure (PMP) algorithms. The IPM does not guarantee identification of the optimal flow path with the minimum threshold pressure, while the PMP algorithm provides an approximate, albeit accurate, solution of the problem. We also compare the computational complexity of the three methods. For large PNs, both the IPM and PMP are much less efficient than the ACO algorithm. Finally, we study the effect of the morphology of the pore space on the minimum threshold pressure.

Solutions for radial flow of a Bingham fluid are analyzed in this paper. It aims to eliminate confusions in the literature concerning the plug flow region in different solutions for analysis and design of grouting in rock fractures. The analyses based on the force balance equation reveal that the plug flow region in Bingham radial flow is independent of the fracture radius, and is not a growth function adapted from the solution of one-dimensional (1D) slit flow according to ‘similarity’. Based on the shear stress distribution, we analytically proposed that a non-uniform plug flow region cannot exist. The Bingham fluid (grout) penetration and flowrate evolution as functions of grouting time are given using the correct expression for the plug flow region. The radius-independent plug flow region and the presented flowrate evolution equation are also verified numerically. For radial flow, the relative penetration length is equal to the relative width of plug flow region, which is the same as that for 1D channel flow. Discrepancies in analytical solutions for grout penetration and flowrate evolution were also illustrated. The clarification of the plug flow region and evaluation of discrepancies in analytical solutions presented in this work could simplify modeling and design of grouting in rock engineering applications.

The current paper is focused on analyzing the flow of viscoplastic fluid in a cavity that is driven by the two walls. The Lattice Boltzmann method (LBM) is used to solve the discrete Boltzmann equation. To represent the stress-strain rate relationship of viscoplastic fluids, the Bingham Papanastasiou constitutive model is considered. Cavity flow filled with Bingham fluids is considered for validating the present LBM code. After successful validation of the code, the analysis is extended for three dissimilar wall motions-simultaneous and opposed movement of the parallel facing walls, and the simultaneous motion of non-facing walls. The flow dynamics of Bingham fluid is influenced by Reynolds and Bingham numbers which can be studied using velocity and streamline plots. Subsequently, the yielded and un-yielded zones in a cavity have been effectively tracked using the limiting condition of yield stress. Further, the effect of wall motion on the variation of those zones inside a cavity has been studied. Finally, the drag coefficient for considered wall motions is presented.

A mathematical study is developed for the electro-osmotic flow of a non-Newtonian fluid in a wavy microchannel in which a Bingham viscoplastic fluid model is considered. For electric potential distributions, a Poisson-Boltzmann equation is employed in the presence of an electrical double layer (EDL). The analytical solutions of dimensionless boundary value problems are obtained with the Debye-Huckel theory, the lubrication theory, and the long wavelength approximations. The effects of the Debyelength parameter, the plug flow width, the Helmholtz-Smoluchowski velocity, and the Joule heating on the normalized temperature, the velocity, the pressure gradient, the volumetric flow rate, and the Nusselt number for heat transfer are evaluated in detail using graphs. The analysis provides important findings regarding heat transfer in electroosmotic flows through a wavy microchannel.

Ice-containing cemented paste backfill (ICPB) is a novel type of backfill material that has both the advantages of traditional CPB and the underground cooling function. The pipeline transport characteristics of ICPB need to be investigated to promote its application, which are influenced by the rheological properties. The objective of this study is to investigate the rheological properties of ICPB by means of experimental and computational fluid dynamics (CFD) methods. The rheological parameters were measured using a rheometer and the slump test was performed in a cylinder slump. A calculation model that considered the slump mold lifting and the ice particle phase change was established to simulate the cylindrical slump. The results indicate that the ICPB behaved as a Bingham fluid. The yield stress, plastic viscosity, and thixotropy increased with an increase in the ice/water ratio and concentration, while the cylindrical slump decreased under the same variation. The simulation results of the cylindrical slump were compatible with the experimental results. Furthermore, the numerical study indicates the slump increased when the test temperature was increased from 285 K to 305 K. The results found in this study can be used as a reference for the pipe design of ICPB.

In the present work, the sloshing phenomenon in a rectangular tank with non-Newtonian Bingham fluid is modeled. The numerical method used is the Incompressible Smoothed particle Hydrodynamics (ISPH...

In an effort to further extend bearing life, the authors have attempted to acquire greater knowledge regarding lubricating grease behavior in a bearing. While conducting experiments, some kinds of difficulties commonly arise when attempting to observe grease behavior directly from the bearing exterior without removing seals and shields. Making a breakthrough such a troubling aspect, X-ray computed tomography (CT), which is one of the non-destructive inspection techniques, was employed and resulted in visualizing remarkable details of grease distribution in a resin ball bearing. Hydrodynamic grease transition from churning to channeling state was well revealed by the mixture distribution of urea and barium-based greases which have different properties of X-ray absorption capability. Furthermore, the three dimensional unsteady liquid-gas multi-phase flows analysis was performed. Hydrodynamic feature of grease was regarded as a non-Newtonian fluid, which shows a highly non-linear flow curve, and the constitutive equation of modified Bingham plastic model proposed by Papanastasiou was applied to rheological property. Through these novel experimental and calculation approaches, several new insights about grease behavior inside a ball bearing were brought out.

Bio-mimicking the natural structure is one of the important aspects which have been developed with various inventions in the fields of engineering. There has been a subsequent effort to mimic the branched structure that can be seen in the leaves of trees and also in the branched structure of veins and arteries in the human body. This is done on the Lifting plate Hele-Shaw Cell (LHSC) setup. In LHSC, Bingham plastic fluid is squeezed between two parallel plates and then one of the plates is lifted up. This leads to some indefinite branched pattern commonly known as viscous fingers or micro fractals. Viscous fingering can be obtained using Newtonian fluid, Bingham plastic fluid. In this paper Bingham plastic fluid is used to study. Saffman-Taylor instability explains the formation of these viscous fingers. The main part of this paper is to control the formation of the viscous fingers and to get it in the desired shape. This can be done by providing some anisotropies on the surface of the plate in the form of pits or holes. The experiments are conducted under wide range of lifting velocities in Hele-Shaw cell. The geometrical dimensions of the viscous fingers are measured and their variation with the lifting velocity is plotted.

The linear stability of channel flows driven by pressure drops is carried out for fluids that exhibit a yield stress, the mechanical response of which is prescribed by the Bingham constitutive relation or two of its regularizations: the model due to Allouche and co-workers that is usually referred to as the “simple model”, and Papanastasiou model. Despite the fact that these two regularized models provide a good approximation of the steady Poiseuille flow of a Bingham fluid, they fail to predict the stability characteristics of the exact Bingham model. The critical thresholds for the onset of turbulence predicted by using the simple and Papanastasiou models are essentially the same, but they differ significantly from that of the exact Bingham model. This discrepancy is shown to be due to the absence of energy dissipation in the rigid core of a Bingham fluid.

A journal bearing is a machine element that is used to keep the shaft rotating about its axis. The increasing demand for journal bearing applications in high-speed machines that are efficient and economical has resulted in the need for improvements to the acoustic and tribological performance of journal bearings. In order to improve the tribological and acoustic performance, this study aims to propose a novel journal bearing design by introducing a roughness condition in a specific zone of the stationary bearing surface. In addition, the impact of the application of engineered roughness on the performance of Bingham-plastic-lubricated bearings is investigated in more detail. Considering the effect of cavitation, the analysis was conducted using a 3D computational fluid dynamics (CFD) model of a journal bearing. In comparison with the Reynolds equation—which is inertialess—for lubrication analysis, the use of a 3D CFD model based on Navier–Stokes equations reflects more detailed flow characteristics. Moreover, in this work, variations in the area of surface roughness were employed, resulting in various roughness patterns on the surface of the journal bearing, so that the acoustic and tribological performances of the journal bearing were anticipated to be enhanced. The findings of this study show that under non-Newtonian lubrication of the bearing, the engineered roughness has a strong effect in altering the tribological performance. Furthermore, the well-chosen roughened surface was proven to be more pronounced in enhancing the load support and reducing the friction force. The simulation results also show that using an engineered surface has little effect on the noise of the bearing.

When drilling an oil well, a viscous fluid is injected to aid drilling. This fluid is also responsible for removing the cuttings and maintaining structural stability of well. The rheology of this drilling fluid has a direct influence on the cleaning of the well, on the dynamics of the fluid in pipe and annular areas. Linear mathematical extrapolations for high pressure and high temperature environments can lead to rheology errors up to 75%. In this study, a finite volume model was developed to simulate the flow of a water-based mud in annular and jetting environments in the drilling environment. Annulars were made by steel pipes and permeable formations. The fluids evaluated were developed empirically with xanthan gum and bentonite clay. The numerical results are consistent with literature and represent characteristics of a Yield Power Law fluid and a Bingham plastic. A comparison was made with water, allowing a correlation between rheological effects and fluid dynamics in annular and high vorticity regions.

Modern concretes such as ultra-high performance concrete (UHPC) show excellent strength properties combined with favorable flow properties. However, the flow properties depend strongly on process parameters during production (temperature, humidity etc.), but also change sensitively even with slight variations in the mixture. In order to ensure desired processing of the fluidlike material and consistent process quality, the flow properties of the concrete must be evaluated quantitatively and objectively. The usual evaluation of measurements from concrete rheometers, for example of the ball probe system type, does not allow the direct determination of the objective material parameters yield stress and plastic viscosity of the sample. We developed a simulation-based method for the evaluation of rheometric measurements of fine grained high performance concretes like self-compacting concrete (SCC) and UHPC. The method is based on a dimensional analysis for ball measuring systems. Through numerical parameter studies we were able to describe the identified relationship between measuring quantities and material parameters quantitatively for two devices of this type. The evaluation method is based on the Bingham model. With this method it is possible to measure both the yield stress and the plastic viscosity of the fresh sample simultaneously. Device independence of the evaluation process is proven and an application to fiber-reinforced UHPC is presented.

On the basis of the laws of conservation and the principles of thermodynamics, a mathematical model of the flow of a two-phase granular fluid is proposed. One of the phases is the viscoplastic granular Bingham fluid; the other phase is a viscous Newtonian fluid. The equations for flows in the Hele–Shaw cell are analyzed asymptotically, i.e., when the flat-channel width is much less than its length. The correlations between the phase flow rates and the pressure gradient leading to equations of filtration for a two-phase granular viscoplastic fluid are constructed. The criterion is found for the initiation of motion of a granular phase in a porous medium. It is established that, depending on the shear-yield stress, such a phase does not flow if either the pressure gradient or the channel thickness is small. The phase flow rates are analyzed numerically at various input parameters such as the phase viscosities, phase resistivities, ultimate shear stress, etc. The factors slowing down the penetrating motion of the solid phase into the porous medium are revealed.

This study examined the magnetorheological (MR) behavior of magnetite-based MR fluids according to the magnetite concentration, magnetic field strength, and different viscosity of silicone oil as the carrier fluid. To this end, an experimental design was chosen to decrease the experimental error and provide a better rheological interpretation. The flow behavior parameters increased remarkably with increasing medium viscosity and nanoparticles (NPs) concentration, whereas the magnetic field strength had a milder effect, exhibiting a saturation value above which its effect became almost negligible. Moreover, both the Bingham plastic and power law models were found to fit their flow behaviors well. Almost all the critical rheological parameters of the aforementioned models experienced remarkable improvement after increasing the magnetic field strength, NPs concentration and carried fluid viscosity. In addition, mathematical correlations were derived to model each of the rheological parameters such as plastic viscosity and yield stress (for Bingham fluid model) and the consistency and power-law indices (for power-law model).

A thin fluid film formed by withdrawing a vertical plate is being studied theoretically. It is the case of a viscoplastic fluid described by the Bingham model. The boundary-value problem is formulated by the lubrication theory. The analysis presented herein required to solve numerically a system of ordinary differential equations at appropriate boundary conditions. For the viscoplastic fluid, we have predicted a film thickness over a wide range of a capillary number (10−4 to 10), with a particular emphasis on finding the free surface location at various conditions: (a) the unyielded region is much larger than the yielded one; (b) without the yielded region; and (c) the unyielded region is smaller but is compared with the yielded one. For the case of Bingham fluid and the limiting case of Newtonian fluid, the agreement was found between our theory and the experiment for the considered capillary number range.A thin fluid film formed by withdrawing a vertical plate is being studied theoretically. It is the case of a viscoplastic fluid described by the Bingham model. The boundary-value problem is formulated by the lubrication theory. The analysis presented herein required to solve numerically a system of ordinary differential equations at appropriate boundary conditions. For the viscoplastic fluid, we have predicted a film thickness over a wide range of a capillary number (10−4 to 10), with a particular emphasis on finding the free surface location at various conditions: (a) the unyielded region is much larger than the yielded one; (b) without the yielded region; and (c) the unyielded region is smaller but is compared with the yielded one. For the case of Bingham fluid and the limiting case of Newtonian fluid, the agreement was found between our theory and the experiment for the considered capillary number range.

The aiding-buoyancy mixed convection heat transfer in Bingham plastic fluids from an isothermal cylinder of elliptical and circular shape in a vertical adiabatic channel is numerically investigated. For a fixed shape of the elliptical cylinder E = 2 (ratio of major to minor axes), the effect of confinement is studied for three values of blockage ratio, B, defined as the ratio of the channel width to the circumference of the cylinder/π, as 6.5, 2.17 and 1.3. In order to delineate the role of cross-section of the cylinder, results are also presented here for a circular cylinder of the same heat transfer area as the elliptical cylinder. The results presented herein span the range of conditions as: Bingham number, 0 ≤ Bn ≤ 100, Reynolds number, 1 ≤ Re ≤ 40, and Prandtl number, 1 ≤ Pr ≤ 100 over the range of Richardson number Ri = 0 (pure forced convection) to Ri = 10. Extensive results on drag coefficient, local and surface averaged values of the Nusselt number and yield surfaces are presented herein to elucidate the combined effects of buoyancy, blockage ratio and fluid yield stress. The morphology of the yield surfaces shows that the unyielded plug regions formed upstream and downstream of the cylinder grow faster at low Reynolds numbers with the increasing yield stress effects under the weak buoyancy forces, i.e., small values of Grashof or Richardson number. The heat transfer enhancement is observed with the increasing channel-confinement due to the sharpening of the temperature gradients near the surface of the cylinder. The average Nusselt number shows a positive dependence on the Reynolds number, Prandtl number and Richardson number irrespective of the shape of the cylinder or the type of fluid. By employing the modified definitions of the dimensionless parameters (based on the two choices of the overall effective fluid velocity), predictive correlations have been established for estimating the value of the average Nusselt number in a new application.