Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transfo...Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.展开更多
A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done ...A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.展开更多
Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is requ...Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid（irregular grid）.To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces.展开更多
This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild slope equation in the orthogonal coordinate system. Two types of coordinate systems are stud...This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the conformal system. The computational results show that the conformal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical conformal grid system.展开更多
随着分布式电源(distributed generation,DG)的容量变化,微电网原有的供电结构发生改变,使得潮流大小、方向和功率结构发生变化,对快速检测和定位微电网中的短路故障区域提出了挑战。在MATLAB/Simulink中搭建低压交流微电网模型;通过高...随着分布式电源(distributed generation,DG)的容量变化,微电网原有的供电结构发生改变,使得潮流大小、方向和功率结构发生变化,对快速检测和定位微电网中的短路故障区域提出了挑战。在MATLAB/Simulink中搭建低压交流微电网模型;通过高尺度小波能量谱算法对微电网与大电网公共连接点(point of common coupling,PCC)处检测到的电流进行分解,提取适应不同容量情况的短路故障特征值,实现了不同容量下微电网短路故障的早期检测;利用小波能量谱特征结合基于正交最小二乘法(orthogonal least square,OLS)的径向基函数(radial basis function,RBF)神经网络算法提出一种适用于不同容量微电网的短路故障区域定位方法,并进行仿真验证;在此基础上设计并网模式微电网短路故障保护硬件系统,并进行实验验证。结果表明,所设计的保护系统能够快速、准确地同时实现并网模式下交流微电网短路故障的早期检测与区域定位。展开更多
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing t...A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.展开更多
为了提高极化敏感阵列中压缩感知类波达方向(Direction Of Arrival,DOA)估计算法的精度,避免网格失配问题,本文使用正交偶极子阵列在原子范数最小化(Atomic Norm Minimization,ANM)的理论基础上提出一种无网格波达方向估计算法.首先,将...为了提高极化敏感阵列中压缩感知类波达方向(Direction Of Arrival,DOA)估计算法的精度,避免网格失配问题,本文使用正交偶极子阵列在原子范数最小化(Atomic Norm Minimization,ANM)的理论基础上提出一种无网格波达方向估计算法.首先,将一维正交偶极子天线接收到的多快拍信号分解为两个子阵再求和,然后通过解决半正定规划问题恢复出一个含有入射信源信息的半正定Toeplitz矩阵,继而对该矩阵进行Vandermonde分解,恢复入射信源的DOA信息.同时结合协方差矩阵的向量化结果和最小二乘法计算得到入射信源的极化辅助角和极化相位角信息.通过仿真实验,在不同快拍数和信噪比下,对比子空间类算法和压缩感知类算法,证明了该算法具有较高的测角精度.展开更多
基金supported by the National Nature Science Foundation of China(Nos.41504102 and 41604037)National Science and Technology Major Project(No.2016ZX05015-006)Yangtze University Youth Found(No.2015cqn32)
文摘Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.
基金financially supported by the fund of the State Key Laboratory of Solidif ication Processing in NWPU(No:SKLSP201006)the fund of the National Basic Research Program of China(No:2011CB610402)
文摘A type of mesh called a body-fi tted Cartesian mesh,very different from the traditional structured body-fi tted mesh,is established.At f irst,the right parallelepiped mesh is generated,then,a feature analysis is done on the cross sections.These cross sections are the intersections of the casting shape with the right parallelepiped grids(under the Cartesian coordinate system).On the basis of the feature analysis,two sorts of body-f itted boundary grids,shape-keeping grids and shape-distortion grids,are def ined.Shape-distortion grids can be removed or weaken by increasing the number of grids or moving the coordinates of the mesh generation region,so actually the body-fitted Cartesian mesh generation is to get shape-keeping grids.A shape-keeping grid mainly consists of two sorts of surfaces(I type face and II type face),and each of them is joined by two types of points(I type point and II type point).If only these two types of points were given,the shape-keeping mesh would be constructed.In this paper,the cases of the above two boundary grids being generated were discussed.An algorithm was put forward to get the shape-keeping grids.Several body-fi tted Cartesian meshes generated on castings show the validity of the algorithm.The mesh generation examples show that the body-fi tted Cartesian mesh is more excellent than the right parallelepiped mesh in aspects of decreasing grids number and being closer to the shape of the casting solid.
基金supported in part by National Basic Research Program of China (No. 2014CB845900)Hubei Provincial Natural Science Foundation of China (No. 2014CFA005)
文摘Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid（irregular grid）.To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces.
文摘This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the conformal system. The computational results show that the conformal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical conformal grid system.
文摘随着分布式电源(distributed generation,DG)的容量变化,微电网原有的供电结构发生改变,使得潮流大小、方向和功率结构发生变化,对快速检测和定位微电网中的短路故障区域提出了挑战。在MATLAB/Simulink中搭建低压交流微电网模型;通过高尺度小波能量谱算法对微电网与大电网公共连接点(point of common coupling,PCC)处检测到的电流进行分解,提取适应不同容量情况的短路故障特征值,实现了不同容量下微电网短路故障的早期检测;利用小波能量谱特征结合基于正交最小二乘法(orthogonal least square,OLS)的径向基函数(radial basis function,RBF)神经网络算法提出一种适用于不同容量微电网的短路故障区域定位方法,并进行仿真验证;在此基础上设计并网模式微电网短路故障保护硬件系统,并进行实验验证。结果表明,所设计的保护系统能够快速、准确地同时实现并网模式下交流微电网短路故障的早期检测与区域定位。
基金financially supported by the Science and Technology Project of the Ministry of Transport (Grant No. 2011329224170)
文摘A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.
文摘为了提高极化敏感阵列中压缩感知类波达方向(Direction Of Arrival,DOA)估计算法的精度,避免网格失配问题,本文使用正交偶极子阵列在原子范数最小化(Atomic Norm Minimization,ANM)的理论基础上提出一种无网格波达方向估计算法.首先,将一维正交偶极子天线接收到的多快拍信号分解为两个子阵再求和,然后通过解决半正定规划问题恢复出一个含有入射信源信息的半正定Toeplitz矩阵,继而对该矩阵进行Vandermonde分解,恢复入射信源的DOA信息.同时结合协方差矩阵的向量化结果和最小二乘法计算得到入射信源的极化辅助角和极化相位角信息.通过仿真实验,在不同快拍数和信噪比下,对比子空间类算法和压缩感知类算法,证明了该算法具有较高的测角精度.
