- 博客(33)
- 收藏
- 关注
原创 typora中显示除号的问题
在latex中“除号(÷)” 通常用\div。但在typora中写数学公式时,却发现 “除号” 如果使用\div并没有显示为 ÷,而是 ∇⋅。
2023-12-10 21:24:42
254
原创 Star CCM+相关资料分享
Star-CCM+ 安装包及相关资料序号资源备注1相关资料序号资源备注1Star CCM+ TutorialFiles 模型文件提取码:454j 。包含软件自带教程案例的相关模型文件和计算文件等。2传热仿真视频教程及实例提取码:6s0m...
2022-11-15 23:19:25
7660
11
原创 有限体积法(13)——SIMPLE算法的一维算例
模型是一个二维平面喷嘴,如图,假设流动是稳态且无粘的,流体不可压。采用后向交错网格,如下图,压力节点5个,速度节点4个,入口给定滞止压强,出口给定静压
2020-11-19 19:06:29
6966
9
原创 有限体积法(12)——SIMPLE算法
承接上篇《交错网格》一文,本文介绍流动方程在交错网格上的离散以及SIMPLE算法。方程离散二维稳态的动量方程和连续性方程如下:∂∂x(ρuu)+∂∂y(ρvu)=∂∂x(μ∂u∂x)+∂∂y(μ∂u∂y)−∂p∂x+Su(1-a)\begin{aligned}\frac{\partial}{\partial x} (\rho u u) + \frac{\partial}{\partial y} (\rho vu) &= \frac{\partial}{\partial x} \left( \
2020-09-29 10:20:51
7715
2
原创 有限体积法(11)——交错网格
前面计算对流扩散问题的算例中都直接假设了计算域的流场是已知的(假设是均匀流场,各地流速的方向和大小都相同),然而,实际情况下流场是未知的,是需要求解的量,而其他输运量(例如ϕ\phiϕ)的对流通量是和当地速度的大小与方向密切相关的。因此,计算出流场是关键性的一步。本文主要介绍求解流场中面对的问题及解决方法。二维的动量方程和连续性方程如下:∂∂x(ρuu)+∂∂y(ρvu)=∂∂x(μ∂u∂x)+∂∂y(μ∂u∂y)−∂p∂x+Su∂∂x(ρuv)+∂∂y(ρvv)=∂∂x(μ∂v∂x)+∂∂y(μ∂v∂
2020-09-14 18:37:37
6330
转载 Star CCM+:RPM要被淘汰了
STAR-CCM+有两种不同的网格划分方法:基于计算域的网格划分(Region Based Meshing,RBM)和基于部件的网格划分(Part Based Meshing, PBM)。基于计算域的网格划分(RBM)是最初为STAR-CCM +开发的网格划分方法。现在已弃用该功能,并计划将其淘汰。 RBM网格与您定义的Regions和Continua关联,并且网格控制参数与Regions和Boundaries关联。另一方面,基于部件的网格划分(PBM)是最新的网格划分方法,它在性能和工作流程方面均具有
2020-08-11 18:16:14
1488
1
原创 有限体积法(10)——格式精度与待定系数法
前面用到的格式及其计算精度总结如下:格式名称公式精度中心差分格式二阶迎风格式一阶QUICK格式三阶
2020-08-11 17:55:40
2658
2
原创 有限体积法(9)——高阶差分格式:QUICK格式
迎风格式和混合格式只有一阶计算精度,虽然迎风格式使用起来非常稳定并且满足输运性要求,但一阶精度容易导致数值扩散的误差,可以通过使用高阶离散格式来降低这些误差。高阶格式一般需要使用更多的节点值,通过考虑更大范围内的影响来降低误差。中心差分格式具有二阶精度,但它的稳定性很差,也不满足输运特性。考虑到向中心差分格式这种不能包含流动方向信息的格式是不稳定的,我们就需要寻找那些包含流动方向信息的高阶差分格式。下面就介绍这类差分格式:QUICK格式。...
2020-07-22 18:53:07
9053
2
原创 有限体积法(8)——混合差分格式
Spalding(1972)提出了混合差分格式,该格式结合了中心差分格式和迎风格式的优点。小Pe数情况下(Pe<2Pe<2Pe<2),使用中心差分格式,它具有二阶计算精度;大Pe数情况下(Pe>2Pe>2Pe>2),使用迎风格式计算控制体界面对流输运量并忽略扩散作用。虽然迎风格式只有一阶精度,但可较好的反应流动的输运特征。混合差分格整合了中心差分格式和迎风格式的计算公式,使用分段线性的计算公式来近似通过网格边界面处的通量。通过左边界单位面积通量的混合差分格式计算公式为
2020-07-13 09:37:37
1820
1
原创 有限体积法(6)——离散格式的特性
为了能够在有限网格数下得到符合物理现实的数值解,那就需要采用的离散格式满足某些特性。最重要的几个有:守恒性、有界性和输运性。守恒性有界性输运性
2020-07-10 17:48:21
3773
原创 有限体积法(5)——对流-扩散方程的离散
关于变量ϕ\phiϕ的输运方程,∂(ρϕ)∂t+∇⋅(ρϕu)=∇⋅(Γ∇ϕ)+Sϕ(1)\frac{\partial (\rho \phi)}{\partial t}+ \nabla \cdot (\rho \phi \bold u) = \nabla \cdot (\Gamma \nabla \phi) + S_\phi \tag{1}∂t∂(ρϕ)+∇⋅(ρϕu)=∇⋅(Γ∇ϕ)+Sϕ(1)省略时间项就是稳态的对流扩散方程,∇⋅(ρϕu)=∇⋅(Γ∇ϕ)+Sϕ(2) \nabla \
2020-07-09 18:08:48
9010
8
原创 有限体积法(3)——一维扩散方程数值求解(第一类边界条件)
例1:无热源导热考虑一根细棒的导热问题,假设截面内温度均匀,问题简化为一维稳态导热问题。棒两端边界是恒温边界条件,TA=100℃T_A=100℃TA=100℃,TB=500℃T_B=500℃TB=500℃。导热系数k=1000W/m.Kk=1000 W/m.Kk=1000W/m.K,长度L=0.5mL=0.5mL=0.5m,截面面积A=0.01m2A=0.01 m^2A=0.01m2,棒内无热源,如图。划分网格首先将计算域均匀划分为5个网格单元,则节点间距δx=0.1m\delta x = 0.
2020-06-18 22:28:01
7342
1
原创 流体控制方程为什么要区分守恒形式和非守恒形式
突然想到这个问题,找了一下资料,发现在John D Anderson的《计算流体力学基础与应用》中对这个详细的阐述了这个问题。书中指出:在20世纪80年代之前,一般的流体力学和空气动力学教科书中都没有提到过这个问题,作者怀疑能否在这些书中找到关于守恒形式和非守恒形式对比的论述。方程放在那儿,但是并没有用特别的名词对它们加以区分。将控制方程分为守恒形式和非守恒形式,同时还关心对于给定的CFD问题,应该使用哪一种形式,这些都来源于现代CFD。简单点儿说,控制方程的守恒和非守恒形式本质上是相同的,数学上也
2020-06-10 19:44:05
8022
7
原创 Star CCM+ Interface踩坑——表面发射率
本文讨论计算辐射传热时多个计算域之间接触面上的发射率问题。以下面这个简单模型为例子,该模型有三个计算域,上、下都是金属铝固体,中间是空气,空气域这里简化为透明固体,忽略其导热(导热系数取1e-5 W/m-K),仅考虑辐射传热。空气层上、下两个接触面在生成Regine时会创建两个Interface。上、下表面为恒温边界,其他表面为绝热壁面。由于存在温差且中间是透明体,上、下两个金属体之间会通过热辐射传递热量。在Regine中可以看到,中间空气域(mid)和上、下金属体的接触面处生成了两个Interfac
2020-06-05 14:39:38
2897
3
原创 Star CCM+ interface踩坑——Interface没有生成的原因
对于存在多个计算域的仿真模型(比如,不同材料固体结构之间的导热、流体与固体之间的对流传热等),Star CCM+会识别体与体之间的接触面,在创建Part时会建立代表这些接触面的Contacts特征,然后在从Part创建Region时建立与Contacts特征对应的Interfaces。进行计算时,不同计算域之间的数据传递和耦合计算都是通过Interface来实现的。如果在创建Interface时失败,则在检查计算结果时就会发现物理量应该连续的地方却出现了间断。Star CMM+在从3D-CAD创建新Pa
2020-06-04 18:46:18
16403
3
原创 湍流模型(2)——雷诺平均方程
本文主要是不可压流动的雷诺平均方程(Reynolds-averaged Navier-Stokes equations,RANS)的推导。在推导雷诺方程之前先来总结一些关于时间平均的运算。对于瞬时量ϕ=Φ+ϕ′\phi=\Phi+\phi^\primeϕ=Φ+ϕ′和 ψ=Ψ+ψ′\psi=\Psi+\psi^\primeψ=Ψ+ψ′,根据时间平均的定义(见湍流的统计描述),有ϕ′‾=ψ′‾=0\overline{\phi^\prime}=\overline{\psi^\prime}=0ϕ′=ψ′=0
2020-06-03 21:50:30
20487
1
原创 湍流模型(1)——湍流的统计描述
湍流是一种混沌的、不规则的流动状态,内部是无数形状和尺度不同的漩涡。对于这种现象,分析数学就显得尴尬了,于是统计理论就登场了。
2020-06-01 23:09:19
3373
3
原创 Star CCM+ :迎浪小船 DFBI
在Star CCM+中,DFBI(Dynamic Fluid Body Interaction)模型可用于模拟刚体与流体的相互作用,可以计算刚体所受的力和力矩,并模拟刚体在与流体相互作用下的运动。下面仿真案例主要翻译自Star CCM+软件的user Guide文件,略有修改。使用的软件版本为 STAR-CCM+ 2019.1(14.02)。模型及材料参数设置1.启动软件,按Ctrl+N快捷键,然后按Enter键确定,新建一个模拟。2.如下图步骤,导入网格文件boat.ccm,然后保存该模拟为
2020-05-15 14:09:01
12496
10
原创 Star CCM+ 如何创建自己的物性参数数据库
Star CCM+提供了常见物质材料的物性参数,在设置物理模型时可以直接选用替换默认物质材料即可,但这些数据基本都是常数参数,对于特殊物质以及变化参数则需要自己根据具体情况进行修改。如果在工作中经常用到的物质则可以建立自己的物性参数数据库,Star CCM+提供了该功能。下面以液态Na为例,说明创建物性参数数据库的过程。首先,找到Tools -> Material Databases...
2020-04-24 11:29:07
4584
4
原创 Star CCM+如何修改默认单位
Star CCM+软件默认单位是国际单位制,例如长度单位默认是m。但建模一般是用mm单位的,可以按如下方法更改默认单位:在Tool节点下找到Units节点; 展开该节点,找到单位“mm”; 单击选中“mm”,在下面Properties中勾选Preferred选项。更改完后默认的长度单位就“mm”了。可以查看一下,在更改之前长度单位“m”的Preferred是默认勾选...
2020-04-05 09:55:19
3704
2
原创 将目录树下所有的某类型文件拷贝到指定目录下
从网上下载一些课程资料,资料已经被整理的很好,每一期的文件都放在一个目录里,包括mp3、mp4和pdf文件,但我想把mp3文件拷到手机上,用手动打开每个文件夹再把mp3文件拷出来太繁琐了,于是写了下面这个小程序。#! pythonimport osimport shutiladdress = '2016年书单' # 需要遍历的目录dest = 'mp3files' # 目标位...
2020-03-07 10:01:50
327
原创 二维平面内的坐标旋转变换
在做数值计算后处理中,需要对计算模型的边界数据进行坐标旋转,网上找到的公式大多是以原点为旋转中心点。因此这里直接考虑旋转中心点为一般坐标点的情况,将该部分内容作简单推导和整理。如上图所示,以点M (a,b)为圆心将点N0(x0,y0)旋转到N1(x1,y1)。点N0与水平轴的夹角为α,从点N0到N1的旋转角度为β,角度均已逆时针方向为正,顺时针方向为负。根据三角形内边角关系有,...
2020-01-05 22:22:35
3914
2
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人