数学

Physics and Maths

濮祖荫(歌手朴树的父亲):“数学背后有物理,物理需要通过数学来表达”。数学方程(表达式)需要你通过物理直觉去思考,去书写。 濮说做研究的过程都是提出一个模型,找到问题的本质,针对围绕本质的几个变量进行研究,来近似替代现实的复杂模型(这个简化的思想对我有帮助)
… 阅读全文>>

TurbSim 8个湍流模型的对比

TurbSim,是NREL提供的一个小的开源工具, 用于产生湍流风,设置的时候可能对于不同模型会有不同的参数,基于范本进行了一番分析。当然测试文件在TurbSim的Test文件下的inp文件夹下。 如果对那些湍流模型做一个对比,那样是否可以更加清晰的看出不同呢? 于是就有了下面的shell分析脚本。 TurbSim产生的风力机一圈的入流截面模型(仅仅一个截面,全场是因为捕捉了很多个截面)
… 阅读全文>>

Incompact3d的Makefile及DNS算法

Incompact3d是一个开源的基于fortran语言编写的DNS求解器,也是我研究生期间使用的开源源代码。 Incompact3d整体框架的介绍可以参考软件的user-guide和文献High-order-compact schemes for incompressible flows: a simple and effcient method with the quasi-spectral accuary,还可以搜索sylvain Laziet 相关的文章 ,只不过他并不是成熟的软件,很多的编译和后处理都可能出现问题,下面是我使用过程遇到的一些问题。 当然主要过程是,通过makefile编译,然后运行incompact3d,最后处理计算结果。
… 阅读全文>>

Incompact3d如何通过结果处理出迹线

Incompact3D获得的数据都是二进制的数据文件,为了获得内部的文件 需要进一步通过编程获得,下面是一个获得pathline的源代码。
… 阅读全文>>

U+-Y+

这个gnuplot脚本是我用于处理Incompact3d产生出来的统计结果(如何统计部分,暂时未写)的画图脚本。
… 阅读全文>>

The greatest weapons

Everything in Haskell has a Type Here are some Type declarations. These are your greatest weapons.
… 阅读全文>>

Haskell几个基本概念--学而无用--无以为用

模式匹配 = +起着一种模式匹配,=不要翻译成等于 最好翻译为如果匹配前面+ 则他的值为后面.-> 也起着一种类型的模式匹配。 分支 |在创建函数的时候 体现着一种分支的效果。 if-else 这其实是一种控制语句,不能说Haskell放掉了控制语句,而应该是循环语句。 一切都是函数。 根据monad的内置结构,函数可以组合生成类似于数学的组合函数(埋藏内部的lambda结构)
… 阅读全文>>

尾递归

尾递归的好处就是快速计算,尾递归实际上是在递归计算的过程中, (印象中递归过程[表示语法形式-调用自己]和递归计算过程[表正线性方式 和非线性方式]是不一样的), 加入了迭代的思想,不断的修改了product和counter的值, 不需要树形展开。
… 阅读全文>>

求定积分


… 阅读全文>>

找出函数的不动点--Navier-Stokes方程

##Beautiful code! 数x为函数的f的不动点(f可以代表N-S方程),如果满足f(x)=x,则称x 为函数f的不动点。 性质 f(x),f(f(x)),f(f(f(x)))…., change the x’s value if f(x)=x,then you found!
… 阅读全文>>