Simulating Two-Phase Seepage in Undisturbed Soil Based on Lattice Boltzmann Method and X-ray Computed Tomography Images.

The two-phase seepage fluid (i.e., air and water) behaviors in undisturbed granite residual soil (U-GRS) have not been comprehensively studied due to a lack of accurate and representative models of its internal pore structure. By leveraging X-ray computed tomography (CT) along with the lattice Boltzmann method (LBM) enhanced by the Shan-Chen model, this study simulates the impact of internal pore characteristics of U-GRS on the water-gas two-phase seepage flow behaviors. Our findings reveal that the fluid demonstrates a preference for larger and straighter channels for seepage, and as seepage progresses, the volume fraction of the water/gas phases exhibits an initial increase/decrease trend, eventually stabilizing. The results show the dependence of two-phase seepage velocity on porosity, while the local seepage velocity is influenced by the distribution and complexity of the pore structure. This emphasizes the need to consider pore distribution and connectivity when studying two-phase flow in undisturbed soil. It is observed that the residual gas phase persists within the pore space, primarily localized at the pore margins and dead spaces. Furthermore, the study identifies that hydrophobic walls repel adjacent fluids, thereby accelerating fluid movement, whereas hydrophilic walls attract fluids, inducing a viscous effect that decelerates fluid flow. Consequently, the two-phase flow rate is found to increase with then-enhanced hydrophobicity. The apex of the water-phase volume fraction is observed under hydrophobic wall conditions, reaching up to 96.40%, with the residual gas-phase constituting 3.60%. The hydrophilic wall retains more residual gas-phase volume fraction than the neutral wall, followed by the hydrophobic wall. Conclusively, the investigations using X-ray CT and LBM demonstrate that the pore structure characteristics and the wettability of the pore walls significantly influence the two-phase seepage process.

Analysis of 'Investigating an extended multiphase flow model that includes specific interfacial area', Computer Methods in Applied Mechanics and Engineering, 418:116594, 2024.

Comments are provided on the recent paper by Ebadi et al. [3], which demonstrates that the formulated model that was solved contains misconceptions or errors that render the work unsuitable for describing the evolution of interfacial areas in two-fluid porous medium systems. The need for kinematic equations is described and components of a theoretically consistent approach are summarized.

Squeezing Droplet Formation in a Flow-Focusing Micro Cross-Junction.

Motivated by the increasing need of optimised micro-devices for droplet production in medical and biological applications, this paper introduces an integrated approach for the study of the liquid-liquid droplet creation in flow-focusing micro cross-junctions. The micro-junction considered is characterised by a restriction of the channels cross-sections in the junction, which has the function of focusing the flow in the region of the droplet formation. The problem is studied numerically in the OpenFOAM environment and validated by a comparison with experimental results obtained by high-speed camera images and micro-PIV measurements. The analysis of the forces acting on the dispersed phase during the droplet formation and the diameter of the droplets obtained numerically are considered for the development of a model of the droplet breakup under the squeezing regime. On the basis of energy balancing during the breakup, a relation between interfacial tension, the size of the cross-sections in the junction, and the time interval needed for droplet creation is obtained, which yields a novel correlation between the dimensionless length of the droplet and the dimensionless flow rate. This research expands our knowledge of the phenomenon of drop creation in micro-junctions with restrictions providing new aid for the optimal design of micro-drop generators.

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.

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.

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.

A numerical study on machine-learning-based ultrasound tomography of bubbly two-phase flows.

Ultrasound tomography (UT) of bubbly two-phase flows using machine learning (ML) was investigated by performing two-dimensional ultrasound numerical simulations using a finite element method simulator. Studies on UT for two-phase flow measurements have been conducted only for some bubbles. However, in an actual bubbly flow, numerous bubbles are complexly distributed in the cross-section of the flow channel. This limitation of previous studies originates from the transmission characteristics of ultrasound waves through a medium. The transmission characteristics of ultrasound waves differ from those of other probe signals, such as radiation, electrical, and optical signals. This study evaluated the feasibility of combining UT with ML for predicting dense bubble distributions with up to 20 bubbles (cross-sectional average void fraction of approximately 0.29). We investigated the effects of the temporal length of the received waveform and the number of sensors to optimize the system on the prediction performance of the bubble distribution. The simultaneous driving of the installed sensors was simulated to reduce the measurement time for the entire cross-section and verify the method's applicability. Thus, it was confirmed that UT using ML has sufficient prediction performance, even for a complex bubble distribution with many bubbles, and that the cross-sectional average void fraction can be predicted with high accuracy.

Numerical study of laser-induced cavitation bubble with consideration of chemical reactions.

Cavitation generated during injector jetting can significantly affect fuel atomization. Laser-induced cavitation bubble is an important phenomenon in laser induced plasma ignition technology. Limited by the difficulties in experimental measurements, numerical simulations have become an important tool in the study of laser-induced cavitation bubble, but most previous numerical models used to study the dynamics of laser-induced cavitation bubble usually ignore the effect of chemical reactions. In this study, the finite volume method is used to solve the compressible two-dimensional reynolds averaged Navier-Stokes equation by considering the heat and mass transfer as well as the chemical reactions within the cavitation bubble. The effects of overall reaction and elementary reactions on the cavitation bubble are evaluated, respectively. It is found that by additionally considering chemical reactions within the numerical model, lower maximum temperatures and higher maximum pressures are predicted within the bubble. And the generated non-condensable gases produced by the chemical reactions enhance the subsequent expansion process of the cavitation bubble. Besides, the effect of the one-sided wall boundary condition on cavitation bubble is compared with the infinite boundary condition. Influenced by the wall boundary, the cavitation bubble forms a localized high pressure on the side of the bubble away from the wall during the collapse process, which causes the bubble to be compressed into a "crescent" shape. The maximum pressure and temperature inside the bubble are lower due to localized losses caused by the wall.

Modeling liquid rate through wellhead chokes using machine learning techniques.

Precise measurement and prediction of the fluid flow rates in production wells are crucial for anticipating the production volume and hydrocarbon recovery and creating a steady and controllable flow regime in such wells. This study suggests two approaches to predict the flow rate through wellhead chokes. The first is a data-driven approach using different methods, namely: Adaptive boosting support vector regression (Adaboost-SVR), multivariate adaptive regression spline (MARS), radial basis function (RBF), and multilayer perceptron (MLP) with three algorithms: Levenberg-Marquardt (LM), bayesian-regularization (BR), and scaled conjugate gradient (SCG). The second is a developed correlation that depends on wellhead pressure (Pwh), gas-to-liquid ratio (GLR), and choke size (Dc). A dataset of 565 data points is available for model development. The performance of the two suggested approaches is compared with earlier correlations. Results revealed that the proposed models outperform the existing ones, with the Adaboost-SVR model showing the best performance with an average absolute percent relative error (AAPRE) of 5.15% and a correlation coefficient of 0.9784. Additionally, the results indicated that the developed correlation resulted in better predictions compared to the earlier ones. Furthermore, a sensitivity analysis of the input variable was also investigated in this study and revealed that the choke size variable had the most significant effect, while the Pwh and GLR showed a slight effect on the liquid rate. Eventually, the leverage approach showed that only 2.1% of the data points were in the suspicious range.

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.

Insights into two-phase flow dynamics in closed-loop pulsating heat pipes utilizing Fe3O4/water: experimental visualization study.

This article discusses a focused study on visualizing the flow patterns in a two-phase pulsating heat pipe (PHP) using Fe3O4/water as the working fluid at 3 V/V% concentration. The research also aims to meticulously examine phase change phenomena in the heating section, particularly focusing on bubble formation and expansion processes. A high-speed video camera was utilized to capture dynamic insights into the behavior of the Fe3O4/water mixture. Based on the findings, a straightforward model was developed to explain bubble generation and growth in the mixture, serving as a useful reference for future PHP designs and optimizations. Visual observations also noted the stable nature of the Fe3O4/water nanofluid over a 4-day period, confirming its consistency throughout the experiments. Moreover, the impact of heat load variation on the evaporator section was assessed using controlled heat inputs ranging from 10 to 80 W. Observations on the arrangement of slugs and plugs at a 50% filling ratio revealed interesting self-adjusting flow patterns in response to increasing heat inputs, providing valuable insights into PHP operational dynamics. Notably, the oscillatory flow behavior of Fe3O4/water, the chosen working fluid, exhibited greater activity in comparison to water. This distinctive flow behavior contributed to achieving heightened thermal performance efficiency for the Fe3O4/water system, attributed to its faster attainment of the annular flow condition.

Precipitation and evaporation affecting landfill gas migration into passive methane oxidation biosystems: Models development and verification.

Passive methane oxidation biosystems (PMOBs) are developed as an innovative and cost-effective solution to reduce methane (CH4) emissions from municipal solid waste landfills. A PMOB consists of a methane oxidation layer (MOL) and an underlying gas distribution layer (GDL). The length of unrestricted gas migration (LUGM) has been recently proposed as the design criterion for PMOBs where the LUGM is calculated as the horizontal length along the MOL-GDL interface with the volumetric gas content (θa) exceeding the threshold volumetric gas content (θa,occ). This paper examined water and gas migration within three PMOBs with different MOL-GDL interfaces subject to precipitation and evaporation using verified numerical models. The results show that the use of a single-phase flow model underestimates the LUGM values of the PMOB for heavy precipitation events, and a two-phase flow model should be used to calculate both the LUGM and the total gas mass flow rate into the MOL when designing PMOBs. Both zig-zag and trapezoidal MOL-GDL interfaces can redistribute the gas mass flow rate at the MOL-GDL interface, while the trapezoidal MOL-GDL interface slightly outperforms the zig-zag MOL-GDL interface for enhancing the total gas mass flow rate into the MOL when comparing with the planar MOL-GDL interface. The zig-zag and trapezoidal MOL-GDL interfaces allow gas migration in the upper part of each PMOB segment even when the lower part of each PMOB segment was filled with water, and thus have a potential to minimize hotspot formation.

The impact of 180° return bend inclination on pressure drop characteristics and phase distribution during oil-water flow.

The present work aims to capture the influence of the inclination of the return bend on flow patterns and pressure drop during oil-water flow. The experiments were carried out for different inclinations (0°, 15°, 30°, and 45°) of return bend for various superficial velocity combinations of oil (kerosene) and water ranging from 0.07 to 0.66 m/s. The experiments showed that pressure drop increases with the increase in inclination. However, the pressure drop at a fixed inclination (say 15°) decreases with the increase in the superficial velocity of the water. Distinct flow patterns observed in the return bend were droplet flow, film inversion, slug flow, plug flow and large slug flow. Droplet flow dominates at the lower range of kerosene (i.e., Usk = 0.07-0.2 m/s) and higher range of water superficial velocity (i.e., Usw = 0.40-0.66 m/s) at all the inclinations considered in this study. Additionally, comparisons between the experimental and numerical simulation results were made. The numerical solution utilized the Euler-Euler approach, considering the different phases as interpenetrating continua. The Volume of Fluid (VOF) model was used within this approach, monitoring the volume fraction of each phase over the domain while calculating one set of momentum equations for each phase. To capture the turbulent effects accurately, the k-ε turbulence model was incorporated. It happened to be found that the numerical findings showed remarkable agreement with the experimental data.

Continuous Sampling of Aerosolized Particles Using Stratified Two-Phase Microfluidics.

Health and security concerns have made it essential to develop integrated, continuous collection and sensing platforms that are compact and capable of real-time detection. In this study, we numerically investigate the flow physics associated with the single-step collection and enrichment of aerosolized polystyrene microparticles into a flowing liquid using a stratified air-water flow in a U-shaped microchannel. We validate our simulation results by comparing them to experimental data from the literature. Additionally, we fabricate an identical microfluidic device using PDMS-based soft lithography and test it to corroborate the previously published experimental data. Diversion and entrapment efficiencies are used as evaluation metrics, both of which increase with increasing particle diameter and superficial air inlet velocity. Overall, our ANSYS Fluent two-dimensional (2D) and three-dimensional (3D) multiphase flow simulations exhibit a good agreement with our experimental data and data in the literature (average deviation of ∼11%) in terms of diversion efficiency. Simulations also found the entrapment efficiency to be lower than the diversion efficiency, indicating discrepancies in the literature in terms of captured particles. The effect of the Dean force on the flow physics was also investigated using 3D simulations. We found that the effect of the Dean flow was more dominant relative to the centrifugal force on the smaller particles (e.g., 0.65 µm) compared to the larger particles (e.g., 2.1 µm). Increasing the superficial air inlet velocity also increases the effect of the centrifugal forces relative to the Dean forces. Overall, this experimentally validated multiphase model decouples and investigates the multiple and simultaneous forces on aerosolized particles flowing through a curved microchannel, which is crucial for designing more efficient capture devices. Once integrated with a microfluidic-based biosensor, this stratified flow-based microfluidic biothreat capture platform should deliver continuous sensor-ready enriched biosamples for real-time sensing.

The Identification of Leidenfrost Phenomenon Formation on TiO2-Coated Surfaces and the Modelling of Heat Transfer Processes.

Experiments on specimen cooling dynamics and possible film boiling around a body are very important in various industrial applications, such as nucleate boiling, to decrease drag reduction or achieve better surface properties in coating technologies. The objective of this study was to investigate the interaction between the heat transfer processes and cooling dynamics of a sample in different boundary conditions. This article presents new experimental data on specimens coated with Al-TiO2 film and Leidenfrost phenomenon (LP) formation on the film's surface. Furthermore, this manuscript presents numerical heat and mass transfer parameter results. The comparative analysis of new experiments on Al-TiO2 film specimens and other coatings such as polished aluminium, Al-MgO, Al-MgH2 and Al-TiH2 provides further detail on oxide and hydride materials. In the experimental cooling dynamics experiments, specimens were heated up to 450 °C, while the sub-cooling water temperatures were 14*‒20 °C (room temperature), 40 °C and 60 °C. The specimens' cooling dynamics were calculated by applying Newton's cooling law, and heat transfer was estimated by calculating the heat flux q transferred from the specimens' surface and the Bi parameter. The metadata results from the performed experiments were used to numerically model the cooling dynamics curves for different material specimens. Approximated polynomial equations are proposed for the polished aluminium, Al-TiO2, Al-MgO, Al-MgH2 and Al-TiH2 materials. The provided comparative analysis makes it possible to see the differences between oxides and hydrides and to choose materials for practical application in the industrial sector. The presented results could also be used in software packages to model heat transfer processes.

Simultaneous carbon, nitrogen and phosphorus removal in sequencing batch membrane aerated biofilm reactor with biofilm thickness control via air scouring aided by computational fluid dynamics.

Membrane aerated biofilm reactor (MABR) is challenged by biofilm thickness control and phosphorus removal. Air scouring aided by computational fluid dynamics (CFD) was employed to detach outer biofilm in sequencing batch MABR treating low C/N wastewater. Biofilm with 177-285 µm thickness in cycle 5-15 achieved over 85 % chemical oxygen demand (COD) and total inorganic nitrogen (TIN) removals at loading rate of 13.2 gCOD/m2/d and 2.64 gNH4+-N/m2/d. Biofilm rheology measurements in cycle 10-25 showed yield stress against detachment of 2.8-7.4 Pa, which were equal to CFD calculated shear stresses under air scouring flowrate of 3-9 L/min. Air scouring reduced effluent NH4+-N by 10 % and biofilm thickness by 78 µm. Intermittent aeration (4h off, 19.5h on) and air scouring (3 L/min, 30 s before settling) in one cycle achieved COD removal over 90 %, TIN and PO43--P removals over 80 %, showing great potential for simultaneous carbon, nitrogen and phosphorus removals.

Review on Bubble Dynamics in Proton Exchange Membrane Water Electrolysis: Towards Optimal Green Hydrogen Yield.

Water electrolysis using a proton exchange membrane (PEM) holds substantial promise to produce green hydrogen with zero carbon discharge. Although various techniques are available to produce hydrogen gas, the water electrolysis process tends to be more cost-effective with greater advantages for energy storage devices. However, one of the challenges associated with PEM water electrolysis is the accumulation of gas bubbles, which can impair cell performance and result in lower hydrogen output. Achieving an in-depth knowledge of bubble dynamics during electrolysis is essential for optimal cell performance. This review paper discusses bubble behaviors, measuring techniques, and other aspects of bubble dynamics in PEM water electrolysis. It also examines bubble behavior under different operating conditions, as well as the system geometry. The current review paper will further improve the understanding of bubble dynamics in PEM water electrolysis, facilitating more competent, inexpensive, and feasible green hydrogen production.