Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of diffe...Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of different components that require additional degrees of freedom for representation in computational fluid dynamics. It is challenging to construct an accurate and efficient numerical algorithm to faithfully represent multiscale flow physics across different regimes. In this paper, a unified gas-kinetic scheme(UGKS) is developed to study non-equilibrium multicomponent gaseous flows. Based on the Boltzmann kinetic equation, an analytical space-time evolving solution is used to construct the discretized equations of gas dynamics directly according to cell size and scales of time steps, i.e., the so-called direct modeling method. With the variation in the ratio of the numerical time step to the local particle collision time(or the cell size to the local particle mean free path), the UGKS automatically recovers all scale-dependent flows over the given domain and provides a continuous spectrum of the gas dynamics. The performance of the proposed unified scheme is fully validated through numerical experiments.The UGKS can be a valuable tool to study multiscale and multicomponent flow physics.展开更多
Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though loo...Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.展开更多
Computational techniques are invaluable to the continued success and development of Magnetic Resonance Imaging (MRI) and to its widespread applications. New processing methods are essential for addressing issues at ea...Computational techniques are invaluable to the continued success and development of Magnetic Resonance Imaging (MRI) and to its widespread applications. New processing methods are essential for addressing issues at each stage of MRI techniques. In this study, we present new sets of non-exponential generating functions representing the NMR transverse magnetizations and signals which are mathematically designed based on the theory and dynamics of the Bloch NMR flow equations. These signals are functions of many spinning nuclei of materials and can be used to obtain information observed in all flow systems. The Bloch NMR flow equations are solved using the Boubaker polynomial expansion scheme (BPES) and analytically connect most of the experimentally valuable NMR parameters in a simplified way for general analyses of magnetic resonance imaging with adiabatic condition.展开更多
For a stochastic differential equation with non-Lipschitz coefficients, we construct, by Euler scheme, a measurable flow of the solution, and we prove the solution is a Markov process.
Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly im...Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.展开更多
In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity a...In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity and linear velocity-density relationship. Macroscopic traffic flow model is first developed by Lighthill Whitham and Richards (LWR) and used to study traffic flow by collective variables such as flow rate, velocity and density. The LWR model is treated as an initial value problem and its numerical simulations are presented using numerical schemes. A variety of numerical schemes are available in literature to solve first order hyperbolic equations. Of these the well-known ones include one-dimensional explicit: Upwind, Downwind, FTCS, and Lax-Friedrichs schemes. Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes, a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. Tolesa numerical scheme is one of the conditionally stable and highest rates of convergence schemes. All the said numerical schemes are applied to solve advection equation pertaining traffic flows. Also the one-dimensional explicit Tolesa numerical scheme is another alternative numerical scheme to solve advection equation and apply to traffic flows model like other well-known one-dimensional explicit schemes. The effect of density of cars on the overall interactions of the vehicles along a given length of the highway and time are investigated. Graphical representations of density profile, velocity profile, flux profile, and in general the fundamental diagrams of vehicles on the highway with different time levels are illustrated. These concepts and results have been arranged systematically in this paper.展开更多
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.展开更多
A finite difference lattice Boltzmann method of second-order accuracy in time is developed based on non-oscillatory scheme with no-free-parameter dissipation （NND） difference scheme in this paper. The NND lattice Bo...A finite difference lattice Boltzmann method of second-order accuracy in time is developed based on non-oscillatory scheme with no-free-parameter dissipation （NND） difference scheme in this paper. The NND lattice Boltzmann method is used to simulate high-speed flows by constructing a new equilibrium distribution function of the lattice Boltzmann method. Compared with a variation of lattice Boltzmann method developed by Qu, et al., the present method can capture shock waves and handle oscillations of high velocity flows accurately in larger time steps and in shorter computing time. Numerical results indicate the correctness and capability of simulating shock wave interactions of the NND lattice Boltzmann method.展开更多
The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence mod...The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.展开更多
A fluid dynamic traffic flow model based on a non-linear velocity-density function is considered. The model provides a quasi-linear first order hyperbolic partial differential equation which is appended with initial a...A fluid dynamic traffic flow model based on a non-linear velocity-density function is considered. The model provides a quasi-linear first order hyperbolic partial differential equation which is appended with initial and boundary data and turns out an initial boundary value problem (IBVP). A first order explicit finite difference scheme of the IBVP known as Lax-Friedrich’s scheme for our model is presented and a well-posedness and stability condition of the scheme is established. The numerical scheme is implemented in order to perform the numerical features of error estimation and rate of convergence. Fundamental diagram, density, velocity and flux profiles are presented.展开更多
This study sought to forecast water flow and sediment flux in the scheme as potential contributions for improved management in the Chókwè Irrigation Scheme (CIS). Fieldwork data was collected during dry (DS)...This study sought to forecast water flow and sediment flux in the scheme as potential contributions for improved management in the Chókwè Irrigation Scheme (CIS). Fieldwork data was collected during dry (DS) and wet (WS) seasons. Flow measurement was performed at 9 stations using a calibrated flow meter OTT C31. Water flow and sediment flux from 2004 to 2019 were used. Hydrodynamic forecast simulations were performed using Mann-Kendall test and ARIMA model for determination of temporal trends. Findings suggest higher values during DS for water discharge and sediment flux. Mann-Kendall test for sediment discharge trends was not significant at 95% significance level, except for the Offtake in WS. ARIMA test for the sediment discharges, at the Intake, for DS and WS, sediments were well described by the ARIMA model and gave a good result for the sediments. Good fit between the observed and the predicted ARIMA model was found. ARIMA model for sediment discharge at CIS based on AIC has a good fit for AR (p = 1), whereby, at the Intake the ARIMA p-value was 0.822 and 0.932, for WS and DS, respectively. Whilst in the Offtake, the ARIMA p-value was 0.877 and 0.893, respectively. These results can be used to improve the CIS management, both for water flow and sediment flux.展开更多
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...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.展开更多
Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics ...Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.展开更多
In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which prov...In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which provides a non-linear second-order parabolic partial differential equation. The analytical solution of the diffusion-type traffic flow model is very complicated to approximate the initial density of the Cauchy problem as a function of x from given data and it may cause a huge error. For the complexity of the analytical solution, the numerical solution is performed by implementing an explicit upwind, explicitly centered, and second-order Lax-Wendroff scheme for the numerical solution. From the comparison of relative error among these three schemes, it is observed that Lax-Wendroff scheme gives less error than the explicit upwind and explicit centered difference scheme. The numerical, analytical analysis and comparative result discussion bring out the fact that the Lax-Wendroff scheme with exponential velocity-density relation of diffusion type traffic flow model is suitable for the congested area and shows a better fit in traffic-congested regions.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11772281,91530319,and 11521091)the Hong Kong Research Grant Council(Nos.16207715 and 16206617)
文摘Compressible flows exhibit a diverse set of behaviors, where individual particle transports and their collective dynamics play different roles at different scales. At the same time, the atmosphere is composed of different components that require additional degrees of freedom for representation in computational fluid dynamics. It is challenging to construct an accurate and efficient numerical algorithm to faithfully represent multiscale flow physics across different regimes. In this paper, a unified gas-kinetic scheme(UGKS) is developed to study non-equilibrium multicomponent gaseous flows. Based on the Boltzmann kinetic equation, an analytical space-time evolving solution is used to construct the discretized equations of gas dynamics directly according to cell size and scales of time steps, i.e., the so-called direct modeling method. With the variation in the ratio of the numerical time step to the local particle collision time(or the cell size to the local particle mean free path), the UGKS automatically recovers all scale-dependent flows over the given domain and provides a continuous spectrum of the gas dynamics. The performance of the proposed unified scheme is fully validated through numerical experiments.The UGKS can be a valuable tool to study multiscale and multicomponent flow physics.
文摘Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.
文摘Computational techniques are invaluable to the continued success and development of Magnetic Resonance Imaging (MRI) and to its widespread applications. New processing methods are essential for addressing issues at each stage of MRI techniques. In this study, we present new sets of non-exponential generating functions representing the NMR transverse magnetizations and signals which are mathematically designed based on the theory and dynamics of the Bloch NMR flow equations. These signals are functions of many spinning nuclei of materials and can be used to obtain information observed in all flow systems. The Bloch NMR flow equations are solved using the Boubaker polynomial expansion scheme (BPES) and analytically connect most of the experimentally valuable NMR parameters in a simplified way for general analyses of magnetic resonance imaging with adiabatic condition.
文摘For a stochastic differential equation with non-Lipschitz coefficients, we construct, by Euler scheme, a measurable flow of the solution, and we prove the solution is a Markov process.
文摘Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.
文摘In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity and linear velocity-density relationship. Macroscopic traffic flow model is first developed by Lighthill Whitham and Richards (LWR) and used to study traffic flow by collective variables such as flow rate, velocity and density. The LWR model is treated as an initial value problem and its numerical simulations are presented using numerical schemes. A variety of numerical schemes are available in literature to solve first order hyperbolic equations. Of these the well-known ones include one-dimensional explicit: Upwind, Downwind, FTCS, and Lax-Friedrichs schemes. Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes, a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. Tolesa numerical scheme is one of the conditionally stable and highest rates of convergence schemes. All the said numerical schemes are applied to solve advection equation pertaining traffic flows. Also the one-dimensional explicit Tolesa numerical scheme is another alternative numerical scheme to solve advection equation and apply to traffic flows model like other well-known one-dimensional explicit schemes. The effect of density of cars on the overall interactions of the vehicles along a given length of the highway and time are investigated. Graphical representations of density profile, velocity profile, flux profile, and in general the fundamental diagrams of vehicles on the highway with different time levels are illustrated. These concepts and results have been arranged systematically in this paper.
文摘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.
基金Project supported by the Aeronautics Science Foundation (Grant No.20061431)the Program for Changjiang Scholars and Innovative Research Team in University (Grant No.IRT0844)
文摘A finite difference lattice Boltzmann method of second-order accuracy in time is developed based on non-oscillatory scheme with no-free-parameter dissipation （NND） difference scheme in this paper. The NND lattice Boltzmann method is used to simulate high-speed flows by constructing a new equilibrium distribution function of the lattice Boltzmann method. Compared with a variation of lattice Boltzmann method developed by Qu, et al., the present method can capture shock waves and handle oscillations of high velocity flows accurately in larger time steps and in shorter computing time. Numerical results indicate the correctness and capability of simulating shock wave interactions of the NND lattice Boltzmann method.
基金Projects(59375211,10771178,10676031) supported by the National Natural Science Foundation of ChinaProject(07A068) supported by the Key Project of Hunan Education CommissionProject(2005CB321702) supported by the National Key Basic Research Program of China
文摘The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.
文摘A fluid dynamic traffic flow model based on a non-linear velocity-density function is considered. The model provides a quasi-linear first order hyperbolic partial differential equation which is appended with initial and boundary data and turns out an initial boundary value problem (IBVP). A first order explicit finite difference scheme of the IBVP known as Lax-Friedrich’s scheme for our model is presented and a well-posedness and stability condition of the scheme is established. The numerical scheme is implemented in order to perform the numerical features of error estimation and rate of convergence. Fundamental diagram, density, velocity and flux profiles are presented.
文摘This study sought to forecast water flow and sediment flux in the scheme as potential contributions for improved management in the Chókwè Irrigation Scheme (CIS). Fieldwork data was collected during dry (DS) and wet (WS) seasons. Flow measurement was performed at 9 stations using a calibrated flow meter OTT C31. Water flow and sediment flux from 2004 to 2019 were used. Hydrodynamic forecast simulations were performed using Mann-Kendall test and ARIMA model for determination of temporal trends. Findings suggest higher values during DS for water discharge and sediment flux. Mann-Kendall test for sediment discharge trends was not significant at 95% significance level, except for the Offtake in WS. ARIMA test for the sediment discharges, at the Intake, for DS and WS, sediments were well described by the ARIMA model and gave a good result for the sediments. Good fit between the observed and the predicted ARIMA model was found. ARIMA model for sediment discharge at CIS based on AIC has a good fit for AR (p = 1), whereby, at the Intake the ARIMA p-value was 0.822 and 0.932, for WS and DS, respectively. Whilst in the Offtake, the ARIMA p-value was 0.877 and 0.893, respectively. These results can be used to improve the CIS management, both for water flow and sediment flux.
文摘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.
基金supported by the National Natural Science Foundation of China (No.10621062)the Research Fund for Next Generation of General Armament Department (No.9140A13050207KG29)
文摘Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.
文摘In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which provides a non-linear second-order parabolic partial differential equation. The analytical solution of the diffusion-type traffic flow model is very complicated to approximate the initial density of the Cauchy problem as a function of x from given data and it may cause a huge error. For the complexity of the analytical solution, the numerical solution is performed by implementing an explicit upwind, explicitly centered, and second-order Lax-Wendroff scheme for the numerical solution. From the comparison of relative error among these three schemes, it is observed that Lax-Wendroff scheme gives less error than the explicit upwind and explicit centered difference scheme. The numerical, analytical analysis and comparative result discussion bring out the fact that the Lax-Wendroff scheme with exponential velocity-density relation of diffusion type traffic flow model is suitable for the congested area and shows a better fit in traffic-congested regions.
