The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central ...The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central schemes to solve these equations. The equilibrium states for each component are coupled in space and time to have a common temperature and velocity. The current schemes can handle strong shocks and the oscillations near the interfaces are negligible, which usually happens in the multi-component flows. The schemes also guarantee the exact mass conservation for each component, the exact conservation of total momentum, and energy in the whole particle system. The central schemes are robust, reliable, compact and easy to implement. Several one- and two-dimensional numerical test cases are included in this paper, which validate the application of these schemes to relativistic multi-component flows.展开更多
In this paper,we propose a second-order moving-water equilibria preserving nonstaggered central scheme to solve the Ripa model via flux globalization.To maintain the moving-water steady states,we use the discrete sour...In this paper,we propose a second-order moving-water equilibria preserving nonstaggered central scheme to solve the Ripa model via flux globalization.To maintain the moving-water steady states,we use the discrete source terms proposed by Britton et al.(J Sci Comput,2020,82(2):Art 30)by incorporating the expression of the source terms as a whole into the flux gradient,which directly avoids the discrete complexity of the source terms in order to maintain the well-balanced properties of the scheme.In addition,since the nonstaggered central scheme requires re-projecting the updated values of the nonstaggered cells onto the staggered cells,we modify the calculation of the global variables by constructing ghost cells and alternating the values of the global variables with the water depths obtained from the solution through the nonlinear relationship between the global flux and the water depth.In order to maintain the second-order accuracy of the scheme on the time scale,we incorporate a new Runge-Kutta type time discretization in the evolution of the numerical solution for the nonstaggered cells.Meanwhile,we introduce the"draining"time step technique to ensure that the water depth is positive and that it satisfies mass conservation.Numerical experiments verify that the scheme is well-balanced,positivity-preserving and robust.展开更多
A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization an...A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization and the optimal third-order Strong Stability Preserving (SSP) Runge-Kutta method for time integration. The third-order compact Central Weighted Essentially Non-Oscillatory (CWENO) reconstruction was adopted to guarantee the non-oscillatory behavior of the presented scheme and improve the resolution. Two kinds of source terms were considered in this work. They were evaluated using different approaches. The resulting scheme does not require Riemann solvers or characteristic decomposition, hence it retains all the attractive features of central schemes such as simplicity and high resolution. To evaluate the performance of the presented scheme, several numerical examples were tested. The results demonstrate that our method is efficient, stable and robust.展开更多
A well-balanced second order finite volume central scheme for the magnetohydrodynamic(MHD)equations with gravitational source term is developed in this paper.The scheme is an unstaggered central scheme that evolves th...A well-balanced second order finite volume central scheme for the magnetohydrodynamic(MHD)equations with gravitational source term is developed in this paper.The scheme is an unstaggered central scheme that evolves the numerical solution on a single grid and avoids solving Riemann problems at the cell interfaces using ghost staggered cells.A subtraction technique is used on the conservative variables with the support of a known steady state in order to manifest the well-balanced property of the scheme.The divergence-free constraint of themagnetic field is satisfied after applying the constrained transport method(CTM)for unstaggered central schemes at the end of each time-step by correcting the components of the magnetic field.The robustness of the proposed scheme is verified on a list of numerical test cases from the literature.展开更多
We study high-resolution central schemes in Lagrangian coordinates for the one-dimensional system of conservation laws describing the evolution of two gases in slab geometry separated by an interface.By using Lagrangi...We study high-resolution central schemes in Lagrangian coordinates for the one-dimensional system of conservation laws describing the evolution of two gases in slab geometry separated by an interface.By using Lagrangian coordinates,the interface is transformed to a fixed coordinate in the computational domain and,as a consequence,the movement of the interface is obtained as a byproduct of the numerical solution.The main contribution is the derivation of a special equation of state to be imposed at the interface in order to avoid non-physical oscillations.Suitable boundary conditions at the piston that guarantee second order convergence are described.We compare the solution of the piston problem to other results available in the literature and to a reference solution obtained within the adiabatic approximation.A shock-interface interaction problem is also treated.The results on these tests are in good agreement with those obtained by other methods.展开更多
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.展开更多
A two-dimensional nonoscillatory central difference scheme was extended to the shallow water equations. A high-resolution numerical method for solving the shallow water equations was presented. In order to prevent osc...A two-dimensional nonoscillatory central difference scheme was extended to the shallow water equations. A high-resolution numerical method for solving the shallow water equations was presented. In order to prevent oscillation, the nonlinear limiter is employed to approximate the discrete slopes. The main advantage of the presented method is simplicity comparable with the upwind schemes. This method does not require Riemann solvers or some form of flux difference splitting methods. Furthermore, the discrete derivatives of flux can be approximated by the component-wise approach and thus the computation of Jacobian can be avoided. The method retains high resolution and high accuracy similar to the upwind results. It is applied to simulating several tests, including circular dam-break problem, shock focusing problem and partial dam-break problem. The results are in good agreement with the numerical results obtained by other methods. The simulated results also demonstrate that the presented method is stable and efficient.展开更多
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 elaborated the bidding scheme design of the Central Recreational Area in Linkong Economic Zone from the perspectives of site analysis, determination of properties, design theme and principle, scheme layout ...This paper elaborated the bidding scheme design of the Central Recreational Area in Linkong Economic Zone from the perspectives of site analysis, determination of properties, design theme and principle, scheme layout and conception, in order to explore the potential enlightenments of the scheme design process, and disclose significance of site analysis and its close relationship with the scheme design.展开更多
The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechan...The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics (English Edition), 2007, 28(7), 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.展开更多
This study examines the evolution of fiscal relations between the central state and local governments in China since the founding of the People's Republic The emphasis is placed on the changing arrangement for t...This study examines the evolution of fiscal relations between the central state and local governments in China since the founding of the People's Republic The emphasis is placed on the changing arrangement for the central state's collection of taxes and allocation of capital among provinces The central local relation has experienced constant reorganization since 1949, but the tax division scheme introduced in the early 1990s represents the most significant change Prior to the 1990s, the central local fiscal relation was shaped by a highly centralized administrative and economic system without rational and scientific design for revenue collection and expenditure The financial responsibility system established after the reforms contributed to regional economic development during the first decade of reforms, but it was a temporary and transitional arrangement that does not meet the requirement of rational resources arrangement according to free market forces The implementation of the tax division scheme has stopped the decline of the ratio of fiscal revenue to gross domestic product, raised the percentage of central revenue in total national income, and strengthened the central function for macro economic control The mechanism for tax refund under the tax division scheme required further improvements Several proposals are made in this study: 1) tax classification should be adjusted according to administration; 2) a scientific and standardized system for regional transfer payment should be developed; and 3) the power for tax legislation should be delineated according to the rational division of administration between the central and local government展开更多
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.展开更多
This paper proposes a negotiation-based TDMA scheme for ad hoc networks, which was modeled as an asynchronous myopic repeated game. Compared to the traditional centralized TDMA schemes, our scheme operates in a decent...This paper proposes a negotiation-based TDMA scheme for ad hoc networks, which was modeled as an asynchronous myopic repeated game. Compared to the traditional centralized TDMA schemes, our scheme operates in a decentralized manner and is scalable to topology changes. Simulation results show that, with respect to the coloring quality, the performance of our scheme is close to that of the classical centralized algorithms with much lower complexity.展开更多
Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisi...Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisit the question of non-oscillatory behavior of high-resolution central schemes in terms of the slope limiter proposed by van Albada et al.(Astron Astrophys 108:76–84,1982).The van Albada(vA)limiter is smoother near extrema,and consequently,in many cases,it outperforms the results obtained using the standard minmod limiter.In particular,we prove that the vA limiter ensures the one-dimensional Total-Variation Diminishing(TVD)stability and demonstrate that it yields noticeable improvement in computation of one-and two-dimensional systems.展开更多
文摘The dynamics of inviscid multi-component relativistic fluids may be modeled by the relativistic Euler equations, augmented by one (or more) additional species equation(s). We use the high-resolution staggered central schemes to solve these equations. The equilibrium states for each component are coupled in space and time to have a common temperature and velocity. The current schemes can handle strong shocks and the oscillations near the interfaces are negligible, which usually happens in the multi-component flows. The schemes also guarantee the exact mass conservation for each component, the exact conservation of total momentum, and energy in the whole particle system. The central schemes are robust, reliable, compact and easy to implement. Several one- and two-dimensional numerical test cases are included in this paper, which validate the application of these schemes to relativistic multi-component flows.
基金supported by the National Natural Science Foundation of China(51879194)。
文摘In this paper,we propose a second-order moving-water equilibria preserving nonstaggered central scheme to solve the Ripa model via flux globalization.To maintain the moving-water steady states,we use the discrete source terms proposed by Britton et al.(J Sci Comput,2020,82(2):Art 30)by incorporating the expression of the source terms as a whole into the flux gradient,which directly avoids the discrete complexity of the source terms in order to maintain the well-balanced properties of the scheme.In addition,since the nonstaggered central scheme requires re-projecting the updated values of the nonstaggered cells onto the staggered cells,we modify the calculation of the global variables by constructing ghost cells and alternating the values of the global variables with the water depths obtained from the solution through the nonlinear relationship between the global flux and the water depth.In order to maintain the second-order accuracy of the scheme on the time scale,we incorporate a new Runge-Kutta type time discretization in the evolution of the numerical solution for the nonstaggered cells.Meanwhile,we introduce the"draining"time step technique to ensure that the water depth is positive and that it satisfies mass conservation.Numerical experiments verify that the scheme is well-balanced,positivity-preserving and robust.
基金Project supported by the National Natural Science Foundation of China (Grant No: 60134010).
文摘A numerical two-dimensional shallow water method was based on method for solving the equations was presented. This the third-order genuinely multidimensional semi-discrete central scheme for spatial discretization and the optimal third-order Strong Stability Preserving (SSP) Runge-Kutta method for time integration. The third-order compact Central Weighted Essentially Non-Oscillatory (CWENO) reconstruction was adopted to guarantee the non-oscillatory behavior of the presented scheme and improve the resolution. Two kinds of source terms were considered in this work. They were evaluated using different approaches. The resulting scheme does not require Riemann solvers or characteristic decomposition, hence it retains all the attractive features of central schemes such as simplicity and high resolution. To evaluate the performance of the presented scheme, several numerical examples were tested. The results demonstrate that our method is efficient, stable and robust.
基金the National Council for Scientific Research of Lebanon(CNRS-L)for granting a doctoral fellowship to Farah Kanbarfunding by theQualification Programof the Julius Maximilians University Wurzburg.
文摘A well-balanced second order finite volume central scheme for the magnetohydrodynamic(MHD)equations with gravitational source term is developed in this paper.The scheme is an unstaggered central scheme that evolves the numerical solution on a single grid and avoids solving Riemann problems at the cell interfaces using ghost staggered cells.A subtraction technique is used on the conservative variables with the support of a known steady state in order to manifest the well-balanced property of the scheme.The divergence-free constraint of themagnetic field is satisfied after applying the constrained transport method(CTM)for unstaggered central schemes at the end of each time-step by correcting the components of the magnetic field.The robustness of the proposed scheme is verified on a list of numerical test cases from the literature.
基金supported by grants of the Catania and Messina Universitiespartially by the Italian“MIUR”.
文摘We study high-resolution central schemes in Lagrangian coordinates for the one-dimensional system of conservation laws describing the evolution of two gases in slab geometry separated by an interface.By using Lagrangian coordinates,the interface is transformed to a fixed coordinate in the computational domain and,as a consequence,the movement of the interface is obtained as a byproduct of the numerical solution.The main contribution is the derivation of a special equation of state to be imposed at the interface in order to avoid non-physical oscillations.Suitable boundary conditions at the piston that guarantee second order convergence are described.We compare the solution of the piston problem to other results available in the literature and to a reference solution obtained within the adiabatic approximation.A shock-interface interaction problem is also treated.The results on these tests are in good agreement with those obtained by other methods.
文摘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.
基金This study was supported by the National Natural Science Foundation of China under contract No.60134010.
文摘A two-dimensional nonoscillatory central difference scheme was extended to the shallow water equations. A high-resolution numerical method for solving the shallow water equations was presented. In order to prevent oscillation, the nonlinear limiter is employed to approximate the discrete slopes. The main advantage of the presented method is simplicity comparable with the upwind schemes. This method does not require Riemann solvers or some form of flux difference splitting methods. Furthermore, the discrete derivatives of flux can be approximated by the component-wise approach and thus the computation of Jacobian can be avoided. The method retains high resolution and high accuracy similar to the upwind results. It is applied to simulating several tests, including circular dam-break problem, shock focusing problem and partial dam-break problem. The results are in good agreement with the numerical results obtained by other methods. The simulated results also demonstrate that the presented method is stable and efficient.
文摘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 elaborated the bidding scheme design of the Central Recreational Area in Linkong Economic Zone from the perspectives of site analysis, determination of properties, design theme and principle, scheme layout and conception, in order to explore the potential enlightenments of the scheme design process, and disclose significance of site analysis and its close relationship with the scheme design.
基金Supported by the National Natural Science Foundation of China (Nos.50876114 and 10602043)the Program for New Century Excellent Talents in University,and the Scientific Research Key Project Fund of Ministry of Education (No.106142)
文摘The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics (English Edition), 2007, 28(7), 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.
文摘This study examines the evolution of fiscal relations between the central state and local governments in China since the founding of the People's Republic The emphasis is placed on the changing arrangement for the central state's collection of taxes and allocation of capital among provinces The central local relation has experienced constant reorganization since 1949, but the tax division scheme introduced in the early 1990s represents the most significant change Prior to the 1990s, the central local fiscal relation was shaped by a highly centralized administrative and economic system without rational and scientific design for revenue collection and expenditure The financial responsibility system established after the reforms contributed to regional economic development during the first decade of reforms, but it was a temporary and transitional arrangement that does not meet the requirement of rational resources arrangement according to free market forces The implementation of the tax division scheme has stopped the decline of the ratio of fiscal revenue to gross domestic product, raised the percentage of central revenue in total national income, and strengthened the central function for macro economic control The mechanism for tax refund under the tax division scheme required further improvements Several proposals are made in this study: 1) tax classification should be adjusted according to administration; 2) a scientific and standardized system for regional transfer payment should be developed; and 3) the power for tax legislation should be delineated according to the rational division of administration between the central and local government
文摘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.
基金supported in part by National Science Fund for Distinguished Young Scholars under Grant No.60725105National Key Basic Research Program of China ( 973 Program ) under Grant No.2009CB320404+2 种基金Program for Changjiang Scholars and Innovative Research Team in University under Grant No.IRT0852National Natural Science Foundation of China under Grants No.60972047, 61072068111 Project under Grant No.B08038
文摘This paper proposes a negotiation-based TDMA scheme for ad hoc networks, which was modeled as an asynchronous myopic repeated game. Compared to the traditional centralized TDMA schemes, our scheme operates in a decentralized manner and is scalable to topology changes. Simulation results show that, with respect to the coloring quality, the performance of our scheme is close to that of the classical centralized algorithms with much lower complexity.
基金Research was supported in part by the ONR Grant N00014-2112773.
文摘Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisit the question of non-oscillatory behavior of high-resolution central schemes in terms of the slope limiter proposed by van Albada et al.(Astron Astrophys 108:76–84,1982).The van Albada(vA)limiter is smoother near extrema,and consequently,in many cases,it outperforms the results obtained using the standard minmod limiter.In particular,we prove that the vA limiter ensures the one-dimensional Total-Variation Diminishing(TVD)stability and demonstrate that it yields noticeable improvement in computation of one-and two-dimensional systems.