学堂的井

家门口有个小学，
去学校的路上有口井，
大家都叫它学堂的井，
它的顶部是石头砌的很宽敞的井台，
平整的石头，已经被挑水的磨得很光滑，
告诉了它的年龄，也告诉了当年工匠的细心；
"学堂井的水，甜"，村里人都知道，
从我记事起，就有很多人吃这口井的水，
于是每次早上出门，门口都会有挑水人留下的水迹，
长长的，伸向大街小巷，
冬天还会结冰；
"挑水啊！"，记得门口的人经常这样打招呼；
后来我上了小学，每天都会经过这口井，
后来通了自来水，仍然有人来这里挑水，因为它甜
电话里，妈妈告诉我
这口井已经被拆了，
腾出地方将要建楼，

村里的房子现在涨的很贵了。
小学的操场都被挪来建房了。

妈妈说嫁到我们村就有这口井，
爸爸说，他小时候就吃这口井的水，

再也喝不到这井的水了

纳新活动小结

昨天，在交大举行了梅花桩协会的纳新表演活动，梅花桩也算是在交大开始生根了，但离发芽、扩枝、散叶还有很多路要走。

这次活动虽然从形式上给我了一个安慰——梅花桩协会成立，表演了，讲课了。但结果上还是比较令人沮丧的，我们一个人也没有招到。

想来想去，算是经验、也是教训吧。

1. 这次活动，通知是下午两点开始，但实际上，各位师兄弟从大老远跑来，时间很难把握，来了之后，又不会是马上表演，还要热身一下，开始表演也是三点的事情了。而我们给交大同学通知也说是两点开始，结果就是一部分有热情的人，来了之后，看没有，等，还没有，走人。都没观众了，还会有人来疯吗？还会有表演热情吗？以后可以给我们的观众通知晚点。
   
2. 宣传，几个武术类版都发过帖，另有300份传单。不过来的人不多。看来以后还要"恶意上十大"，这才是王道。
   
3. 活动场地，南区体育场是最好的活动场地，同时也是最好的踢球场地，结果来的人都忙着踢球呢，没人看。以后可以试着光明体育场，或思源湖边（虽然校动力科的老师说没办法连电，可是昨天中午看了有那么长的拖线板，其实那个老师还是骗人的）。
   
4. 没有组织好我们的会员，以后尽量将任务细化，分工做好。
   
5. 还有些自身的不足，嗯，不说了。
   

这次是我在交大组织的第一个活动，对于本科没有在此生活过，对这社团活动陌生的人来说，能做到这些还是不错的。呵呵，自我安慰下。

-----------------------------------------------------------

学校的一些资源，及如何获得。

如果是周一至周五，可以到学联添张单子，按照流程跑，可以跑的妥妥当当。

但是到了周末，就不好说了。

电：去游泳池附件的动力科找人，不过他们也不白跑，要"人工费"50块，比方说我们这次，说好了50块，结果人家放了我们鸽子，活动时没有人过来，还是自己连的电（在此感谢上午足球活动，及给我们帮助的那位老先生）

桌椅：跟体育馆借的，没要钱。原来还打算去租个。

条幅：图书馆下边

传单：据说印刷厂5分6一张。不过这次在小店里打的。

电脑音箱那些的都是周围同学帮着借来的，呵呵。

最后把活动花费贴出来，备忘

项目	花费（元）	说明
晚餐	200	四十张餐票
条幅	48	6米长条幅
饮用水	30	杯子、矿泉水，及饮料
传单	30	300张打印传单
总计	308

fluent材料库备份

;;;这是我做工作需要的材料，及其属性。

(window

     solid

    (chemical-formula . #f)

    (density (constant . 2500))

    (specific-heat (constant . 840))

    (thermal-conductivity (constant . 0.69999999))

     (formation-entropy (constant . 164448.08))

    (electric-conductivity (constant . 35410000))

     (magnetic-permeability (constant . 1.257e-06)))

 

(wall

    solid

    (chemical-formula . #f)

    (density (constant . 34))

    (specific-heat (constant . 2010))

    (thermal-conductivity (constant . 0.027000001))

    (formation-entropy (constant . 164448.08))

    (electric-conductivity (constant . 35410000))

    (magnetic-permeability (constant . 1.257e-06)))

    

(floor

    solid

    (chemical-formula . #f)

    (density (constant . 700))

    (specific-heat (constant . 2310))

    (thermal-conductivity (constant . 0.15000001))

    (formation-entropy (constant . 164448.08))

    (electric-conductivity (constant . 35410000))

    (magnetic-permeability (constant . 1.257e-06)))

    

(door

    solid

(chemical-formula . #f)

    (density (constant . 34))

    (specific-heat (constant . 2010))

    (thermal-conductivity (constant . 0.027000001))

    (formation-entropy (constant . 164448.08))

    (electric-conductivity (constant . 35410000))

    (magnetic-permeability (constant . 1.257e-06)))

    

(ceiling

    solid

    (chemical-formula . #f)

    (density (constant . 34))

    (specific-heat (constant . 2010))

    (thermal-conductivity (constant . 0.027000001))

    (formation-entropy (constant . 164448.08))

    (electric-conductivity (constant . 35410000))

    (magnetic-permeability (constant . 1.257e-06)))

；；；；；；；；；；；；；；；；；；；；

；；；其实至定义了density; specific-heat; formation-entropy三个参数，不知为什么在cas文件中，会出现这么多的内容，估计应该是程序动默认生成的。

去徐汇修电脑

这个blog真是让人又爱又恨，总也打不开，确还可以让你写。
实验室的电脑实在烦死人了，昨天去修，晚上就坏了，今天去修，回来看看还是没有好。我都快忙死了，还得忙这种事。不过那个老板给了句狠话：再坏，给你换全套。
这点，毫无悬念。
今天得知燕老师要过来。大概23号左右吧。看来要抓紧时间成立协会了。

网格划分

1. 面网格划分
   

Type

Map	产生规则的结构化网格
Submap	把一个非mappable面分成几个mappable面，从而在每个区域产生结构化网格这是个需要技巧的方法的，必须要对点的类型进行调整。
Pave	产生非结构化网格
Tri Primitive	把一个三边形面分成三个四边形部分，在每个部分生成结构化网格
Wedge Primitive	在楔形面的顶点产生三角形网格单元，从顶点往外生成发散性的网格

组合

方法	Quad	Tri	Quad/Tri
Map	*		*
Submap	*
Pave	*	*	*
Tri Primifive	*
Wedge Primitive			*

map-quad

用于边大于或者等于四的面，将其划为非常规整的网格，必要时需要对点的类型进行设置。

Quad/Tri-Map

网格划分格式仅仅适用于包含一个狭窄的由两个侧边组成的长条。（长条的两条自然是三角形，其余的就被划为四边形了，顶点类型：条形的尖部——Trielement，所有其它顶点——Side）

Quad-Pave & Tri-Pave

这是最简单的划分网格的方法。几乎没有任何限制。

Quad/Tri-Pave

当用户将Quad/Tri-Pave网格划分格式应用于一个面时，GAMBIT将生成一个平面网格，主要包含四边形网格单元但是在边界之间交角很小的角落采用三角形网格。用户可以设定某个点的类型决定该点附近是否要做三边形。

Wedge Primitive

这个方法画出来就像扇形，以某点为中心向外散射，注意要有三个点，trielement/end/end（组成了抽象的三角形），其他点设为side好了。

1. 体网格划分
   

Element

选项	详细说明
Hex	指定网格仅仅包含六面体网格单元
Hex/Wedge	指定网格主要有六面体网格单元组成但是也包括在适当地位置的楔形网格
Tet/Hybrid	指定网格主要由四面体网格构成但是在适当的位置可以包含六面体、锥形和楔形网格单元

Type

选项	详细说明
Map	生成一般六面体结构化网格单元
Submap	将一个不可图示化体积分割成可图示化区域并在每个区域生成六面体结构化网格单元
Tet Primitive	将一个四个侧面的体积分成四个六面体区域并在每个区域生成可图示化网格
Cooper	扫描整个体积的指定的源面的网格节点类型
Tet/Hybrid	指定该网格主要包含四面体网格单元但是在合适的位置也可以包含六面体、锥体和楔形单元
Stairstep	生成普通六面体网格和一个与原是提及形状近似的平滑的体积

组合

Type 选项	Elements 选项
	Hex	Hex/Wedge	Tet/Hybrid
Map	*
Submap	*
Tet Primitive	*
Cooper	*	*
TGrid			*
Stairstep	*

可以看出，Hex就是不可以与TGrid搭配，看名字就能看出来，至于Hex/Wedge和Tet/Hybrid，明明只有唯一的方法，偏要出来冒充可以和别人搭配的样子。

说明

在划分网格时，我们总是每种方法都尝试后，抱怨大部分的方法都不可以用，实际上不是这样的，而是在对体网格划分时，每种方法对顶点的类型都有很高的要求。顶点不去设置，或者类型不对，当然没有办法用了。所以当你想采用一种方法先去看看帮助吧。

另外，体网格的划分方法可选方案很多，不像面那样具有很强的针对性。所以体网格划分的确是个技术含量很高的划分阶段。

结构网格和非结构网格

在fluent教程中，我们常看到这两个名词，什么意思呢？

1.什么是结构化网格和非结构化网格：
结构化网格的拓扑结构具有严格的有序性。网格的定位能够用空间上的三个指标i，j，k识 别，且网格单元之间的拓扑连接关系是简单的i，j，k递增或递减的关系，在计算过程中不 需要存储它的拓扑结构。当流动区域易于被结构化网格所剖分、流动结构不需要作自适应 处理时，结构化网格仍被研究者乐于采用。 非结构化网格是一种无规则随机的网格结构，网格的定位只能用一维变量识别，网格的拓 扑连接关系是无规则的，需要在网格生成过程中存储网格的拓扑结构，这意味着非结构化 网格对计算存储需求量大。但利用非结构化网格易于剖分具有复杂边界的流动区域，且在 计算过程中可以在流场变化剧烈的区域内随意加密网格，故非结构化网格成为近十年以来 计算流体动力学研究的热点之一。（http://edu.stuccess.com/knowcenter/NumComp/algorithm/00000050.htm）
1.1结构化网格从严格意义上讲，结构化网格是指网格区域内所有的内部点都具有相同的毗邻单元。它可以很容易地实现区域的边界拟合，适于流体和表面应力集中等方面的计算。
它的主要优点是：
网格生成的速度快。
网格生成的质量好。
数据结构简单。
对曲面或空间的拟合大多数采用参数化或样条插值的方法得到，区域光滑，与实际的模型更容易接近。
它的最典型的缺点是适用的范围比较窄，只适用于形状规则的图形。尤其随着近几年的计算机和数值方法的快速发展，人们对求解区域的几何形状的复杂性的要求越来越高，在这种情况下，结构化网格生成技术就显得力不从心了。
1.2非结构化网格同结构化网格的定义相对应，非结构化网格是指网格区域内的内部点不具有相同的毗邻单元。即与网格剖分区域内的不同内点相连的网格数目不同。从定义上可以看出，结构化网格和非结构化网格有相互重叠的部分，即非结构化网格中可能会包含结构化网格的部分。

2.如果一个几何造型中既有结构化网格，也有非结构化网格，分块完成的，分别生成网格后，也可以直接就调入fluent中计算。

3.在fluent中，对同一个几何造型，如果既可以生成结构化网格，也可生成非结构化网格，当然前者要比后者的生成复杂的多，那么应该选择哪种网格，两者计算结果是否相同，哪个的计算结果更好些呢？
一般来说，结构网格的计算结果比非结构网格更容易收敛，也更准确。但后者容易做。影响精度主要是网格质量，和你是用那种网格形式关系并不是很大，如果结构话网格的质量很差，结果同样不可靠，相对而言，结构化网格更有利于计算机存储数据和加快计算速度。 结构化网格据说计算速度快一些，但是网格划分需要技巧和耐心。非结构化网格容易生成，但相对来说速度要差一些。

4.在gambit中，只有map和submap生成的是结构化网格，其余均为非结构化网格。

http://hi.baidu.com/Äã˵ÎÒÊÇÁ÷Ã¥ÍÃ/blog/item/b67925d92b9c032e11df9b9f.html

Intensity and Hydraulic Diameter的设置

在边界层设置，例如velocity inlet 和 pressure outlet时，往往需要设置 turbulence specification method，这里主要讲解其中Intensity and Hydraulic Diameter的设置方法。

Turbulence Intensity[%]，表示湍流的强度，该值越大，湍流就越厉害。它的定义式为I=u'/u=sqrt(u'^2+v'^2+w'^2)/u_avg，分子各量u'、v'、w'表示脉动值，分母u_avg就是平均速度。这个公式明显不适合用来计算湍流强度。对于管道流，可按以下公式计算I=u'/u_avg=0.16*(Re)^(-1/8)，Re的特征长度为水力直径。

http://219.232.54.3/cgi-bin/LB5000/topic.cgi?forum=6&topic=6112

对于入口边界条件建议使用k-e方法，k=0.05V^2,e=0.09k^2或是0，V为实际入流速度。

现在有很多人都是人为给出turbulence
intensity，这是很不合理的。对于出口边界条件建议采用intensity and Hydraulic
Diameter方法，尤其是当出现回流时候易于收敛，对于低Re建议使用intensity 用1%，高Re可用4-5%等。

Hydraulic Diameter[m],就是水力直径。

http://combust.hit.edu.cn/cgi-bin/topic.cgi?forum=109&topic=315

水力直径，当量直径有些时候，我们会涉及到水力直径这个量的计算。水力直径是水力半径的二倍，水力半径是总流过流断面面积与湿周之比。当量直径是总流过流断面面积的四倍与湿周之比。因此当量直径是水力直径的二倍。所谓湿周，就是在总流的过流断面上与流体相接触的固体边壁周长

举例说明一下，例如某风口的尺寸为0.49*0.07，各个数如下显示

尺寸	水利半径	水力直径
0.07*0.49	0.030625	0.6125

速度V	水力直径	粘性系数	Re	湍流强度
4.3	0.6125	1.49e-5	1.77e+4	0.047