论文:2015,Vol:33,Issue(2):265-270
引用本文:
刘雪梅, 邓子辰, 胡伟鹏. 饱和多孔弹性杆热传导的广义多辛方法及其数值实现[J]. 西北工业大学学报
Liu Xuemei, Deng Zichen, Hu Weipeng. Generalized Multi-Symplectic Method and Numerical Experiment for Thermal Conduction of Saturated Poroelastic Rod[J]. Northwestern polytechnical university

饱和多孔弹性杆热传导的广义多辛方法及其数值实现
刘雪梅1,2, 邓子辰1, 胡伟鹏1
1. 西北工业大学力学与土木建筑学院, 陕西西安 710072;
2. 长安大学工程力学系, 陕西西安 710064
摘要:
首先根据多孔介质理论,利用饱和多孔介质的能量方程和本构关系,推导出饱和多孔弹性杆局部热平衡的热传导方程;继而引入正交变量,将热传导方程导入Hamilton系统,得到饱和多孔弹性杆热传导方程的广义多辛形式和多种局部守恒律形式;接着采用中点离散方法对热传导方程的广义多辛形式进行数值离散;最后利用计算机数值实现了饱和多孔弹性杆的热传导过程,并且讨论了参数取值的不同对热传导过程的影响,同时在数值模拟过程中记录了广义多辛格式的局部动量误差。研究结果表明,构造的广义多辛方法能够很好地模拟系统的热传导过程和耗散效应,同时也可长时间保持系统的固有几何性质。
关键词:    多孔介质    广义多辛    耗散    热传导   
Generalized Multi-Symplectic Method and Numerical Experiment for Thermal Conduction of Saturated Poroelastic Rod
Liu Xuemei1,2, Deng Zichen1, Hu Weipeng1
1. Department of Engineering Mechanics, Northwestern Polytechnical University, Xi'an 710072, China;
2. Department of Engineering Mechanics, Chang'an University, Xi'an 710064, China
Abstract:
Based on the theory of porous media, the thermal conduction equation of fluid saturated poroelastic rod is established by using the energy equation and constitutive relations of the two constitutes firstly in this paper. Then introducing orthogonal variables, we use the generalized multi-symplectic method to derive a first-order generalized multi-symplectic form for thermal conduction equation and several errors of conservation laws illustrating the local properties of the system. Thirdly, a midpoint box generalized multi-symplectic scheme is constructed; furthermore, discrete errors of generalized multi-symplectic conservation law and generalized local momentum conservation law are also obtained. Finally, the dissipation effect in thermal conduction process of saturated poroelastic rod and generalized local momentum conversation law are investigated numerically; moreover, the influence of parameter values for thermal conduction process is established later. From results of the numerical experiments, it can be preliminarily concluded that the generalized multi-symplectic scheme constructed in this paper has excellent accuracy, long-time numerical behavior and good conservation properties.
Key words:    boundary conditions    constitutive equations    errors    matrix model    momentum    numerical methods    porosity    temperature distribution    thermal conductivity    dissipation    generalized multi-symplectic    saturated poroelastic rod    thermal conduction   
收稿日期: 2014-09-10     修回日期:
DOI:
基金项目: 国家自然科学基金(11372252、11172239和11372253)与中央高校基金(2014G1121096)资助
通讯作者:     Email:
作者简介: 刘雪梅(1980-),女,西北工业大学博士研究生,主要从事多孔介质问题的保结构算法研究。
相关功能
PDF(1114KB) Free
打印本文
把本文推荐给朋友
作者相关文章
刘雪梅  在本刊中的所有文章
邓子辰  在本刊中的所有文章
胡伟鹏  在本刊中的所有文章

参考文献:
[1] 徐曾和,徐小荷.饱和多孔地层中定量抽放的流固耦合问题[J].岩土工程学报, 1999, 21:737-741 Xu Zenghe, Xu Xiaohe. Fluid-Solid Coupling Problem in the Liquid Extraction at Fixed Flux[J]. Chinese Journal of Geotechnical Engineering, 1999, 21:737-741(in Chinese)
[2] De Boer R, Svanadze M. Fundamental Solution of the System of Equations of Steady Oscillations in the Theory of Fluid-Saturated Porous Media[J]. Transport in Porous Media, 2004, 56:39-50
[3] 秦冰,陈正汉,方振东,孙树国,方祥位,王驹.基于混合物理论的非饱和土的热-水-力耦合分析模型Ⅰ[J].应用数学与力学,2010, 31:1476-1488 Qin Bing, Chen Zhenghan, Fang Zhendong, Sun Shuguo, Fang Xiangwei, Wang Ju. Analysis of Coupled Thermo-Hydro-Mechanical Behavior of Unsaturated Soils Based on Theory of MixturesⅠ[J]. Applied Mathematics and Mechanics, 2010, 31:1476-1488(in Chinese)
[4] Feng K. On Difference Schemes and Symplectic Geometry[C] //Proceeding of the 1984 Beijing Symposium on D D, 1984:42-58
[5] 赵长印,廖新浩,刘林.辛算法在动力天文中的应用[J].计算物理, 1992, 4:812-812 Zhao Changyin, Liao Xinhao, Liu Ling. Application of Symplectic Algorithm in Dynamic Astronomy[J]. Chinese Journal of Computational Physics, 1992, 9(4):812-812(in Chinese)
[6] Bridge T J. Multi-Symplectic Structures and Wave Propagation[J]. Mathematical Proceedings of the Cambridge Philosophical Society, 1997, 121:147-190
[7] Hu W P, Han S M, Deng Z C, Fan W. Analyzing Dynamic Response of Non-Homogeneous String Fixed at Both Ends[J]. International Journal of Non-Linear Mechanics, 2012, 47:1111-1115
[8] Hu W P, Deng Z C, Han S M,Zhang W R. Generalized Multi-Sympletic Integrators for a Class of Hamiltonian Nonlinear Wave PDEs[J]. Journal of Computational Physics, 2013, 235:394-406
[9] 胡伟鹏,邓子辰.大阻尼杆振动的广义多辛算法[J].动力学与控制学报, 2013, 11:1-4 Hu Weipeng, Deng Zichen. Generalized Multi-Sympletic Method for Vibration of Big Damping Rod[J]. Journal of Dynamics and Control, 2013, 11:1-4(in Chinese)
[10] Yang X. Gurtin-Type Variational Principles for Dynamics of a Non-Local Thermal Equilibrium Saturated Porous Medium[J]. Acta Mechanica Solida Sinica, 2005, 18:37-45
[11] 杨骁,程昌钧.流体饱和多孔介质的动力学Gurtin型变分原理和有限元模拟[J].固体力学学报, 2003, 24:267-276 Yang Xiao, Cheng Changjun. Gurtin Variational Principle and Finite Element Simulation for Dynamical Problems of Fluid-Saturated Porous Media[J]. Acta Mechanica Solida Sinica, 2003, 24:267-276(in Chinese)