Two-phase flow interface evolution simulation method and system

Abstract

The embodiment of the invention provides a two-phase flow interface evolution simulation method and system. The method comprises the following steps: building a conservation type phase field Allen-Cahn equation about the interface evolution of a two-phase flow under a periodic boundary condition, introducing a stable term, and obtaining a semi-discrete system of a two-phase flow phase field model after spatial discretization; performing time discretization by adopting an explicit integral factor method to obtain a space-time full discretization format; performing reasonable approximation on an exponential function in an explicit integral factor method to obtain an explicit single-step space-time full discrete format for unconditionally maintaining the extreme value principle and mass conservation of the two-phase flow phase field model; and obtaining an approximate value of the interface on the next time layer, and sequentially iterating to obtain a simulation result of the phase field interface of the two-phase flow at the termination moment in the calculation region. Based on improvement processing of a time discretization method, especially selection of stabilization parameters and approximation of an exponential function, a time discretization format can achieve a fourth-order unconditional extreme value maintaining principle and mass conservation.

Description

technical field [0001] The invention relates to the field of two-phase flow simulation, in particular to a two-phase flow interface evolution simulation method and system. Background technique [0002] Two-phase flow exists widely in nature and in modern industrial production processes, and is closely related to human production and life. The rapid development of industrial simulation based on computational fluid dynamics requires complex two-phase flow simulation, and puts forward higher and higher requirements for the calculation effect of two-phase flow models and algorithms. For the Allen-Cahn equation of the phase field and its correction equation for maintaining the mass conservation law, most of the current schemes have strict requirements on the time step when maintaining the extremum principle, while the high-order integral factor Runge-Kutta (IFRK) scheme is in When the time step is large, the exponential decay phenomenon will appear again, which makes it difficul...

AI Technical Summary

Problems solved by technology

For the Allen-Cahn equation of the phase field and its correction equation for maintaining the mass conservation law, most of the current schemes have strict requirements on the time step when maintaining the extremum principle, while the high-order integral factor Runge-Kutta (IFRK) scheme is in Exponential decay occurs when the time step is large, which makes it difficult to quickly solve the Allen-Cahn equation and efficiently capture the two-phase interface

Images