CN105808820B - 一种求解区间热对流扩散问题的高精度数值方法 - Google Patents
一种求解区间热对流扩散问题的高精度数值方法 Download PDFInfo
- Publication number
- CN105808820B CN105808820B CN201610101396.8A CN201610101396A CN105808820B CN 105808820 B CN105808820 B CN 105808820B CN 201610101396 A CN201610101396 A CN 201610101396A CN 105808820 B CN105808820 B CN 105808820B
- Authority
- CN
- China
- Prior art keywords
- point
- thermal convection
- temperature
- diffusion problem
- interval
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000009792 diffusion process Methods 0.000 title claims abstract description 32
- 238000012546 transfer Methods 0.000 claims abstract description 25
- HBMJWWWQQXIZIP-UHFFFAOYSA-N silicon carbide Chemical compound [Si+]#[C-] HBMJWWWQQXIZIP-UHFFFAOYSA-N 0.000 claims abstract description 18
- 229910010271 silicon carbide Inorganic materials 0.000 claims abstract description 18
- 230000004044 response Effects 0.000 claims abstract description 17
- 239000000463 material Substances 0.000 claims description 5
- 239000000758 substrate Substances 0.000 claims description 3
- 239000013529 heat transfer fluid Substances 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013076 uncertainty analysis Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000004448 titration Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种求解区间热对流扩散问题的高精度数值方法,步骤如下:建立热对流扩散问题的传热控制方程;引入区间变量对传热模型中输入参数的不确定性进行表征,建立热对流扩散问题的区间控制方程;利用勒让德多项式对区间控制方程中的温度响应进行近似表示;根据张量积运算,确定高维空间中的配点集合;利用有限元程序计算所有配点处的温度响应,建立关于温度响应近似表达式中展开系数的线性方程组,并采用最小二乘法进行求解;基于多项式函数的光滑性,确定温度响应近似表达式的极值点,得到区间温度响应的上下界。本发明可***化解决含有区间不确定参数的热对流扩散问题,有效提高了区间数值方法的计算精度。
Description
技术领域
本发明属于机械工程领域,具体涉及一种求解区间热对流扩散问题的高精度数值方法。
背景技术
在自然界和各种生产技术领域中,通过流体流动来进行热量的传递是一种极其普遍的物理现象。尤其是在航空航天、能源化工等工业装备以及电子器件等精密产品的生产制造过程中,如何有效的实现热量传递,已成为产品设计的一个重要方面。现有关于热分析的许多研究都是针对确定性模型进行的,没有考虑模型输入参数的不确定性。实际工程中,由于制造工艺的限制、测量误差以及认知的局限,结构的材料属性、外部载荷和边界条件等物理参数不可避免的受到多种不确定因素的影响,使得***的温度响应也表现出一定的不确定性。
基于概率理论和数理统计的随机建模与计算方法在不确定***的分析、试验和设计方面已经取得了不少研究成果,成功解决了许多工程实际问题。但用随机理论求解问题时,往往假定不确定因素为随机变量或随机过程,这就事先需要大量的试验信息确定其概率分布规律。然而,多数实际情况下,获取充足的试验数据往往代价昂贵。如此一来,信息的缺乏使得概率模型不能真实反映客观实际,也就使得随机不确定分析方法失去意义。实际工程中的多数情况,设计者往往只关注某些响应的变化幅度,而获得不确定模型输入参数的取值范围相对概率密度函数来说要容易的多,且所需的不确定信息也大大减少,这就有效拓宽了不确定分析的应用范围。近些年来,利用区间理论和分析方法来处理工程中的不确定因素受到学者们越来越多的重视。目前,将区间理论与有限元计算方法相结合衍生出来的区间有限元法在不确定结构的静、动力特性分析方面已经取得了不少研究成果,但对具有区间参数的传热问题分析的文献却十分少见。另外,传统区间分析方法因区间运算所引起的区间扩张问题还比较严重,计算精度还亟待提高。因此,如何建立高精度的区间分析方法对不确定传热问题进行数值求解,是目前学术领域的一个研究热点,对于弥补现有传热数值计算方法的不足,具有重要的工程应用价值。
发明内容
本发明所要解决的技术问题为:克服现有技术在热对流扩散问题求解中存在的不足,充分考虑传热问题中的区间不确定因素,基于多项式逼近思想和配点分析理论,提出了一种预测温度响应区间变化范围的高精度数值计算方法,可***化解决含有区间不确定参数的温度场预测问题,在保证计算精度的同时,进一步降低了传统抽样方法的计算耗费。
本发明为解决上述技术问题采用的技术方案为:一种求解区间热对流扩散问题的高精度数值方法,包括以下步骤:
步骤一:根据传热模型建立热对流扩散问题的传热控制方程;
步骤二:引入区间变量对传热模型中输入参数的不确定性进行表征,根据步骤一中的传热控制方程建立热对流扩散问题的区间控制方程;
步骤三:利用勒让德多项式对步骤二区间控制方程中的温度响应进行逼近,得到区间温度响应的近似表达式;
步骤四:根据张量积运算,由一维空间的配点集合确定整个高维空间的配点集合;所述高维是指多于一维,维数等于区间变量的个数;
步骤五:利用有限元程序计算步骤四配点集合中所有配点处的温度响应,建立关于步骤三区间温度响应近似表达式中展开系数的线性方程组,并采用最小二乘法对此线性方程组进行求解,得到展开系数的一组值;
步骤六:将步骤五中得到的展开系数的一组值代回到步骤三温度响应的近似表达式中,基于多项式函数的光滑性,确定此近似表达式的极值点,进而得到区间温度响应的上下界。
其中,所述步骤三中利用勒让德多项式对温度响应进行近似表示,多项式的截断阶数并不是固定不变的,根据逼近精度要求选取,截断阶数越高,逼近精度就越高。
其中,所述步骤四中配点方案的选取并不是固定不变的,根据计算耗费和计算精度的要求来选取配点数量,配点数量越多,计算精度就越高,而计算耗费就越大。
其中,所述步骤四中的高维是指多于一维,维数等于区间变量的个数。
上述各步骤具体包括以下过程:
步骤一:根据传热模型建立热对流扩散问题的传热控制方程:
其中x表示物理坐标,T(x)表示温度响应,ρ,c,k分别表示材料的密度、比热容和热传导系数,u为传热流体的流动速度,Q(x)表示***的热源强度。
步骤二:引入n个区间变量对传热模型中输入参数的不确定性进行表征,并记为向量的形式其中上标I是区间符号,和表示区间变量的下界和上界,和称作区间变量的中点和半径,为标准区间变量根据步骤一中的控制方程建立热对流扩散问题的区间控制方程:
步骤三:利用有限阶的勒让德多项式对步骤二区间控制方程中的温度响应T(x,αI)进行逼近,得到区间温度响应的近似表达式:
其中Φi(αI)为事先选定的勒让德正交多项式基底,wi(x)为对应的展开系数,i=(i1,i2,...,in)表示多维指标,且满足|i|=i1+i2+...+in,N为此多项式的截断阶数。根据多项式理论,上述近似表达式中展开项的个数可用变量数n和截断阶数N表示为
步骤四:根据张量积运算,由一维空间的配点集合确定整个高维空间的配点集合。首先,对于一维区间变量来说,设定其配点数量为mi,则各个配点的具***置为:
其次,用点集表示一维区间变量内所有配点组成的集合,那么对于n个变量组成的高维空间而言,直接利用张量积运算,可得其配点集合Θ
而配点总数M为:
因此将Θ改写为的形式,用来表示高维空间所有的配点
步骤五:利用有限元程序计算步骤四配点集合中所有配点处的温度响应,建立关于步骤三区间温度响应近似表达式中展开系数的线性方程组,并采用最小二乘法对此线性方程组进行求解,得到展开系数的一组值。首先,步骤二中所建立的区间控制方程在配点处可改写为:
利用有限元程序对上述方程进行求解,可以得到所有配点处的温度响应
其次,基于步骤三中区间温度响应近似表达式,可以建立关于所有展开系数wi(x)的线性方程组:
其中表示多项式基底函数Φi(αI)在配点处的取值。
然后,利用最小二乘法求解上述方程组,得到展开系数wi(x)的一组值。
步骤六:将步骤五中计算得到的展开系数wi(x)的一组值代回到步骤三温度响应的近似表达式中,基于多项式函数TN(x,αI)的光滑性,令其一阶导数为零,易确定其极值点,连同边界点一起记为其中r为极值点和边界点的总体数量。比较这r个点处的温度值大小,最终确定区间温度响应T(x,αI)的下界和上界
本发明与现有技术相比的优点在于:
(1)与传统的热对流扩散问题分析方法相比,所提出的数值计算方法充分考虑实际工程中材料属性、外部载荷和边界条件的区间不确定性,计算结果对温度场分析具有更重要的指导意义。
(2)利用高阶勒让德多项式对温度响应进行近似表示,可有效提高逼近精度。同时,利用多项式函数的光滑性,可快速确定其极值点,进而得到温度响应的上下界。
(3)基于配点理论,可以充分利用原有确定性模型的有限元计算程序而无需对其做进一步的修改,保证了计算的可移植性。
(4)本发明操作简单,实施方便,有效提高了计算精度。
附图说明
图1为本发明的一种求解区间热对流扩散问题的高精度数值方法流程图;
图2为本发明的三维热交换器模型示意图;
图3为上部面板中心线区间温度响应示意图;
图4为流管中心线区间温度响应示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
本发明可***化解决含有区间不确定参数的热对流扩散问题,有效提高了区间数值方法的计算精度。本发明适用于含有区间不确定参数的热对流扩散问题的温度场预测。本发明实施方式以某三维热交换器模型为例,具体说明所述的一种求解区间热对流扩散问题的高精度数值方法。另外,此三维热交换器模型的区间温度响应数值计算方法可以推广到其他含有区间参数的热对流扩散问题温度场预测中。
一种求解区间热对流扩散问题的高精度数值方法的计算过程如图1所示,根据具体的传热模型建立热对流扩散问题的传热控制方程,引入区间变量对传热模型中输入参数的不确定性进行表征,建立热对流扩散问题的区间控制方程,利用勒让德多项式对温度响应进行近似表示,同时根据张量积运算,确定配点集合,利用有限元程序计算所有配点处的温度响应,对温度响应近似表达式中展开系数进行求解,并基于函数的光滑性,快速得到区间温度响应的上下界。可分为如下几个步骤进行:
步骤一:考虑图2所示的长为300mm的热交换器传热模型,截面为40mm×20mm的矩形,中间有30mm×15mm的方孔通过速度为u的冷却空气,入口处8的空气温度为Ts,结构的上面板7承受密度为qs=20000×sin(ω)W/m2的热流载荷,在上面板中心线上选定编号为1~3的三个点,在流管中心线上选定编号为4~6的三个点作为温度场的观测点。根据传热模型建立热对流扩散问题的传热控制方程:
其中x,y,z表示三个空间方向上的物理坐标,T(x,y,z)表示温度响应,ρ,c,k分别表示空气的密度、比热容和热传导系数,u为空气的流动速度。
步骤二:由于材料制造工艺的限制以及测量的误差,所有模型输入参数均含有一定的区间不确定性,引入六个区间变量对不确定性进行表征ρI=[1.3,1.5]kg/m3,cI=[900,1100]J/(kg·℃),kI=[0.023,0.029]W/(m·℃),uI=[1.8,2.2]m/s,Ts I=[18,22]℃,ωI=[1.5,2.5]。将所有区间变量统一表示为向量形式其中上标I是区间符号,和表示区间变量的下界和上界,和称作区间变量的中点和半径,为标准区间变量根据步骤一中的控制方程建立热对流扩散问题的区间控制方程:
步骤三:根据区间变量的分布特点,选用勒让德多项式对步骤二区间控制方程中的温度响应T(x,y,z,αI)进行逼近,截断阶数设定为N=3,得到区间温度响应的近似表达式:
其中Φi(αI)为事先选定的勒让德正交多项式基底,wi(x,y,z)为对应的展开系数,i=(i1,i2,...,i6)表示多维指标,且满足|i|=i1+i2+...+i6。此时上述近似表达式中展开项的个数为
步骤四:根据张量积运算,由一维空间的配点集合确定整个高维空间的配点集合。首先,对于一维区间变量来说,设定其配点数量为mi=5,则各个配点的具***置为:
其次,用点集表示一维区间变量内所有配点组成的集合,那么对于6个变量组成的六维空间而言,直接利用张量积运算,可得其配点集合Θ:
而配点总数为将Θ改写为的形式,用来表示六维空间所有的配点
步骤五:利用有限元程序计算步骤四配点集合中所有配点处的温度响应,建立关于步骤三区间温度响应近似表达式中展开系数的线性方程组,并采用最小二乘法对此线性方程组进行求解,得到展开系数的一组值。首先,步骤二中所建立的区间控制方程在配点处可改写为:
利用软件Nastran中的有限元程序对上述传热问题进行求解,可以得到所有配点处的温度响应
其次,基于步骤三中区间温度响应近似表达式,建立关于所有展开系数wi(x,y,z)的线性方程组:
其中表示多项式基底函数Φi(αI)在配点处的取值。
然后,利用最小二乘法求解上述方程组,得到展开系数wi(x,y,z)的一组值。
步骤六:将步骤五中计算得到的展开系数wi(x,y,z)的一组值代回到步骤三温度响应的近似表达式中,利用多项式函数TN(x,y,z,αI)的光滑性,令其一阶导数为零,确定其极值点,连同边界点一起记为其中r为极值点和边界点的总体数量。比较这r个点处的温度值大小,最终确定区间温度响应T(x,y,z,αI)的下界和上界
六个观测点处温度响应的计算结果如表1所示。与样本数为106的传统蒙特卡洛抽样方法对比可以看出,本发明方法的计算误差小于1%,计算精度完全满足工程需求。另外,从样本数量上看,本发明方法的样本数为15625,计算耗费远远小于蒙特卡洛方法。
表1观测点处区间温度响应上下界
除了上述六个观测点外,沿x轴方向,上部面板和流管中心线区间温度响应如图3和图4所示,横坐标表示沿x轴方向的空间位置,纵坐标表示空间位置处的温度值,实线和虚线分别表示蒙特卡洛抽样方法和本发明方法计算得到的结果。可以看出,本发明方法计算得到的温度响应上下界曲线与传统蒙特卡洛抽样得到的参考值吻合程度很好,计算结果真实可信。用本发明方法可以解决含有区间不确定输入参数的热对流扩散问题,计算精度高,此功能是一般商用软件所不能实现的。
以上所述的仅为本发明的较佳实施例而已,本发明不仅仅局限于上述实施例,凡在本发明的精神和原则之内所作的局部改动、等同替换、改进等均应包含在本发明的保护范围之内。
Claims (5)
1.一种求解区间热对流扩散问题的高精度数值方法,其特征在于包括以下步骤:
步骤一:根据传热模型建立热对流扩散问题的传热控制方程:
其中x表示物理坐标,T(x)表示温度响应,ρ,c,k分别表示材料的密度、比热容和热传导系数,u为传热流体的流动速度,Q(x)表示***的热源强度;
步骤二:引入n个区间变量对传热模型中输入参数的不确定性进行表征,并记为向量的形式其中上标I是区间符号,α i和表示区间变量的下界和上界,和称作区间变量的中点和半径,为标准区间变量进而根据步骤一中的传热控制方程建立热对流扩散问题的区间控制方程:
步骤三:利用勒让德多项式对步骤二区间控制方程中涉及的温度响应T(x,αI)进行逼近,得到区间温度响应的近似表达式:
其中Φi(αI)为事先选定的勒让德正交多项式基底,wi(x)为对应的展开系数,i=(i1,i2,...,in)表示多维指标,且满足|i|=i1+i2+...+in,N为此多项式的截断阶数;
步骤四:用点集表示一维区间变量内所有配点组成的集合,其中mi表示配点数量,表示一维空间中的配点;根据张量积运算,由一维空间的配点集合确定由n个变量组成的n维空间中的配点集合Θ
其中配点总数M为
进而将n维空间中的配点集合Θ改写为的形式,其中表示n维空间中的配点;
步骤五:基于步骤四中所建立的n维空间配点集合Θ,步骤二中所建立的区间控制方程在其所有配点处改写为:
利用有限元程序对上述方程进行求解,计算步骤四n维空间配点集合Θ中所有配点处的温度响应,建立关于步骤三区间温度响应近似表达式中展开系数的线性方程组,并采用最小二乘法对此线性方程组进行求解,得到展开系数的一组值;
步骤六:将步骤五中得到的展开系数的一组值代回到步骤三温度响应的近似表达式中,基于多项式函数的光滑性,确定此近似表达式的极值点,进而得到区间温度响应的上下界。
2.根据权利要求1所述的一种求解区间热对流扩散问题的高精度数值方法,其特征在于:所述步骤三中利用勒让德多项式对区间控制方程中的温度响应进行逼近时,其中的多项式的截断阶数并不是固定不变的,根据逼近精度要求选取,截断阶数越高,逼近精度就越高。
3.根据权利要求1所述的一种求解区间热对流扩散问题的高精度数值方法,其特征在于:所述步骤四中根据张量积运算,由一维空间的配点集合确定整个高维空间的配点集合中的配点方案的选取并不是固定不变的,根据计算耗费和计算精度的要求来选取配点数量,配点数量越多,计算精度就越高,而计算耗费就越大。
4.根据权利要求1所述的一种求解区间热对流扩散问题的高精度数值方法,其特征在于:所述步骤四中的n维是指多于一维,维数等于区间变量的个数。
5.根据权利要求1所述的一种求解区间热对流扩散问题的高精度数值方法,其特征在于:所述步骤五中利用有限元程序计算步骤四配点集合中所有配点处的温度响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610101396.8A CN105808820B (zh) | 2016-02-24 | 2016-02-24 | 一种求解区间热对流扩散问题的高精度数值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610101396.8A CN105808820B (zh) | 2016-02-24 | 2016-02-24 | 一种求解区间热对流扩散问题的高精度数值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105808820A CN105808820A (zh) | 2016-07-27 |
CN105808820B true CN105808820B (zh) | 2018-12-21 |
Family
ID=56466286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610101396.8A Expired - Fee Related CN105808820B (zh) | 2016-02-24 | 2016-02-24 | 一种求解区间热对流扩散问题的高精度数值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105808820B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846149B (zh) * | 2018-04-20 | 2020-02-07 | 北京航空航天大学 | 一种基于多源不确定性的结构分布式动态载荷识别的方法 |
CN112036039B (zh) * | 2020-09-01 | 2022-12-23 | 内蒙古科技大学 | 一种共轭传热材料热性能的高精度数值分析方法 |
CN112434447B (zh) * | 2020-12-17 | 2022-04-26 | 湖南大学 | 一种丝杠加工的时变可靠性分析***及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100094603A1 (en) * | 2001-12-31 | 2010-04-15 | George Danko | Multiphase physical transport modeling method and modeling system |
CN103366065A (zh) * | 2013-07-17 | 2013-10-23 | 北京航空航天大学 | 一种基于区间可靠性的飞行器热防护***尺寸优化设计方法 |
-
2016
- 2016-02-24 CN CN201610101396.8A patent/CN105808820B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100094603A1 (en) * | 2001-12-31 | 2010-04-15 | George Danko | Multiphase physical transport modeling method and modeling system |
CN103366065A (zh) * | 2013-07-17 | 2013-10-23 | 北京航空航天大学 | 一种基于区间可靠性的飞行器热防护***尺寸优化设计方法 |
Non-Patent Citations (3)
Title |
---|
Chong Wang等.Uncertainty analysis for heat convection-diffusion problem with large uncertain-but-bounded parameters.《Acta Mech》.2015, * |
Interval analysis of steady-state heat convection-diffusion problem with uncertain-but-bounded parameters;Chong Wang等;《International Journal of Heat and Mass Transfer》;20150831;全文 * |
含区间参数的结构-声耦合***可靠性优化设计;王冲等;《振动工程学报》;20141031;第27卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105808820A (zh) | 2016-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ozden et al. | Shell side CFD analysis of a small shell-and-tube heat exchanger | |
Wu | An approach combining body-fitted grid generation and conjugate gradient methods for shape design in heat conduction problems | |
Sivasankaran et al. | Numerical simulation on mixed convection in a porous lid-driven cavity with nonuniform heating on both side walls | |
Sivasankaran et al. | Natural convection in a porous cavity with sinusoidal heating on both sidewalls | |
Ismail et al. | Performance of a thermal storage system of the vertical tube type | |
CN108345741A (zh) | 基于无网格rkpm的各向异性材料结构二维热变形和热应力分析方法 | |
Karimi et al. | Numerical simulation of unsteady natural convection from heated horizontal circular cylinders in a square enclosure | |
Koo et al. | Shape optimization of inlet part of a printed circuit heat exchanger using surrogate modeling | |
CN106897520B (zh) | 一种含有模糊参数的传热***可靠性分析方法 | |
CN105808820B (zh) | 一种求解区间热对流扩散问题的高精度数值方法 | |
Wołoszyn et al. | Modelling of a borehole heat exchanger using a finite element with multiple degrees of freedom | |
Dalal et al. | Natural convection in a rectangular cavity heated from below and uniformly cooled from the top and both sides | |
Gao et al. | Element differential method with the simplest quadrilateral and hexahedron quadratic elements for solving heat conduction problems | |
Malinowski et al. | Inverse method implementation to heat transfer coefficient determination over the plate cooled by water spray | |
Wu et al. | Numerical simulation of subcooled nucleate boiling by coupling level-set method with moving-mesh method | |
Dash et al. | Natural convection from an eccentric square cylinder using a novel flexible forcing IB-LBM method | |
CN105677993B (zh) | 一种热传导热源位置识别反问题的数值通解方法 | |
Lan et al. | Shape design for heat conduction problems using curvilinear grid generation, conjugate gradient, and redistribution methods | |
CN107291969A (zh) | 一种变压器绕组温升计算的体网耦合分析法 | |
CN105760586A (zh) | 一种基于配点理论的模糊温度响应隶属度函数求解方法 | |
CN105302964A (zh) | 一种用于芯片结构的热分布分析方法 | |
Shah et al. | Correlation of experimental and theoretical heat transfer in mantle tanks used in low flow SDHW systems | |
CN105677995A (zh) | 一种基于全网格配点理论的模糊稳态热传导问题数值求解方法 | |
Wu et al. | Numerical analysis of natural convection in a porous cavity with the sinusoidal thermal boundary condition using a thermal nonequilibrium model | |
D'Avignon et al. | Modeling horizontal storage tanks with encapsulated phase change materials for building performance simulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181221 |
|
CF01 | Termination of patent right due to non-payment of annual fee |