The data obtained from a high resolution seismic refraction profile, which was carded out in Jiashi, Xinjiang, strong earthquake swarm area, were processed with both finite difference inversion and Hagedoorn refractor...The data obtained from a high resolution seismic refraction profile, which was carded out in Jiashi, Xinjiang, strong earthquake swarm area, were processed with both finite difference inversion and Hagedoorn refractor wavefront imaging technique and the fine upper crustal structure was determined. The results show that the upper crustal structure is relatively well-distributed in laterally and obviously by layers vertically.From surface to 11.0 km depth, there are about four layers. The P wave velocity of top two layers range from 1.65 to 4.5 km/s and their bottom boundaries, the buried depths of which are 0.4, 2.96-3.0 km respectively, are almost horizontal; The third layer is comparatively complicated and its P wave velocity presents inhomogeneous in both laterally and vertically. The bottom boundary of third layer is crystalline basement and shows a little uplift, which seemly suggest that the upper crust had been resisted while the hard Tarim block inserting into Tianshan Mountain; The forth layer is relatively even and its P wave velocity is about 6.3 km/s. There are a lateral velocity variation at the depth of about 4.0 km, and suggest that it has something to do with the hidden Meigaiti fault and Meigaiti-Xiasuhong fault but there are no the structure features about these faults stretching to the surface and passing through the crystalline basement. The seismogenic tectonic of Jiashi strong earthquake swarm at least lies in middle or lower crust beneath 11.0 km depth.展开更多
This paper presents a class of high resolution local time step schemes for nonlinear hyperbolic conservation laws and the closely related convection-diffusion equations, by projecting the solution increments of the un...This paper presents a class of high resolution local time step schemes for nonlinear hyperbolic conservation laws and the closely related convection-diffusion equations, by projecting the solution increments of the underlying partial differential equations (PDE) at each local time step. The main advantages are that they are of good consistency, and it is convenient to implement them. The schemes are L^∞ stable, satisfy a cell entropy inequality, and may be extended to the initial boundary value problem of general unsteady PDEs with higher-order spatial derivatives. The high resolution schemes are given by combining the reconstruction technique with a second order TVD Runge-Kutta scheme or a Lax-Wendroff type method, respectively. The schemes are used to solve a linear convection-diffusion equation, the nonlinear inviscid Burgers' equation, the one- and two-dimensional compressible Euler equations, and the two-dimensional incompressible Navier-Stokes equations. The numerical results show that the schemes are of higher-order accuracy, and efficient in saving computational cost, especially, for the case of combining the present schemes with the adaptive mesh method [15]. The correct locations of the slow moving or stronger discontinuities are also obtained, although the schemes are slightly nonconservative.展开更多
This paper is concerned with a new version of the Osher-Solomon Riemann solver and is based on a numerical integration of the path-dependent dissipation matrix.The resulting scheme is much simpler than the original on...This paper is concerned with a new version of the Osher-Solomon Riemann solver and is based on a numerical integration of the path-dependent dissipation matrix.The resulting scheme is much simpler than the original one and is applicable to general hyperbolic conservation laws,while retaining the attractive features of the original solver:the method is entropy-satisfying,differentiable and complete in the sense that it attributes a different numerical viscosity to each characteristic field,in particular to the intermediate ones,since the full eigenstructure of the underlying hyperbolic system is used.To illustrate the potential of the proposed scheme we show applications to the following hyperbolic conservation laws:Euler equations of compressible gasdynamics with ideal gas and real gas equation of state,classical and relativistic MHD equations as well as the equations of nonlinear elasticity.To the knowledge of the authors,apart from the Euler equations with ideal gas,an Osher-type scheme has never been devised before for any of these complicated PDE systems.Since our new general Riemann solver can be directly used as a building block of high order finite volume and discontinuous Galerkin schemes we also show the extension to higher order of accuracy and multiple space dimensions in the new framework of PNPM schemes on unstructured meshes recently proposed in[9].展开更多
基金National Natural Science Foundation of China (40334040) and Joint Seismological Foundation (106076).
文摘The data obtained from a high resolution seismic refraction profile, which was carded out in Jiashi, Xinjiang, strong earthquake swarm area, were processed with both finite difference inversion and Hagedoorn refractor wavefront imaging technique and the fine upper crustal structure was determined. The results show that the upper crustal structure is relatively well-distributed in laterally and obviously by layers vertically.From surface to 11.0 km depth, there are about four layers. The P wave velocity of top two layers range from 1.65 to 4.5 km/s and their bottom boundaries, the buried depths of which are 0.4, 2.96-3.0 km respectively, are almost horizontal; The third layer is comparatively complicated and its P wave velocity presents inhomogeneous in both laterally and vertically. The bottom boundary of third layer is crystalline basement and shows a little uplift, which seemly suggest that the upper crust had been resisted while the hard Tarim block inserting into Tianshan Mountain; The forth layer is relatively even and its P wave velocity is about 6.3 km/s. There are a lateral velocity variation at the depth of about 4.0 km, and suggest that it has something to do with the hidden Meigaiti fault and Meigaiti-Xiasuhong fault but there are no the structure features about these faults stretching to the surface and passing through the crystalline basement. The seismogenic tectonic of Jiashi strong earthquake swarm at least lies in middle or lower crust beneath 11.0 km depth.
基金This research was partially sponsored by the National Basic Research Program under the Grant 2005CB321703, National Natural Science Foundation of China (No. 10431050, 10576001), SRF for R0CS, SEM, the Alexander von Humboldt foundation, and the Deutsche Forschungsgemeinschaft (DFG Wa 633/10-3).Acknowledgments. The authors thank Professor Tao Tang for numerous discussions during the preparation of this work, and also thank the referees for many helpful suggestions.
文摘This paper presents a class of high resolution local time step schemes for nonlinear hyperbolic conservation laws and the closely related convection-diffusion equations, by projecting the solution increments of the underlying partial differential equations (PDE) at each local time step. The main advantages are that they are of good consistency, and it is convenient to implement them. The schemes are L^∞ stable, satisfy a cell entropy inequality, and may be extended to the initial boundary value problem of general unsteady PDEs with higher-order spatial derivatives. The high resolution schemes are given by combining the reconstruction technique with a second order TVD Runge-Kutta scheme or a Lax-Wendroff type method, respectively. The schemes are used to solve a linear convection-diffusion equation, the nonlinear inviscid Burgers' equation, the one- and two-dimensional compressible Euler equations, and the two-dimensional incompressible Navier-Stokes equations. The numerical results show that the schemes are of higher-order accuracy, and efficient in saving computational cost, especially, for the case of combining the present schemes with the adaptive mesh method [15]. The correct locations of the slow moving or stronger discontinuities are also obtained, although the schemes are slightly nonconservative.
基金financed by the Italian Ministry of Research(MIUR)under the project PRIN 2007 and by MIUR and the British Council under the project British-Italian Partnership Programme for young researchers 2008-2009。
文摘This paper is concerned with a new version of the Osher-Solomon Riemann solver and is based on a numerical integration of the path-dependent dissipation matrix.The resulting scheme is much simpler than the original one and is applicable to general hyperbolic conservation laws,while retaining the attractive features of the original solver:the method is entropy-satisfying,differentiable and complete in the sense that it attributes a different numerical viscosity to each characteristic field,in particular to the intermediate ones,since the full eigenstructure of the underlying hyperbolic system is used.To illustrate the potential of the proposed scheme we show applications to the following hyperbolic conservation laws:Euler equations of compressible gasdynamics with ideal gas and real gas equation of state,classical and relativistic MHD equations as well as the equations of nonlinear elasticity.To the knowledge of the authors,apart from the Euler equations with ideal gas,an Osher-type scheme has never been devised before for any of these complicated PDE systems.Since our new general Riemann solver can be directly used as a building block of high order finite volume and discontinuous Galerkin schemes we also show the extension to higher order of accuracy and multiple space dimensions in the new framework of PNPM schemes on unstructured meshes recently proposed in[9].