构造热演化的数值模拟方法

METHOD OF NUMERICAL SIMULATION ON TECTONICALLY THERMAL EVOLUTION

对模拟构造热演化的二维和三维非稳态热传导-平流微分方程的数值求解进行了研究.在空间和时间上分别应用有限单元法和差分法进行离散化处理;推导了非稳态热传导-平流有限元方程.在求解技术上发展了“迎风”权函数方法,将有限元法、有限差分法和“迎风”权函数法结合起来,使用变网格、变结点方法对地质体的运动介质、运动边界、多热源的非线性、非稳态的构造热演化进行了模拟.并对各种热源体的处理、初始条件和边界条件的处理进行了分析.

Using numerical solutions of 2-D and 3-D nonsteady-state heat advective-conducrion equations to simulate tectonically thermal evolution have been studied in this paper.The finite element method and differential method have been used for discreteness in space and time domains respectively, and the advection finite element equations have been derived.The authors developed an upwind scheme and proposed a method which combined finite element method, differential method, upwind scheme with net-changeable and nod-changeable methods to simulate tectonically thermal evolution, including advection.The paper also disscussed the processes of initial conditions, boundary conditions and various heat sources.