A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are ...A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.展开更多
In complex networks,identifying influential spreader is of great significance for improving the reliability of networks and ensuring the safe and effective operation of networks.Nowadays,it is widely used in power net...In complex networks,identifying influential spreader is of great significance for improving the reliability of networks and ensuring the safe and effective operation of networks.Nowadays,it is widely used in power networks,aviation networks,computer networks,and social networks,and so on.Traditional centrality methods mainly include degree centrality,closeness centrality,betweenness centrality,eigenvector centrality,k-shell,etc.However,single centrality method is onesided and inaccurate,and sometimes many nodes have the same centrality value,namely the same ranking result,which makes it difficult to distinguish between nodes.According to several classical methods of identifying influential nodes,in this paper we propose a novel method that is more full-scaled and universally applicable.Taken into account in this method are several aspects of node’s properties,including local topological characteristics,central location of nodes,propagation characteristics,and properties of neighbor nodes.In view of the idea of the multi-attribute decision-making,we regard the basic centrality method as node’s attribute and use the entropy weight method to weigh different attributes,and obtain node’s combined centrality.Then,the combined centrality is applied to the gravity law to comprehensively identify influential nodes in networks.Finally,the classical susceptible-infected-recovered(SIR)model is used to simulate the epidemic spreading in six real-society networks.Our proposed method not only considers the four topological properties of nodes,but also emphasizes the influence of neighbor nodes from the aspect of gravity.It is proved that the new method can effectively overcome the disadvantages of single centrality method and increase the accuracy of identifying influential nodes,which is of great significance for monitoring and controlling the complex networks.展开更多
An investigation is made of the magnetic Rayleigh problem where a semi_infinite plate is given an impulsive motion and thereafter moves with constant velocity in a non_Newtonian power law fluid of infinite extent. The...An investigation is made of the magnetic Rayleigh problem where a semi_infinite plate is given an impulsive motion and thereafter moves with constant velocity in a non_Newtonian power law fluid of infinite extent. The solution of this highly non_linear problem is obtained by means of the transformation group theoretic approach. The one_parameter group transformation reduces the number of independent variables by one and the governing partial differential equation with the boundary conditions reduce to an ordinary differential equation with the appropriate boundary conditions. Effect of the some parameters on the velocity u(y,t) has been studied and the results are plotted.展开更多
The hierarchical reconstruction (HR) [Liu, Shu, Tadmor and Zhang, SINUM '07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectra...The hierarchical reconstruction (HR) [Liu, Shu, Tadmor and Zhang, SINUM '07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.展开更多
This paper is concerned with Godunov's scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov's scheme, which satisfies the boundary entropy condition, was given. By use of the ...This paper is concerned with Godunov's scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov's scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.展开更多
The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux...The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.展开更多
This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At fi...This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.展开更多
This paper summarizes a Riemann-solver-free spacetime discontinuous Galerkin method developed for general conservation laws. The method integrates the best features of the spacetime Conservation Element/Solution Eleme...This paper summarizes a Riemann-solver-free spacetime discontinuous Galerkin method developed for general conservation laws. The method integrates the best features of the spacetime Conservation Element/Solution Element (CE/SE) method and the discontinuous Galerkin (DG) method. The core idea is to construct a staggered spacetime mesh through alternate cell-centered CEs and vertex-centered CEs within each time step. Inside each SE, the solution is approximated using high-order spacetime DG basis polynomials. The spacetime flux conservation is enforced inside each CE using the DG concept. The unknowns are stored at both vertices and cell centroids of the spatial mesh. However, the solutions at vertices and cell centroids are updated at different time levels within each time step in an alternate fashion. Thanks to the staggered spacetime formulation, there are no left and right states for the solution at the spacetime interface. Instead, the solution available to evaluate the flux is continuous across the interface. Therefore, no (approximate) Riemann solvers are needed to provide a unique numerical flux. The current method can be used to solve arbitrary conservation laws including the compressible Euler equations, shallow water equations and magnetohydrodynamics (MHD) equations without the need of any form of Riemann solvers. A set of benchmark problems of various conservation laws are presented to demonstrate the accuracy of the method.展开更多
Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to sol...Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.展开更多
The Henry’s law constant of volatiles in polymer systems is a crucial parameter reflecting the gas-liquid equilibrium,which is very important for devolatilization.In this research,polyolefin elastomer(POE)-cyclohexan...The Henry’s law constant of volatiles in polymer systems is a crucial parameter reflecting the gas-liquid equilibrium,which is very important for devolatilization.In this research,polyolefin elastomer(POE)-cyclohexane and polydimethylsiloxane(PDMS)-hexane systems were studied,and the Henry’s law constant was obtained by measuring the gas phase equilibrium partial pressure when polymer solutions containing different mass fractions of volatiles reached a saturated state.The effects of temperature,type of volatiles,and polymer viscosity on the gas phase equilibrium partial pressure and Henry’s law constant of the volatiles were investigated.The results indicate that,with the increase of temperature and polymer viscosity,the gas phase equilibrium partial pressure and Henry’s law constant of volatiles increase.As temperature increases,the solubility of gas in liquid decreases.The relationship between the Henry’s law constant and temperature is consistent with the Arrhenius law.In the PDMS-hexane system,the gas phase equilibrium partial pressure and Henry’s law constant of n-hexane are higher than those of cyclohexane.The obtained Henry’s law constants can be used as a reference for perfecting the devolatilization process and improving the devolatilization effect.展开更多
In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified cons...In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.展开更多
基金supported by the National Natural Science Foundation of China(11390363 and 11172041)Beijing Higher Education Young Elite Teacher Project(YETP1190)
文摘A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.
基金Project support by the National Key Research and Development Program of China(Grant No.2018YFF0301000)the National Natural Science Foundation of China(Grant Nos.71673161 and 71790613)。
文摘In complex networks,identifying influential spreader is of great significance for improving the reliability of networks and ensuring the safe and effective operation of networks.Nowadays,it is widely used in power networks,aviation networks,computer networks,and social networks,and so on.Traditional centrality methods mainly include degree centrality,closeness centrality,betweenness centrality,eigenvector centrality,k-shell,etc.However,single centrality method is onesided and inaccurate,and sometimes many nodes have the same centrality value,namely the same ranking result,which makes it difficult to distinguish between nodes.According to several classical methods of identifying influential nodes,in this paper we propose a novel method that is more full-scaled and universally applicable.Taken into account in this method are several aspects of node’s properties,including local topological characteristics,central location of nodes,propagation characteristics,and properties of neighbor nodes.In view of the idea of the multi-attribute decision-making,we regard the basic centrality method as node’s attribute and use the entropy weight method to weigh different attributes,and obtain node’s combined centrality.Then,the combined centrality is applied to the gravity law to comprehensively identify influential nodes in networks.Finally,the classical susceptible-infected-recovered(SIR)model is used to simulate the epidemic spreading in six real-society networks.Our proposed method not only considers the four topological properties of nodes,but also emphasizes the influence of neighbor nodes from the aspect of gravity.It is proved that the new method can effectively overcome the disadvantages of single centrality method and increase the accuracy of identifying influential nodes,which is of great significance for monitoring and controlling the complex networks.
文摘An investigation is made of the magnetic Rayleigh problem where a semi_infinite plate is given an impulsive motion and thereafter moves with constant velocity in a non_Newtonian power law fluid of infinite extent. The solution of this highly non_linear problem is obtained by means of the transformation group theoretic approach. The one_parameter group transformation reduces the number of independent variables by one and the governing partial differential equation with the boundary conditions reduce to an ordinary differential equation with the appropriate boundary conditions. Effect of the some parameters on the velocity u(y,t) has been studied and the results are plotted.
基金Research was supported in part by NSF grant DMS-0800612Research was supported by Applied Mathematics program of the US DOE Office of Advanced Scientific Computing ResearchThe Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830
文摘The hierarchical reconstruction (HR) [Liu, Shu, Tadmor and Zhang, SINUM '07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10671120)
文摘This paper is concerned with Godunov's scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov's scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.
文摘The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.
文摘This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.
文摘This paper summarizes a Riemann-solver-free spacetime discontinuous Galerkin method developed for general conservation laws. The method integrates the best features of the spacetime Conservation Element/Solution Element (CE/SE) method and the discontinuous Galerkin (DG) method. The core idea is to construct a staggered spacetime mesh through alternate cell-centered CEs and vertex-centered CEs within each time step. Inside each SE, the solution is approximated using high-order spacetime DG basis polynomials. The spacetime flux conservation is enforced inside each CE using the DG concept. The unknowns are stored at both vertices and cell centroids of the spatial mesh. However, the solutions at vertices and cell centroids are updated at different time levels within each time step in an alternate fashion. Thanks to the staggered spacetime formulation, there are no left and right states for the solution at the spacetime interface. Instead, the solution available to evaluate the flux is continuous across the interface. Therefore, no (approximate) Riemann solvers are needed to provide a unique numerical flux. The current method can be used to solve arbitrary conservation laws including the compressible Euler equations, shallow water equations and magnetohydrodynamics (MHD) equations without the need of any form of Riemann solvers. A set of benchmark problems of various conservation laws are presented to demonstrate the accuracy of the method.
文摘Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.
基金Polyolefin Elastomer Technology Development project(2020B-2619).
文摘The Henry’s law constant of volatiles in polymer systems is a crucial parameter reflecting the gas-liquid equilibrium,which is very important for devolatilization.In this research,polyolefin elastomer(POE)-cyclohexane and polydimethylsiloxane(PDMS)-hexane systems were studied,and the Henry’s law constant was obtained by measuring the gas phase equilibrium partial pressure when polymer solutions containing different mass fractions of volatiles reached a saturated state.The effects of temperature,type of volatiles,and polymer viscosity on the gas phase equilibrium partial pressure and Henry’s law constant of the volatiles were investigated.The results indicate that,with the increase of temperature and polymer viscosity,the gas phase equilibrium partial pressure and Henry’s law constant of volatiles increase.As temperature increases,the solubility of gas in liquid decreases.The relationship between the Henry’s law constant and temperature is consistent with the Arrhenius law.In the PDMS-hexane system,the gas phase equilibrium partial pressure and Henry’s law constant of n-hexane are higher than those of cyclohexane.The obtained Henry’s law constants can be used as a reference for perfecting the devolatilization process and improving the devolatilization effect.
文摘In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.