Effect of forced lamina flow on microsegregation simulated by phase field method quantitatively

The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.

RESEARCH ON METHOD TO CALCULATE VELOCITIES OF SOLID PHASE AND LIQUID PHASE IN DEBRIS FLOW

Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two-phase liquid composed of solid phase with the same diameter particles and liquid phase with the same mechanical features. Assume debris flow was one-dimension two-phase liquid moving to one direction, then general equations of velocities of solid phase and liquid phase were founded in two-phase theory. Methods to calculate average pressures, volume forces and surface forces of debris flow control volume were established. Specially, surface forces were ascertained using Bingham＇s rheology equation of liquid phase and Bagnold＇s testing results about interaction between particles of solid phase. Proportional coefficient of velocities between liquid phase and solid phase was put forward, meanwhile, divergent coefficient between theoretical velocity and real velocity of solid phase was provided too. To state succinctly before, method to calculate velocities of solid phase and liquid phase was obtained through solution to general equations. The method is suitable for both viscous debris flow and thin debris flow. Additionally, velocities every phase can be identified through analyzing deposits in-situ after occurring of debris flow. It is obvious from engineering case the result in the method is consistent to that in real-time field observation.

Simulation of faceted dendrite growth of non-isothermal alloy in forced flow by phase field method

Numerical simulation based on a new regularized phase field model was presented to simulate the dendritic shape of a non-isothermal alloy with strong anisotropy in a forced flow. The simulation results show that a crystal nucleus grows into a symmetric dendrite in a free flow and into an asymmetry dendrite in a forced flow. As the forced flow velocity is increased, both of the promoting effect on the upstream arm and the inhibiting effects on the downstream and perpendicular arms are intensified, and the perpendicular arm tilts to the upstream direction. With increasing the anisotropy value to 0.14, all of the dendrite arms tip velocities are gradually stabilized and finally reach their relative saturation values. In addition, the effects of an undercooling parameter and a forced compound flow on the faceted dendrite growth were also investigated.

A method and device for online magnetic resonance multiphase flow detection

Most multiphase flow separation detection methods used commonly in oilfields are low in efficiency and accuracy,and have data delay.An online multiphase flow detection method is proposed based on magnetic resonance technology,and its supporting device has been made and tested in lab and field.The detection technology works in two parts:measure phase holdup in static state and measure flow rate in flowing state.Oil-water ratio is first measured and then gas holdup.The device is composed of a segmented magnet structure and a dual antenna structure for measuring flowing fluid.A highly compact magnetic resonance spectrometer system and intelligent software are developed.Lab experiments and field application show that the online detection system has the following merits:it can measure flow rate and phase holdup only based on magnetic resonance technology;it can detect in-place transient fluid production at high frequency and thus monitor transient fluid production in real time;it can detect oil,gas and water in a full range at high precision,the detection isn’t affected by salinity and emulsification.It is a green,safe and energy-saving system.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

Effect of fluid flow on solidified equiaxed dendrite morphology evolution based on phase field-lattice Boltzmann method

Fluid flow can significantly change the evolution of microstructural morphology. However, relatively little is known how the fluid flow, concentration and microstructure affect each other quantitatively, which is essential to optimize processing parameters. A quantitative simulation study of Al-Cu solidified equiaxed dendrite evolution under forced flow based on phase field-lattice Boltzmann method(PF-LBM) is carried out. Results obtained are validated by Gibbs-Thomson relation at the dendrite tip. Compared with the equiaxed dendrite evolution without flow, the upstream dendrite arm is enhanced while the downstream arm is inhibited. Besides, as the inlet flow rate increases, the secondary arms attached onto the upstream primary arm and the upstream side of the primary arm normal to the inflow velocity has been well developed. Results show that sidewise instabilities of the primary dendrite arm and onset of secondary arm is caused by the local concentration perturbation and will be enhanced or inhabited by the flow. It is believed that the coupled PF-LBM method is able to handle dendrite evolution under forced flow quantitatively, which helps in investigating the solidified dendrite morphology evolution.

Numerical Simulation of Unsteady Friction in Transient Two-Phase Flow with Godunov Method

Most numerical transient flow models that consider dynamic friction employ a finite differences approach or the method of characteristics. These models assume a single fluid (water only) with constant density and pressure wave velocity. But when transient flow modeling attempts to integrate the presence of air, which produces a variable density and pressure-wave velocity, the resolution scheme becomes increasingly complex. Techniques such as finite volumes are often used to improve the quality of results because of their conservative form. This paper focuses on a resolution technique for unsteady friction using the Godunov approach in a finite volume method employing single-equivalent twophase flow equations. The unsteady friction component is determined by taking into account local and convective instantaneous accelerations and the sign of both convective acceleration and velocity values. The approach is used to reproduce a set of transient flow experiments reported in the literature, and good agreement between simulated and experimental results is found.

Velocity Projection with Upwind Scheme Based on the Discontinuous Galerkin Methods for the Two Phase Flow Problem

The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

Gradient-augmented hybrid interface capturing method for incompressible two-phase flow

Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives（normals and curvatures） and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.

THE SOLUTION FOR TRANSIENT TWO-PHASE FLOW BY SPLIT FLUX VECTOR METHOD

In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.

Three-dimensional multi-relaxation-time lattice Boltzmann front-tracking method for two-phase flow

We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.

Multi-relaxation-time lattice Boltzmann front tracking method for two-phase flow with surface tension

In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.

Liquid-Gas Two-Phase Flow Simulation in Mold Filling Process with Level Set Method

Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.

VARIABLE GRID FINITE DIFFERENCE METHOD FOR TWO-DIMENSIONAL TWO-PHASEIMMISCIBLE FLOW

An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.

Review of gas-solid two phase flow rate-concentration detection technology

The indirect detection method basic principle of rate and concentration,application range and research results on gassolid two phase flow were discussed.The present development situation and the existing problems of rate and concentration detection technology were analyzed and summarized.Emphatically analyzed the existing problems in the industrial application and research status of electrostatic method in measuring phase concentration.Design criterion of electrostatic phase concentration sensor is given,the superiority and wide industrial application prospect of the sensor used for phase concentration measurement are clarified.

Phase transition model of water flow irradiated by high-energy laser in a chamber

In the absorption chamber of a high-energy laser energy meter, water is directly used as an absorbing medium and the interaction of the high-power laser and the water flow can produce a variety of physical phenomena such as phase transitions. The unit difference method is adopted to deduce the phase transition model for water flow irradiated by a high-energy laser. In addition, the model is simulated and verified through experiments. Among them, the experimental verification uses the photographic method, shooting the distribution and the form of the air mass of water flow in different operating conditions, which are compared with the simulation results. The research shows that it is achievable to reduce the intensity of the phase transition by increasing the water flow, reducing the power intensity of the beam, shortening the distance the beam covers, reducing the initial water temperature or adopting a shorter wavelength laser. The study＇s results will provide the reference for the design of a water-direct-absorption-type high-energy laser energy meter as well as an analysis of the interaction processes of other similar high-power lasers and water flow.

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.

An Improved Optical Flow Method for Image Registration with Large-scale Movements

在这份报纸，为图象登记的一个改进光流动方法被建议。它以它与一个起始的运动评估者改进光流动方法的方法是新奇的：扩大了阶段关联技术(EPCT ) ，用补偿前者的缺乏的后者的优点。以一种更多的详细方式，光流动方法能到达亚象素精确性并且计算象啁啾并且倾斜一样的复杂失真模式，这能被说但是与大规模运动弱。因为 EPCT 用象素水平精确性盖住大翻译和旋转的大小并且在精明的负担是有效的，它能为光流动方法被当作一个好起始的运动评估者。测试证明了特别，这个改进方法将显著地与大规模运动为图象提高登记表演并且对随机的噪音柔韧。

TWO-PHASE FLOW FOR A HORIZONTAL WELL PENETRATING A NATURALLY FRACTURED RESERVOIR WITH EDGE WATER INJECTION

This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamlinemodel,The mathematical model of the verical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established ,and whose accuratesolutions are obtained by using the characteristic method .The saturation distributionsin the fractured system and the matrix system as well as the formula of the time of water free production are presented .All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fracturedreservirs.