The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Lamer. Since that many numerical methods on the subject have been developed, such as pure numerical techniques, ray m...The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Lamer. Since that many numerical methods on the subject have been developed, such as pure numerical techniques, ray method and boundary method. The boundary method based on boundary integral equation is a semi-analytical method which is suitable to modeling wave field induced by irregular border. According to the property of the applied Green's function the boundary methods can be sorted into space domain boundary method and wavenumber domain boundary method. For both of them it is necessary to solve a large equation, which means much computation is needed. Thus, it is difficult for the boundary methods to be applied in simulating wave field with high frequency or in large range. To develop a new method with less computation is meaningful. For this purpose, localized boundary integral equation, i.e., discrete wavenumber method is proposed. It is rooted in the Bouchon-Campillo method, an important wavenumber domain boundary method. Firstly the force on interface is separated into two parts: one is on flat part and the other on irregular part of the interface. Then Fourier transform is applied to identify their relation, the unknown distributes only on irregular part. Consequently computation efficiency is dramatically improved. Importantly its accuracy is the same as that of Bouchon-Campillo.展开更多
A parallel algorithm of circulation numerical model based on message passing interface(MPI) is developed using serialization and an irregular rectangle decomposition scheme. Neighboring point exchange strategy(NPES...A parallel algorithm of circulation numerical model based on message passing interface(MPI) is developed using serialization and an irregular rectangle decomposition scheme. Neighboring point exchange strategy(NPES) is adopted to further enhance the computational efficiency. Two experiments are conducted on HP C7000 Blade System, the numerical results show that the parallel version with NPES(PVN) produces higher efficiency than the original parallel version(PV). The PVN achieves parallel efficiency in excess of 0.9 in the second experiment when the number of processors increases to 100, while the efficiency of PV decreases to 0.39 rapidly. The PVN of ocean circulation model is used in a fine-resolution regional simulation, which produces better results. The capability of universal implementation of this algorithm makes it applicable in many other ocean models potentially.展开更多
基金supported by National Natural Science Foundation of China (Nos.40874027,90715020 and 90915012)IGPCEA(DQJB07B06)Special Public Welfare Industry (Nos.20070804 and 200808008)
文摘The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Lamer. Since that many numerical methods on the subject have been developed, such as pure numerical techniques, ray method and boundary method. The boundary method based on boundary integral equation is a semi-analytical method which is suitable to modeling wave field induced by irregular border. According to the property of the applied Green's function the boundary methods can be sorted into space domain boundary method and wavenumber domain boundary method. For both of them it is necessary to solve a large equation, which means much computation is needed. Thus, it is difficult for the boundary methods to be applied in simulating wave field with high frequency or in large range. To develop a new method with less computation is meaningful. For this purpose, localized boundary integral equation, i.e., discrete wavenumber method is proposed. It is rooted in the Bouchon-Campillo method, an important wavenumber domain boundary method. Firstly the force on interface is separated into two parts: one is on flat part and the other on irregular part of the interface. Then Fourier transform is applied to identify their relation, the unknown distributes only on irregular part. Consequently computation efficiency is dramatically improved. Importantly its accuracy is the same as that of Bouchon-Campillo.
基金The National High Technology Research and Development Program(863 Program)of China under contract No.2013AA09A505
文摘A parallel algorithm of circulation numerical model based on message passing interface(MPI) is developed using serialization and an irregular rectangle decomposition scheme. Neighboring point exchange strategy(NPES) is adopted to further enhance the computational efficiency. Two experiments are conducted on HP C7000 Blade System, the numerical results show that the parallel version with NPES(PVN) produces higher efficiency than the original parallel version(PV). The PVN achieves parallel efficiency in excess of 0.9 in the second experiment when the number of processors increases to 100, while the efficiency of PV decreases to 0.39 rapidly. The PVN of ocean circulation model is used in a fine-resolution regional simulation, which produces better results. The capability of universal implementation of this algorithm makes it applicable in many other ocean models potentially.