In this paper,we investigate the performance of the exponential time differencing(ETD)method applied to the rotating shallow water equations.Comparing with explicit time stepping of the same order accuracy in time,the...In this paper,we investigate the performance of the exponential time differencing(ETD)method applied to the rotating shallow water equations.Comparing with explicit time stepping of the same order accuracy in time,the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability.To accelerate the ETD simulations,we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition.By dividing the original problem into many subdomain problems of smaller sizes and solving them locally,the proposed approach could speed up the calculation of matrix exponential vector products.Several standard test cases for shallow water equations of one or multiple layers are considered.The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.展开更多
This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the frac...This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the fracture is represented as an interface between subdomains and is assumed to have larger permeability than the surrounding area.We develop three global-in-time domain decomposition methods coupled with operator splitting for the reduced fracture model,where the advection and the diffusion are treated separately by different numerical schemes and with different time steps.Importantly,smaller time steps can be used in the fracture-interface than in the subdomains.The first two methods are based on the physical transmission conditions,while the third one is based on the optimized Schwarz waveform relaxation approach with Ventcel-Robin transmission conditions.A discrete space-time interface system is formulated for each method and is solved iteratively and globally in time.Numerical results for two-dimensional problems with various P′eclet numbers and different types of fracture are presented to illustrate and compare the convergence and accuracy in time of the proposed methods with local time stepping.展开更多
基金supported by U.S.Department of Energy through the grants DE-SC0016540,DE-SC0020270U.S.National Science Foundation through the grant DMS-1912626,Office of the Vice President for Research at the University of South Carolina through an ASPIRE grantNatural Science Foundation of China through the grant 11871454.
文摘In this paper,we investigate the performance of the exponential time differencing(ETD)method applied to the rotating shallow water equations.Comparing with explicit time stepping of the same order accuracy in time,the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability.To accelerate the ETD simulations,we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition.By dividing the original problem into many subdomain problems of smaller sizes and solving them locally,the proposed approach could speed up the calculation of matrix exponential vector products.Several standard test cases for shallow water equations of one or multiple layers are considered.The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.
基金partially supported by the US National Science Foundation under grant number DMS-1912626.
文摘This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the fracture is represented as an interface between subdomains and is assumed to have larger permeability than the surrounding area.We develop three global-in-time domain decomposition methods coupled with operator splitting for the reduced fracture model,where the advection and the diffusion are treated separately by different numerical schemes and with different time steps.Importantly,smaller time steps can be used in the fracture-interface than in the subdomains.The first two methods are based on the physical transmission conditions,while the third one is based on the optimized Schwarz waveform relaxation approach with Ventcel-Robin transmission conditions.A discrete space-time interface system is formulated for each method and is solved iteratively and globally in time.Numerical results for two-dimensional problems with various P′eclet numbers and different types of fracture are presented to illustrate and compare the convergence and accuracy in time of the proposed methods with local time stepping.
