力学分析开源软件简介.doc

Palabos (/) PalaBos的是一款高效的流体模拟及其建模库，开发基于C＋＋的STL(标准模板库)，有极强的拓展性！尽管其源代码是开放，但是基于PalaBos的FlowKit公司已于2011年9月开始运营（/），主要为流体力学相关领域提供解决方案，并定制软件。主要的开发者为我的日内瓦朋友Jonas Latt博士，另外一个重要开发成员Orestis博士也是我的合作者和好朋友，其主要的贡献在于湍流模型和多块加密的代码的开发。在版本1.0中，目前二维的多块加密是可用的，三维的曲面边界可用，需要提供stl几何文件（参：examples/showCases/aneurysm）。PalaBos的主要特点在于，其在并行结构上采取并行机制与模型分离的方式，使得应用建模与并行机制不相关。这也使得PalaBos的易于扩展。下面举例来说明其代码特点：对于二维计算下面两个基本的文件必须包括#include palabos2D.h#include palabos2D.hh#include vector#include cmath#include iostream#include fstream#include iomanip基本的名字空间using namespace plb;using namespace plb::descriptors;using namespace std;typedef double T;基本模型的描述。对于PalaBos，众多模型的应用，都是通过DnQmDescriptor来描述的。用户可自定义！ #define DESCRIPTOR D2Q9Descriptor初场的定义，建议使用这种方法T poiseuilleVelocity(plint iY, IncomprFlowParamT const parameters) {T y = (T)iY / parameters.getResolution();return 4.*parameters.getLatticeU() * (y-y*y);}压力的定义T poiseuillePressure(plint iX, IncomprFlowParamT const parameters) {T Lx = parameters.getNx()-1;T Ly = parameters.getNy()-1;return 8.*parameters.getLatticeNu()*parameters.getLatticeU() / (Ly*Ly) * (Lx/(T)2-(T)iX);}密度场的定义T poiseuilleDensity(plint iX, IncomprFlowParamT const parameters) {return poiseuillePressure(iX,parameters)*DESCRIPTORT::invCs2 + (T)1;}根据坐标进行速度初始化templatetypename Tclass PoiseuilleVelocity {public:PoiseuilleVelocity(IncomprFlowParamT parameters_): parameters(parameters_){ }void operator()(plint iX, plint iY, ArrayT,2 u) const {u[0] = poiseuilleVelocity(iY, parameters);u[1] = T();}private:IncomprFlowParamT parameters;};根据坐标进行密度初始化templatetypename Tclass PoiseuilleDensity {public:PoiseuilleDensity(IncomprFlowParamT parameters_): parameters(parameters_){ }T operator()(plint iX, plint iY) const {return poiseuilleDensity(iX,parameters);}private:IncomprFlowParamT parameters;};零速度场初始化templatetypename Tclass PoiseuilleDensityAndZeroVelocity {public:PoiseuilleDensityAndZeroVelocity(IncomprFlowParamT parameters_): parameters(parameters_