Pore-Scale Investigation of Coupled Two-Phase and Reactive Transport in the Cathode Electrode of Proton Exchange Membrane Fuel Cells

A three-dimensional multicomponent multiphase lattice Boltzmann model(LBM)is established to model the coupled two-phase and reactive transport phenomena in the cathode electrode of proton exchange membrane fuel cells.The gas diff usion layer(GDL)and microporous layer(MPL)are stochastically reconstructed with the inside dynamic distribution of oxygen and liquid water resolved,and the catalyst layer is simplifi ed as a superthin layer to address the electrochemical reaction,which provides a clear description of the fl ooding eff ect on mass transport and performance.Diff erent kinds of electrodes are reconstructed to determine the optimum porosity and structure design of the GDL and MPL by comparing the transport resistance and per-formance under the fl ooding condition.The simulation results show that gradient porosity GDL helps to increase the reactive area and average concentration under fl ooding.The presence of the MPL ensures the oxygen transport space and reaction area because liquid water cannot transport through micropores.Moreover,the MPL helps in the uniform distribution of oxygen for an effi cient in-plane transport capacity.Crack and perforation structures can accelerate the water transport in the assembly.The systematic perforation design yields the best performance under fl ooding by separating the transport of liquid water and oxygen.

A USM-Θ two-phase turbulence model for simulating dense gas-particle flows

A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.

Coupled Model of Two-phase Debris Flow,Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.

REVIEW ON MATHEMATICAL ANALYSIS OF SOME TWO-PHASE FLOW MODELS

The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.

IMPROVED SUBGRID SCALE MODEL FOR DENSE TURBULENT SOLID-LIQUID TWO-PHASE FLOWS

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.

GLOBAL EXISTENCE OF CLASSICAL SOLUTION FOR A VISCOUS LIQUID-GAS TWO-PHASE MODEL WITH MASS-DEPENDENT VISCOSITY AND VACUUM

In this work, we obtain the global existence and uniqueness of classical solu-tions to a viscous liquid-gas two-phase model with mass-dependent viscosity and vacuum in one dimension, where the initial vacuum is allowed. We get the upper and lower bounds of gas and liquid masses n and m by the continuity methods which we use to study the compressible Navier-Stokes equations.

Transient pressure analysis of polymer flooding fractured wells with oil-water two-phase flow

The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.

Two-phase nonlocal integral models with a bi-Helmholtz averaging kernel for nanorods

In this work,the static tensile and free vibration of nanorods are studied via both the strain-driven(Strain D)and stress-driven(Stress D)two-phase nonlocal models with a bi-Helmholtz averaging kernel.Merely adjusting the limits of integration,the integral constitutive equation of the Fredholm type is converted to that of the Volterra type and then solved directly via the Laplace transform technique.The unknown constants can be uniquely determined through the standard boundary conditions and two constrained conditions accompanying the Laplace transform process.In the numerical examples,the bi-Helmholtz kernel-based Strain D(or Stress D)two-phase model shows consistently softening(or stiffening)effects on both the tension and the free vibration of nanorods with different boundary edges.The effects of the two nonlocal parameters of the bi-Helmholtz kernel-based two-phase nonlocal models are studied and compared with those of the Helmholtz kernel-based models.

Comparative numerical study of single and two-phase models of nanofluid heat transfer in wavy channel

The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.

THE STUDY ON THREE-DIMENSIONAL MATHEMATICAL MODEL OF RIVER BED EROSION FOR WATER-SEDIMENT TWO-PHASE FLOW

Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.

LARGE TIME BEHAVIOR OF GLOBAL STRONG SOLUTIONS TO A TWO-PHASE MODEL WITH A MAGNETIC FIELD

In this paper,the Cauchy problem for a two-phase model with a magnetic field in three dimensions is considered.Based on a new linearized system with respect to(c−c_(∞),P−P_(∞),u,H)for constants c_(∞)≥0 and P_(∞)>0,the existence theory of global strong solution is established when the initial data is close to its equilibrium in three dimensions for the small H^(2) initial data.We improve the existence results obtained by Wen and Zhu in[40]where an additional assumption that the initial perturbations are bounded in L^(1)-norm was needed.The energy method combined with the low-frequency and high-frequency decomposition is used to derive the decay of the solution and hence the global existence.As a by-product,the time decay estimates of the solution and its derivatives in the L^(2)-norm are obtained.

CFD studies of scraper built in SPA clarifying tank based on mixture solid-liquid two-phase flow model

In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.

Modelling of a Two-Phase Thermosyphon Loop for Passive Air-Conditioning of a House in Hot and Dry Climate Countries

The two-phase thermosyphon loop is an efficient solution for space cooling. This paper presents the simulation results of numerical studies on the heat transfer and thermal performance of a two-phase thermosiphon loop for passive air-conditioning of a house. The fluid considered in this study is methanol, which is compatible with copper and is environmentally friendly. These numerical results show that the temperature at the evaporator wall drops from 23°C to 13°C and increases at the condenser. The solar flux density has a strong influence on the condenser temperature. The mass flow rates and masses at the evaporator and condenser increase with temperature. The variation of evaporating and condensing temperature affects the performance of the system. For a constant evaporating and condensing temperature of 2°C and 29°C, the COP is 0.77 and 0.84 respectively. With these results, the use of the two-phase thermosyphon loop in air conditioning is possible to obtain a thermal comfort of the occupants acceptable by the standards but with a large exchange surface of the evaporator.

A K-εTWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.

On well-posedness of two-phase nonlocal integral models for higher-order refined shear deformation beams

Due to the conflict between equilibrium and constitutive requirements,Eringen’s strain-driven nonlocal integral model is not applicable to nanostructures of engineering interest.As an alternative,the stress-driven model has been recently developed.In this paper,for higher-order shear deformation beams,the ill-posed issue(i.e.,excessive mandatory boundary conditions(BCs)cannot be met simultaneously)exists not only in strain-driven nonlocal models but also in stress-driven ones.The well-posedness of both the strain-and stress-driven two-phase nonlocal(TPN-Strain D and TPN-Stress D)models is pertinently evidenced by formulating the static bending of curved beams made of functionally graded(FG)materials.The two-phase nonlocal integral constitutive relation is equivalent to a differential law equipped with two restriction conditions.By using the generalized differential quadrature method(GDQM),the coupling governing equations are solved numerically.The results show that the two-phase models can predict consistent scale-effects under different supported and loading conditions.

Prediction of shear bands in sand based on granular flow model and two-phase equilibrium

In contrast to the traditional interpretation of shear bands in sand as a bifurcation problem in continuum mechanics,shear bands in sand are considered as high-strain phase(plastic phase) of sand and the materials outside the bands are still in low-strain phase(elastic phase),namely,the two phases of sand can coexist under certain condition.As a one-dimensional example,the results show that,for materials with strain-softening behavior,the two-phase solution is a stable branch of solutions,but the method to find two-phase solutions is very different from the one for bifurcation analysis.The theory of multi-phase equilibrium and the slow plastic flow model are applied to predict the formation and patterns of shear bands in sand specimens,discontinuity of deformation gradient and stress across interfaces between shear bands and other regions is considered,the continuity of displacements and traction across interfaces is imposed,and the Maxwell relation is satisfied.The governing equations are deduced.The critical stress for the formation of a shear band,both the stresses and strains inside the band and outside the band,and the inclination angle of the band can all be predicted.The predicted results are consistent with experimental measurements.

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall,gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions.In this paper by modification of a friction model available in the literature,an improved semiempirical model is proposed.The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved.Pressure gradient and velocity profiles are validated against experimental data.Using the improved model,the pressure gradient deviation from experiments diminishes by about 3%;the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.

Three-dimensional CFD simulation of inlet structure flow in pumping station based on Eulerian solid- liquid two-phase flow model

Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- liquid two-phase flow model. The numerical results of the preliminary scheme show that sediment deposition occurs in the forebay of pumping station because of poor flow pattern therein. In order to improve hydraulic configuration in the forebay,one modified measure reconstructs water diversion weir shape,and another measure sets a water retaining sill in the approach channel. The simulation results of the modified scheme prove that back flow in the forebay has been eliminated and the sediment deposition region has also been reduced.

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.

Numerical Study of Heat Transfer and Flow Bifurcation of CuO Nanofluid in Sudden Expansion Microchannel Using Two-Phase Model

In this paper, laminar forced convection of CuO nanofluid is numerically investigated in sudden expansion microchannel with isotherm walls and different expansion ratios (ER). An Eulerian two-fluid model is considered to simulate the nanofluid flow inside the microchannel and the governing mass, momentum and energy equations for both phases are solved using the finite volume method. Eulerian-Eulerian two-phase model is very efficient because of considering the relative velocity and temperature of the phases and the nanoparticle concentration distribution. In solving the flow equations for both phases, the SIMPLE algorithm is modified for the coupling of the velocity and pressure and the continuity equations for both phases are combined in order to create the pressure correction equations. However, the Eulerian-Eulerian modeling results show higher heat transfer enhancement in comparison to pure water, so that for a 2% copper-water nanofluid, it has been observed a 35% increase of the heat transfer. The heat transfer enhancement increases with increase in Reynolds number and nanoparticle volume concentration, while the pressure drop increases only slightly. An investigation of the expansion ratio of microchannel shows that the average Nusselt number increases with decrease in expansion ratio as well as with increase in Reynolds number. Also, the Bifurcation has been occurred in higher Reynolds number that is different for each expansion ratio of the microchannel.