The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm t...The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件...以Aspen Open Solver接口集中的非线性代数方程组(NLA)部分作为研究对象,在对接口集进行系统地分析之后,利用AspenTech提供的接口代码将分别基于梯度和非基于梯度的四种求解算法嵌入生成solver组件,并实现用Aspen Plus调用该solver组件观察各种算法嵌入的结果。展开更多
In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open ...In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open FOAM.Reynolds-averaged Navier？Stokes（RANS） equations are chosen as governing equations and the volume of fluid（VOF） is employed to capture the two phases interface.Incoming wave group is generated by imposing the boundary conditions of the tank inlet.A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum.This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank.After 3 to 4 iterations,the measured spectrum agrees well with the target one.In order to validate this method,several wave spectra are chosen and validated in the numerical wave tank,with comparison between the final measured and target spectra.In order to investigate a practical situation,a modified Wigley hull is placed in the wave tank with incoming irregular waves.The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses,showing good agreements with the measurements.展开更多
基金partial financial support from Gazpromneft Science and Technology Center。
文摘The aim of this study is to create a fast and stable iterative technique for numerical solution of a quasi-linear elliptic pressure equation. We developed a modified version of the Anderson acceleration(AA)algorithm to fixed-point(FP) iteration method. It computes the approximation to the solutions at each iteration based on the history of vectors in extended space, which includes the vector of unknowns, the discrete form of the operator, and the equation's right-hand side. Several constraints are applied to AA algorithm, including a limitation of the time step variation during the iteration process, which allows switching to the base FP iterations to maintain convergence. Compared to the base FP algorithm, the improved version of the AA algorithm enables a reliable and rapid convergence of the iterative solution for the quasi-linear elliptic pressure equation describing the flow of particle-laden yield-stress fluids in a narrow channel during hydraulic fracturing, a key technology for stimulating hydrocarbon-bearing reservoirs. In particular, the proposed AA algorithm allows for faster computations and resolution of unyielding zones in hydraulic fractures that cannot be calculated using the FP algorithm. The quasi-linear elliptic pressure equation under consideration describes various physical processes, such as the displacement of fluids with viscoplastic rheology in a narrow cylindrical annulus during well cementing,the displacement of cross-linked gel in a proppant pack filling hydraulic fractures during the early stage of well production(fracture flowback), and multiphase filtration in a rock formation. We estimate computational complexity of the developed algorithm as compared to Jacobian-based algorithms and show that the performance of the former one is higher in modelling of flows of viscoplastic fluids. We believe that the developed algorithm is a useful numerical tool that can be implemented in commercial simulators to obtain fast and converged solutions to the non-linear problems described above.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.
基金financially supported by National Natural Science Foundation of China(Grant Nos.51379125,51411130131,11432009,and 51490675)the Chang Jiang Scholars Program(Grant No.T2014099)+3 种基金the Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(Grant No.2016-23)the Foundation of State key Laboratory of Ocean Engineering(Grant No.GKZD010065)the Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)center for HPC at Shanghai Jiao Tong University,and Lloyd’s Register Foundation(LRF)
文摘In this paper,a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented.The naoe-FOAM-SJTU solver is developed using an open source tool kit,Open FOAM.Reynolds-averaged Navier？Stokes（RANS） equations are chosen as governing equations and the volume of fluid（VOF） is employed to capture the two phases interface.Incoming wave group is generated by imposing the boundary conditions of the tank inlet.A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum.This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank.After 3 to 4 iterations,the measured spectrum agrees well with the target one.In order to validate this method,several wave spectra are chosen and validated in the numerical wave tank,with comparison between the final measured and target spectra.In order to investigate a practical situation,a modified Wigley hull is placed in the wave tank with incoming irregular waves.The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses,showing good agreements with the measurements.
