Dimensionless Control Parameters Research Articles

Non-similar analysis of bio-convective micropolar nanofluid flow including gyrotactic microorganisms across a stretched geometry

This study explores the novel interaction of gyrotactic microorganisms within a two-dimensional bioconvective micropolar nanofluid flow along a stretched surface, employing the Tiwari-Das nanofluid model. Uniquely, blood is modeled as a micropolar fluid incorporating titanium dioxide ( Ti O 2 ) and silver ( Ag ) nanoparticles. The analysis incorporates influential factors such as viscous dissipation, magnetic field, and heat source effects. By converting the control system into dimensionless nonlinear coupled differential equations, a complete understanding of the system is attained through the application of the local non-similarity approach and bvp4c (MATLAB tool). Particularly, the study provides valuable insight into the influence of dimensionless control parameters on key profiles such as velocity, microrotation, temperature, concentration, and microbial profiles. The findings show that temperature, microrotation, and velocity profiles are all improved by raising the values of micropolar parameters. Furthermore, the microbe profile shows an inverse connection with the Lewis number and a positive association with the Peclet number, whereas the concentration profile declines with increasing Schmidt number values. Furthermore, the study provides a quantitative analysis through tables, revealing variations in the drag coefficient and Nusselt number. This innovative exploration sheds light on the unique behavior of micropolar bioconvective nanofluids, offering valuable contributions to the field of fluid dynamics and bioengineering.

Scaling laws for planetary sediment transport from DEM-RANS numerical simulations

We use an established discrete element method (DEM) Reynolds-averaged Navier–Stokes (RANS)-based numerical model to simulate non-suspended sediment transport across conditions encompassing almost seven orders of magnitude in the particle–fluid density ratio$s$, ranging from subaqueous transport ($s=2.65$) to aeolian transport in the highly rarefied atmosphere of Pluto ($s=10^7$), whereas previous DEM-based sediment transport studies did not exceed terrestrial aeolian conditions ($s\approx 2000$). Guided by these simulations and by experiments, we semi-empirically derive simple scaling laws for the cessation threshold and rate of equilibrium aeolian transport, both exhibiting a rather unusual$s^{1/3}$-dependence. They constitute a simple means to make predictions of aeolian processes across a large range of planetary conditions. The derivation consists of a first-principle-based proof of the statement that, under relatively mild assumptions, the cessation threshold physics is controlled by only one dimensionless control parameter, rather than two expected from dimensional analysis. Crucially, unlike existing models, this proof does not resort to coarse-graining the particle phase of the aeolian transport layer above the bed surface. From the pool of existing models, only that by Pähtzet al.(J. Geophys. Res.: Earth, vol. 126, 2021, e2020JF005859) is somewhat consistent with the combined numerical and experimental data. It captures the scaling of the cessation threshold and the$s^{1/3}$-dependence of the transport rate, but fails to capture the latter's superimposed grain size dependence. This hints at a lack of understanding of the transport rate physics and calls for future studies on this issue.

General relationships of dimensionless parameters for thermocapillary-buoyancy flow in side-heated annular pools

It is well-known that thermocapillary-buoyancy flow instability can be affected by several dimensionless control parameters. The general relationships among them are helpful for actively controlling flow instability but were not clear. This work reported some experiments in the annular pool with a radius ratio of 0.25, and proposed the general relationships among these dimensionless control parameters, especially for the critical Marangoni number of the flow destabilization and the dimensionless oscillatory frequency of pulsating pattern. Combined with previous studies, these general relationships are achieved within no more than 15% relative deviations, whose applicability is for the thermocapillary-buoyancy flow of medium Prandtl number fluids in side-heated annular pools. Based on the experimental findings, we can confirm that the flow destabilization can be suppressed by buoyancy. The steady multicellular flow appears before flow destabilization. When the flow is destabilized, the pulsating pattern first appears, and then a superimposition of the pulsating pattern and the flower-like wave.

Stability of interlinked neutron vortex and proton flux-tube arrays in a neutron star – III. Proton feedback

ABSTRACT The coupled, time-dependent Gross–Pitaevskii and Ginzburg–Landau equations are solved simultaneously in three dimensions to investigate the equilibrium state and far-from-equilibrium, spin-down dynamics of an interpenetrating neutron superfluid and proton type-II superconductor, as an idealized description of the outer core of a neutron star. The simulations generalize previous calculations without the time-dependent Ginzburg–Landau equation, where proton feedback is absent. If the angle θ between the rotation and magnetic axes does not equal zero, the equilibrium state consists of geometrically complicated neutron vortex and proton flux-tube tangles, as the topological defects pin to one another locally but align with different axes globally. During spin-down, new types of motion are observed. For θ = 0, entire vortices pair rectilinearly with flux tubes and move together while pinned. For θ ≠ 0, vortex segments pair with segments from one or more flux tubes, and the paired segments move together while pinned. The degree to which proton feedback impedes the deceleration of the crust is evaluated as a function of θ and the pinning strength, η. Key geometric properties of vortex-flux-tube tangles, such as filament length, mean curvature, and polarity are analysed. It is found that proton feedback smooths the deceleration of the crust, reduces the rotational glitch sizes, and stabilizes the vortex tangle dynamics. The dimensionless control parameters in the simulations are mutually ordered to match what is expected in a real neutron star, but their central values and dynamic ranges differ from reality by many orders of magnitude due to computational limitations.

Quasi-steady-state modelling and characterization of diffusion-controlled dissolution from monodisperse prolate and oblate spheroidal particles

A quasi-steady-state model of the dissolution of a single prolate or oblate spheroidal particle has been developed based on the exact solution of the steady-state diffusion equation for mass transfer in an unconfined media. With appropriate treatment of bulk concentration, the model can predict the detailed dissolution process of a single particle in a container of finite size. The dimensionless governing equations suggest that the dissolution process is determined by three dimensionless control parameters, initial solid particle concentration, particle aspect ratio and the product of specific volume of solid particles and saturation concentration of the dissolved substance. Using this model, the dissolution processes of felodipine particles are analysed in a broad range of space of the three control parameters and some characteristics are identified. The effects of material properties indicated by the product of specific volume and saturation concentration are also analysed. The model and the analysis are applicable to the system of monodisperse spheroidal particles of the same shape.

Analytical and computational study of cascade reaction processes in catalytic fibrous membranes

Multistep catalytic reactions use two different catalysts for the A→B and the subsequent B→C reaction, respectively. Often the employed catalysts are chemically incompatible, such as acid–base systems, which prohibits simple mixing in one solution. In this work, we study the efficiency of reactors where the incompatible catalytic sites are immobilized on fibrous membranes. We compare a lattice Boltzmann based solver for the advection–diffusion–reaction equation, a random walk particle tracking method and a simple theoretical model to investigate the reaction efficiency as a function of two dimensionless control parameters: the Péclet and the Damköhler number. We find that, while the efficiency decreases with higher flow speed (due to the reduced reaction time), the total production nevertheless increases due to the higher mass flux in most cases. Our results further show that, even at high flow speeds, spatial proximity of the two catalysts increases reaction efficiency, which supports recent experimental efforts to locate both catalysts on a single fiber in a side-by-side geometry.

Jammed disks of two sizes and weights in a channel: Alternating sequences

Disks of two sizes and weights in alternating sequence are confined to a long and narrow channel. The axis of the channel is horizontal and its plane vertical. The channel is closed off by pistons that freeze jammed microstates out of loose disk configurations subject to moderate pressure, gravity, and random agitations. Disk sizes and channel width are such that under jamming no disk remains loose and all disks touch one wall. We present exact results for the characterization of jammed macrostates including volume and entropy. The rigorous analysis divides the disk sequences of jammed microstates into overlapping tiles from which we construct a small number of species of statistically interacting particles. Jammed macrostates depend on dimensionless control parameters inferred from ratios between measures of expansion work against the pistons, gravitational potential energy, and intensity of random agitations. These control parameters enter the configurational statistics via the activation energies prior to jamming of the particles. The range of disk weights naturally divides into regimes where qualitatively different features come into play. We sketch a path toward generalizations that include random sequences under a modified jamming protocol.

Trapping patterns during capillary displacements in disordered media

We investigate the impact of wettability distribution, pore size distribution and pore geometry on the statistical behaviour of trapping in pore-throat networks during capillary displacement. Through theoretical analyses and numerical simulations, we propose and prove that the trapping patterns, defined as the percentage and distribution of trapped elements, are determined by four dimensionless control parameters. The range of all possible trapping patterns and how the patterns are dependent on the four parameters are obtained. The results help us to understand the impact of wettability and structure on trapping behaviour in disordered media.

Unsteady forcing of turbulence by a randomly actuated impeller array

We investigate the unsteady forcing of turbulent flow in a well-stirred reactor using opposing arrays of pitched-blade impellers which randomly and independently reverse rotation. We systematically explore the dependence of the large-scale motions and the homogeneity and isotropy of the turbulence upon the forcing. We identify three dimensionless control parameters: the source fraction (the fraction of time spent in clockwise motion), the dimensionless forcing period and an impeller Reynolds number. We find the timescale of unsteady motion corresponds to the forcing period T, the average period of impeller reversal, independently of the impeller angular speed varOmega and source fraction. As in jet-stirred tanks, unsteady forcing substantially increases the unsteady kinetic energy, energy dissipation, integral length scale and Taylor microscale Reynolds number (R_lambda) and improves the homogeneity and isotropy of the flow, provided the source fraction is chosen optimally and the forcing period is sufficiently large (varOmega T > 10^3); impeller Reynolds number has a relatively small influence. The forcing period must be matched to angular speed: decreasing the forcing period below this threshold results in a less intense, more inhomogeneous turbulent flow. Spectra of two-point velocity increments demonstrate that unsteady energy injection is dominated by axial shear generated across impellers and becomes less prominent at smaller scales. However, even at R_lambda approx 354, the signature of this unsteady forcing can still be detected in near-dissipation-range statistics. These observations provide insight into optimisation of forcing and the mechanism of energy transfer when using unsteady forcing to generate turbulence in confined vessels.Graphical abstract

Time scales and rheology of visco-cohesive granular flows

In the presence of viscous and cohesive interactions between particles, a granular flow is governed by several characteristic time and stress scales that determine its rheological properties (shear stress, packing fraction, effective viscosities). In this paper, we revisit and extend the scaling arguments used previously for dry cohesionless granular flows and suspensions. We show that the rheology can be in principle described by a single dimensionless control parameter that includes all characteristic times. We also briefly present simulation results for 2D sheared suspensions and 3D wet granular flows where the effective friction coefficient and packing fraction are consistently described as functions of this unique control parameter.

Tidal locking and the gravitational fold catastrophe

The purpose of this work is to study the phenomenon of tidal locking in a pedagogical framework by analyzing the effective gravitational potential of a two-body system with two spinning objects. It is shown that the effective potential of such a system is an example of a fold catastrophe. In fact, the existence of a local minimum and saddle point, corresponding to tidally locked circular orbits, is regulated by a single dimensionless control parameter that depends on the properties of the two bodies and on the total angular momentum of the system. The method described in this work results in compact expressions for the radius of the circular orbit and the tidally locked spin/orbital frequency. The limiting case in which one of the two orbiting objects is point-like is studied in detail. An analysis of the effective potential, which in this limit depends on only two parameters, allows one to clearly visualize the properties of the system. The notorious case of the Mars' moon Phobos is presented as an example of a satellite that is past the no return point and, therefore, will not reach a stable or unstable tidally locked orbit.

Vortices of electro-osmotic flow in heterogeneous porous media

Traditional models of electrokinetic transport in porous media are based on homogenized material properties, which neglect any macroscopic effects of microscopic fluctuations. This perspective is taken not only for convenience, but also motivated by the expectation of irrotational electro-osmotic flow, proportional to the electric field, for uniformly charged surfaces (or constant zeta potential) in the limit of thin double layers. Here, we show that the inherent heterogeneity of porous media generally leads to macroscopic vortex patterns, which have important implications for convective transport and mixing. These vortical flows originate due to competition between pressure-driven and electro-osmotic flows, and their size are characterized by the correlation length of heterogeneity in permeability or surface charge. The appearance of vortices is controlled by a single dimensionless control parameter, defined as the ratio of a typical electro-osmotic velocity to the total mean velocity.

Biomathematical model for gyrotactic free-forced bioconvection with oxygen diffusion in near-wall transport within a porous medium fuel cell

Bioconvection has shown significant promise for environmentally friendly, sustainable “green” fuel cell technologies. The improved design of such systems requires continuous refinements in biomathematical modeling in conjunction with laboratory and field testing. Motivated by exploring deeper the near-wall transport phenomena involved in bio-inspired fuel cells, in the present paper, we examine analytically and numerically the combined free-forced convective steady boundary layer flow from a solid vertical flat plate embedded in a Darcian porous medium containing gyrotactic microorganisms. Gyrotaxis is one of the many taxes exhibited in biological microscale transport, and other examples include magneto-taxis, photo-taxis, chemotaxis and geo-taxis (reflecting the response of microorganisms to magnetic field, light, chemical concentration or gravity, respectively). The bioconvection fuel cell also contains diffusing oxygen species which mimics the cathodic behavior in a proton exchange membrane (PEM) system. The vertical wall is maintained at iso-solutal (constant oxygen volume fraction and motile microorganism density) and iso-thermal conditions. Wall values of these quantities are sustained at higher values than the ambient temperature and concentration of oxygen and biological microorganism species. Similarity transformations are applied to render the governing partial differential equations for mass, momentum, energy, oxygen species and microorganism species density into a system of ordinary differential equations. The emerging eight order nonlinear coupled, ordinary differential boundary value problem features several important dimensionless control parameters, namely Lewis number (Le), buoyancy ratio parameter i.e. ratio of oxygen species buoyancy force to thermal buoyancy force (Nr), bioconvection Rayleigh number (Rb), bioconvection Lewis number (Lb), bioconvection Péclet number (Pe) and the mixed convection parameter ([Formula: see text] spanning the entire range of free and forced convection. The transformed nonlinear system of equations with boundary conditions is solved numerically by a finite difference method with central differencing, tridiagonal matrix manipulation and an iterative procedure. Computations are validated with the symbolic Maple 14.0 software. The influence of buoyancy and bioconvection parameters on the dimensionless temperature, velocity, oxygen concentration and motile microorganism density distribution, Nusselt, Sherwood and gradient of motile microorganism density are studied. The work clearly shows the benefit of utilizing biological organisms in fuel cell design and presents a logical biomathematical modeling framework for simulating such systems. In particular, the deployment of gyrotactic microorganisms is shown to stimulate improved transport characteristics in heat and momentum at the fuel cell wall.

From droplets to waves: periodic instability patterns in highly viscous microfluidic flows

We experimentally study the transition from droplet to wave regimes in microfluidic liquid–liquid multiphase flows having large differences in viscosity. A unified approach based on periodic pattern analysis is employed to study relationships between dispersed and separated flow regimes, including dripping, jetting, capillary waves, inertial waves and core–annular flows over a wide range of flow rates and viscosity contrasts. We examine the morphology and dynamics of each flow regime based on wavelength, frequency and velocity of repeating unit cells to elucidate their connections and to develop predictive capabilities based on dimensionless control parameters. We demonstrate in particular that pattern selection is contingent upon the propagation velocity of droplets and waves at the transition. We also investigate microfluidic wave breaking phenomena with the formation of ligaments and droplets from wave crests in both capillary and inertial wave regimes. This work expands conventional multiphase flow regimes observed in microchannels and shows new routes to disperse highly viscous materials using interfacial waves dynamics in confined microsystems.

Numerical Treatment of MHD Flow of Casson Nanofluid via Convectively Heated Non-Linear Extending Surface with Viscous Dissipation and Suction/Injection Effects

This paper introduces the effect of heat absorption (generation) and suction (injection) on magnetohydrodynamic (MHD) boundary-layer flow of Casson nanofluid (CNF) via a non-linear stretching surface with the viscous dissipation in two dimensions. By utilizing the similarity transformations, the leading PDEs are transformed into a set of ODEs with adequate boundary conditions and then resolved numerically by (4–5)th-order Runge-Kutta Fehlberg procedure based on the shooting technique. Numerical computations are carried out by Maple 15 software. With the support of graphs, the impact of dimensionless control parameters on the nanoparticle concentration profiles, the temperature, and the flow velocity are studied. Other parameters of interest, such as the skin friction coeffi- cient, heat, and mass transport at the diverse situation and dependency of various parameters are inspected through tables and graphs. Additionally, it is verified that the numerical computations with the reported earlier studies are in an excellent approval. It is found that the heat and mass transmit rates are enhanced with the increasing values of the power-index and the suction (blowing) parameter, whilst are reduced with the boosting Casson and the heat absorption (generation) parameters. Also, the drag force coefficient is an increasing function of the powerindex and a reduction function of Casson parameter.

Boundary condition control of fluvial obstacle mark formation – framework from a geoscientific perspective

ABSTRACTObstacle marks are sedimentary bedforms, typically composed of an upstream local scour hole and a downstream sediment accumulation in the vicinity of an obstruction that is exposed to a current. However, specific morphologies are variable in fluvial, coastal and submarine environments. Although obstacle marks and the phenomenon of local scouring are subject to different scientific disciplines, the objectives of investigations are rather incoherent and no systematic framework for analysing and evaluating boundary condition control exists yet, especially concerning limited knowledge of the cause and effect relationship of obstacle mark formation at instream boulders or vegetation elements in variable environmental conditions. Thus, a parameter framework is developed which identifies a spectrum of extrinsic and intrinsic boundary conditions that control the major process dynamics of obstacle mark formation. The framework is composed of dimensionless control parameters that are separated by a hierarchical order regarding their significance for obstacle mark formation. Primary control parameters determine the geometrical scale of flow field at the obstacle, and therefore control the potential maximum size of the obstacle. Secondary control parameters affect the dynamics of the flow field in geometrical scale and limit the potential maximum size of the emerging sedimentary structure if thresholds are crossed. The framework is supposed to be a foundation for subsequent quantification and determination of thresholds by systematic laboratory studies. To elucidate this, flume‐based research is presented, evaluating the influence of different flow levels at boulder‐like obstacles of different shapes. The results show that obstacle mark dimensions were maximized at shallow flow depths compared to obstacle dimensions, while deep flows at submerged boulder‐like obstructions caused considerably smaller obstacle marks. In interdependency with a rounded and more streamlined obstacle shape, deep flows even cause a deviation of morphology if the flow depth above an obstacle exceeds 1.6 times the obstacle's dimensions. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.

Particle-based modeling of the mechanical behavior of porous fluid-saturated viscoelastic solids

In the paper, we developed a macroscopic discrete element model of permeable fluid-saturated materials with solid skeleton characterized by viscoelastic rheological properties. The Biot’s linear model of poroelasticity was used as a mathematical basis for describing the mechanical interrelation between the solid skeleton and interstitial fluid. Using this model, we numerically studied the dependences of the effective Young’s modulus and strength of fluid-saturated viscoelastic materials on the loading rate, sample size and the mechanical parameters, which determine the relaxation time of the solid-phase skeleton and the time scale of redistribution of fluid in the pore space. We revealed two dimensionless control parameters that determine the dynamic values of the effective mechanical characteristics of the samples under compression loading. We obtained the general relations that describe the above-mentioned dependences in terms of the two proposed control parameters. These relations have a logistic nature and are described by sigmoid functions. The importance of the proposed empirical expressions is determined by the possibility of their application for predicting the mechanical response of fluid-saturated materials of different nature (bone tissue, rocks, porous materials with polymeric skeleton, including elastomers, etc.) under dynamic loading.

Phase transition in time-reversible Navier-Stokes equations.

We present a comprehensive study of the statistical features of a three-dimensional (3D) time-reversible truncated Navier-Stokes (RNS) system, wherein the standard viscosity ν is replaced by a fluctuating thermostat that dynamically compensates for fluctuations in the total energy. We analyze the statistical features of the RNS steady states in terms of a non-negative dimensionless control parameter R_{r}, which quantifies the balance between the fluctuations of kinetic energy at the forcing length scale ℓ_{f} and the total energy E_{0}. For small R_{r}, the RNS equations are found to produce "warm" stationary statistics, e.g., characterized by the partial thermalization of the small scales. For large R_{r}, the stationary solutions have features akin to standard hydrodynamic ones: they have compact energy support in k space and are essentially insensitive to the truncation scale k_{max}. The transition between the two statistical regimes is observed to be smooth but rather sharp. Using insights from a diffusion model of turbulence (Leith model), we argue that the transition is in fact akin to a continuous second-order phase transition, where R_{r} indeed behaves as a thermodynamic control parameter, e.g., a temperature. A relevant order parameter can be suitably defined in terms of a (normalized) enstrophy, while the symmetry-breaking parameter h is identified as (one over) the truncation scale k_{max}. We find that the signatures of the phase transition close to the critical point R_{r}^{★} can essentially be deduced from a heuristic mean-field Landau free energy. This point of view allows us to reinterpret the relevant asymptotics in which the dynamical ensemble equivalence conjectured by Gallavotti [Phys. Lett. A 223, 91 (1996)PYLAAG0375-960110.1016/S0375-9601(96)00729-3] could hold true. We argue that Gallavotti's limit is precisely the joint limit R_{r}→[over >]R_{r}^{★} and h→[over >]0, with the overset symbol ">" indicating that those limits are approached from above. The limit therefore relates to the statistical features at the critical point. In this regime, our numerics indicate that the low-order statistics of the 3D RNS are indeed qualitatively similar to those observed in direct numerical simulations of the standard Navier-Stokes equations with viscosity chosen so as to match the average value of the reversible thermostat. This result suggests that Gallavotti's equivalence conjecture could indeed be of relevance to model 3D turbulent statistics, and provides a clear guideline for further numerical investigations involving higher resolutions.

Self-similar dynamics of two-phase flows injected into a confined porous layer

We study the dynamics of two-phase flows injected into a confined porous layer. A model is derived to describe the evolution of the fluid–fluid interface, where the effective saturation of the injected fluid is zero. The flow is driven by pressure gradients due to injection, the buoyancy due to density contrasts and the interfacial tension between the injected and ambient fluids. The saturation field is then computed after the interface evolution is obtained. The results demonstrate that the flow behaviour evolves from early-time unconfined to late-time confined behaviours. In particular, at early times, the influence of capillary forces drives fluid flow and produces a new self-similar spreading behaviour in the unconfined limit that is distinct from the gravity current solution. At late times, we obtain two new similarity solutions, a modified shock solution and a compound wave solution, in addition to the rarefaction and shock solutions in the sharp-interface limit. A schematic regime diagram is also provided, which summarises all possible similarity solutions and the time transitions between them for the partially saturating flows resulting from fluid injection into a confined porous layer. Three dimensionless control parameters are identified and their influence on the fluid flow is also discussed, including the viscosity ratio, the pore-size distribution and the relative contributions of capillary and buoyancy forces. To underline the relevance of our results, we also briefly describe the implications of the two-phase flow model to the geological storage of $\text{CO}_{2}$, using representative geological parameters from the Sleipner and In Salah sites.

Buckling of a compliant hollow cylinder attached to a rotating rigid shaft

Bifurcations in the equilibrium shape of a thick elastic layer attached to a circular rigid cylinder rotating about its axis are investigated analytically and numerically. The centrifugal force breaks the symmetry of the system leading to deformations invariant along the axis as the result of an instability. The instability threshold depends at linear order on the relative thickness of the compliant layer, and on a dimensionless control parameter based on the elastic modulus, the angular velocity and the outer radius. A weakly non linear analysis, carried out for layers following the Mooney–Rivlin constitutive law, points out the discontinuous (sub-critical) features of the bifurcation, except for relative thickness laying in a very narrow range in which the bifurcation is super-critical. Numerical simulations in the fully post-buckled regime yield the absolute instability threshold, and the order in the rotational symmetry of the developed equilibrium shape.