- 博客(7)
- 资源 (3)
- 收藏
- 关注
原创 数据处理器的执行和集成机制:applyProcessingFunctional
源码位置:/src/atomicBlock/dataProcessorWrapper3D.hh首先就为我们展示了applyProcessingFunctional。template<typename T, template<typename U> class Descriptor>void applyProcessingFunctional ( LatticeBoxProcessingFunctional3D<T,Descriptor>* func
2020-12-27 03:57:02 534 3
原创 addNaiveForce和addGuoForce的源码位置
/src/latticeBoltzmann/externalForceTemplates3D.h以addNaiveForce为例,我们可以看到palabos是如何将外力场的3个分量分别施加到相应的18的方向的分布函数f’s上。static void addNaiveForce( Array<T,descriptors::ForcedD3Q19Descriptor<T>::q>& f, T* externalS
2020-12-27 03:43:54 288
原创 Palabos源码:momentTemplates3D
这是比较底层的代码,基本上我们在palabos的各种操作都是基于数据处理器上的,而数据处理器则是通过各种方式调用这些底层的代码进行基于分布函数f的运算。以D3Q19为例,介绍一些常见的功能。typedef descriptors::D3Q19DescriptorBase Descriptor;get_rhoBarT rhoBar = f[0] + f[1] + f[2] + f[3] + f[4] + f[5] + f[6] + f[7] + f[8]
2020-12-22 23:33:23 225
原创 Palabos的数据处理器
Palabos算例涉及大量数据处理器,我稍作整理,因博客不便于ctrl+f检索,我将其解读写在pdf文件上,上传至github供大家参考。我以movingWall为例演示一下如何使用这本工具书,在算例中有这两个重点数据处理器:std::vector<MultiBlock3D*> rhoBarJarg;rhoBarJarg.push_back(lattice);rhoBarJarg.push_back(rhoBar);rhoBarJarg.push_back(j);integrateP
2020-12-21 16:50:53 447 2
原创 Galerkin method 热传导公式推导过程
以下公式推导来源于我2020 fall学期的有限元课程作业,作业并没有反馈推导过程和结果是否百分之百正确,在这里供各位同学参考,如果有错误也欢迎指出。代码实现可以参考deal.ii的step 26。最后分享一下如何用Matlab显示vtk格式的数据,因为课程要求3D展示2d的算例数据,z轴为温度,而matlab不是很好直接处理vtk的数据。这里我们可以用paraview先打开vtk的数据,然后导出为csv格式文件,此时用matlab打开,便易于处理了。...
2020-12-19 23:29:42 2154 6
原创 step26 heat conduction theta scheme
在这个算例中,我们需要注意的是,当theta大于等于0.5时,算例才会无条件收敛,误差不会线性增加,当theta为0显式法对时间步长有一定的要求。
2020-12-19 23:16:17 197 1
原创 测试了一下Palabos访问cell的population
之前在论坛和其他地方,得知了一些访问f’s的方法,但是不够确信,现在跑多次测试看看效果。结论:cell[iPop],lattice.get(ix,iy,iz)[iPop],cell.getRawPopulations()三者得到的结果相同。
2020-12-14 12:36:52 279
Tutorial_ the open-source library Palabos in your daily work
2020-10-08
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人