We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high...We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.展开更多
Stress response of a tension leg platform (TLP) in extreme environments was investigated in this paper. A location on one of the gussets was selected as the object point, where directional stresses were numerically ...Stress response of a tension leg platform (TLP) in extreme environments was investigated in this paper. A location on one of the gussets was selected as the object point, where directional stresses were numerically simulated and also experimentally verified by a strain gage. Environmental loading and the platform's structural strength were analyzed in accordance with industrial standards, utilizing linear wave theory and the finite element method (FEM). The fast Fourier transform technique was used to calculate the stress response amplitude operators (RAO) from the records of measurements. A comparison was performed between the stress RAO of the numerical simulation and that of the actual measurements. The results indicated that the stress RAO of the numerical simulation fitted well with measured data at specified wave headings with different periods.展开更多
The discrete element method is used to simulate specimens under three different loading conditions(conventional triaxial compression,plane strain,and direct shear)with different initial conditions to explore the und...The discrete element method is used to simulate specimens under three different loading conditions(conventional triaxial compression,plane strain,and direct shear)with different initial conditions to explore the underlying mechanics of the specimen deformation from a microscale perspective.Deformations of specimens with different initial void ratios at different confining stresses under different loading conditions are studied.Results show that the discrete element models successfully capture the specimen deformation and the strain localization.Particle behaviors including particle rotation and displacement and the mesoscale void ratio distributions are used to explain the strain localization and specimen deformation.It is found that the loading condition is one of the most important factors controlling the specimen deformation mode.Microscale behavior of the granular soil is the driving mechanics of the macroscale deformation of the granular assembly.展开更多
时域有限差分(finite difference time domain,FDTD)方法广泛应用于电磁场仿真领域,并与量子力学理论相结合来求解时域薛定谔方程,然而数值计算中的稳定性研究缺少理论方面的探讨。基于冯.诺依曼的稳定性分析方法得到了时域薛定谔方程...时域有限差分(finite difference time domain,FDTD)方法广泛应用于电磁场仿真领域,并与量子力学理论相结合来求解时域薛定谔方程,然而数值计算中的稳定性研究缺少理论方面的探讨。基于冯.诺依曼的稳定性分析方法得到了时域薛定谔方程的一维以及多维的稳定性条件,并且讨论了在不同势能情况下该稳定性条件的表现形式。数值结果充分证明了结论的正确性。展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
Rainfall is an important factor to trigger the debris flow.Numerical simulation on the responses of slopes and the initiation of debris flow under rainfall was processed by using the software FLAC2D based on the soil ...Rainfall is an important factor to trigger the debris flow.Numerical simulation on the responses of slopes and the initiation of debris flow under rainfall was processed by using the software FLAC2D based on the soil parameters in Weijia Gully,Beichuan County,Sichuan Province,China.The effects of the slope angle,rainfall intensity,soil parameters on the developments of the stress and pore pressure and deformation of the slope were studied.It indicates that large displacements of the slope are mainly located near the slope toe.With the increase of the rainfall intensity the stability of the slope decreases and so the debris-flow is easy to occur.展开更多
Homogeneous mixtures of CH4/air under moderate or intense low-oxygen dilution(MILD) combustion conditions were numerically studied to clarify the fundamental effects of exhaust gas recirculation(EGR),espe-cially C...Homogeneous mixtures of CH4/air under moderate or intense low-oxygen dilution(MILD) combustion conditions were numerically studied to clarify the fundamental effects of exhaust gas recirculation(EGR),espe-cially CO2 in EGR gases,on ignition characteristics.Specifically,effects of CO2 addition on autoignition delay time were emphasized at temperature between 1200 K and 1600 K for a wide range of the lean-to-rich equivalence ratio(0.2~2).The results showed that the ignition delay time increased with equivalence ratio or CO2 dilution ratio.Fur-thermore,ignition delay time was seen to be exponentially related with the reciprocal of initial temperature.Special concern was given to the chemical effects of CO2 on the ignition delay time.The enhancement of ignition delay time with CO2 addition can be mainly ascribed to the decrease of H,O and OH radicals.The predictions of tem-perature profiles and mole fractions of CO and CO2 were strongly related to the chemical effects of CO2.A single ignition time correlation was obtained in form of Arrhenius-type for the entire range of conditions as a function of temperature,CH4 mole fraction and O2 mole fraction.This correlation could successfully capture the complex be-haviors of ignition of CH4/air/CO2 mixture.The results can be applied to MILD combustion as "reference time",for example,to predict ignition delay time in turbulent reacting flow.展开更多
This paper studies the perturbed nonlinear diffusion-convection equation with source term via the approximate generalized conditional symmetry (A GCS). Complete classification of those perturbed equations which admi...This paper studies the perturbed nonlinear diffusion-convection equation with source term via the approximate generalized conditional symmetry (A GCS). Complete classification of those perturbed equations which admit certain types of AGCSs is derived. Some approximate invariant solutions to the resulting equations can also be obtained.展开更多
Damage caused by underground coal mining is a serious problem in mining areas in China; therefore, studying and obtaining the rules of ground movement and deformation under different geological conditions is of great ...Damage caused by underground coal mining is a serious problem in mining areas in China; therefore, studying and obtaining the rules of ground movement and deformation under different geological conditions is of great importance. The numerical software ANSYS was used in this study to simulate mining processes under two special geological conditions: (1) thick unconsolidated soil layer and thin bedrock; (2) thin soil layer and thick bedrock. The rules for ground movement and deformation for different soil layer to bedrock ratios were obtained. On the basis of these rules, a prediction parameter modified model of the influence function was proposed, which is suitable for different values of unconsolidated soil layer thickness. The prediction results were verified using two sets of typical field data.展开更多
According to geological conditions of No. 3 and No. 4 coal seams (namely A3 and B4) of the Pan'er coal mine and the parameters of panels 11223, 11224, and 11124 with fully-mechanical coal mining, we built 2D simila...According to geological conditions of No. 3 and No. 4 coal seams (namely A3 and B4) of the Pan'er coal mine and the parameters of panels 11223, 11224, and 11124 with fully-mechanical coal mining, we built 2D similar material simulation and FLAC3D numerical simulation models to investigate the development of mining-induced stress and the extraction effect of pressure-relief gas with large height and upward mining. Based on a comprehensive analysis of experimental data and observations, we obtained the deformation and breakage characteristics of strata overlying the coal seam, the development patterns of the mining-induced stress and fracture, and the size of the stress-relief area. The stress-relief effect was investigated and analyzed in consideration with mining height and three thick hard strata. Because of the group of three hard thick strata located in the main roof and the residual stress of mined panel 11124, the deformation, breakage, mining-induced stress and fracture development, and the stress-relief coefficient were discontinuous and asymmetrical. The breakage angle of the overlying strata, and the compressive and expansive zones of coal deformation were mainly controlled by the number, thickness, and strength of the hard stratum. Compared with the value of breakage angle derived by the traditional empirical method, the experimental value was lower than the traditional results by 3°-4°below the hard thick strata group, and by 13°-19° above the hard thick strata group. The amount of gas extracted from floor drainage roadway of B4 over 17 months was variable and the amount of gas per month differed considerably, being much smaller when panel 11223 influenced the area of the three hard thick strata. Generally, the stress-relief zone of No. 4 coal seam was small under the influence of the hard thick strata located in the main roof, which played an important role in delaying the breakage time and increasing the breakage space. In this study we gained understanding of the stress-relief mechanism influenced by the hard thick roof. The research results and engineering practice show that the main roof of the multiple hard thick strata is a critical factor in the design of panel layout and roadways for integrated coal exploitation and gas extraction, provides a theoretical basis for safe and high-efficient mining of coal resources.展开更多
In industrial applications involving metal and polymer sheets, the flow situation is strongly unsteady and the sheet temperature is a mixture of prescribed surface temperature and heat flux. Further, a proper choice o...In industrial applications involving metal and polymer sheets, the flow situation is strongly unsteady and the sheet temperature is a mixture of prescribed surface temperature and heat flux. Further, a proper choice of cooling liquid is also an important component of the analysis to achieve better outputs. In this paper, we numerically investigate Darcy-Forchheimer nanoliquid flows past an unsteady stretching surface by incorporating various effects, such as the Brownian and thermophoresis effects, Navier’s slip condition and convective thermal boundary conditions. To solve the governing equations, using suitable similarity transformations, the nonlinear ordinary differential equations are derived and the resulting coupled momentum and energy equations are numerically solved using the spectral relaxation method. Through the systematically numerical investigation, the important physical parameters of the present model are analyzed. We find that the presence of unsteadiness parameter has significant effects on velocity, temperature, concentration fields, the associated heat and mass transport rates. Also, an increase in inertia coefficient and porosity parameter causes an increase in the velocity at the boundary.展开更多
The influence of an upper,mined coal seam on the stability of rock surrounding a roadway in a lower coal seam is examined.The technical problems of roadway control are discussed based on the geological conditions exis...The influence of an upper,mined coal seam on the stability of rock surrounding a roadway in a lower coal seam is examined.The technical problems of roadway control are discussed based on the geological conditions existing in the Liyazhuang Mine No.2 coal seam.The stress distribution and floor failure in the lower works after mining the upper coal is studied through numerical simulations.The failure mechanism of the roof and walls of a roadway located in the lower coal seam is described.The predicted deformation and failure of the roadway for different distances between the two coal seams are used to design two ways of supporting the lower structure.One is a combined support consisting of anchors with a joist steel tent and a combined anchor truss.A field test of the design was performed to good effect.The results have significance for the design of supports for roadways located in similar conditions.展开更多
This paper takes No.52 return uphill roadway of Yangquhe coal mine as a research project. Based on the research, especially its geological condition, indoor experiments, numerical simulation and theoretical analysis w...This paper takes No.52 return uphill roadway of Yangquhe coal mine as a research project. Based on the research, especially its geological condition, indoor experiments, numerical simulation and theoretical analysis were employed to determine the difficult coefficients of Yangquhe project. By using these means,the difficult coefficients of the deep rock engineering were determined. From a study of the effects of crustal stress and the roof mechanism on roadway stability, the transformation mechanism in Yangquhe coal mine has been determined. As a result of this research, the interactive support technology of prestressed cable mesh was developed and the technology tested in mining engineering, which proved to be feasible.展开更多
Containment booms are commonly used in collecting and containing spilled oil on the sea surface and in protecting specific sea areas against oil slick spreading.In the present study,a numerical model is proposed based...Containment booms are commonly used in collecting and containing spilled oil on the sea surface and in protecting specific sea areas against oil slick spreading.In the present study,a numerical model is proposed based on the N-S equations in a mesh frame.The proposed model tracks the outline of the floating boom in motion by using the fractional area/volume obstacle representation technique.The boom motion is then simulated by the technique of general moving object.The simulated results of the rigid oil boom motions are validated against the experimental results.Then,the failure mechanism of the boom is investigated through numerical experiments.Based on the numerical results,the effects of boom parameters and dynamic factors on the oil containment performance are also assessed.展开更多
To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat s...To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.展开更多
基金supported by the National Nature Science Foundation of China(Grant No.U1262208)the Important National Science & Technology Specific Projects(Grant No.2011ZX05019-008)
文摘We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.
基金supported by the Fund of "111 Project" (Grant No.B07019) from the State Administration of Foreign Experts Affairs and the Ministry of Education of China
文摘Stress response of a tension leg platform (TLP) in extreme environments was investigated in this paper. A location on one of the gussets was selected as the object point, where directional stresses were numerically simulated and also experimentally verified by a strain gage. Environmental loading and the platform's structural strength were analyzed in accordance with industrial standards, utilizing linear wave theory and the finite element method (FEM). The fast Fourier transform technique was used to calculate the stress response amplitude operators (RAO) from the records of measurements. A comparison was performed between the stress RAO of the numerical simulation and that of the actual measurements. The results indicated that the stress RAO of the numerical simulation fitted well with measured data at specified wave headings with different periods.
基金The National Natural Science Foundation of China(No.51079030)
文摘The discrete element method is used to simulate specimens under three different loading conditions(conventional triaxial compression,plane strain,and direct shear)with different initial conditions to explore the underlying mechanics of the specimen deformation from a microscale perspective.Deformations of specimens with different initial void ratios at different confining stresses under different loading conditions are studied.Results show that the discrete element models successfully capture the specimen deformation and the strain localization.Particle behaviors including particle rotation and displacement and the mesoscale void ratio distributions are used to explain the strain localization and specimen deformation.It is found that the loading condition is one of the most important factors controlling the specimen deformation mode.Microscale behavior of the granular soil is the driving mechanics of the macroscale deformation of the granular assembly.
文摘时域有限差分(finite difference time domain,FDTD)方法广泛应用于电磁场仿真领域,并与量子力学理论相结合来求解时域薛定谔方程,然而数值计算中的稳定性研究缺少理论方面的探讨。基于冯.诺依曼的稳定性分析方法得到了时域薛定谔方程的一维以及多维的稳定性条件,并且讨论了在不同势能情况下该稳定性条件的表现形式。数值结果充分证明了结论的正确性。
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
文摘Rainfall is an important factor to trigger the debris flow.Numerical simulation on the responses of slopes and the initiation of debris flow under rainfall was processed by using the software FLAC2D based on the soil parameters in Weijia Gully,Beichuan County,Sichuan Province,China.The effects of the slope angle,rainfall intensity,soil parameters on the developments of the stress and pore pressure and deformation of the slope were studied.It indicates that large displacements of the slope are mainly located near the slope toe.With the increase of the rainfall intensity the stability of the slope decreases and so the debris-flow is easy to occur.
基金Supported by the National Natural Science Foundation of China (50206014)the Shuguang Scholar Program of Shanghai Education Development Foundation (05SG23)
文摘Homogeneous mixtures of CH4/air under moderate or intense low-oxygen dilution(MILD) combustion conditions were numerically studied to clarify the fundamental effects of exhaust gas recirculation(EGR),espe-cially CO2 in EGR gases,on ignition characteristics.Specifically,effects of CO2 addition on autoignition delay time were emphasized at temperature between 1200 K and 1600 K for a wide range of the lean-to-rich equivalence ratio(0.2~2).The results showed that the ignition delay time increased with equivalence ratio or CO2 dilution ratio.Fur-thermore,ignition delay time was seen to be exponentially related with the reciprocal of initial temperature.Special concern was given to the chemical effects of CO2 on the ignition delay time.The enhancement of ignition delay time with CO2 addition can be mainly ascribed to the decrease of H,O and OH radicals.The predictions of tem-perature profiles and mole fractions of CO and CO2 were strongly related to the chemical effects of CO2.A single ignition time correlation was obtained in form of Arrhenius-type for the entire range of conditions as a function of temperature,CH4 mole fraction and O2 mole fraction.This correlation could successfully capture the complex be-haviors of ignition of CH4/air/CO2 mixture.The results can be applied to MILD combustion as "reference time",for example,to predict ignition delay time in turbulent reacting flow.
基金The project supported by National Natural Science Foundation of China under Grant Nos.10371098 and 10447007the Natural Science Foundation of Shanxi Province of China under Grant No.2005A13
文摘This paper studies the perturbed nonlinear diffusion-convection equation with source term via the approximate generalized conditional symmetry (A GCS). Complete classification of those perturbed equations which admit certain types of AGCSs is derived. Some approximate invariant solutions to the resulting equations can also be obtained.
文摘Damage caused by underground coal mining is a serious problem in mining areas in China; therefore, studying and obtaining the rules of ground movement and deformation under different geological conditions is of great importance. The numerical software ANSYS was used in this study to simulate mining processes under two special geological conditions: (1) thick unconsolidated soil layer and thin bedrock; (2) thin soil layer and thick bedrock. The rules for ground movement and deformation for different soil layer to bedrock ratios were obtained. On the basis of these rules, a prediction parameter modified model of the influence function was proposed, which is suitable for different values of unconsolidated soil layer thickness. The prediction results were verified using two sets of typical field data.
基金Acknowledgments This work is supported by the National Nature Science Foundation of China (51374011).
文摘According to geological conditions of No. 3 and No. 4 coal seams (namely A3 and B4) of the Pan'er coal mine and the parameters of panels 11223, 11224, and 11124 with fully-mechanical coal mining, we built 2D similar material simulation and FLAC3D numerical simulation models to investigate the development of mining-induced stress and the extraction effect of pressure-relief gas with large height and upward mining. Based on a comprehensive analysis of experimental data and observations, we obtained the deformation and breakage characteristics of strata overlying the coal seam, the development patterns of the mining-induced stress and fracture, and the size of the stress-relief area. The stress-relief effect was investigated and analyzed in consideration with mining height and three thick hard strata. Because of the group of three hard thick strata located in the main roof and the residual stress of mined panel 11124, the deformation, breakage, mining-induced stress and fracture development, and the stress-relief coefficient were discontinuous and asymmetrical. The breakage angle of the overlying strata, and the compressive and expansive zones of coal deformation were mainly controlled by the number, thickness, and strength of the hard stratum. Compared with the value of breakage angle derived by the traditional empirical method, the experimental value was lower than the traditional results by 3°-4°below the hard thick strata group, and by 13°-19° above the hard thick strata group. The amount of gas extracted from floor drainage roadway of B4 over 17 months was variable and the amount of gas per month differed considerably, being much smaller when panel 11223 influenced the area of the three hard thick strata. Generally, the stress-relief zone of No. 4 coal seam was small under the influence of the hard thick strata located in the main roof, which played an important role in delaying the breakage time and increasing the breakage space. In this study we gained understanding of the stress-relief mechanism influenced by the hard thick roof. The research results and engineering practice show that the main roof of the multiple hard thick strata is a critical factor in the design of panel layout and roadways for integrated coal exploitation and gas extraction, provides a theoretical basis for safe and high-efficient mining of coal resources.
基金Project(NRF-2016R1A2B4011009)supported by National Research Foundation of KoreaProject(KSTePS/VGST-KFIST(L1)/2017)supported by Vision Group of Science and Technology,Government of Karnataka,India
文摘In industrial applications involving metal and polymer sheets, the flow situation is strongly unsteady and the sheet temperature is a mixture of prescribed surface temperature and heat flux. Further, a proper choice of cooling liquid is also an important component of the analysis to achieve better outputs. In this paper, we numerically investigate Darcy-Forchheimer nanoliquid flows past an unsteady stretching surface by incorporating various effects, such as the Brownian and thermophoresis effects, Navier’s slip condition and convective thermal boundary conditions. To solve the governing equations, using suitable similarity transformations, the nonlinear ordinary differential equations are derived and the resulting coupled momentum and energy equations are numerically solved using the spectral relaxation method. Through the systematically numerical investigation, the important physical parameters of the present model are analyzed. We find that the presence of unsteadiness parameter has significant effects on velocity, temperature, concentration fields, the associated heat and mass transport rates. Also, an increase in inertia coefficient and porosity parameter causes an increase in the velocity at the boundary.
基金supported by the National Natural Science Foundation of China (No.50874103)the National Basic Research Program of China (No.2010CB226805)+1 种基金the Natural Science Foundation of Jiangsu Province (No.BK2008135)by the Open Foundation of State Key Laboratory of Geomechanics and Deep Underground Engineering (No.SKLGDUEK0905)
文摘The influence of an upper,mined coal seam on the stability of rock surrounding a roadway in a lower coal seam is examined.The technical problems of roadway control are discussed based on the geological conditions existing in the Liyazhuang Mine No.2 coal seam.The stress distribution and floor failure in the lower works after mining the upper coal is studied through numerical simulations.The failure mechanism of the roof and walls of a roadway located in the lower coal seam is described.The predicted deformation and failure of the roadway for different distances between the two coal seams are used to design two ways of supporting the lower structure.One is a combined support consisting of anchors with a joist steel tent and a combined anchor truss.A field test of the design was performed to good effect.The results have significance for the design of supports for roadways located in similar conditions.
文摘This paper takes No.52 return uphill roadway of Yangquhe coal mine as a research project. Based on the research, especially its geological condition, indoor experiments, numerical simulation and theoretical analysis were employed to determine the difficult coefficients of Yangquhe project. By using these means,the difficult coefficients of the deep rock engineering were determined. From a study of the effects of crustal stress and the roof mechanism on roadway stability, the transformation mechanism in Yangquhe coal mine has been determined. As a result of this research, the interactive support technology of prestressed cable mesh was developed and the technology tested in mining engineering, which proved to be feasible.
基金supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(No.51321065)the Program of International S&T Cooperation(No.S2015ZR1030)
文摘Containment booms are commonly used in collecting and containing spilled oil on the sea surface and in protecting specific sea areas against oil slick spreading.In the present study,a numerical model is proposed based on the N-S equations in a mesh frame.The proposed model tracks the outline of the floating boom in motion by using the fractional area/volume obstacle representation technique.The boom motion is then simulated by the technique of general moving object.The simulated results of the rigid oil boom motions are validated against the experimental results.Then,the failure mechanism of the boom is investigated through numerical experiments.Based on the numerical results,the effects of boom parameters and dynamic factors on the oil containment performance are also assessed.
基金Supported by China Postdoctoral Science Foundation under Grant No.20090460102
文摘To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.