中级 地球科学-热传递.ppt

COMSOL Multiphysics“地球科学模块”中级培训——热传递 热传递 “Conduction”模式 “Convection and Conduction”模式 偏微分方程– 传导和对流 边界设定 Example: Free Convection in Porous Media 背景：这个算例展示了多孔介质中对流和传导的数值模拟，随着温度流体的密度发生变化，由于密度不同产生了浮力成为流体流动的驱动力（浮力驱动流动），模拟结果与发表的研究论文进行了对比。 模型定义： 右侧的图片给出了几何模型和边界条件。流体在计算域内流动，但不能流出。温度在外边界由高到底变化。初始谁是静止的，但温度改变是流体密度发生变化继而发生流动。模拟稳态的状况。 浮力驱动流动用布里克曼方程定义，从布里克曼方程计算的流体的速度被应用到热传递模式中。 控制方程 布里克曼方程——方程右项是由于热扩散导致的力 对流、传导方程 Implementation: Initial Conditions for Boussinesq Approximation 求解非线性问题的策略 为了减轻非线性方程的数值收敛困难，用户可以在Boussinesq项前加一个松弛因子damp，来逐步找到一个求解真值的初始条件。当damp=0的时候，相当于流场和温度场是不耦合的，此时收敛是很容易的。然后逐步增大damp值，用之前的求解作为下一个参数步的初值，直到damp=1。这种迭代方式用COMSOL Multiphysics 的参数求解器是很容易实现的。 建模步骤： Options and settings OptionsExpressionsScalar Expressions 创建几何模型 求解域控制方程、边界条件设定 Multiphysics - Brinkman Equations (chns). initial conditions will at their default values: u = v = p = 0 边界设定 默认的边界条件是 No slip ，在本算例中不需要被改变。 点设定 Multiphysics - Convection and Conduction in Porous Media (eshcc). “Conduction”模式 偏微分方程 应用模式属性 求解域控制方程设定 边界设定 点和边设定 标量变量 偏微分方程—— 传导 T = temperature, the dependent variable Ceq = volumetric heat capacity – calc optional Keq = thermal conductivity – calc optional QG = geothermal heat source QH = general heat source QC = convective heat source QR = radiative heat source δts is an optional coefficient added mostly for user convenience 地球科学模块自动定义 Ceq 和Keq Θ-volume fraction; ρ-density; Cp -specific heat capacity; K -thermal conductivity, 下标“L” 和 “P” 指代流体和固体 地热源 热传导应用模式预定义了 均匀分布和指数分布的热源，用户也可以自定义。. 指数型 均布型 ρgeo - the geothermal density, qgeo - the radiogenic heating per unit mass zgeo - a model coordinate x, y, or z hgeo - length scale 对于 1D 和 2D 模型，热源项还包括下述两项 QC –对流热源 h - a heat transfer coefficient, Tinf - the external temperature QR – 地表辐射的影响 Cconst often is defined as the product of the Stefan-Boltzmann constant, σ, and the emissivity of the surface, ε. 应用模式属性 求解域控制方程设定 求解域方程设定对话框包含五类选项： (1) Materials, (2) Thermal Properties, (3) Sources/Sinks, (4) Init, (5) Element. (