kelie

导航

2014年2月27日 #

稳态与瞬态的区别

摘要: 稳态与瞬态的区别主要体现在控制方程是否存在时间项上。换句话说,其区别在于计算结果是否是与时间相关。但是我们观察现实生活,似乎找不到什么现象是与时间无关的。于是我们可以这样理解:稳态是一种近似。还是不太好理解,我们来举个例子。假设雨滴从高空落下,其阻力与运动速度的平方成正比,比例系数为1。雨滴质量为1,重力加速度为g,假设雨滴运动初速度为0,则依据牛顿定律很容易得出当阻力与重力平衡时,该雨滴将获得最大速度。学过物理的人都知道在0.32s时雨滴达到最大速度3.13m/s,阻力9.8与重力平衡后其将保持匀速运动。好了,我们可以将运动状态分为两部分,以t=0.32s为界,在此之前,运动速度与时间有关, 阅读全文

posted @ 2014-02-27 15:39 kelie 阅读(24291) 评论(0) 推荐(5) 编辑

CFD计算

摘要: 47 求解器为flunet5/6在设置边界条件时,specify boundary types下的types中有三项关于interior,interface,internal设置,在什么情况下设置相应的条件?它们之间的区别是什么?interior好像是把边界设置为内容默认的一部分;interface是两个不同区域的边界区,比如说离心泵的叶轮旋转区和叶轮出口的交界面;internal;请问以上三种每个的功能?最好能举一两个例子说明一下,因为这三个都是内部条件吧,好像用的很多。interface,interior,internal boundary区别?在Fluent中,Interface意思为“ 阅读全文

posted @ 2014-02-27 15:37 kelie 阅读(3208) 评论(0) 推荐(3) 编辑

边界层网格的一些理解

摘要: 我大概说一下关于边界层这里的东西,大都是自己理解后的,希望能帮到在边界层上遇到困难的朋友,写的不好,轻喷。首先,为什么会出现边界层、Y+这些东西,追根究底是因为湍流的数值求解的需要。先回顾一下湍流数值模拟的方法,我们只针对Euler方法,其下分为直接数值模拟(DNS)和非直接数值模拟。为什么不用直接模拟,第一,湍流中充满了各宗尺寸的涡,第二,即使计算机硬件条件达到了DNS模拟要求,要精确给出满足最小尺度量合理的边界条件和初始条件是不可能的,第三,为了减少耗散和色散,DNS中常采用高阶方案离散方案,由此生成边界条件和处理复杂集合外形的流动很困难。那我们只能采用非直接模拟,现在多采用Reynold 阅读全文

posted @ 2014-02-27 15:24 kelie 阅读(9580) 评论(0) 推荐(1) 编辑

网格无关性验证方法

摘要: 弟:“上次说到了在进行计算结果评估的时候需要做无关性评价,这个无关性的概念应该怎么去理解呢?”哥:“这里的无关性验证主要是指网格无关性,在一些特殊在场合中可能包括有时间步长无关性检验。但是稍微有点数值计算常识的人都知道,计算结果不可能与网格大小无关的。我们这里的无关是一种近似的概念。”弟:“求真相。”哥:“我们先讨论网格无关的概念,步长无关的概念与这个相似。数值计算中之所以需要网格,是由所采取的算法密切相关的。当前的主流偏微分方程数值离散方法都是先计算节点上的物理量,然后通过插值在方式求得节点间的值。因此,从理论上讲,网格点布置得越密集,所得到的计算结果也越精确。”哥:“但是网格不可能无限制的 阅读全文

posted @ 2014-02-27 14:55 kelie 阅读(32312) 评论(0) 推荐(0) 编辑

ICEM相关

摘要: 1,几何体建模不用讲(可以不学,因为通常是其他软件导入)在初始分块前,建立part,为建立边界条件使用(这是部分的定义最重要的作用,所以你可以按照不同的情况来定义,划分网格只是块的工作),所以对于三维,点part不一定要建立,可以直接缺省的geom,线也不必,2,块划分:首先初始分块,然后spilt 块,然后删除多余块。构造几何的外形。3,建立曲线映射对应关系,移动点来更好表现。(在这先前可以群组曲线,也可不用),另一种是先映射点,(对于以后要移动的点先别映射),再映射边。(为了方便以后的点移动,可以先映射边,再移动点。为了方便观察,可以先映射点,这样自动移动点到指定位置,然后映射边)也就是说 阅读全文

posted @ 2014-02-27 14:50 kelie 阅读(2764) 评论(0) 推荐(0) 编辑

Y+的查看及FLUENT壁面函数的选择

摘要: y+的查看其实,我们关心的应该是壁面y+值。那么我们看云图的话,是可以直接看到的,但是个人感觉,如果case大的话,也不是很方便。此外,你要是看云图的话,要用filled的方式,而且把node value选上,不然显示的是插值结果,那样不对。推荐你用plot图看,在Plot里面,有xy-plot和histogram两个。这两个都要选择统计位置,请把所有的壁面选上,注意里面有内部体的boundary名称,如那些interior,就不要选。y+的值在turbulence/Wall Y plus里面。然后进行historam或plot。histogram里面能告诉你不同y+的网格个数都有几个,其中最 阅读全文

posted @ 2014-02-27 14:45 kelie 阅读(54867) 评论(0) 推荐(0) 编辑

Y+的一些讨论

摘要: 一、关于fluent计算时壁面函数法和网格的关系,还有一个小问题1:各位用fluent的同仁和高手们,我想要比较好的使用fluent软件,最重要的就是要学好理论,在这里我想请教各位一个问题,在使用标准k-eplison和一些其他的封闭模型时,对于近壁区的流动要使用壁面函数法求解。那么在划分网格时,是不是一定要把把第一个内节点布置在湍流充分发展的区域内呢?我们如果自动生成网格时,如果说第一个节点在壁面的粘性底层内,是不是对计算有一定的影响呢?还有一个问题就是在gambit中设置的wall壁面,怎么到fluent设置为内部表面interior,好像在边界条件设置时没有这个边界呀。2:为什么要用壁面 阅读全文

posted @ 2014-02-27 14:43 kelie 阅读(12080) 评论(0) 推荐(1) 编辑

MATLAB中FFT的使用方法

摘要: MATLAB中FFT的使用方法说明:以下资源来源于《数字信号处理的MATLAB实现》万永革主编一.调用方法X=FFT(x);X=FFT(x,N);x=IFFT(X);x=IFFT(X,N)用MATLAB进行谱分析时注意:(1)函数FFT返回值的数据结构具有对称性。例:N=8;n=0:N-1;xn=[4 3 2 6 7 8 9 0];Xk=fft(xn)→Xk =39.0000 -10.7782 + 6.2929i 0 - 5.0000i 4.7782 - 7.7071i 5.0000 4.7782 + 7.7071i 0 + 5.0000i -10.7782 - 6.2929iXk与xn的维数 阅读全文

posted @ 2014-02-27 13:56 kelie 阅读(1289) 评论(0) 推荐(0) 编辑