Recent advances in deep learning have expanded new possibilities for fluid flow simulation in petroleum reservoirs.However,the predominant approach in existing research is to train neural networks using high-fidelity ...Recent advances in deep learning have expanded new possibilities for fluid flow simulation in petroleum reservoirs.However,the predominant approach in existing research is to train neural networks using high-fidelity numerical simulation data.This presents a significant challenge because the sole source of authentic wellbore production data for training is sparse.In response to this challenge,this work introduces a novel architecture called physics-informed neural network based on domain decomposition(PINN-DD),aiming to effectively utilize the sparse production data of wells for reservoir simulation with large-scale systems.To harness the capabilities of physics-informed neural networks(PINNs)in handling small-scale spatial-temporal domain while addressing the challenges of large-scale systems with sparse labeled data,the computational domain is divided into two distinct sub-domains:the well-containing and the well-free sub-domain.Moreover,the two sub-domains and the interface are rigorously constrained by the governing equations,data matching,and boundary conditions.The accuracy of the proposed method is evaluated on two problems,and its performance is compared against state-of-the-art PINNs through numerical analysis as a benchmark.The results demonstrate the superiority of PINN-DD in handling large-scale reservoir simulation with limited data and show its potential to outperform conventional PINNs in such scenarios.展开更多
We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. ...We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge.展开更多
To achieve high parallel efficiency for the global MASNUM surface wave model, the algorithm of an irregular quasirectangular domain decomposition and related serializing of calculating points and data exchanging schem...To achieve high parallel efficiency for the global MASNUM surface wave model, the algorithm of an irregular quasirectangular domain decomposition and related serializing of calculating points and data exchanging schemes are developed and conducted, based on the environment of Message Passing Interface(MPI). The new parallel version of the surface wave model is tested for parallel computing on the platform of the Sunway BlueLight supercomputer in the National Supercomputing Center in Jinan. The testing involves four horizontal resolutions, which are 1°×1°,(1/2)°×(1/2)°,(1/4)°×(1/4)°, and(1/8)°×(1/8)°. These tests are performed without data Input/Output(IO) and the maximum amount of processors used in these tests reaches to 131072. The testing results show that the computing speeds of the model with different resolutions are all increased with the increasing of numbers of processors. When the number of processors is four times that of the base processor number, the parallel efficiencies of all resolutions are greater than 80%. When the number of processors is eight times that of the base processor number, the parallel efficiency of tests with resolutions of 1°×1°,(1/2)°×(1/2)° and(1/4)°×(1/4)° is greater than 80%, and it is 62% for the test with a resolution of(1/8)°×(1/8)° using 131072 processors, which is the nearly all processors of Sunway BlueLight. When the processor's number is 24 times that of the base processor number, the parallel efficiencies for tests with resolutions of 1°×1°,(1/2)°×(1/2)°, and(1/4)°×(1/4)° are 72%, 62%, and 38%, respectively. The speedup and parallel efficiency indicate that the irregular quasi-rectangular domain decomposition and serialization schemes lead to high parallel efficiency and good scalability for a global numerical wave model.展开更多
The near-wall domain decomposition method(NDD)has proved to be very efficient for modeling near-wall fully turbulent flows.In this paper the NDD is extended to non-equilibrium regimeswith laminar-turbulent transition(...The near-wall domain decomposition method(NDD)has proved to be very efficient for modeling near-wall fully turbulent flows.In this paper the NDD is extended to non-equilibrium regimeswith laminar-turbulent transition(LTT)for the first time.The LTT is identified with the use of the e^(N)-method which is applied to both incompressible and compressible flows.TheNDD ismodified to take into account LTT in an efficientway.In addition,implementation of the intermittency expands the capabilities of NDD to model non-equilibrium turbulent flows with transition.Performance of the modified NDD approach is demonstrated on various test problems of subsonic and supersonic flows past a flat plate,a supersonic flow over a compression corner and a planar shock wave impinging on a turbulent boundary layer.The results of modeling with and without decomposition are compared in terms of wall friction and show good agreement with each other while NDD significantly reducing computational resources needed.It turns out that the NDD can reduce the computational time as much as three times while retaining practically the same accuracy of prediction.展开更多
In this paper the parallel computing of a grid-point nine-level atmospheric general circulation model on the Dawn 1000 is introduced. The model was developed by the Institute of Atmospheric Physics (IAP), Chinese Acad...In this paper the parallel computing of a grid-point nine-level atmospheric general circulation model on the Dawn 1000 is introduced. The model was developed by the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS). The Dawn 1000 is a MIMD massive parallel computer made by National Research Center for Intelligent Computer (NCIC), CAS. A two-dimensional domain decomposition method is adopted to perform the parallel computing. The potential ways to increase the speed-up ratio and exploit more resources of future massively parallel supercomputation are also discussed.展开更多
Partition property plays a central role in domain decomposition methods. Existing theory essentially assumes certain ellipticity. We prove the partition property for prod lems without ellipticity which are of practica...Partition property plays a central role in domain decomposition methods. Existing theory essentially assumes certain ellipticity. We prove the partition property for prod lems without ellipticity which are of practical importance. Example applications include implicit schemes applied to degenerate parabolic partial differential equations arising from superconductors, superfluids and liquid crystals. With this partition property, Schwarz algorithms can be applied to general non-elliptic problems with an h-independent optimal convergence rate. Application to the time-dependent Ginzburg-Landau model of superconductivity is illustrated and numerical results are presented.展开更多
Accurate displacement prediction is critical for the early warning of landslides.The complexity of the coupling relationship between multiple influencing factors and displacement makes the accurate prediction of displ...Accurate displacement prediction is critical for the early warning of landslides.The complexity of the coupling relationship between multiple influencing factors and displacement makes the accurate prediction of displacement difficult.Moreover,in engineering practice,insufficient monitoring data limit the performance of prediction models.To alleviate this problem,a displacement prediction method based on multisource domain transfer learning,which helps accurately predict data in the target domain through the knowledge of one or more source domains,is proposed.First,an optimized variational mode decomposition model based on the minimum sample entropy is used to decompose the cumulative displacement into the trend,periodic,and stochastic components.The trend component is predicted by an autoregressive model,and the periodic component is predicted by the long short-term memory.For the stochastic component,because it is affected by uncertainties,it is predicted by a combination of a Wasserstein generative adversarial network and multisource domain transfer learning for improved prediction accuracy.Considering a real mine slope as a case study,the proposed prediction method was validated.Therefore,this study provides new insights that can be applied to scenarios lacking sample data.展开更多
In the Windows XP 64 bit operating system environment, several common PC were used to build a cluster system, establishing the distributed memory parallel (DMP) computing system. A finite element model of whole aircra...In the Windows XP 64 bit operating system environment, several common PC were used to build a cluster system, establishing the distributed memory parallel (DMP) computing system. A finite element model of whole aircraft with about 260 million degrees of freedom (DOF) was developed using three-node and four-node thin shell element and two-node beam element. With the large commercial finite element software MSC.MARC and employing two kinds of domain decomposition method (DDM) respectively, realized the parallel solving for the static strength analysis of the whole aircraft model, which offered a high cost-effective solution for solving large-scale and complex finite element models.展开更多
The building sector significantly contributes to climate change.To improve its carbon footprint,applications like model predictive control and predictive maintenance rely on system models.However,the high modeling eff...The building sector significantly contributes to climate change.To improve its carbon footprint,applications like model predictive control and predictive maintenance rely on system models.However,the high modeling effort hinders practical application.Machine learning models can significantly reduce this modeling effort.To ensure a machine learning model’s reliability in all operating states,it is essential to know its validity domain.Operating states outside the validity domain might lead to extrapolation,resulting in unpredictable behavior.This paper addresses the challenge of identifying extrapolation in data-driven building energy system models and aims to raise knowledge about it.For that,a novel approach is proposed that calibrates novelty detection algorithms towards the machine learning model.Suitable novelty detection algorithms are identified through a literature review and a benchmark test with 15 candidates.A subset of five algorithms is then evaluated on building energy systems.First,on two-dimensional data,displaying the results with a novel visualization scheme.Then on more complex multi-dimensional use cases.The methodology performs well,and the validity domain could be approximated.The visualization allows for a profound analysis and an improved understanding of the fundamental effects behind a machine learning model’s validity domain and the extrapolation regimes.展开更多
When carrying out calculations for turbulent flow simulation,one inevitably has to face the choice between accuracy and speed of calculations.In order to simultaneously obtain both a computationally efficient and more...When carrying out calculations for turbulent flow simulation,one inevitably has to face the choice between accuracy and speed of calculations.In order to simultaneously obtain both a computationally efficient and more accurate model,a surrogate model can be built on the basis of some fast special model and knowledge of previous calculations obtained by more accurate base models from various test bases or some results of serial calculations.The objective of this work is to construct a surrogate model which allows to improve the accuracy of turbulent calculations obtained by a special model on unstructured meshes.For this purpose,we use 1D Convolutional Neural Network(CNN)of the encoder-decoder architecture and reduce the problem to a single dimension by applying space-filling curves.Such an approach would have the benefit of being applicable to solutions obtained on unstructured meshes.In this work,a non-local approach is applied where entire flow fields obtained by the special and base models are used as input and ground truth output respectively.Spalart-Allmaras(SA)model and Near-wall Domain Decomposition(NDD)method for SA are taken as the base and special models respectively.The efficiency and accuracy of the obtained surrogate model are demonstrated in a case of supersonic flow over a compression corner with different values for angleαand Reynolds number Re.We conducted an investigation into interpolation and extrapolation by Re and also into interpolation byα.展开更多
We provide a concise review of the exponentially convergent multiscale finite element method(ExpMsFEM)for efficient model reduction of PDEs in heterogeneous media without scale separation and in high-frequency wave pr...We provide a concise review of the exponentially convergent multiscale finite element method(ExpMsFEM)for efficient model reduction of PDEs in heterogeneous media without scale separation and in high-frequency wave propagation.The ExpMsFEM is built on the non-overlapped domain decomposition in the classical MsFEM while enriching the approximation space systematically to achieve a nearly exponential convergence rate regarding the number of basis functions.Unlike most generalizations of the MsFEM in the literature,the ExpMsFEM does not rely on any partition of unity functions.In general,it is necessary to use function representations dependent on the right-hand side to break the algebraic Kolmogorov n-width barrier to achieve exponential convergence.Indeed,there are online and offline parts in the function representation provided by the ExpMsFEM.The online part depends on the right-hand side locally and can be computed in parallel efficiently.The offline part contains basis functions that are used in the Galerkin method to assemble the stiffness matrix;they are all independent of the right-hand side,so the stiffness matrix can be used repeatedly in multi-query scenarios.展开更多
A new penalty-free neural network method,PFNN-2,is presented for solving partial differential equations,which is a subsequent improvement of our previously proposed PFNN method[1].PFNN-2 inherits all advantages of PFN...A new penalty-free neural network method,PFNN-2,is presented for solving partial differential equations,which is a subsequent improvement of our previously proposed PFNN method[1].PFNN-2 inherits all advantages of PFNN in handling the smoothness constraints and essential boundary conditions of self-adjoint problems with complex geometries,and extends the application to a broader range of non-self-adjoint time-dependent differential equations.In addition,PFNN-2 introduces an overlapping domain decomposition strategy to substantially improve the training efficiency without sacrificing accuracy.Experiments results on a series of partial differential equations are reported,which demonstrate that PFNN-2 can outperform state-of-the-art neural network methods in various aspects such as numerical accuracy,convergence speed,and parallel scalability.展开更多
基金funded by the National Natural Science Foundation of China(Grant No.52274048)Beijing Natural Science Foundation(Grant No.3222037)+1 种基金the CNPC 14th Five-Year Perspective Fundamental Research Project(Grant No.2021DJ2104)the Science Foundation of China University of Petroleum-Beijing(No.2462021YXZZ010).
文摘Recent advances in deep learning have expanded new possibilities for fluid flow simulation in petroleum reservoirs.However,the predominant approach in existing research is to train neural networks using high-fidelity numerical simulation data.This presents a significant challenge because the sole source of authentic wellbore production data for training is sparse.In response to this challenge,this work introduces a novel architecture called physics-informed neural network based on domain decomposition(PINN-DD),aiming to effectively utilize the sparse production data of wells for reservoir simulation with large-scale systems.To harness the capabilities of physics-informed neural networks(PINNs)in handling small-scale spatial-temporal domain while addressing the challenges of large-scale systems with sparse labeled data,the computational domain is divided into two distinct sub-domains:the well-containing and the well-free sub-domain.Moreover,the two sub-domains and the interface are rigorously constrained by the governing equations,data matching,and boundary conditions.The accuracy of the proposed method is evaluated on two problems,and its performance is compared against state-of-the-art PINNs through numerical analysis as a benchmark.The results demonstrate the superiority of PINN-DD in handling large-scale reservoir simulation with limited data and show its potential to outperform conventional PINNs in such scenarios.
文摘We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge.
基金supported by National Basic Research Program of China (Grant Nos. 2010CB950300, 2010CB950500)Public Science and Technology Research Funds Projects of Ocean (Grant No. 201105019)+1 种基金Key Supercomputing Science-Technology Project of Shandong Province of China (Grant No. 2011YD01107)Scientific Research Foundation of the First Institute of Oceanography, State Oceanic Administration (Grant No. GY02-2010G22)
文摘To achieve high parallel efficiency for the global MASNUM surface wave model, the algorithm of an irregular quasirectangular domain decomposition and related serializing of calculating points and data exchanging schemes are developed and conducted, based on the environment of Message Passing Interface(MPI). The new parallel version of the surface wave model is tested for parallel computing on the platform of the Sunway BlueLight supercomputer in the National Supercomputing Center in Jinan. The testing involves four horizontal resolutions, which are 1°×1°,(1/2)°×(1/2)°,(1/4)°×(1/4)°, and(1/8)°×(1/8)°. These tests are performed without data Input/Output(IO) and the maximum amount of processors used in these tests reaches to 131072. The testing results show that the computing speeds of the model with different resolutions are all increased with the increasing of numbers of processors. When the number of processors is four times that of the base processor number, the parallel efficiencies of all resolutions are greater than 80%. When the number of processors is eight times that of the base processor number, the parallel efficiency of tests with resolutions of 1°×1°,(1/2)°×(1/2)° and(1/4)°×(1/4)° is greater than 80%, and it is 62% for the test with a resolution of(1/8)°×(1/8)° using 131072 processors, which is the nearly all processors of Sunway BlueLight. When the processor's number is 24 times that of the base processor number, the parallel efficiencies for tests with resolutions of 1°×1°,(1/2)°×(1/2)°, and(1/4)°×(1/4)° are 72%, 62%, and 38%, respectively. The speedup and parallel efficiency indicate that the irregular quasi-rectangular domain decomposition and serialization schemes lead to high parallel efficiency and good scalability for a global numerical wave model.
文摘The near-wall domain decomposition method(NDD)has proved to be very efficient for modeling near-wall fully turbulent flows.In this paper the NDD is extended to non-equilibrium regimeswith laminar-turbulent transition(LTT)for the first time.The LTT is identified with the use of the e^(N)-method which is applied to both incompressible and compressible flows.TheNDD ismodified to take into account LTT in an efficientway.In addition,implementation of the intermittency expands the capabilities of NDD to model non-equilibrium turbulent flows with transition.Performance of the modified NDD approach is demonstrated on various test problems of subsonic and supersonic flows past a flat plate,a supersonic flow over a compression corner and a planar shock wave impinging on a turbulent boundary layer.The results of modeling with and without decomposition are compared in terms of wall friction and show good agreement with each other while NDD significantly reducing computational resources needed.It turns out that the NDD can reduce the computational time as much as three times while retaining practically the same accuracy of prediction.
文摘In this paper the parallel computing of a grid-point nine-level atmospheric general circulation model on the Dawn 1000 is introduced. The model was developed by the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS). The Dawn 1000 is a MIMD massive parallel computer made by National Research Center for Intelligent Computer (NCIC), CAS. A two-dimensional domain decomposition method is adopted to perform the parallel computing. The potential ways to increase the speed-up ratio and exploit more resources of future massively parallel supercomputation are also discussed.
基金This work was supported in part by Hong Kong RGC DAG93/94 SC10, Competitive Earmarked ResearchGrant HKUST593/94E and the speci
文摘Partition property plays a central role in domain decomposition methods. Existing theory essentially assumes certain ellipticity. We prove the partition property for prod lems without ellipticity which are of practical importance. Example applications include implicit schemes applied to degenerate parabolic partial differential equations arising from superconductors, superfluids and liquid crystals. With this partition property, Schwarz algorithms can be applied to general non-elliptic problems with an h-independent optimal convergence rate. Application to the time-dependent Ginzburg-Landau model of superconductivity is illustrated and numerical results are presented.
基金supported by the National Natural Science Foundation of China(Grant No.51674169)Department of Education of Hebei Province of China(Grant No.ZD2019140)+1 种基金Natural Science Foundation of Hebei Province of China(Grant No.F2019210243)S&T Program of Hebei(Grant No.22375413D)School of Electrical and Electronics Engineering。
文摘Accurate displacement prediction is critical for the early warning of landslides.The complexity of the coupling relationship between multiple influencing factors and displacement makes the accurate prediction of displacement difficult.Moreover,in engineering practice,insufficient monitoring data limit the performance of prediction models.To alleviate this problem,a displacement prediction method based on multisource domain transfer learning,which helps accurately predict data in the target domain through the knowledge of one or more source domains,is proposed.First,an optimized variational mode decomposition model based on the minimum sample entropy is used to decompose the cumulative displacement into the trend,periodic,and stochastic components.The trend component is predicted by an autoregressive model,and the periodic component is predicted by the long short-term memory.For the stochastic component,because it is affected by uncertainties,it is predicted by a combination of a Wasserstein generative adversarial network and multisource domain transfer learning for improved prediction accuracy.Considering a real mine slope as a case study,the proposed prediction method was validated.Therefore,this study provides new insights that can be applied to scenarios lacking sample data.
文摘In the Windows XP 64 bit operating system environment, several common PC were used to build a cluster system, establishing the distributed memory parallel (DMP) computing system. A finite element model of whole aircraft with about 260 million degrees of freedom (DOF) was developed using three-node and four-node thin shell element and two-node beam element. With the large commercial finite element software MSC.MARC and employing two kinds of domain decomposition method (DDM) respectively, realized the parallel solving for the static strength analysis of the whole aircraft model, which offered a high cost-effective solution for solving large-scale and complex finite element models.
基金the financial support by the Federal Ministry for Economic Affairs and Climate Action(BMWK),promotional reference 03EN1066A and 03EN3060Dfunding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No.101023666.
文摘The building sector significantly contributes to climate change.To improve its carbon footprint,applications like model predictive control and predictive maintenance rely on system models.However,the high modeling effort hinders practical application.Machine learning models can significantly reduce this modeling effort.To ensure a machine learning model’s reliability in all operating states,it is essential to know its validity domain.Operating states outside the validity domain might lead to extrapolation,resulting in unpredictable behavior.This paper addresses the challenge of identifying extrapolation in data-driven building energy system models and aims to raise knowledge about it.For that,a novel approach is proposed that calibrates novelty detection algorithms towards the machine learning model.Suitable novelty detection algorithms are identified through a literature review and a benchmark test with 15 candidates.A subset of five algorithms is then evaluated on building energy systems.First,on two-dimensional data,displaying the results with a novel visualization scheme.Then on more complex multi-dimensional use cases.The methodology performs well,and the validity domain could be approximated.The visualization allows for a profound analysis and an improved understanding of the fundamental effects behind a machine learning model’s validity domain and the extrapolation regimes.
文摘When carrying out calculations for turbulent flow simulation,one inevitably has to face the choice between accuracy and speed of calculations.In order to simultaneously obtain both a computationally efficient and more accurate model,a surrogate model can be built on the basis of some fast special model and knowledge of previous calculations obtained by more accurate base models from various test bases or some results of serial calculations.The objective of this work is to construct a surrogate model which allows to improve the accuracy of turbulent calculations obtained by a special model on unstructured meshes.For this purpose,we use 1D Convolutional Neural Network(CNN)of the encoder-decoder architecture and reduce the problem to a single dimension by applying space-filling curves.Such an approach would have the benefit of being applicable to solutions obtained on unstructured meshes.In this work,a non-local approach is applied where entire flow fields obtained by the special and base models are used as input and ground truth output respectively.Spalart-Allmaras(SA)model and Near-wall Domain Decomposition(NDD)method for SA are taken as the base and special models respectively.The efficiency and accuracy of the obtained surrogate model are demonstrated in a case of supersonic flow over a compression corner with different values for angleαand Reynolds number Re.We conducted an investigation into interpolation and extrapolation by Re and also into interpolation byα.
基金part supported by the NSF Grants DMS-1912654 and DMS 2205590。
文摘We provide a concise review of the exponentially convergent multiscale finite element method(ExpMsFEM)for efficient model reduction of PDEs in heterogeneous media without scale separation and in high-frequency wave propagation.The ExpMsFEM is built on the non-overlapped domain decomposition in the classical MsFEM while enriching the approximation space systematically to achieve a nearly exponential convergence rate regarding the number of basis functions.Unlike most generalizations of the MsFEM in the literature,the ExpMsFEM does not rely on any partition of unity functions.In general,it is necessary to use function representations dependent on the right-hand side to break the algebraic Kolmogorov n-width barrier to achieve exponential convergence.Indeed,there are online and offline parts in the function representation provided by the ExpMsFEM.The online part depends on the right-hand side locally and can be computed in parallel efficiently.The offline part contains basis functions that are used in the Galerkin method to assemble the stiffness matrix;they are all independent of the right-hand side,so the stiffness matrix can be used repeatedly in multi-query scenarios.
基金supported in part by National Natural Science Foundation of China(grants#12131002 and#12288101)Huawei Technologies Co.,Ltd.
文摘A new penalty-free neural network method,PFNN-2,is presented for solving partial differential equations,which is a subsequent improvement of our previously proposed PFNN method[1].PFNN-2 inherits all advantages of PFNN in handling the smoothness constraints and essential boundary conditions of self-adjoint problems with complex geometries,and extends the application to a broader range of non-self-adjoint time-dependent differential equations.In addition,PFNN-2 introduces an overlapping domain decomposition strategy to substantially improve the training efficiency without sacrificing accuracy.Experiments results on a series of partial differential equations are reported,which demonstrate that PFNN-2 can outperform state-of-the-art neural network methods in various aspects such as numerical accuracy,convergence speed,and parallel scalability.