This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is com...This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for the water depth is implemented to resolve the stationary or wet/dry fronts. The bed slope source term is discretized using a central difference method to capture the static flow state over the irregular topography. Second-order accuracy in space is achieved by using the slope limited linear reconstruction method. With the proposed method, the model can avoid the partially wetting/drying cell problem and maintain the mass conservation. The proposed model is tested and verified against three theoretical benchmark tests and two experimental dam break flows. Further, the model is applied to predict the maximum water level and the flood arrival time at different gauge points for the Malpasset dam break event. The predictions agree well with the numerical results and the measurement data published in literature, which demonstrates that with the present model, a well-balanced state can be achieved and the water depth can be nonnegative when the Courant number is kept less than 0.25.展开更多
A kinetic flux-vector splitting (KFVS) scheme is applied for solving a reduced six-equation two-phase flow model of Saurel et al. [1]. The model incorporates single velocity, two pressures and relaxation terms. An add...A kinetic flux-vector splitting (KFVS) scheme is applied for solving a reduced six-equation two-phase flow model of Saurel et al. [1]. The model incorporates single velocity, two pressures and relaxation terms. An additional seventh equation, describing the total mixture energy, is added to the model to guarantee the correct treatment of shocks in the single phase limit. Some salient features of the model are that it is hyperbolic with only three wave propagation speeds and the volume fraction remains positive. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge-Kutta time stepping method. Moreover, a pressure relaxation procedure is used to fulfill the interface conditions. For validation, the results of suggested scheme are compared with those from the high resolution central upwind and HLLC schemes. The central upwind scheme is also applied for the first time to this model. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.展开更多
Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-var...Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-variable cell population balance models (PBMs) offer the most general way to describe the complicated phenomena associated with cell growth, substrate consumption and product formation. For that reason, solving and understanding of such models are essential to predict and control cell growth in the processes of biotechnological interest. Such models typically consist of a partial integro-differential equation for describing cell growth and an ordinary integro-differential equation for representing substrate consumption. However, the involved mathematical complexities make their numerical solutions challenging for the given numerical scheme. In this article, the central upwind scheme is applied to solve the single-variate and bivariate cell population balance models considering equal and unequal partitioning of cellular materials. The validity of the developed algorithms is verified through several case studies. It was found that the suggested scheme is more reliable and effective.展开更多
This paper presents a well-balanced two-dimensional (2D) finite volume model to simulate the propagation, runup and rundown of long wave. Non-staggered grid is adopted to discretize the governing equation and the inte...This paper presents a well-balanced two-dimensional (2D) finite volume model to simulate the propagation, runup and rundown of long wave. Non-staggered grid is adopted to discretize the governing equation and the intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for water depth is implemented in the present model to treat the appearance of wet/dry fronts, and the friction term is solved by a semi-implicit scheme to ensure the stability of the model. The Euler method is applied to update flow variable to the new time level. The model is verified against two experimental cases and good agreements are observed between numerical results and observed data.展开更多
A numerical model for simulating the dambreak problems was presented. The model was based on a high-resolution semi-discrete central-upwind difference scheme. In order to reduce spurious oscillation, the uniformly non...A numerical model for simulating the dambreak problems was presented. The model was based on a high-resolution semi-discrete central-upwind difference scheme. In order to reduce spurious oscillation, the uniformly non-oscillatory limiter was employed. A third-order total variation diminishing Runge-Kutta method is used for time integration. The main feature of the presented method is its simplicity. It requires no Riemann solvers, no flux splitting and no flux limiter. It is explicit and does not require dimensional splitting for two dimensions. The Simpson quadrature rule was employed to compute the source term. To verify the effectiveness and accuracy of the proposed method, the 1D dambreak, circular dam-break and partial dam-break problems were simulated. The results are shown to be in good agreement with analytical solution and numerical results obtained by other methods.展开更多
基金Project supported by Natural Science Foundation of Zhejiang Province(Grant No.LR16E090001)the Research Funding of Shenzhen City(Grant No.JCYJ20160425164642646)the Zhejiang Province Science and Technology Research Funding(Grant No.2015C03015)
文摘This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for the water depth is implemented to resolve the stationary or wet/dry fronts. The bed slope source term is discretized using a central difference method to capture the static flow state over the irregular topography. Second-order accuracy in space is achieved by using the slope limited linear reconstruction method. With the proposed method, the model can avoid the partially wetting/drying cell problem and maintain the mass conservation. The proposed model is tested and verified against three theoretical benchmark tests and two experimental dam break flows. Further, the model is applied to predict the maximum water level and the flood arrival time at different gauge points for the Malpasset dam break event. The predictions agree well with the numerical results and the measurement data published in literature, which demonstrates that with the present model, a well-balanced state can be achieved and the water depth can be nonnegative when the Courant number is kept less than 0.25.
文摘A kinetic flux-vector splitting (KFVS) scheme is applied for solving a reduced six-equation two-phase flow model of Saurel et al. [1]. The model incorporates single velocity, two pressures and relaxation terms. An additional seventh equation, describing the total mixture energy, is added to the model to guarantee the correct treatment of shocks in the single phase limit. Some salient features of the model are that it is hyperbolic with only three wave propagation speeds and the volume fraction remains positive. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge-Kutta time stepping method. Moreover, a pressure relaxation procedure is used to fulfill the interface conditions. For validation, the results of suggested scheme are compared with those from the high resolution central upwind and HLLC schemes. The central upwind scheme is also applied for the first time to this model. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.
文摘Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-variable cell population balance models (PBMs) offer the most general way to describe the complicated phenomena associated with cell growth, substrate consumption and product formation. For that reason, solving and understanding of such models are essential to predict and control cell growth in the processes of biotechnological interest. Such models typically consist of a partial integro-differential equation for describing cell growth and an ordinary integro-differential equation for representing substrate consumption. However, the involved mathematical complexities make their numerical solutions challenging for the given numerical scheme. In this article, the central upwind scheme is applied to solve the single-variate and bivariate cell population balance models considering equal and unequal partitioning of cellular materials. The validity of the developed algorithms is verified through several case studies. It was found that the suggested scheme is more reliable and effective.
文摘This paper presents a well-balanced two-dimensional (2D) finite volume model to simulate the propagation, runup and rundown of long wave. Non-staggered grid is adopted to discretize the governing equation and the intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for water depth is implemented in the present model to treat the appearance of wet/dry fronts, and the friction term is solved by a semi-implicit scheme to ensure the stability of the model. The Euler method is applied to update flow variable to the new time level. The model is verified against two experimental cases and good agreements are observed between numerical results and observed data.
文摘A numerical model for simulating the dambreak problems was presented. The model was based on a high-resolution semi-discrete central-upwind difference scheme. In order to reduce spurious oscillation, the uniformly non-oscillatory limiter was employed. A third-order total variation diminishing Runge-Kutta method is used for time integration. The main feature of the presented method is its simplicity. It requires no Riemann solvers, no flux splitting and no flux limiter. It is explicit and does not require dimensional splitting for two dimensions. The Simpson quadrature rule was employed to compute the source term. To verify the effectiveness and accuracy of the proposed method, the 1D dambreak, circular dam-break and partial dam-break problems were simulated. The results are shown to be in good agreement with analytical solution and numerical results obtained by other methods.
