力学分析开源软件简介

《力学分析开源软件简介》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、Palabos(http://www.lbmethod.org/)PalaBos的是一款高效的流体模拟及其建模库，开发基于C＋＋的STL(标准模板库)，有极强的拓展性！尽管其源代码是开放，但是基于PalaBos的FlowKit公司已于2011年9月开始运营（http://www.flowkit.com/），主要为流体力学相关领域提供解决方案，并定制软件。主要的开发者为我的日内瓦朋友JonasLatt博士，另外一个重要开发成员Orestis博士也是我的合作者和好朋友，其主要的贡献在于湍流模型和多块加密的代码的开发。在版本1.0中，目前二维的多块

2、加密是可用的，三维的曲面边界可用，需要提供stl几何文件（参：examples/showCases/aneurysm）。PalaBos的主要特点在于，其在并行结构上采取并行机制与模型分离的方式，使得应用建模与并行机制不相关。这也使得PalaBos的易于扩展。下面举例来说明其代码特点：对于二维计算下面两个基本的文件必须包括#include"palabos2D.h"#include"palabos2D.hh"#include#include#include#include#inc

3、lude基本的名字空间usingnamespaceplb;usingnamespaceplb::descriptors;usingnamespacestd;typedefdoubleT;基本模型的描述。对于PalaBos，众多模型的应用，都是通过DnQmDescriptor来描述的。用户可自定义！#defineDESCRIPTORD2Q9Descriptor初场的定义，建议使用这种方法TpoiseuilleVelocity(plintiY,IncomprFlowParamconst¶meters){Ty=(T)

4、iY/parameters.getResolution();return4.*parameters.getLatticeU()*(y-y*y);}压力的定义TpoiseuillePressure(plintiX,IncomprFlowParamconst¶meters){TLx=parameters.getNx()-1;TLy=parameters.getNy()-1;return8.*parameters.getLatticeNu()*parameters.getLatticeU()/(Ly*Ly)*(Lx/(T)2-(T)i

5、X);}密度场的定义TpoiseuilleDensity(plintiX,IncomprFlowParamconst¶meters){returnpoiseuillePressure(iX,parameters)*DESCRIPTOR::invCs2+(T)1;}根据坐标进行速度初始化templateclassPoiseuilleVelocity{public:PoiseuilleVelocity(IncomprFlowParamparameters_):parameters(parameter

6、s_){}voidoperator()(plintiX,plintiY,Array&u)const{u[0]=poiseuilleVelocity(iY,parameters);u[1]=T();}private:IncomprFlowParamparameters;};根据坐标进行密度初始化templateclassPoiseuilleDensity{public:PoiseuilleDensity(IncomprFlowParamparameters_):parameters(parameter

7、s_){}Toperator()(plintiX,plintiY)const{returnpoiseuilleDensity(iX,parameters);}private:IncomprFlowParamparameters;};零速度场初始化templateclassPoiseuilleDensityAndZeroVelocity{public:PoiseuilleDensityAndZeroVelocity(IncomprFlowParamparameters_):parameters(paramete

8、rs_){}voidoperator()(plintiX,plintiY,T&rho,Array&u)const{rho=poiseuilleDensit