Discrete Droplets Research Articles

Photofabricated channel-digital microfluidics (pCDMF): A promising lab-on-a-chip platform for fully integrated digital PCR

Microfluidics-based digital PCR (dPCR) is a potent technique capable of absolute quantification of nucleic acids with ultrahigh sensitivity. However, an indispensable process, nucleic acid extraction (NAE) remains hard to be integrated with dPCR into a monolithic device, because both of them involve complicated processes and exquisite fluid manipulation. An ideal solution to this issue is to integrate the two main microfluidic platforms, channel-based microfluidics (CMF) and digital microfluidics (DMF), which, however, remains very challenging. To address these challenges, we herein develop a photofabricated channel-digital microfluidic chip (pCDMF) with significant advantages including: the fabrication of pCDMF is very convenient due to the flexible material choices of photocuring, without tedious process of hydrophobic coating; a parallel-plate configuration can be easily formed, enabling all basic functions of DMF; various channels and microstructures can be freely designed, allowing for fluidic functions of CMF. It is therefore the first example of full integration of DMF and CMF as a monolithic device, with capabilities to precisely handle both discrete droplets and continuous fluids. Finally, an all-in-one device for sample-to-answer digital PCR has been successfully demonstrated. We also envisioned that the pCDMF may open an alternative route toward the fully integrated lab-on-a-chip for various biological applications.

Discrete quantum droplets in one-dimensional binary Bose–Einstein condensates

We discuss the propagation and interaction of discrete quantum droplets in one-dimensional binary Bose mixtures by variationally and numerically solving the underlying generalized discrete Gross–Pitaevskii equation. Useful analytical expressions are obtained for the energy, the chemical potential, the density profile, and the droplet separation in the regime of a small droplet. We show that for weak lattice coupling, the stationary state of small droplets remains intact and does not display any special structures. An irregular dynamics is observed in both single and dual droplets with flat-top plateau in the limit of relatively strong coupling. Our results reveal also that two droplets interact weakly and cannot form a bound-state due to the intriguing interplay of the discreteness and quantum fluctuations.

Characterizing infection risk in a restaurant environment due to airborne diseases using discrete droplet dispersion simulations

The use of masks as a measure to control the spread of respiratory viruses has been widely acknowledged. However, there are instances where wearing a mask is not possible, making these environments potential vectors for virus transmission. Such environments can contain multiple sources of infection and are challenging to characterize in terms of infection risk. To address this issue, we have developed a methodology to investigate the role of ventilation in reducing the infection risk in such environments. We use a restaurant setting as a representative scenario to demonstrate the methodology. Using implicit large eddy simulations along with discrete droplet dispersion modeling we investigate the impact of ventilation and physical distance on the spread of respiratory viruses and the risk of infection. Our findings show that operating ventilation systems, such as mechanical mixing and increasing physical distance between subjects, can significantly reduce the average room infection risk and number of newly infected subjects. However, this observation is subject to the transmissibility of the airborne viruses. In the case of a highly transmissible virus, the use of mechanical mixing may be inconsequential when compared to only fresh air ventilation. These findings provide valuable insights into the mitigation of infection risk in situations where the use of masks is not possible.

Simulation of Slag Droplet Entrainment by Volume of Fluid and Lagrangian Particle Tracking Coupled Modeling.

The slag droplet entrainment is a common phenomenon in steel refining processes, which may lead to inclusions and defects. In the multiphase flow system, the distinct interface and tiny blobs possess a wide range of spatial and temporal scales and make it hard to be simulated. In numerical methods, the volume of fluid (VOF) approach is appropriate for capturing the interface, but for the unresolvable tiny blobs, the Lagrangian particle tracking (LPT) is preferable. This work newly implements a bidirectional VOF-LPT transformation algorithm for developing a multiscale solver in OpenFOAM to simulate the slag droplet entrainment. The interIsoFoam solver is selected as the main solver to resolve the interface, and the resolution is improved with using the geometric reconstruction and the adaptive mesh refinement (AMR). For capturing tiny droplets, a connected component labeling (CCL) method is adopted for detecting discrete droplets in the VOF field, and then the VOF-to-LPT transition takes place for saving computational costs. Conversely, the LPT-to-VOF transformation for droplets touching the interface is also incorporated to achieve the bidirectional transition. The solver is first validated by a simple case, indicating that the two-way transition algorithm and the Eulerian-Lagrangian momentum coupling are accurate. Then the solver is applied to simulate the slag layer behavior for revealing the mechanisms of slag droplet formation and entrainment. Two main mechanisms of slag droplet formation are identified, and it is found that fewer discrete droplets are generated when the surface tension increases.

Measurement of a Three-Dimensional Rotating Flow Field and Analysis of the Internal Oil Droplet Migration

Investigating the motion of discrete oil droplets in a rotating flow field can provide a theoretical basis for optimizing the flow field and structural parameters of hydrocyclones and centrifugal separation equipment. In the present work, the particle image velocimetry (PIV) method was applied to study the velocity distribution of a three-dimensional axial-rotor-driven rotating flow field and the influence of the velocity distribution of different rotor speeds on the flow field. The radial migration of oil droplets with different particle sizes in the rotating flow field was visually analyzed using high-speed video (HSV). The results showed that the oil droplets with the same radial position had diameters of 2.677 and 4.391 mm, whereas the movement times to the axis were 0.902 and 0.752 s. The larger the oil droplet size, the shorter the time to move to the axial center of the rotating flow field. The radial velocities of oil droplets with diameters of 2.677 and 2.714 mm were 0.0221 and 0.02 m·s−1, respectively. In addition, a mathematical expression was established between the radial migration time and the oil droplet size in the rotating flow field. The accuracy of the proposed expression was verified using experiments.

Prediction of carbon capture and sequestration (CCS) technology in a 125 MW tangentially coal-fired subcritical thermal power plant for retrofitting in Bangladesh

The main objective of this study is to assess the use of carbon capturing and sequestration (CCS) technology in the existing plant to reduce the emissions leading to retrofitting. A typical 125 MW tangentially coal-fired Barapukuria Thermal Power Plant (BTPP), Bangladesh has been scrutitinized considering air-firing (existing) and five different oxy-firing (OF) cases have been numerically investigated. For predicting the particle flow, radiative heat transfer and turbulence effects, discrete droplet method (DDM), discrete transfer radiation method (DTRM) and k–ϵ model is considered respectively. For homogeneous reaction for the devolatilization and heterogeneous reactions for the char oxidation process, a two steps reaction mechanism is considered. Model has been validated by comparing the predicted flue gas temperature against the plant data at different positions of the furnace and reasonable agreement have been observed. This study assessed the flow characteristics, thermal behavior and species generations of the selected plant under 80% fuel loading. It is predicted that with the increase of cases from air-firing, level of oxiding elements are increasing in the main combustion zone due to enriched O2 level leading to complete combustion of the particles. It is also found that for OF 31 case, level of CO2 is close to the air-firing case. it has been found that wall temperature is comparatively higher for 100% loading and comparatively lower for 50% loading for almost all similar type of air or oxy-firing conditions. For 100% load, the wall temperature has been reached around 2900 K for OF31 case. It is predicted that air operation equivalent to oxy-firing have been found at < OF 30+ % (Flame temp), at = OF 30+ % (Convec Heat transfer) and at > OF 30+ % (radiative Heat transfer) respectively.

A numerical study on the interaction of droplet collisions and air flow impact in cross-impinging spray

The spray impingement method has been widely used to effectively achieve better atomization of liquid fuel. The factors influencing the formation of impinged spray are analysed using computational fluid dynamics (CFD) modelling. Free spray and impinging sprays with three impinging distances are simulated using the Eulerian-Lagrangian multiphase method. The results are validated with experimental data. Meanwhile, some advanced methods are employed to obtain better results of simulation such as the AMR technique and adaptive collision mesh method. The spray widths in different impingement distances are compared to analyse the distribution of spray droplets. Upon the binary collision theory, the distribution of discrete droplets is attributed to the collisions of huge amounts of particles. In the process of spray impingement, the impact of air flow, generated by gas entrainment of spray, also influences the spray. As a supplement to the collisions of discrete droplets, the impingement of continuous air flow and their interaction is mainly discussed. The impact waves observed in the flow field demonstrate that the impingement of entrained air flow is similar to jet impingement.

Experimental analysis of drop size distribution and nucleation site density during dropwise condensation from humid air flow

Condensation of the water vapor present in the air is a heat and mass transfer process encountered in many applications as humid air dehumidification and water harvesting. Depending on the wettability characteristics of the surface, condensation can take place in filmwise mode or in dropwise mode with the formation of discrete liquid droplets over the condensing surface. While dropwise condensation (DWC) of pure steam was found to promote a considerable enhancement of the heat transfer compared to filmwise condensation, when dealing with humid air DWC more investigation is needed. Modeling of DWC from humid air requires the calculation of the heat flow rate through a single droplet and the determination of the drop-size distribution. The heat exchanged through a single droplet depends on the heat and mass transfer resistances, while the drop-size distribution is also affected by nucleation site density and droplets mobility. Therefore, to better understand the DWC phenomenon with humid air and for the validation of the models, it is necessary to measure the heat flux (total and latent), droplet population and nucleation site density. In the present work, condensation tests from humid air are performed over two square (40 mm x 40 mm) aluminum samples that display different wettability. The experimental apparatus consists of a closed air loop with two main components: the environmental chamber and the test chamber. The air is conditioned in the environmental chamber and then it flows inside the test section where the vapor present in the humid air is condensed over the vertical metallic sample. Two variable speed fans are used to circulate the air. The test section is designed for heat and mass transfer measurements and for simultaneous visualization of the condensation process. As a peculiar characteristic of the present experimental technique, all the test section assembly is suspended on a high precision balance allowing a precise measurement of the mass of condensate. The effect of surface wettability on the heat and mass transfer during DWC is investigated. Time-lapse videos of the condensation process are acquired at different magnifications. By using a homemade MATLAB® program for droplet detection, recorded images are analysed allowing the determination of both the drop size density distribution (small and large droplet population) and the nucleation sites density.

Settling behavior of polydisperse droplets in homogeneous isotropic turbulence

The settling behavior of polydisperse droplets in homogeneous and isotropic turbulence was measured by an ultra-high-resolution two-dimensional Particle Image Velocimetry. The aim of the present study is to provide new insight on the dependence of multi-scale particle settling behavior on characteristic parameters of two-phase turbulent flow via a sophisticate conditional analysis. The relative settling strength (defined as the ratio of mean droplet settling velocity to root mean square velocity of turbulence), whose effect on droplet settling behavior is of the primary interest, ranges as SvL=0.5–2.0. The turbulence Taylor Reynolds number is Reλ=200–300, and the droplet Stokes number is Stp=0.1–10. Voronoï analysis is performed to obtain the concentration field of discrete droplets from particle images. Particle structures including clusters or voids are detected, and the droplet settling velocities corresponding to various probing conditions, such as Stp, local particle concentration, and size of particle structures, were then analyzed. For the present configuration (droplet net sedimentation), there is a non-monotonic dependency of the settling velocity on local particle concentration. The negative correlation between them occurs in the moderate-concentration sub-regime and is insensitive to the variation of SvL, in which individual droplets interact with turbulent flow independently. It can be well explained by the commonly invoked preferential sweeping mechanisms. On the other hand, the dense-concentration regime, in which droplets prefer to accumulate into clusters, presents a positive correlation; namely, the conditional-averaged settling velocity decreases with the increase in local particle concentration. In this sub-regime, it is not the scale of single particles but the scale of particle clusters and the relative strength of turbulence (measured by SvL) that jointly determines the droplet settling behavior. Such a process, to our knowledge, is consistent with the so-called multi-scale preferential sweeping effect.

Automated and label-free detection of HIV DNA via digital microfluidics-chemiluminescence analysis

Automated and sensitive detection of nucleic acid is crucial for reliable and safe disease diagnostics because the minimum manual operation can effectively avoid sample contamination and infection risk. Herein, a sensitive, label-free and automated method was developed for detecting trace DNA by integrating digital microfluidics (DMF) with chemiluminescence (CL) detection. DMF can automatically manipulate discrete droplets. Taking advantage of sensitive and simple CL detection, DMF-CL integration can automatically and sensitively detect HIV DNA via the preset pathway. To further improve the sensitivity, on-chip exponential amplification reaction (EXPAR) was explored to produce large quantities of double-stranded DNA (dsDNA) in the presence of HIV DNA. Then, SYBR Green I (SG) intercalated into dsDNA to generate singlet oxygen upon light irradiation, which in situ oxidized potassium ferrocyanide into potassium ferricyanide. Finally, the CL reaction of potassium ferricyanide and luminol ensured the label-free, reliable and specific detection of HIV DNA in the range of 10 pM to 2 nM. Saving time-consuming and laborious manual operations, the DMF-CL integration greatly reduced the reagent consumption and avoided the contamination of DNA sample as well. Therefore, DMF-CL offers a promising and useful technology for developing automated and reliable biosensing and bioanalysis.

High-throughput screening of microbial strains in large-scale microfluidic droplets.

The transformation of engineered microbial cells is a pivotal link in green biomanufacturing. Its distinctive research application involves genetic modification of microbial chassis to impart targeted traits and functions for effective synthesis of the desired products. Microfluidics, as an emerging complementary solution, focuses on controlling and manipulating fluid in channels at the microscopic scale. One of its subcategories is droplet-based microfluidics (DMF), which can generate discrete droplets using immiscible multiphase fluids at kHz frequencies. To date, droplet microfluidics has been successfully applied to a variety of microbes, including bacteria, yeast, and filamentous fungi, and the detection of massive metabolites of strain products, such as polypeptides, enzymes, and lipids, has been realized. In summary, we firmly believe that droplet microfluidics has evolved into a powerful technology that will pave the way for high-throughput screening of engineered microbial strains in the green biomanufacturing industry.

Fine structural aspects on the web glue production in the golden orb-web spider Trichonephila clavata

ABSTRACT The water-soluble glue substance of the capture threads in Trichonephila clavata is solely produced from two pairs of aggregate silk glands. During the web glue production, secretory vesicles were synthesized via the extensive rough endoplasmic reticulum of epithelial cells. Unlike the clearly described fibrous web production in spiders, the process of aqueous web glue production appears to involve either a condensing or a packaging step by the Golgi complex. In particular, the fine structure of secretory vesicles varies from cell to cell and may represent the secretory cycle. The electron-dense multivesicular bodies were clearly visible as discrete droplets, and the mature secretory product in the glandular epithelium appeared as a spherical vacuole grown by fusion with surrounding small vesicles. Our fine structural observation reveals that the secretion occurs when the release of secreted material involves the loss of part of the cytoplasm. The bleb along the luminal surface of the secretory cells and membrane-bound extracellular vesicles which pinched off from the cell suggests that the secretory product is released by the mechanism of apocrine secretion.

Condensation heat transfer on slippery lubricant-impregnated nanohole-arrayed stainless-steel surface

The condensation heat transfer performance of lubricant-impregnated surfaces has a competitive advantage owing to the improved mobility of droplets and self-healing characteristics. In this study, we fabricated the lubricant-impregnated surfaces featured with Teflon-coated anodized nanoporous oxide layer on stainless steel infused with Krytox GPL103 oil. The characterization shows that the droplet state on the surface changes from the pinning to slippery states (sliding angle α of 1.8 ± 0.3°) enabled by the presence of the lubricant oil. The condensation heat transfer performances of the lubricant-impregnated Teflon-coated anodized stainless-steel surface (L-T-AS) are investigated from the aspects of the surface wettability and droplet mobility, compared to that of Teflon-coated stainless-steel surface (T-S) and Teflon-coated anodized stainless-steel surface (T-AS). Whereas coalescence-induced droplet shedding and droplet sweeping are found on the T-S and T-AS, the shedding of discrete droplets and their sweeping without further coalescence are noticed on the L-T-AS. The heat transfer coefficient of the L-T-AS is 398.4% higher than that of the T-AS at the low subcooling temperature of 0.2 K. The heat transfer performance of the dropwise condensation of the L-T-AS has further been elucidated by the theoretical model considering the heat transfer rate through a single droplet associated with the density distribution function of the droplet.

Plasma Actuation for the Turbulent Mixing of Fuel Droplets and Oxidant Air in an Aerospace Combustor

In order to explore plasma-assisted turbulent mixing in aerospace engines, the dielectric barrier discharge plasma actuation for the turbulent mixing of fuel droplets and oxidant air in a ramjet combustor was studied using computational fluid dynamics. A two-way coupling of turbulent air and discrete droplets was realized by Eulerian–Lagrangian simulation, and the dielectric barrier discharge plasma action on flow was modeled by body force. The results show that the plasma actuation can rearrange the recirculation zone behind the evaporative V-groove flameholder, and the main mechanism of actuation is to increase the local momentum of the fluid; the actuation dimension, actuation intensity, and actuation position of the dielectric barrier discharge plasma have strong effects on the turbulent mixing of fuel droplets and oxidant air; and a relatively optimal turbulent mixing can be achieved by adjusting the actuation parameters.

SPRAY IMPINGEMENT FILM ANALYSIS: CHARACTERISTICS EVALUATION AND CORRESPONDING SIMULATION METHOD

Spray impingement is a process during which discrete spray droplets contact a solid surface and form a continuous liquid film. To thoroughly understand such a process is challenging due to the complex gas-liquid-solid interaction and coupling. The liquid converts from a continuous phase into discrete droplets, and finally back to the continuous phase again on the plate. On the basis of experimental analysis, this work investigates a computational fluid dynamics model in the Lagrangian-Eulerian system to focus more on the film dynamics during the impingement. The impingement criterion is modified to incorporate a more accurate momentum transfer within the liquid film. Furthermore, a submodel of droplet slide is coupled into the numerical model to analyze the effect where droplets with a high horizontal velocity will glide away from the plate without contacting it. The numerical model is validated by free spray experiments in the aspects of spray morphology, droplet size, and droplet velocity. Afterward, film dynamics are validated from experimental film thickness measurement with the high-speed laser-induced fluorescence technique. The results show that the modified Bai-Gosman model in the Lagrangian-Eulerian framework can well predict the motion and characteristics of the impingement film, and incorporating the glide model further improves the approximation in many aspects.

Construction of an approach for evaluation of droplet collision models under DDM framework based on an improved Beer's law

In spray simulations, uncertainties in the collision incidence predictions always exist for droplet collision models under the discrete droplet model framework. In this study, Beer's law, which is used for electromagnetic radiation and light absorption, was improved to evaluate the accuracy of droplet collision models by eliminating the shadow effect among the light-absorbing (stationary) parcels. In addition, the range of parcel numbers in the computational domain was carefully calibrated to ensure the appropriate implementation of Beer's law. Based on the improved Beer's law, a theoretical evaluation method of the computational efficiency of droplet collision models was proposed to conduct a qualitative analysis. Using the above methods, a complete approach for evaluating the accuracy and efficiency of droplet collision models was established. Three representative collision models, including the O'Rourke model, the Nordin model, and the no-time-counter (NTC) model, as well as a new hybrid stochastic/trajectory (HST) collision model, were evaluated comprehensively. The convergence and dependence of the predicted collision incidence on various factors, including parcel radius, grid size, reference frame velocity, and droplet number represented by one parcel, were investigated. It was found that the HST model shows relatively better performance on the collision incidence compared with the other collision models in most cases, especially when one parcel represents multiple droplets. The NTC model achieves the highest computational efficiency, and the computational cost of the construction of the adaptive collision cells in the NTC model and the HST model constitutes a major proportion of the total time consumption.

Investigation of flow characteristics on porous gas diffusion layer microstructure that generated with binder and polytetrafluoroethylene distribution

The drainage properties of a gas diffusion layer (GDL) are essential factors in the performance of proton exchange membrane fuel cells. The GDL consisting of a three-dimensional (3D) carbon paper microstructure was developed and meshed with pore-scale reconstruction models in this paper. Localized binder and polytetrafluoroethylene (PTFE) structures were added to the carbon paper microstructure through 3D morphological imaging processing. The monitoring data of 1000 planes were multi-peaky fitted as a function of gas permeability and height to amend the macroscopic porous medium model. We analyzed drainage properties under different contact angles (θ) for the carbon paper with binder and PTFE. We described the mutual intrusion of moisture and air in GDL under different pressure differences. The results show that the pore-scale reconstruction model has the advantages of describing the flow in GDL accurately and with details, detecting low-flow resistance channels that spontaneously formed in GDL, and describing the variation of permeability as a function of location. In a hydrophobic environment, the liquid film connected to a GDL is challenging to split spontaneously. At the same time, the splitting motion of discrete droplets is more prominent than that of the liquid film. The pressure that enables complete water intrusion into the GDL is between 1 and 10 MPa.

Automated study on kinetics and biosensing of glow-type luminescence reaction via digital microfluidics-chemiluminescence.

Automated manipulation of discrete droplets by digital microfluidics (DMF) combined with chemiluminescence (CL) is promising to achieve automated and sensitive biosensing and bioanalysis. Herein, a DMF-CL device was built to automatically study CL kinetics and biosensing of a glow-type CL reaction. Copper-cysteine nanoparticles (Cu/CysNP) were synthesized as a new CL catalyst to extend the CL reaction of luminol-H2O2 to more than 10 min. The automated manipulation of droplets reduced reagent costs and manual errors, leading to real-time, automated, and reliable study of CL kinetics. The CL kinetics curve collected by the DMF-CL integration device is in accordance with that of a commercial CL analyser. The long-lasting luminescence ensured automated, sensitive, and reliable detection of H2O2 as a direct or indirect analyte of the cascade catalytic reaction. Moreover, an innovative asymmetrical splitting method is proposed to quickly and precisely generate daughter droplets to ensure uniformity of the droplets and good repeatability of the DMF-CL measurements. Therefore, the DMF-CL analysis holds great potential for achieving online and automatic analysis and biosensing.

Modeling and LES of high-pressure liquid injection under evaporating and non-evaporating conditions by a real fluid model and surface density approach

Numerical modeling of high-pressure liquid fuel injection remains a challenge in various applications. Indeed, experimental observations have shown that injected liquid fuel jet undergoes a continuous change of state from classical two-phase atomization and spray droplets evaporation to a dense-fluid mixing phenomenon depending on the ambient pressure, temperature, and fuel properties. Accordingly, a predictive and efficient computational fluid dynamics (CFD) model that can represent the possible coexistence of subcritical and supercritical regimes during the fuel injection event is required. The widely used Lagrangian Discrete Droplet Method (DDM) requires parameter tuning of model constants and cannot model the dense near-nozzle region. Meanwhile, the high computational cost of Interface Capturing Methods (ICM) has prohibited their application to industrial cases. Thus, another alternative is an Eulerian Diffuse Interface Model (DIM), where the unresolved interface features are modeled instead of being tracked. Accordingly, the current work proposes a fully compressible multi-component two-phase real-fluid model (RFM) with a diffused interface and closed by a thermodynamic equilibrium tabulation method based on a real-fluid equation of state. The RFM model is complemented with a postulated surface density equation for fuel atomization modeling within the Large Eddy Simulation (LES) framework. The Engine Combustion Network (ECN) Spray A injector non-evaporating and nominal evaporating conditions are used as a reference for the proposed model validation. Simulations are performed using the proposed RFM model that has been implemented in the CONVERGE CFD solver. Under the non-evaporating condition, the RFM model can capture well the fuel mass distribution in the near-nozzle field, but also the interfacial surface area. Besides, the predicted drop size from simulations falls within the experimental data range. On the other hand, under the evaporating condition, spray liquid and vapor penetrations and fuel mixture fraction distribution are also accurately predicted. The vaporization effect on the surface area density is revealed to enhance surface generation in the dense spray region while reducing the surface density in the dilute spray region. The mean droplet size is also relatively reduced under the evaporating condition in the diluted spray region. Overall, the accuracy and computationally efficiency of the proposed RFM model coupled with the surface density equation for high-pressure fuel injection modeling are confirmed, allowing its use for high pressure industrial configurations in future studies.

Multiscale multiphase flow simulations using interface capturing and Lagrangian particle tracking

Numerical simulations of multiphase flows with both interfaces and discrete particles are challenging because they possess a wide range of length and time scales. Meanwhile, the volume of fluid (VOF) method is suitable for resolving the interface, while the discrete particle model (DPM) under the Lagrangian frame better simulates unresolvable particles; a multiscale VOF–DPM combined model is urgently needed for multiscale multiphase flows. The present work implements a VOF–DPM solver that includes a two-way transition algorithm to model the transformation between discrete and continuous phases for bubbles or droplets using OpenFOAM. The interface-capturing scheme in the solver is based on the interIsoFoam solver, which supports the geometric reconstruction of the interface and adaptive mesh refinement. A connected component labeling approach is used for particle detection and VOF-to-DPM transition for discrete bubbles or droplets produced by interface breakup. Conversely, a DPM-to-VOF transition algorithm for particles touching the interface is incorporated to achieve a two-way transition. In addition, phase change modeling between continuous phases and bubble dynamic modeling for cavitating flow cases are also implemented in the solver. Test simulations are performed for validation, including the gas–liquid two-phase dam break and cavitating flow in a convergent–divergent test section. The results demonstrate that the solver is reasonably accurate and can adequately represent the complex phase structure, including the interface and discrete particles.