MHD Maxwell dusty fluid in thermally stratified radiative flow with temperature-dependent thermal conductivity and Cattaneo-Christov model

It might be very important for the polymer processing industries to comprehend how Maxwell fluids behave on a stretched cylinder. Optimizing the extrusion and drawing processes can ensure the desired product qualities while avoiding faults. The objective of this study is heat transfer analysis on a Maxwell dusty fluid flow cylindrical surface with the Cattaneo-Christov concept. We immerse the cylinder in porous media, with a two-dimensional fluid regulating the flow. Our mathematical model further considers the effects of variable thermal conductivity, radiation, viscous and joule heating, magnetic field, thermal stratification, and slip velocity. Based on the presumptions, partial differential equations (PDE's) have been used to evolve the mathematical model. Using similarity transformations, the PDE's for heat and momentum for both phases are transformed into highly nonlinear ODE's.The numerical results have been obtained on these ordinary differential equations by using the RKF-45 method. This issue's main characteristic is that it examines the scenario's liquid and dust phases throughout. Results are given both visually and tabularly for the major parameters over a velocity, temperature, skin friction coefficient, and Nusselt number. When we compared our method to a previously published paper, we discovered a decent match. The findings, which were obtained for our system, show that the velocity and thermal gradient of both the phases of fluid and dust behave in an opposite trend in favor of rising Maxwell parameter, where the curvature parameter makes the rise in the same manner. Furthermore, the thermal transport profiles for both phases decline for the rising thermal time relaxation parameter.

A homogenized two-phase computational framework for meso- and macroscale blood flow simulations

BACKGROUND AND OBJECTIVE

Owing to the complexity of physics linked with blood flow and its associated phenomena, appropriate modeling of the multi-constituent rheology of blood is of primary importance. To this effect, various kinds of computational fluid dynamic models have been developed, each with merits and limitations. However, when additional physics like thrombosis and embolization is included within the framework of these models, computationally efficient scalable translation becomes very difficult. Therefore, this paper presents a homogenized two-phase blood flow framework with similar characteristics to a single fluid model but retains the flow resolution of a classical two-fluid model. The presented framework is validated against four different sets of experiments.

METHODS

The two-phase model of blood presented here is based on the classical diffusion-flux framework. Diffusion flux models are known to be less computationally expensive than two-fluid multiphase models since the numerical implementation resembles single-phase flow models. Diffusion flux models typically use empirical slip velocity correlations to resolve the motion between phases. However, such correlations do not exist for blood. Therefore, a modified slip velocity equation is proposed, derived rigorously from the two-fluid governing equations. An additional drag law for red blood cells (RBCs) as a function of volume fraction is evaluated using a previously published cell-resolved solver. A new hematocrit-dependent expression for lift force on RBCs is proposed. The final governing equations are discretized and solved using the open-source software OpenFOAM.

RESULTS

The framework is validated against four sets of experiments: (i) flow through a rectangular microchannel to validate RBC velocity profiles against experimental measurements and compare computed hematocrit distributions against previously reported simulation results (ii) flow through a sudden expansion microchannel for comparing experimentally obtained contours of hematocrit distributions and normalized cell-free region length obtained at different flowrates and inlet hematocrits, (iii) flow through two hyperbolic channels to evaluate model predictions of cell-free layer thickness, and (iv) flow through a microchannel that mimics crevices of a left ventricular assist device to predict hematocrit distributions observed experimentally. The simulation results exhibit good agreement with the results of all four experiments.

CONCLUSION

The computational framework presented in this paper has the advantage of resolving the multiscale physics of blood flow while still leveraging numerical techniques used for solving single-phase flows. Therefore, it becomes an excellent candidate for addressing more complicated problems related to blood flow, such as modeling mechanical entrapment of RBCs within blood clots, predicting thrombus composition, and visualizing clot embolization.

3D pore-scale characterization of colloid aggregation and retention by confocal microscopy: Effects of fluid structure and ionic strength

Understanding the mechanisms of colloid transport and retention as well as the spatial distribution of colloids in porous media is an important topic for contamination transport and remediation in subsurface environments. Utilizing advanced three-dimensional visualization experiments, we effectively capture the intricate distribution characteristics of colloids in the 3D pore space and quantify the size of colloid clusters that aggregate at fluid-fluid interfaces and solid surfaces during two-phase flow. Our experimental results reveal the influence of pore-scale events, such as Haines jumps and pinch-off, on colloid retention. Our results also indicate that large drainage rates can facilitate colloid retention on solid surfaces, especially under the condition of high ionic strength. This can be attributed to the migration of colloids from the fluid-fluid interface to the solid surface, propelled by transients in the local fluid structure. The findings reveal a synergistic effect of the ionic strength and hydrodynamic conditions on colloid transport and retention during two-phase flow and provide important insights for predicting the fate and transport of contaminants in soil and groundwater environments involving multiple fluid phases.

Numerical and experimental investigation on the performance of rapid ultrasonic-assisted nucleic acid extraction based on dispersive two-phase flow

BACKGROUND

Nucleic acid extraction (NAE) is an essential step in the whole process of nucleic acid detection (NAT). Traditional manual extraction methods are time-consuming and laborious, unfavorable to the point-of-care testing of nucleic acids. Ultrasound has been emphasized due to its noncontact and easy-to-manipulate characteristics, and integration with microfluidic chip can realize rapid NAE through acoustic streaming effect. The uniformity of magnetic bead mixing in this process is a critical factor affecting the extraction effect. In this study, we developed an ultrasound-assisted NAE technique based on the magnetic bead method and optimized the chip structure to achieve rapid NAE.

RESULT

We use ultrasonic-assisted coupled with magnetic bead method for ultra-fast NAE. The mixing process of magnetic beads driven by acoustic streaming is simulated by a dispersive two-phase flow model, and the ultrasonic incidence angle (θin), cone structure aspect ratio (Dc/Hc) and sheet structure thickness (Hp) are optimized to enhance the mixing performance. Furthermore, the effectiveness of NAE is validated by utilizing quantitative real-time PCR (qPCR) detection. The findings reveal that a θin value of 10° yields superior mixing performance compared to other incidence angles, resulting in a maximum increase of 84 % in mixing intensity. When Dc/Hc = 0.5 and Hp = 0.5 mm, the maximum mixing index in the localized region of the chamber after 1 s of ultrasound action can reach 83.6 % and 92.5 %, respectively. Compared to the original chamber, the CT values extracted after 5 s of ultrasound action shifted forward by up to 1.9 ct and 4.1 ct, respectively.

SIGNIFICANCE

The dispersed two-phase flow model can effectively simulate the mixing process of magnetic beads, which plays an important role in assisting the structural design of chip extraction chambers. The single-step mixing of ultrasound-assisted NAE takes only 15s to achieve an extraction performance comparable to manual extraction. The extraction process can be completed within 7 min after integrating this technology with microfluidic chips and automated equipment, providing a solution for automated and efficient NAE.

Comprehensive experimental and theoretical investigations on the effect of microbubble two-phase flow on the performance of direct-contact membrane distillation

This study provides a comprehensive and systematic overview of the application of gas-liquid two-phase flow with microbubbles in the feed stream to improve heat and mass transfer in direct-contact membrane distillation (DCMD) processes for seawater desalination. A swirl-flow-type microbubble generator (MBG) was installed at the feed-side inlet of the DCMD module to investigate its effect on transmembrane flux. The maximum improvement in the MBG-assisted DCMD permeation flux was found to be approximately 18% at a lower feed temperature (40 °C) and optimal air flow rate (50 cc/min), and an optimal MBG geometry comprising a swirler, a nozzle tip of diameter 2 mm, and a diffuser at an angle of 30°. The results were observed to be related to the number density of microbubbles less than 100 µm in size, which plays an important role in improving heat and mass transfer in two-phase flow. In addition, the simulation results based on conventional heat transfer correlations of bubbly flow underestimated the experimental results. Therefore, this study also aims to propose and verify a new two-phase flow heat transfer correlation. The proposed correlation considers the effects of bubble size distribution to accurately predict the performance of MBG-assisted DCMD processes.

Numerical simulations of high viscosity DNAPL recovery in highly permeable porous media under isothermal and non-isothermal conditions

We developed a decimetric size model based on coupling generalized Darcy's law and heat-transfer equations to model viscous dense non-aqueous phase liquid (DNAPL) pumping through highly permeable porous media under non-isothermal conditions. The presence of fingering and non-wetting phase ganglia was modeled through an unsteady capillary diffusion coefficient and an arbitrary heterogeneous permeability field. The model was validated using existing experimental data of a simple case, an oil injection in a 2D tank packed with glass beads. Next, we compared the results of this model against a DNAPL extracting situation in the 2D tank to better understand the two-phase flow behavior in highly permeable porous media. We found that natural convection during heating plays an essential role in heat transfer, especially in the wetting phase zone. By adding the dynamic effect (unsteady conditions) we were better able to describe the presence of the ganglia in porous media. We observed good agreement between modeled and experimental oil saturation curves until the breakthrough point, with a mean relative error of about 10% for low and high flow rates, and 8% and 16% after breakthrough for low and high flow rates, respectively. Extracting viscous oil at low flow rates and high temperature generates less fingering and is well described by the generalized Darcy's law. The remobilization of residual non-wetting ganglia after the breakthrough point at the outlet is, however, difficult to simulate using the generalized Darcy's law. In the end, we treated this issue by using a perturbed permeability field to simulate the observed fingering in the 2D tank.

Influence of the injection of densified polymer suspension on the efficiency of DNAPL displacement in contaminated saturated soils

Nowadays the remediation of DNAPL contaminated zones near groundwater has gained great prominence in environmental fields due to the high importance of water resources. In this work, we suggest injecting a densified polymer suspension by adding barite particles to displace DNAPL. To evaluate the efficiency of the densification of polymer suspensions on the displacement of DNAPL, various densities of barite-polymer suspension; lower, equal, and higher than the density of DNAPL were prepared and their rheological behavior was analyzed. Then flow experiments were performed using a decimetric-scale 2D tank. The displacement procedure was monitored with an imaging technique and the production and injection process data were recorded by mass balance interpretation. It was shown that the densification of the polymer suspension could improve the displacement efficiency of DNAPL up to four times. The clogging behavior of barite-polymer suspension was assessed in a 1D column. Generalized Darcy's law and the continuity equation were used to numerically simulate the experimental two-phase flow. To take into account the clogging behavior of the suspension, the transport equation of diluted species was implemented into the model. The simulation results show that the model can properly predicts the experimental consequences.

Experimental study of thermally enhanced recovery of high-viscosity DNAPL in saturated porous media under non-isothermal conditions

Thermal enhancement is known to be an efficient way to decrease the residual saturation of some common dense non-aqueous phase liquids (DNAPLs) after pumping. However, the effect of transient heat transfer during the recovery of a high-viscosity contaminant, such as coal tar, in highly permeable porous media is still unknown. A 2D tank experimental setup allowing monitoring of temperature and saturation fields during DNAPL pumping has been developed. Experiments were run under isothermal and non-isothermal conditions, at low and high flow rates. We investigated the presence of viscous fingering and how that influences the shape of the cone of depression, as well as the residual saturation. The saturation fields show that less viscous fingering occurs in pre-heated cases and that heating increases the recovery efficiency. Increasing the temperature increases the critical velocity and the viscosity ratio and helps to stabilize the interface between the non-wetting and wetting phase. Observations were first made on an oil and ethanol fluid pair because its properties were known, before extending the experiments to a coal tar and water fluid pair. Residual oil saturation after pumping was decreased by 6-16% in all pre-heated conditions. Pumping at low flow rate in these conditions leaves the smallest oil residual saturation (20%) after pumping. A low flow rate increases the recovery efficiency by reducing viscous fingering and by spreading the generated heat to a larger part of the tank. Finally, results on coal tar pumping show that the high thermal conductivity of water helps in keeping the temperature high during pumping. The residual coal tar saturation was reduced from 40% at 20 °C to 28% when pre-heating the tank. Operating at a low flow rate and with a uniform temperature is the key to recovering the highest amount of a viscous DNAPL such as coal tar from the soil and satisfying cleanup goals when using thermally enhanced pumping.

Aerosol deposition and airflow dynamics in healthy and asthmatic human airways during inhalation

Inhalation of aerosols such as pharmaceutical aerosols or virus aerosol uptake is of great concern to the human population. To elucidate the underlying aerosol dynamics, the deposition fractions (DFs) of aerosols in healthy and asthmatic human airways of generations 13-15 are predicted. The Navier-stokes equations governing the gaseous phase and the discrete phase model for particles' motion are solved using numerical methods. The main forces responsible for deposition are inertial impaction forces and complex secondary flow velocities. The curvatures and sinusoidal folds in the asthmatic geometry lead to the formation of complex secondary flows and hence higher DFs. The intensities of complex secondary flows are strongest at the generations affected by asthma. The DF in the healthy airways is 0%, and it ranges from 1.69% to 52.93% in the asthmatic ones. From this study, the effects of the pharmaceutical aerosol particle diameters in the treatment of asthma patients can be established, which is conducive to inhibiting the inflammation of asthma airways. Furthermore, with the recent development of COVID-19 which causes pneumonia, the predicted physics and effective simulation methods of bioaerosols delivery to asthma patients are vital to prevent the exacerbation of the chronic ailment and the epidemic.

Local dynamics of pharmaceutical powder fluidization using high speed long distance microscopy and particle image velocimetry

The local dynamics of fluidized pharmaceutical carrier powders in a turbulent channel flow was studied using particle image velocimetry (PIV) and High-speed, long-distance microscopy (HS-LDM). Four different lactose powders which have been used as a drug carrier in dry powder inhalers were used in this study. These powders have median powder particle diameters ranging between 61 and 121 μm. Air flow velocities ranging between 13.3 m/s and 66.7 m/s were examined. In addition, the effect of grid blockage ratio (ranging from ~25% to ~40% of the area of channel cross-section) was also investigated. Results show that the high-speed, long-distance microscopy (HS-LDM) technique was able to capture the mean velocity of the particles, and the results corresponded well with the PIV measurements. Results from the high-speed, long-distance microscopy (HS-LDM) method also demonstrate that the span of particle velocity closely follows that of the particle size distribution both for cohesive and non-cohesive powders. This study contributes towards an improved understanding of pharmaceutical carrier dynamics in turbulent channel flows and demonstrates how advanced image processing can be used to capture local particle dynamics.

Benchmark cases for a multi-component Lattice-Boltzmann method in hydrostatic conditions

Hydrostatic properties of partially saturated granular materials at the pore scale are evaluated by the lattice Boltzmann method (LBM) using Palabos implementation of the multi-component multiphase Shan-Chen model. Benchmark cases are presented to quantify the discretization errors and the sensitivity to geometrical and physical properties. This work offers practical guidelines to design LBM simulations of multiphase problems in porous media. Namely, a solid walls retraction procedure is proposed to reduce discretization errors significantly, leading to quadratic convergence. On this basis the equilibrium shapes of pendular bridges simulated numerically are in good agreement with the Young-Laplace equation. Likewise, entry capillary pressure and meniscus profiles in tubes of various cross-sectional shapes are in agreement with analytical predictions. The main points of this article are summarized as:•Benchmark cases for a multi-component Lattice-Boltzmann method are illustrated to be a guideline to calibrate the method in hydrostatic conditions.•A wall retraction procedure is introduced to minimize discretization errors.

Evaluation of methods using topology and integral geometry to assess wettability

HYPOTHESIS

The development of high-resolution in situ imaging has allowed contact angles to be measured directly inside porous materials. We evaluate the use of concepts in integral geometry to determine contact angle. Specifically, we test the hypothesis that it is possible to determine an average contact angle from measurements of the Gaussian curvature of the fluid/fluid meniscus using the Gauss-Bonnet theorem.

THEORY AND SIMULATION

We show that it is not possible to unambiguously determine an average contact angle from the Gauss-Bonnet theorem. We instead present an approximate relationship: 2πn(1-cosθ)=4π-∫κ_G12dS₁₂, where n is the number of closed loops of the three-phase contact line where phases 1 and 2 contact the surface, θ is the average contact angle, while κ_G12 is the Gaussian curvature of the fluid meniscus which is integrated over its surface S₁₂. We then use the results of pore-scale lattice Boltzmann simulations to assess the accuracy of this approach to determine a representative contact angle for two-phase flow in porous media.

FINDINGS

We show that in simple cases with a flat solid surface, the approximate expression works well. When applied to simulations on pore space images, the equation provides a robust estimate of contact angle, accurate to within 3°, when averaged over many fluid clusters, although individual values can have significant errors because of the approximations used in the calculation.

Pore-scale study of water salinity effect on thin-film stability for a moving oil droplet

HYPOTHESES

The interfacial dynamics in natural porous media are affected not only by the interplay between viscous and capillary forces but also the solid surface wettability. It has been hypothesized that the wettability alteration induced by changes in the water salinity is primarily caused by electric double-layer force expansion, which strongly affects the multiphase flow dynamics.

SIMULATIONS

We investigate the effect of water ionic composition and surface roughness on pore-scale wettability alteration. Multiphase hydrodynamics is numerically captured by a lubrication approximation describing the evolution of thin-films coupled with a multiscale level-set approach. An oil blob mobilized by water within a single pore is considered as a case study. The effect of brine ionic composition is accounted for by an electric double-layer through the water ionic strength and zeta-potential parameters.

FINDINGS

We demonstrate that high-salinity water thin-films collapse to an adsorbed nanometer layer, leading to a large pressure drop during mobilization of the blob induced by the attractive surface forces. However, low-salinity water thin-films are stable due to the repulsive electric double-layer forces, leading to less pressure drop during mobilization of the blob. The novelty of this work lies in efficiently capturing the nanoscale effects of the electric double-layer in pore-scale multiphase flow displacements. Our quantitative investigations provide fundamental insights into the efficiency of low-salinity waterflooding.

Embedded delta shocks

In 1977 Korchinski presented a new type of shock discontinuity in conservation laws. These singular solutions were coined δ-shocks since there is a time dependent Dirac delta involved. A naive description is that such δ-shock is of the overcompressive type: a single shock wave belonging to both families, the four characteristic lines of which impinge into the shock itself. In this work, we open the fan of solutions by studying two-family waves without intermediate constant states but possessing central rarefactions or comprising δ-shocks.

3D-Printed micro-optofluidic device for chemical fluids and cells detection

In this work, it is presented a micro-optofluidic flow detector used for on-chip biological and chemical samples investigation. It is made in Poly-dimethyl-siloxane using a master-slave approach based on the 3D-Printing techniques. The micro-optofluidic device is made by assembling a microfluidic T-junction with a micro-optical section that consists of two optical fiber insertions and a PDMS gold-spattered micro-waveguide. The working principle in the detection is based on a different light transmission correlated to the fluid interfering with the laser beam in a micro-channel section. The proposed solution allows to realize a PDMS micro-device taking the advantage of 3D- Printing and goes beyond the restriction in the material selection. The device's performances were tested in the fluids detection and in the evaluation of the cell concentrations. Additionally, the micro-device was used as a real-time two-phase fluids flow detector. The two-phases flows were successfully monitored in different experimental conditions, varying both hydrodynamic and optical external stimuli.

Two-phase flow in microfluidic-chip design of hydrodynamic filtration for cell particle sorting

As one of the flow-based passive sorting, the hydrodynamic filtration using a microfluidic-chip has shown to effectively separate into different sizes of subpopulations from cell or particle suspensions. Its model framework involving two-phase Newtonian or generalized Newtonian fluid (GNF) was developed, by performing the complete analysis of laminar flow and complicated networks of main and multiple branch channels. To predict rigorously what occurs in flow fields, we estimated pressure drop, velocity profile, and the ratio of the flow fraction at each branch point, in which the analytical model was validated with numerical flow simulations. As a model fluid of the GNF, polysaccharide solution based on Carreau type was examined. The objective parameters aiming practical channel design include the number of the branches and the length of narrow section of each branch for arbitrary conditions. The flow fraction and the number of branches are distinctly affected by the viscosity ratio between feed and side flows. As the side flow becomes more viscous, the flow fraction increases but the number of branches decreases, which enables a compact chip designed with fewer branches being operated under the same throughput. Hence, our rational design analysis indicates the significance of constitutive properties of each stream.

Micro-bubble flow simulation of dissolved air flotation process for water treatment using computational fluid dynamics technique

A dissolved air flotation (DAF) system is one of the water treatment processes that purifies contaminants through a buoyancy effect by attaching the moiety of micro-bubbles on their free surface. Since the DAF system was first used in the drinking water treatment in the 1960s, it has been recognized as an effective treatment for the water purification process. Most previous works laid great emphasis on the internal flow behaviors of fluid to improve the purification efficiency of the DAF system. Nevertheless, the practical implementation with a pilot plant indeed revealed some technical incompleteness for the DAF system. To circumvent for the technical incompleteness, numerical simulation based on computational fluid dynamics (CFD) has been carried out to understand the in-depth knowledge on internal flow phenomena in the DAF system. However, the standard k-ε turbulence model has been conventionally used in the most studies without any proper consideration process. Accordingly, the objectives of this study were to investigate the major effects on the internal flow behaviors for an efficient numerical simulation of DAF when a different turbulence model and micro-bubble parameters are used. As a result, the present study found that the standard k-ε model would be not proper for the internal flow simulation of the DAF process and a careful consideration would be required for a more accurate prediction. In addition, the present study examined a desirable internal flow pattern with various operating conditions of the micro-bubble. Consequently, the main findings of this study are expected to provide realistic information to related researchers for designing the DAF system with the optimal operating parameters.

Governing factors for preparation of silver nanoparticles using droplet-based microfluidic device

Droplet-based microfluidic devices are now more than ever used for the synthesis of nanoparticles with low polydispersity and well-defined properties suitable for various industrial applications. Very small reaction volumes (microlitre to femtolitre) and short diffusion lengths, provide superior mixing efficiency and heat transport. Both play the dominant role in case of ultra-fast chemical reactions triggered upon reactant mixing, e.g. preparation of colloidal silver by reduction of silver salt. The high sensitivity of these systems to process variables makes otherwise more straightforward batch-wise production prone to suffer from inconsistency and poor reproducibility, which has an adverse effect on the reliability of production and further particle utilisation. This work presents a rigorous description of microfluidic droplet formation, reactant mixing, and nanoparticle synthesis using CFD simulations and experimental methods. The reaction mixture inside of droplets was homogenized in less than 40 milliseconds, which has been confirmed by simulations. Silver nanoparticles produced by droplet-based microfluidic chip showed superior to batch-wise preparation in terms of both particle uniformity and polydispersity.

Local Concentrating, Not Shear Stress, That May Lead to Possible Instability of Protein Molecules During Syringe Injection: A Fluid Dynamic Study with Two-Phase Flow Model

A two-phase flow computational fluid dynamic (CFD) model was developed to study the hydrodynamic forces and the protein concentration changes of a protein solution in a syringe injector. Proteins were assumed to be nanosized solid spheres commensurate with their molecular weight and suspended in an aqueous environment, passing through the rapidly constricted sections of the syringe. Interaction between the solid and the liquid phase was taken into account, and four laminar flow cases were studied under the extensional flow. Profiles of pressure, velocity, and shear stress of the different cases were examined and compared. Hydrodynamic forces on a single protein particle were further analyzed. Our results indicate that the hydrodynamic forces are too small to affect significant conformational changes in proteins. The plunger rate showed limited impact on the distribution of protein particles inside the syringe. Nonetheless, the larger velocity gradient at the connection section of the hub toward the needle resulted in considerable accumulation of proteins. Such a concentrating effect may lead to protein aggregation and subsequent structural changes and will be examined in future studies.LAY ABSTRACT: Concentrated protein in the liquid form has become a preferred formulation strategy for delivering protein products, but it suffers from a high possibility of aggregation and precipitation, which may trigger a structural change and denaturation of the protein molecules and eventually cause the loss of the therapeutic functions of the protein products. To understand the effect of hydrodynamic forces on the change in local protein concentration in a syringe injector, we developed a two-phase flow computational fluid dynamic (CFD) model in this work. It was found that the local concentration of protein strongly depends on the velocity gradient of the fluid. A higher amount of proteins accumulates at the connection section of the hub toward the needle where the maximum velocity gradient is exhibited. The model results also imply a limited effect by hydrodynamic forces on possible unfolding of protein molecules. The local concentrating effect may enhance the possibility of molecular collision, resulting in aggregation and structural change.

Characterising the two-phase flow and mixing performance in a gas-mixed anaerobic digester: Importance for scaled-up applications

This study aimed to characterise the gas-liquid flow and mixing behaviour in a gas-mixed anaerobic digester by improving phase interaction modelling using Computational Fluid Dynamics (CFD). A 2D axisymmetric model validated with experimental data was set up using an Eulerian-Eulerian method. Uncertainty factors, including bubble size, phase interaction forces and liquid rheology were found to significantly influence the flow field. A more reliable and complete validation was obtained by critical comparison and assessment of the referred experimental data, compared to the models reported in other studies. Additionally, justifiable corrections and predictions in detail were obtained. Mixing was evaluated by trajectory tracking of a large number of particles based on an Euler-Lagrange method. The mixing performance approximated to a laminar-flow reactor (LFR) that distinctly deviated from expected continuous stirred tank reactor (CSTR) design, indicating limited enhancement from the applied gas-sparging strategy in the studied digester. The study shows the importance of a proper phase-interaction description for a reliable hydrodynamic characterisation and mixing evaluation in gas-mixed digesters. Validations, bend to experimental data without a critical assessment, may lead to an inaccurate model for further scaled-up applications.

800 fps neutron radiography of air-water two-phase flow

We have demonstrated dynamic cold neutron imaging of air-water two-phase flows up to 800 frames per second imaging rates. This has been achieved by using a high-efficiency (relatively thick) scintillator screen in combination with the highest available flux on a continuous spallation source and a high-speed sCMOS camera. This combination renders the spatial resolution to relatively modest value of about 0.5 mm, which is nevertheless sufficient for resolution of bubbles of the size down to about 1.0 mm in motion with unprecedented framerate using neutron imaging. We show the feasibility of the technique on the two-phase flow at ambient temperature and atmospheric pressure conditions, with the foreseen aim of measurements of two phase flows at high-temperatures and high pressures. It is also foreseen that the technique will be further utilized for quantification of the time-resolved instantaneous gas fraction and the gas phase velocity.

•

Demonstration of up to 800 frames per second dynamic cold neutron radiography.

•

Application of such technique for non-periodic (transient) process of bubbly flow in water.

•

Potential for quantification of (i) instantaneous gas volume fraction in dynamic two-phase flow and (ii) instantaneous gas phase velocimetry.

Simulation of immiscible liquid-liquid flows in complex microchannel geometries using a front-tracking scheme

The three-dimensional two-phase flow dynamics inside a microfluidic device of complex geometry is simulated using a parallel, hybrid front-tracking/level-set solver. The numerical framework employed circumvents numerous meshing issues normally associated with constructing complex geometries within typical computational fluid dynamics packages. The device considered in the present work is constructed via a module that defines solid objects by means of a static distance function. The construction combines primitive objects, such as a cylinder, a plane, and a torus, for instance, using simple geometrical operations. The numerical solutions predicted encompass dripping and jetting, and transitions in flow patterns are observed featuring the formation of drops, 'pancakes', plugs, and jets, over a wide range of flow rate ratios. We demonstrate the fact that vortex formation accompanies the development of certain flow patterns, and elucidate its role in their underlying mechanisms. Experimental visualisation with a high-speed imaging are also carried out. The numerical predictions are in excellent agreement with the experimental data.

Manufacturing a Micro-model with Integrated Fibre Optic Pressure Sensors

The measurement of fluid pressure inside pores is a major challenge in experimental studies of two-phase flow in porous media. In this paper, we describe the manufacturing procedure of a micro-model with integrated fibre optic pressure sensors. They have a circular measurement window with a diameter of 260  μ m , which enables the measurement of pressure at the pore scale. As a porous medium, we used a PDMS micro-model with known physical and surface properties. A given pore geometry was produced following a procedure we had developed earlier. We explain the technology behind fibre optic pressure sensors and the procedure for integrating these sensors into a micro-model and demonstrate their utility for the measurement of pore pressure under transient two-phase flow conditions. Finally, we present and analyse results of single and two-phase flow experiments performed in the micro-model and discuss the link between small-scale fast pressure changes with pore-scale events.

Numerical investigation of mist/air impingement cooling on ribbed blade leading-edge surface

The working gas turbine blades are exposed to the environment of high temperature, especially in the leading-edge region. The mist/air two-phase impingement cooling has been adopted to enhance the heat transfer on blade surfaces and investigate the leading-edge cooling effectiveness. An Euler-Lagrange particle tracking method is used to simulate the two-phase impingement cooling on the blade leading-edge. The mesh dependency test has been carried out and the numerical method is validated based on the available experimental data of mist/air cooling with jet impingement on a concave surface. The cooling effectiveness on three target surfaces is investigated, including the smooth and the ribbed surface with convex/concave columnar ribs. The results show that the cooling effectiveness of the mist/air two-phase flow is better than that of the single-phase flow. When the ribbed surfaces are used, the heat transfer enhancement is significant, the surface cooling effectiveness becomes higher and the convex ribbed surface presents a better performance. With the enhancement of the surface heat transfer, the pressure drop in the impingement zone increases, but the incremental factor of the flow friction is smaller than that of the heat transfer enhancement.

CO2 saturated water as two-phase flow for fouling control in reverse electrodialysis

When natural feed waters are used in the operation of a reverse electrodialysis (RED) stack, severe fouling on the ion exchange membranes and spacers occurs. Fouling of the RED stack has a strong influence on the gross power density output; which can decrease up to 50%. Moreover, an increase in the pressure loss occurs between the feed water inlet and outlet, increasing the pumping energy and thus decreasing the net power density that can be obtained. In this work, we extensively investigated the use of CO₂ saturated water as two-phase flow cleaning for fouling mitigation in RED using natural feed waters. Experiments were performed in the REDstack research facility located at the Afsluitdijk (the Netherlands) using natural feed waters for a period of 60 days. Two different gas combinations were experimentally investigated, water/air sparging and water/CO₂ (saturated) injection. Air is an inert gas mixture and induces air sparging in the stack. In the case of CO₂, nucleation, i.e. the spontaneous formation of bubbles, occurs at the spacer filaments due to depressurization of CO₂ saturated water, inducing cleaning. Results showed that stacks equipped with CO₂ saturated water can produce an average net power density of 0.18 W/m² under real fouling conditions with minimal pre-treatment and at a low outside temperature of only 8 °C, whereas the stacks equipped with air sparging could only produce an average net power density of 0.04 W/m². Electrochemical impedance spectroscopy measurements showed that the stacks equipped with air sparging increased in stack resistance due to the presence of stagnant bubbles remaining in the stack after every air injection. Furthermore, the introduction of CO₂ gas in the feed water introduces a pH decrease in the system (carbonated solution) adding an additional cleaning effect in the system, thus avoiding the use of environmentally unwanted cleaning chemicals.

Review-Physicochemical hydrodynamics of gas bubbles in two phase electrochemical systems

Electrochemical systems suffer from poor management of evolving gas bubbles. Improved understanding of bubbles behavior helps to reduce overpotential, save energy and enhance the mass transfer during chemical reactions. This work investigates and reviews the gas bubbles hydrodynamics, behavior, and management in electrochemical cells. Although the rate of bubble growth over the electrode surface is well understood, there is no reliable prediction of bubbles break-off diameter from the electrode surface because of the complexity of bubbles motion near the electrode surface. Particle Image Velocimetry (PIV) and Laser Doppler Anemometry (LDA) are the most common experimental techniques to measure bubble dynamics. Although the PIV is faster than LDA, both techniques are considered expensive and time-consuming. This encourages adapting Computational Fluid Dynamics (CFD) methods as an alternative to study bubbles behavior. However, further development of CFD methods is required to include coalescence and break-up of bubbles for better understanding and accuracy. The disadvantages of CFD methods can be overcome by using hybrid methods. The behavior of bubbles in electrochemical systems is still a complex challenging topic which requires a better understanding of the gas bubbles hydrodynamics and their interactions with the electrode surface and bulk liquid, as well as between the bubbles itself.

Pressure drop of two-phase helium along long cryogenic flexible transfer lines to support a superconducting RF operation at its cryogenic test stand

Background

Establishing a stand-alone cryogenic test stand is of vital importance to ensure the highly reliable and available operation of superconducting radio-frequency module in a synchrotron light source. Operating a cryogenic test stand relies strongly on a capability to deliver two-phase helium along long cryogenic transfer lines. A newly constructed cryogenic test stand with flexible cryogenic transfer lines of length 220 m at National Synchrotron Radiation Research Center is required to support a superconducting radio-frequency module operated at 126.0 kPa with a 40-W dynamic load for a long-term reliability test over weeks. It is designed based on a simple analytical approach with the introduction of a so-called tolerance factor that serves to estimate the pressure drops in transferring a two-phase helium flow with a substantial transfer cryogenic heat load. Tolerance factor 1.5 is adopted based on safety factor 1.5 commonly applied in cryogenic designs to estimate the total mass flow rate of liquid helium demanded. A maximum 60-W dynamic load is verified with experiment measured with heater power 60 W instead after the cryogenic test stand has been installed.

Results

Aligning the modeled cryogenic accumulated static heat load with the results measured in situ, actual tolerance factor 1.287 is obtained. The feasibility and validity of our simple analytical approach with actual tolerance factor 1.287 have been scrutinized by using five test cases with varied operating conditions. Calculated results show the discrepancies of the pressure drops between the estimated and measured values for both liquid helium and cold gaseous helium transfer lines have an underestimate 0.11 kPa and an overestimate 0.09 kPa, respectively. A discrepancy is foreseen, but remains acceptable for engineering applications from a practical point of view.

Conclusions

The simple analytical approach with the introduction of a tolerance factor can provide not only insight into optimizing the choice of each lossy cryogenic piping element of the transfer lines in the design phase but also firm guidance for upgrading the present cryogenic transfer lines for its subsequent application.

Numerical study of the faithful replication of micro/nanostructures on curved surfaces by the electrohydrodynamic instability process

This paper reports the numerical study of the one-step faithful replication of micro/nano-scale structures on a fiber surface by using the electrohydrodynamic instability patterning (EHDIP) process. By employing a rigorous numerical analysis method, conditions are revealed under which the faithful replication of a pattern can be achieved from a curved master electrode. It is found that the radius of curvature of the fiber plays an important role in determining the final morphology of the pattern when the destabilizing electric field is dominant in both the flat and patterned template cases. In general, stronger electric fields and larger radii of curvature of the substrate are favorable for the faithful replication of the pattern. In addition, theoretical analysis shows that higher aspect ratio of micro/nanostructures can be obtained on curved surfaces by using a master with a much lower aspect ratio. The results demonstrated in this study aims to provide guidelines for the faithful fabrication of micro/nanostructures on curved surfaces by the EHDIP process.

Droplet dynamics in rotating flows

This paper deals with investigations of droplet dynamics in rotating flows. In many previous studies droplet dynamics was analyzed in simple unidirectional flows. To fill this gap, the focus of this study is an overview on investigations of droplet dynamics in a complex rotating flow. A Lattice Boltzmann Method with high potential in simulation of two-phase unsteady flows is applied to simulate the physics of the problem in a lid-driven cavity. In spite of its simple geometry, there is a complex rotating flow field containing different vortices and shear regions. The Reynolds number based on the cavity length scale and the upper wall velocity, ReL, is considered to be 1000. We discuss here effects of different parameters such as: density ratios (1, 5, 10, 100, and 1000), droplet sizes (D/L=0.097, 0.114, 0.131 and 0.2), and droplet initial positions (1/8, 2/8, and 3/8 of the cavity length, L, out of center). The results are discussed in terms of global flow physics and its interaction with the droplet, drop deformation during its motion along with the main flow, and droplet trajectories. It is shown that there are strong interactions between the droplet and the main carrying flow. During motion, the droplets pass through different flow regions containing acceleration/deceleration zones. Consequently, the droplets experience different shear forces resulting in stretching, shrinking, rotating and dilatation which all contribute to the dynamics of the droplet.

Single- and two-phase flow simulation based on equivalent pore network extracted from micro-CT images of sandstone core

Due to the intricate structure of porous rocks, relationships between porosity or saturation and petrophysical transport properties classically used for reservoir evaluation and recovery strategies are either very complex or nonexistent. Thus, the pore network model extracted from the natural porous media is emphasized as a breakthrough to predict the fluid transport properties in the complex micro pore structure. This paper presents a modified method of extracting the equivalent pore network model from the three-dimensional micro computed tomography images based on the maximum ball algorithm. The partition of pore and throat are improved to avoid tremendous memory usage when extracting the equivalent pore network model. The porosity calculated by the extracted pore network model agrees well with the original sandstone sample. Instead of the Poiseuille's law used in the original work, the Lattice-Boltzmann method is employed to simulate the single- and two- phase flow in the extracted pore network. Good agreements are acquired on relative permeability saturation curves of the simulation against the experiment results.

A global sensitivity analysis of two-phase flow between fractured crystalline rock and bentonite with application to spent nuclear fuel disposal

Geological disposal of spent nuclear fuel in deep crystalline rock is investigated as a possible long term solution in Sweden and Finland. The fuel rods would be cased in copper canisters and deposited in vertical holes in the floor of deep underground tunnels, embedded within an engineered bentonite buffer. Recent experiments at the Äspö Hard Rock Laboratory (Sweden) showed that the high suction of unsaturated bentonite causes a de-saturation of the adjacent rock at the time of installation, which was also independently predicted in model experiments. Remaining air can affect the flow patterns and alter bio-geochemical conditions, influencing for instance the transport of radionuclides in the case of canister failure. However, thus far, observations and model realizations are limited in number and do not capture the conceivable range and combination of parameter values and boundary conditions that are relevant for the thousands of deposition holes envisioned in an operational final repository. In order to decrease this knowledge gap, we introduce here a formalized, systematic and fully integrated approach to study the combined impact of multiple factors on air saturation and dissolution predictions, investigating the impact of variability in parameter values, geometry and boundary conditions on bentonite buffer saturation times and on occurrences of rock de-saturation. Results showed that four parameters consistently appear in the top six influential factors for all considered output (target) variables: the position of the fracture intersecting the deposition hole, the background rock permeability, the suction representing the relative humidity in the open tunnel and the far field pressure value. The combined influence of these compared to the other parameters increases as one targets a larger fraction of the buffer reaching near-saturation. Strong interaction effects were found, which means that some parameter combinations yielded results (e.g., time to saturation) far outside the range of results obtained by the rest of the scenarios. This study also addresses potential air trapping by dissolution of part of the initial air content of the bentonite, showing that neglecting gas flow effects and trapping could lead to significant underestimation of the remaining air content and the duration of the initial aerobic phase of the repository.

Hydrogel-coated feed spacers in two-phase flow cleaning in spiral wound membrane elements: a novel platform for eco-friendly biofouling mitigation

Biofouling is still a major challenge in the application of nanofiltration and reverse osmosis membranes. Here we present a platform approach for environmentally friendly biofouling control using a combination of a hydrogel-coated feed spacer and two-phase flow cleaning. Neutral (polyHEMA-co-PEG10MA), cationic (polyDMAEMA) and anionic (polySPMA) hydrogels have been successfully grafted onto polypropylene (PP) feed spacers via plasma-mediated UV-polymerization. These coatings maintained their chemical stability after 7 days incubation in neutral (pH 7), acidic (pH 5) and basic (pH 9) environments. Anti-biofouling properties of these coatings were evaluated by Escherichia coli attachment assay and nanofiltration experiments at a TMP of 600 kPag using tap water with additional nutrients as feed and by using optical coherence tomography. Especially the anionic polySPMA-coated PP feed spacer shows reduced attachment of E. coli and biofouling in the spacer-filled narrow channels resulting in delayed biofilm growth. Employing this highly hydrophilic coating during removal of biofouling by two-phase flow cleaning also showed enhanced cleaning efficiency, feed channel pressure drop and flux recoveries. The strong hydrophilic nature and the presence of negative charge on polySPMA are most probably responsible for the improved antifouling behavior. A combination of polySPMA-coated PP feed spacers and two-phase flow cleaning therefore is promising and an environmentally friendly approach to control biofouling in NF/RO systems employing spiral-wound membrane modules.

Modeling blast waves, gas and particles dispersion in urban and hilly ground areas

The numerical simulation of shock and blast waves as well as particles dispersion in highly heterogeneous media such as cities, urban places, industrial plants and part of countries is addressed. Examples of phenomena under study are chemical gas products dispersion from damaged vessels, gas dispersion in urban places under explosion conditions, shock wave propagation in urban environment. A three-dimensional simulation multiphase flow code (HI2LO) is developed in this aim. To simplify the consideration of complex geometries, a heterogeneous discrete formulation is developed. When dealing with large scale domains, such as countries, the topography is considered with the help of elevation data. Meteorological conditions are also considered, in particular regarding complex temperature and wind profiles. Heat and mass transfers on sub-scale objects, such as buildings, trees and other obstacles are considered as well. Particles motion is addressed through a new turbulence model involving a single parameter to describe accurately plumes. Validations against experiments in basic situations are presented as well as examples of industrial and environmental computations.

High-frame rate imaging of two-phase flow in a thin rectangular channel using fast neutrons

We have demonstrated the feasibility of performing high-frame-rate, fast neutron radiography of air-water two-phase flows in a thin channel with rectangular cross section. The experiments have been carried out at the accelerator facility of the Physikalisch-Technische Bundesanstalt. A polychromatic, high-intensity fast neutron beam with average energy of 6 MeV was produced by 11.5 MeV deuterons hitting a thick Be target. Image sequences down to 10 ms exposure times were obtained using a fast-neutron imaging detector developed in the context of fast-neutron resonance imaging. Different two-phase flow regimes such as bubbly slug and churn flows have been examined. Two phase flow parameters like the volumetric gas fraction, bubble size and mean bubble velocities have been measured. The first results are promising, improvements for future experiments are also discussed.