Increase In Density Ratio Research Articles

Lens Evaporation on Immiscible Liquid Surface with an Interfacial Cooling Effect.

A theoretical heat and mass transfer model of volatile liquid lens evaporation on the surface of an immiscible liquid substrate is established in toroidal coordinates. According to the coupled boundary conditions of heat and mass transfer at the lens surface as well as the interfacial cooling effect, the analytical solutions of the temperature field inside the lens and the vapor concentration field around the lens are derived for the first time. Compared with the isothermal model, the change of contact radius calculated by the present model agrees well with the experimental data, especially when the liquid substrate reaches a relatively high temperature. It also reveals that the temperature distribution inside the lens is not uniform, which is similar to the sessile droplet evaporation on a solid substrate surface. In addition, the excess temperature, heat flux, and evaporation flux of the lens–air interface increase monotonically from the lens center to the contact line. Finally, the influences of density ratio and evaporative cooling number E0 on lens mass evaporation rate are analyzed, which shows that the lens mass evaporation rate decreases with increasing density ratio and evaporative cooling number.

Numerical simulation of local and global mixing/segregation characteristics in a gas–solid fluidized bed

Researches on solids mixing and segregation are of great significance for the operation and design of fluidized bed reactors. In this paper, the local and global mixing and segregation characteristics of binary mixtures were investigated in a gas–solid fluidized bed by computational fluid dynamics-discrete element method (CFD-DEM) coupled approach. A methodology based on solids mixing entropy was developed to quantitatively calculate the mixing degree and time of the bed. The mixing curves of global mixing entropy were acquired, and the distribution maps of local mixing entropy and mixing time were also obtained. By comparing different operating conditions, the effects of superficial gas velocity, particle density ratio and size ratio on mixing/segregation behavior were discussed. Results showed that for the partial mixing state, the fluidized bed can be divided into three parts along the bed height: complete segregation area, transition area and stable mixing area. These areas showed different mixing/segregation processes. Increasing gas velocity promoted the local and global mixing of binary mixtures. The increase in particle density ratio and size ratio enlarged the complete segregation area, reduced the mixing degree and increased the mixing time in the stable mixing area.

Influences of groove configuration and density ratio on grooved leading-edge showerhead film cooling using the pressure sensitive paint measurement technique

In response to the cooling requirements of high adiabatic cooling effectiveness and low convective heat transfer coefficient on the outer wall of the showerhead leading-edge of the turbine vane, a groove configuration is set near-stagnation region, and a grooved leading-edge showerhead cooling structure is designed. Using the steady-state pressure-sensitive paint technology as the measurement method, the interference effect of the groove configuration on the leading-edge showerhead film cooling was experimentally investigated. The influences of the depth and width of the groove on the adiabatic cooling effectiveness of the showerhead leading-edge were compared and analyzed in detail. The present study investigated the adiabatic cooling effectiveness of the showerhead leading-edge with counter-inclined film-holes under varied momentum flux ratio MR varies in the wide range of 0.09–6.04. The influence of coolant mass flow rate on the showerhead cooling effectiveness for the original leading-edge with counter-inclined holes was studied thoroughly; that is, the area-averaged adiabatic cooling effectiveness increases first, then decreases, and then increases again as MR increases. The arrangement of the grooves also enables the grooved structure to ensure the trend mentioned above. Compared with the original case, the effects of the groove on the showerhead cooling effectiveness are mainly reflected in the coolant retention within the groove, the jet-out and superimposition of the downstream of the outer rows. As the MR continues to increase, the grooved case is more likely to have a counter-inclined film distribution in advance, and the turning point of the adiabatic cooling effectiveness drops first and then rises earlier than the original case at a higher MR. The film-coverage of the shallow-grooved case is slightly better than that of the deep-grooved case. The widening of the groove width is not conducive to improving the film-coverage of the showerhead leading-edge. On the contrary, the interference of flow within the groove on the outflow causes a negative effect under lower coolant injection rates. The influence of the density ratio on the showerhead cooling effectiveness of the grooved leading-edge is similar to the original case. At the constant MR, the cooling effectiveness increases as the density ratio increases, and the core region of film-coverage downstream is broader and more extended.

Transverse kink oscillations of inhomogeneous prominence threads: Numerical analysis and Hα forward modelling

Context. Prominence threads are very long and thin flux tubes that are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Aims. The properties of oscillations in threads are greatly affected by the density distribution along the flux tube. Here, we aim to study how the density inhomogeneities in the longitudinal and radial directions modify the periods and damping times of kink oscillations and how this effect would be reflected in the observations. Methods. We solved the ideal magnetohydrodynamics equations using two different methods: (a) performing 3D numerical simulations and (b) solving a 2D generalised eigenvalue problem. We studied the dependence of the periods, damping times, and amplitudes of transverse kink oscillations on the ratio between the densities at the centre and at the ends of the tube, and also on the average density. We applied forward modelling to our 3D simulations to compute synthetic Hα profiles. Results. We confirm that the ratio of the period of the fundamental oscillation mode to the period of the first harmonic increases as the ratio of the central density to the footpoint density is increased, or as the averaged density of the tube is decreased. We find that the damping times due to resonant absorption decrease as the central-to-footpoint density ratio increases. Contrary to the case of longitudinally homogeneous tubes, we find that the damping-time-to-period ratio also increases as the density ratio is increased or the average density is reduced. We present snapshots and time-distance diagrams of the emission in the Hα line. Conclusions. The results presented here have implications for the field of prominence seismology. While the Hα emission can be used to detect the fundamental mode, the first harmonic is barely detectable in Hα. This may explain the lack of detections of the first harmonic. A combination of different spectral lines is required to obtain information about the period ratio and to use it to infer physical properties of the threads.

Oxygen transport in proton exchange membrane fuel cells with metal foam flow fields

Metal foam flow fields have shown great potential in improving the performance of proton exchange membrane (PEM) fuel cells, while their effect on the oxygen transport process remains inadequately understood. In this study, oxygen transport in metal foam flow fields (under zero-humidity operating conditions) is simulated by using a three-dimensional multi-species lattice Boltzmann model. Comparison is done between the metal foam flow field and the conventional channel-rib flow field, and parametric studies are conducted on the metal foam porosity, pore density, and compression ratio. Results show that the metal foam flow field enhances mass transfer of oxygen to the catalyst layer and improves the oxygen distribution homogeneity. Within the range of parameters considered, decrease in the metal foam porosity yields nonmonotonic variation of the mass transfer rate of oxygen to the catalyst layer, which increases at high inlet velocities (higher than 2 m/s) but decreases at low inlet velocities (less than 2 m/s). The increase in metal foam pore density and compression ratio leads to enhanced mass transfer of oxygen, which becomes increasingly prominent at high inlet velocity. The results of this study could be insightful for the implementation of metal foam flow fields in PEM fuel cells. • Oxygen transport in PEM ofuel cells is studied by using lattice Boltzmann method. • Metal foam flow field has stronger oxygen transport than channel-rib flow field. • Increasing metal foam pore density/compression ratio facilitates oxygen transport. • Reducing porosity suppresses/enhances oxygen transport at low/high inlet velocity.

Hydrodynamic behavior of a circular floating solar pond with an entrapped two-layer fluid

The resonant behavior in a moonpool of a floating circular solar pond, with an entrapped two-layer fluid, is studied. The problem is solved by applying a domain-decomposition method using eigenfunction matching. The surface- and internal-wave elevations and the hydrodynamic coefficients of a typical floating solar pond under forced heave or surge motion are computed. The effects of density stratification on surface-wave elevation, added mass, and damping coefficients are analyzed. A collection of resonance frequencies of surface and internal waves is examined, together with the corresponding variations of modal shapes. For heave resonance, the surface and internal waves are characterized by axisymmetric sloshing modes, as opposed to antisymmetric sloshing modes under surge resonances. A frozen-mode approximation method that treats the moonpool fluid as a density-stratified solid is developed to estimate piston-mode frequencies. Non-dimensional resonance frequencies corresponding to antisymmetric and axisymmetric sloshing modes are estimated based on the standing-wave approximation and reciprocity relations between surface and internal wavenumbers. Satisfactory agreement between the estimated resonance frequencies and those computed by eigenfunction matching method is achieved. It is found that the first resonance of the internal wave, rather than higher-order resonances, is more likely to affect the surface-wave behavior, whereas resonances of the surface-wave modes have significant effects on the internal waves. Parametric analyses are performed to study the effects of geometry configurations of the pond. It is found that the resonance frequencies of internal waves under forced heave or surge motion decrease with an increasing density ratio.

Study of shrinkage effect of aluminium based binary alloys as phase change materials for latent heat thermal energy storage applications

In this paper, the effect of shrinkage voids during the solidification of binary alloy as phase change material (PCM) on the thermal performance of a latent heat thermal energy storage is studied. Al-Cu 4.1% wt., Al-Si 7% wt., and Al-Mg 30% wt. are the three binary alloys considered as the PCMs. A generalized numerical model considering differences in density, thermal conductivity, and specific heat between liquid and solid phases is developed. A novel method of fluid fraction update scheme is employed to account for the change in volume fraction due to concentration variation induced due to shrinkage and buoyancy flows. The proposed numerical model is validated using reported numerical and experimental results. The density ratio is found to have a pronounced effect on the solidification of binary alloys. With an increasing density ratio, a faster solidification rate along with a thicker mushy region is observed. This effect of density ratio is found to be much more prominent for higher cold wall temperatures. The performance of three binary alloys, namely, Al-Cu 4.1% wt., Al-Si 7% wt., and Al-Mg 30% wt. is examined. Al-Cu 4.1% wt. alloy is found to solidify much earlier with the release of higher energy, followed by Al-Si and Al-Mg, respectively. The model is further used to predict the size and shape of shrinkage void produced for the binary alloys.

Performance of bioethanol and diesel fuel by thermodynamic simulation of the miller cycle in the diesel engine

The rapid growth in automotive industry and increased fuel consumption and environmental pollution, along with decrease in oil and fossil fuel resources, necessitate the adoption of new technologies to improve engine thermal efficiency and the use of alternative fuels. So, the present study proposes the thermodynamic cycle model of engine based on the Miller cycle and the use of bioethanol and diesel fuels. To analyze the cycle performance after examining the thermodynamic processes of the Miller cycle in the pressure diagram in terms of volume and calculating the amount of the operating fluid (fuel) and its thermodynamic properties based on bioethanol and diesel fuels and a homogeneous mixture of air, temperature was determined at the key points of engine operation through thermodynamic relationships. Then, the changes of engine performance parameters including output power and thermal efficiency were compared in terms of compression ratio based on Miller cycle using bioethanol and diesel fuels and they were also compared to the Otto cycle. The results showed that by increasing density ratio, the output power and thermal efficiency first increase and then decrease. Assuming to have the same volumetric efficiency and use bioethanol fuel, the engine output power and thermal efficiency would increase. Under optimum conditions of engine performance, power at the maximum point of efficiency and efficiency at the maximum point of power in the Miller cycle were greater than those of the Otto cycle. The results can be generalized to Atkinson cycle and other alternative fuels and used in a more practical way in the design and performance evaluation of internal combustion engines.

Effects of liquid and air swirls on the three-dimensional instability of an annular viscous liquid sheet

A three-dimensional temporal linear instability model is proposed for an annular viscous liquid sheet moving between two inviscid gas streams. The model is numerically solved by Chebyshev collocation method. Effects of liquid swirl, inner and outer gas swirls, liquid viscosity, radius ratio, and density ratio on the liquid sheet instability are systematically investigated. The energy budget is calculated to provide insight into the instability driving mechanism. Results show that without gas swirl the liquid swirl has a dual effect on the sheet instability mainly through the axisymmetric mode. For a non-swirling liquid sheet, inner gas swirl induces the maximum wave growth rates of non-axisymmetric modes about 2–3 times higher than those induced by outer gas swirl, suggesting that inner gas swirl is more effective than outer one in promoting the sheet instability. For a swirling liquid sheet, simultaneous inner and outer gas swirls have the greatest destabilizing effect; outer gas swirl shows a slightly greater destabilizing effect than the inner one. Regarding the relative rotating direction, the combination of both co-rotating inner and outer gases is more favored to producing finer droplets. Liquid viscosity is found to have a evident stabilizing effect when Reynolds number is smaller than 500, it dampens the disturbance by dissipating the disturbance energy. Increasing radius ratio decreases the maximum wave growth rates of modes with n=0 and 5⩽n⩽7 but increases those of modes with 1⩽n⩽4. Maximum wave growth rate and optimal wave number are found to first decrease and then increase with increasing density ratio.

DEM simulation of binary mixing particles with different density in an intensive mixer

DEM simulations of binary mixing of particles with different densities were conducted to study the influence of density ratio, impeller speed, and vessel speed on the particle dynamics and mixing performance in an intensive mixer. Six particles at different locations were tagged to discuss the influence of three factors on the particle trajectory and residence time in the impeller influence area. Velocity profile at a fixed point was analyzed to explore the local phenomena. Mean velocity fluctuation, coordination number and relative standard deviation (RSD) were used to quantify the mixing performance in the mixing process. It was found that the residence time increases with the increase of particle density, and the mean velocity fluctuation increases with the increase of impeller speed. The trajectory radius in impeller influence area and the vertical movement increase with the increase of impeller speed. Results obtained indicate that the mixing performance of binary particles with different density increases with the increase of impeller speed and vessel speed, while it decreases with the increase of density ratio.

Development of ecosystem‐based flood mitigation approach – investigations by experiments and numerical simulation

AbstractThe purpose of this research was to understand how to minimise the destructive effects of floods by establishing a bio‐shield of natural vegetation. Experiments were conducted in a water flume to determine the flow properties under different vegetation conditions whereas; volume of fluid (VOF) multiphase approach with a standard k‐ɛ turbulence model was developed for a simulation. The study involved comparative analysis of flow structures and loss of energy in a single patch of finite length, layered (single or double layer with gaps) and continuous (full width) vegetation models. The results showed that, for a single patch, the increase in vegetation density and aspect ratio increased the backwater rise in front of vegetation, reduced velocity and enhanced energy dissipation. Furthermore, in the case of horizontally layered vegetation models (single and double horizontal layers), the loss of energy was higher in double horizontally layered vegetation layouts. Gaps between vegetation patches of a single horizontally layered vegetation model possessed high flow velocity, which was countered by providing a second horizontal layer of vegetation with an alternate gap. The findings from this research are critical in designing optimum vegetation layouts for flood mitigation.

Visualization study on the condensation in a propylene loop heat pipe operating at condenser temperatures between 153 and 283 K

A cryogenic loop heat pipe (CLHP) is an effective means of temperature control at low temperatures, but the physical properties of cryogenic working fluids are quite different from those of room-temperature working fluids, and their phase change and heat transfer characteristics need to be further studied. In this paper, a test system for evaluating the loop heat pipe (LHP) heat transfer performance and visualizing the condensation phenomenon in the condenser was fabricated. The performance of a propylene LHP with 20 W of heat exchange was tested as the condenser temperature was varied from 283 to 153 K, and the flow patterns and heat transfer characteristics during condensation were obtained. With the decreasing condenser temperature, slug flow disappeared, and pseudo-slug flow occurred and transitioned into wavy-annular flow, which was mainly caused by the increase of density ratio and viscosity ratio of the liquid to vapor. These change of flow pattern further affected the two-phase heat transfer and the performance of the CLHP. Pseudo-slug flow was observed at a condenser temperature of 273 to 203 K and it had a notable influence on the condensation pressure drop. The heat transfer coefficient of the two-phase region ranged from 2500 to 5500 W/m2K. When the condenser temperature ranged from 203 to 173 K, the transition of flow pattern and the fluctuation in the system caused a large deviation between the obtained and predicted heat transfer coefficients, and at other condenser temperatures, the obtained values were consistent with the predicted values for the corresponding flow patterns.

Lattice Inclusion of Copper Ions in Ammonium Perchlorate through Co-crystallization: Impact on Lattice, Physical, and Thermal Characteristics

Ammonium perchlorate (AP) is an important oxidizer extensively used in composite solid propellants. Any alterations in the lattice configuration of AP could bring in a dramatic change in its physical, thermal, and ballistic characteristics though the basic thermodynamic properties could remain unaltered. In this work, we attempt to dope AP with copper nitrate through co-crystallization and examine its impact on the lattice, physical, and thermal characteristics. The effect of copper ion on the crystal morphology, bulk density, friability, moisture content, and the decomposition behavior is compared with normal propellant-grade AP. The incorporation of copper ion into AP resulted in an increase in bulk density and aspect ratio and a marginal decrease in the average particle size. The shape factor remained intact. The presence of copper ion remarkably decreased both the low- and high-temperature decomposition and reduced the activation energy for both stages, confirming the catalytic nature of copper-doped AP.

Evolution and modulational instability of interfacial waves in a two-layer fluid with arbitrary layer depths

A nonlinear Schrödinger equation (NLSE) describing the evolution of interfacial waves in a gravitationally stable, inviscid, incompressible, and irrotational two-layer fluid with arbitrary constant layer depths is derived using the multiple scale analysis method. The modulational instability (MI) of the interfacial waves is then analyzed using this NLSE. It is shown that the unstable region shrinks as the density ratio of the two layers increases and as each layer gets thinner. A requirement for unstable waves is that both the upper and lower layers are thicker than the critical depths for those layers. The critical depth of each layer as a function of the density ratio of two layers is obtained by curve fitting and used as a criterion for MI. Moreover, nine cases with various upper- and lower-layer depths are investigated. The relationships of the dark soliton to modulational stability and the bright soliton to MI are discussed in each case. In the unstable regions of the nine cases, it is found that the steepness of the perturbed interface amplitude increases, and the perturbed interface elevation decays more rapidly as the depth of each layer increases. Both the height and the steepness of the perturbed interface elevation increase with increasing density ratio of the two layers.

Nonlinear interaction phenomenon of dust acoustic solitary and shock waves in dusty plasma

AbstractThe effects of head‐on collision on dust acoustic (DA) solitary and shock waves in dusty plasma are investigated considering positively charged inertial dust, Boltzmann distributed negatively charged heavy ions, positively charged light ions, and superthermal electrons in the plasma system. The nonlinear Korteweg‐de‐Vries (KdV) Burger equations are derived taking the extended Poincaré‐Lighthill‐Kuo method into account to study the characteristic properties of nonlinearity and production of solitary shock due to collisions. The study reveals that the amplitudes and widths of the DA shock waves are decreasing with increasing viscosity, electron to dust density ratio, and dust to ion temperature ratio, while they are increasing due to the presence of superthermal electrons. The nonlinearity of DA waves are enhanced with increasing density ratio of electron to dust and temperature ratio of dust to ion and electron, respectively, but it is reducing with superthermal electrons. The phase shifts of DA solitary waves are found to decrease with rising superthermality of electrons and increase with the density ratio of electron to dust.

Oblique collisional effects of dust acoustic waves in unmagnetized dusty plasma

Effects of oblique collisions of the dust acoustic (DA) waves in dusty plasma are studied by considering unmagnetized fully ionized plasma. The plasma consists of inertial warm negatively charged massive dusts, positively charged dusts, superthermal kappa distributed electrons, and isothermal ions. The extended Poincaré–Lighthill–Kuo (ePLK) method is employed for the drivation of two-sided Korteweg–de Vries (KdV) equations (KdVEs). The KdV soliton solutions are derived by using the hyperbolic secant method. The effects of superthermality index of electrons, temperature ratio of isothermal ion to electron, and the density ratio of isothermal ions to negatively charged massive dusts on nonlinear coefficients are investigated. The effects of oblique collision on amplitude, phase shift, and potential profile of right traveling solitons of DA waves are also studied. The study reveals that the new nonlinear wave structures are produced in the colliding region due to head-on collision of the two counter propagating DA waves. The nonlinearity is found to decrease with the increasing density ratio of ion to negative dust in the critical region. The phase shifts decrease (increase) with increasing the temperature ratio of ion to electron (κe). The hump (compressive, κe < κec) and dipshaped (rarefactive, κe > κec) solitons are produced depending on the angle (θ) of oblique collision between the two waves.

Effect of Crystal Chemistry Properties on the Distribution Characteristics of H2O and Na+ in Na-Montmorillonite Interlayer Space: Molecular Dynamics Simulation Study

At monolayer hydration state, the spatial distribution of H2O and Na+ in the interlayer of Na-montmorillonite (Na-MT) with different crystal chemistry properties was investigated by the molecular dynamics simulation method. The simulation results show that when layer charge density increases, H2O will move and form hydrogen bonds with O in tetrahedral surfaces (Ot) at a distance of 1.676 ± 0.043 Å. The impact of isomorphic substitution on the relative concentration of H2O depends largely on the layer charge density of Na-MT, when layer charge density is high, H2O move obviously to both sides of Na-MT sheets with the increase of octahedral substitution ratio. Nevertheless, Na+ coordinate with Ot at a distance of 2.38 Å, and the effect of isomorphic substitution ratio on the diffusion of Na+ is opposite to that of H2O. The mobility of both H2O and Na+ decreases with the increase of layer charge density or tetrahedral substitution ratio. The radial distribution function of Na-Ow (O in H2O) shows that the coordination strength between Na+ and Ow decreases with the increase of layer charge density or tetrahedral substitution ratio, and Na+ are hydrated by four H2O at a Van der Waals radius of 2.386 ± 0.004 Å. The research results can provide a theoretical basis for the efficient application of Na-MT at the molecular and atomic levels.

Interface-resolved numerical simulations of particle-laden turbulent channel flows with spanwise rotation

Interface-resolved simulations of particle-laden turbulent channel flows with spanwise rotation at a Reynolds number of 180 and different rotation numbers ranging from 0.1 to 1.0 are performed with a fictitious domain method. The difficulty of the centrifugal force on the particles not satisfying the periodic boundary condition is circumvented by the feature of the fictitious domain formulation for the neutrally buoyant case, where the centrifugal force in the particle motion equation vanishes, and by only considering a low rotation number of 0.1 and setting the rotation center to be far away from the channel for the non-unity density ratio case. Our results show that the heavy particles (i.e., the particle density being larger than the fluid density) migrate towards the pressure wall, whereas the light particles migrate towards the suction wall. For the density ratio being unity, the particle concentration is higher near the pressure wall than near the suction wall, and we attribute the reason to the effects of the mean secondary flow structure (i.e., the Taylor–Görtler vortices), since similar particle concentration distribution and secondary flow structure are observed in a rotating laminar channel flow. The mean velocities of heavy particles are smaller in the pressure-side half channel except the near-wall region, and larger in the suction-side half channel, compared to the fluid mean velocity; the opposite occurs for the light particle case. The addition of the finite-size particles increases the flow drag. The flow drag is not sensitive to the density ratio for the light particles and increases with increasing density ratio for the heavy particles. The effects of the particles on the fluid root-mean-square velocities of the rotating turbulent channel flow are generally similar to the non-rotating channel case, but become more complicated because of the asymmetric turbulence intensity and particle concentration distribution near two walls caused by the channel rotation.

Numerical study of the particles segregation phenomenon in the fluidized beds

Purpose The purpose of this paper is to present a computational fluid dynamic simulation for the investigation of the particles segregation phenomenon in the gas–solid fluidized beds. Design/methodology/approach These particles have the same size and different densities. The k–ε model and multiphase particle-in-cell method have been utilized for modeling the turbulent fluid flow and solid particles behaviors, respectively. The coupled equations of the velocity and pressure have been solved by using a combination of SIMPLE and PISO algorithms. After validating the simulation, different mixing indices, with different calculation bases, have been investigated, and it has been found that the Lacey mixing index, which was defined based on statistical concepts, is suitable to investigate the segregation/mixing phenomena of this bed in different conditions. Finally, the effects of parameters such as velocity, particle density ratio, jetsam concentration, and initial arrangement on the segregation/mixing behaviors of the bed have been studied. Findings The results show that the increase in the superficial gas velocity decreases the mixing index to a minimum value and then increases this index in the beds with mixed initial condition, unlike the beds with separated initial condition. Moreover, an increase in the particle density ratio increases the minimum fluidization velocity of the bed, and also the amount of segregation, and increase in the jetsam concentration increases the value of the mixing index. Originality/value A computational fluid dynamics simulation has been presented for the particles segregation phenomenon in the gas–solid fluidized beds.

Simulation on mass transfer at immiscible liquid interface entrained by single bubble using particle method

As a Lagrangian particle method, Moving Particle Semi-implicit (MPS) method has great capability to capture interface/surface. In recent years, the multiphase flow simulation using MPS method has become one of the important directions of its developments. In this study, some key methods for multiphase flow have been introduced. The interface tension model in multiphase flow is modified to maintain the smooth of the interface and suitable for the three-phase flow. The mass transfer at immiscible liquid interface entrained by single bubble which could occur in Molten Core-Concrete Interaction (MCCI) has been investigated using this particle method. With the increase of bubble size, the height of entrainment column also increases, but the time of film rupture is slightly different. With the increase of density ratio between the two liquids, the height of entrained column decreases significantly due to the decreasing buoyancy of the denser liquid in the lighter liquid. In addition, the larger the interface tension coefficient is, the more rapidly the entrained denser liquid falls. This study validates that the MPS method has shown great performance for multiphase flow simulation. Besides, the influence of physical parameters on the mass transfer at immiscible interface has also been investigated in this study.