Continental subduction and collision normally follows oceanic subduction,with the remarkable event of formation and exhumation of high-to ultra-high-pressure(HP-UHP)metamorphic rocks.Based on the summary of numerical ...Continental subduction and collision normally follows oceanic subduction,with the remarkable event of formation and exhumation of high-to ultra-high-pressure(HP-UHP)metamorphic rocks.Based on the summary of numerical geodynamic models,six modes of continental convergence have been identified:pure shear thickening,folding and buckling,one-sided steep subduction,flat subduction,two-sided subduction,and subducting slab break-off.In addition,the exhumation of HP-UHP rocks can be formulated into eight modes:thrust fault exhumation,buckling exhumation,material circulation,overpressure model,exhumation of a coherent crustal slice,episodic ductile extrusion,slab break-off induced eduction,and exhumation through fractured overriding lithosphere.During the transition from subduction to exhumation,the weakening and detachment of subducted continental crust are prerequisites.However,the dominant weakening mechanisms and their roles in the subduction channel are poorly constrained.To a first degree approximation,the mechanism of continental subduction and exhumation can be treated as a subduction channel flow model,which incorporates the competing effects of downward Couette(subduction)flow and upward Poiseuille(exhumation)flow in the subduction channel.However,the(de-)hydration effect plays significant roles in the deformation of subduction channel and overriding lithosphere,which thereby result in very different modes from the simple subduction channel flow.Three-dimensionality is another important issue with highlighting the along-strike differential modes of continental subduction,collision and exhumation in the same continental convergence belt.展开更多
All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equ...All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equations. The current computational fluid dynamics （CFD） focuses on the numerical solution of partial differential equations （PDEs）, and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numer- ical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require fur- ther expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional dis- tinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of con- structing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm develop- ment. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be mod- eled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct construction of dis- crete numeric展开更多
In this paper, we develop a novel mathematical model to estimate the probability distribution function of the number of tags discovered after a certain number of interrogation rounds. In addition, the pdfs of the numb...In this paper, we develop a novel mathematical model to estimate the probability distribution function of the number of tags discovered after a certain number of interrogation rounds. In addition, the pdfs of the number of rounds needed to discover all the tags are also calculated. The estimation of such pdfs will be helpful in estimating the number of interrogation rounds and the optimal parameter configuration of the RFID system which in turn will be helpful in estimating the time needed to discover all tags. Our results show that the proposed model accurately predicts the tags detection probability. We then use the proposed model to optimally configure the reader parameters (i.e. the frame size and the number of interrogation rounds).展开更多
2D thermo-mechanical models are constructed to investigate the effects of oceanic and continental crustal eclogitization on plate dynamics at three successive stages of oceanic subduction, slab breakoff, and continent...2D thermo-mechanical models are constructed to investigate the effects of oceanic and continental crustal eclogitization on plate dynamics at three successive stages of oceanic subduction, slab breakoff, and continental subduction. Crustal eclogitization directly increases the average slab density and accordingly the slab pull force, which makes the slab subduct deeply and steeply. Numerical results demonstrate that the duration time from initial continental collision to slab breakoff largely depends on the slab pull force. Specifically, eclogitization of subducted crust can greatly decrease the duration time, but increase the breakoff depth. The detachment of oceanic slab from the pro-continental lithosphere is accompanied with obvious exhumation of the subducted continental crust and a sharp uplift of the collision zone in response to the disappearance of downward drag force and the induced asthenospheric upwelling, especially under the condition of no or incomplete crustal eclogitization. During continental subduction, the slab dip angle is strongly correlated with eclogitization of subducted continental lower crust, which regulates the slab buoyancy nature. Our model results can provide several important implications for the Himalayan-Tibetan collision zone. For example, it is possible that the lateral variations in the degree of eclogitization of the subducted Indian crust might to some extent contribute to the lateral variations of subduction angle along the Himalayan orogenic belt. Moreover, the accumulation of highly radiogenic sediments and upper continental crustal materials at the active margin in combination with the strong shear heating due to continuous continental subduction together cause rising of isotherms in the accretionary wedge, which facilitate the development of crustal partial melting and metamorphism.展开更多
基金supported by the Start-up Research Fund from Institute of Geology of CAGS(Grant No.J1219)National Natural Science Foundation of China(Grant No.40921001)China Geological Survey Projects(Grant Nos.1212011121275,1212011120161)and Sinoprobe Project
文摘Continental subduction and collision normally follows oceanic subduction,with the remarkable event of formation and exhumation of high-to ultra-high-pressure(HP-UHP)metamorphic rocks.Based on the summary of numerical geodynamic models,six modes of continental convergence have been identified:pure shear thickening,folding and buckling,one-sided steep subduction,flat subduction,two-sided subduction,and subducting slab break-off.In addition,the exhumation of HP-UHP rocks can be formulated into eight modes:thrust fault exhumation,buckling exhumation,material circulation,overpressure model,exhumation of a coherent crustal slice,episodic ductile extrusion,slab break-off induced eduction,and exhumation through fractured overriding lithosphere.During the transition from subduction to exhumation,the weakening and detachment of subducted continental crust are prerequisites.However,the dominant weakening mechanisms and their roles in the subduction channel are poorly constrained.To a first degree approximation,the mechanism of continental subduction and exhumation can be treated as a subduction channel flow model,which incorporates the competing effects of downward Couette(subduction)flow and upward Poiseuille(exhumation)flow in the subduction channel.However,the(de-)hydration effect plays significant roles in the deformation of subduction channel and overriding lithosphere,which thereby result in very different modes from the simple subduction channel flow.Three-dimensionality is another important issue with highlighting the along-strike differential modes of continental subduction,collision and exhumation in the same continental convergence belt.
基金supported by Hong Kong Research Grant Council (Grants 621011,620813 and 16211014)HKUST(IRS15SC29 and SBI14SC11)
文摘All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equations. The current computational fluid dynamics （CFD） focuses on the numerical solution of partial differential equations （PDEs）, and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numer- ical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require fur- ther expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional dis- tinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of con- structing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm develop- ment. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be mod- eled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct construction of dis- crete numeric
文摘In this paper, we develop a novel mathematical model to estimate the probability distribution function of the number of tags discovered after a certain number of interrogation rounds. In addition, the pdfs of the number of rounds needed to discover all the tags are also calculated. The estimation of such pdfs will be helpful in estimating the number of interrogation rounds and the optimal parameter configuration of the RFID system which in turn will be helpful in estimating the time needed to discover all tags. Our results show that the proposed model accurately predicts the tags detection probability. We then use the proposed model to optimally configure the reader parameters (i.e. the frame size and the number of interrogation rounds).
基金financial supports by the National Natural Science Foundation of China(Nos.41490613,41190073 and 41304071)the National Basic Research Program of China(Nos.2014CB440901 and 2015CB856106)
文摘2D thermo-mechanical models are constructed to investigate the effects of oceanic and continental crustal eclogitization on plate dynamics at three successive stages of oceanic subduction, slab breakoff, and continental subduction. Crustal eclogitization directly increases the average slab density and accordingly the slab pull force, which makes the slab subduct deeply and steeply. Numerical results demonstrate that the duration time from initial continental collision to slab breakoff largely depends on the slab pull force. Specifically, eclogitization of subducted crust can greatly decrease the duration time, but increase the breakoff depth. The detachment of oceanic slab from the pro-continental lithosphere is accompanied with obvious exhumation of the subducted continental crust and a sharp uplift of the collision zone in response to the disappearance of downward drag force and the induced asthenospheric upwelling, especially under the condition of no or incomplete crustal eclogitization. During continental subduction, the slab dip angle is strongly correlated with eclogitization of subducted continental lower crust, which regulates the slab buoyancy nature. Our model results can provide several important implications for the Himalayan-Tibetan collision zone. For example, it is possible that the lateral variations in the degree of eclogitization of the subducted Indian crust might to some extent contribute to the lateral variations of subduction angle along the Himalayan orogenic belt. Moreover, the accumulation of highly radiogenic sediments and upper continental crustal materials at the active margin in combination with the strong shear heating due to continuous continental subduction together cause rising of isotherms in the accretionary wedge, which facilitate the development of crustal partial melting and metamorphism.
