位于pdf第25页公式(3.7)的上边:
"接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
应改为:
“接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x2, y2, z2) 点的时候”
位于pdf第25页公式(3.7)的上边:
"接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
应改为:
“接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x2, y2, z2) 点的时候”
厉害了!
解决了,因为频率太小,看不出来波动
研究物体入水的问题
网格动了,gnuplot测力也有,但是就是alphawater不动,我hull的boundary也是wall,怎么能不影响液体呢,我调下6Dof看下吧,有遇到相同情况的吗?
只是猜测不一定正确
你先检查你背景网格和移动网格是否很准确的重叠好了,检查方法:paraview里的后处理检查zoneID(0为背景网格,1为移动网格)
然后去重叠网格里看boundary力是否有oveset
或者oversetPatch
边界,没有就createPatch
添加一个
这些都ok就检查你的constant/dynamicMeshDict
的dynamicFvMesh
是不是dynamicOversetFvMesh
,solver
应该是multiSolidBodyMotionSolver
。
这些都没问题你再看看你的6DoF
是不是定义正确,还有可能就是你的alpha.water
中的边界定义问题
你都检查下吧
@Joann 这个报错是没找到你的c0域,你检查一下constant/polymesh/set
中是否有c0,没有的话就是你的topoSet
出了问题
今天空闲时间又看了一遍笔记发现有些笔误:
P9-第三章第一节第二段第五行“将手竖直的防止在水中” 应为“将手竖直的放置在水中”
前两天刚重做了0.0 输入了密码不进入桌面 我直接回档了
重叠网格 setFields设定水和水滴的区域就好了,只是想法
对于重叠网格运动网格方面做了新的尝试,在算例tut/multiphase/interFoam/laminar/sloshingTank3D6Dof中拷贝gen6Dof文件夹进入自己的background模型文件夹,/constant中的dynamicMeshDict也要拷贝到自己的/constant中,如果只想重叠网格中Zone1(模型)网格运动 需要把dynamicMeshDict改成如下,当然之前的topoSetDict也要增加movingZone的部分
dynamicFvMesh dynamicOversetFvMesh;
solver multiSolidBodyMotionSolver;
movingZone
{
solidBodyMotionFunction tabulated6DoFMotion;
CofG (0 0 0);
timeDataFileName "<constant>/6Dof.dat";
}
下面说下gen6Dof文件问题:
先进入/Make中
gen6Dof.C
EXE = ./gen6Dof
回到上阶文件夹,更改gen6Dof.C文件 ,更改所需运动,6个自由度随便改改,然后终端输入wmake
编译,完成之后终端输入 ./gen6Dof
,全部没有问题之后把生成的6Dof.dat拷贝到constant文件夹中,然后回到上阶文件夹运行 movingDynamicMesh
(别忘了改controDict文件的时常/6Dof时常,然后就可以观察 paraFoam
了。
大致思路就这样,我也试了的确可行。
如果有相似重叠网格中运动网格类问题可以参考下
https://www.cfd-china.com/topic/3126/overset来处理6dof的网格变形问题
这个兄弟的案例很好的解决了我的问题,添加topoSet就好了,背景网格设定1个hex就好
这两天改算例floatingBody时候卡住了,我在background文件夹中改了blockmesh文件,多加了两个hex,结果发现生成网格时候自动生成了2个region,所以在做topoSet时候变成了4个region,我就想知道怎么把background的3个region合并成1个region啊,有相关算例或者操作命令吗
@fireztw
Examples刚刚没注意格式 不好意思
wing motion
$FOAM_TUTORIALS/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam
floating object
$FOAM_TUTORIALS/multiphase/interDyMFoam/ras/floatingObject
DTC hull
$FOAM_TUTORIALS/multiphase/interDyMFoam/ras/DTCHull
pdf是步骤图。。。 效果图手机拍的没截图
厉害了!顶下!
@李东岳 我可是关注公众号好久了哈哈哈 估计老师猜到我是谁了
推导算法就像大学微积分或者物理力学每个公式老师都要给你推导一下怎么来的类似
另:NS笔记我看了才两遍 但是收益很大,现在当作字典来用,哪忘了就翻一下0.0
完结帖子,后来找到了,我的openmpi不是装在系统文件夹的,换一下bash文件里的export WM_MPLIB=OPENMPI就好了(如果openmpi-x.xx.x装在系统文件夹就是SYSTEMOPENMPI)
@Samuel-Tu 什么v? 你说v1906吗,那个没关系啊
按照官网的做的,没缺什么啊,openmpi也安了啊,咋就不行呢
source /home/chitie/OpenFOAM/OpenFOAM-v1906/etc/bashrc
mpicc: error while loading shared libraries: libopen-pal.so.13: cannot open shared object file: No such file or directory
是怎么回事啊
就是研究楔形物入水的附加质量啊
@东岳 东岳老师,你们都没算过附加质量吗。。
最后留两个Hess-smith方法求added mass的文献
附加质量-lamada
合力Fsum=(mass+lamada)*acceleration
lamada=(P*S)/g-m
想了下应该不是simpleFoam,因为并不是稳态,附加质量是和惯性有关,即加速度,那么就不是稳态,所以应该是pimple或者piso求解器
找到了个文献 http://html.rhhz.net/BJHKHTDXXBZRB/20160423.htm#R-19
差不多懂了
摘自CFDonline的原文:
I want to use potentialFoam to find added mass of a cylinder. It is well known this geometry has an added mass equal to the displaced fluid (rhopir^2*L).
Frank White, Fluid Mechanics, writes: "According to potential theory, mh depends only on the shape and direction of motion and can be computed by summing the total kinetic energy of the fluid relative to the body and setting this equal to an equivalent body energy: KEfluid=integral(0.5dmVrel^2)=0.5mhU^2)", where mh is the added mass. I interpret this for a numeric solver so that mh=sum(dmVrel^2)/U^2=sum(dVolumerho*Vrel^2)/U^2.
I have set up a model in OpenFOAM similarily to the example "Flow around a cylinder". The biggest difference being that the uppper and lower boundaries are set to walls with slip boundary condition, to enable uniform flow over the geometry. I have used 10 non-orthogonal corrector steps, to get correct results. I use a fluid volume that is about 4 diameters in the flow direction and 3 diameters across, use snappyHexMesh and extrudeMesh to get a nice mesh with boundary layers around the cylinder.
After solving with potentialFoam I extract the cell volume of each cell with "postProcess -func writeCellVolumes" and use an octave script to extract the cellvolume and velocity vectors from the OpenFOAM U and V files and calculate the added mass through a loop:
Code:
for i=1:L
Urel=U(i,:)-Uinit;
Urelsq=Urel.^2;
dm=rhoV(i);
mVsqtot=mVsqtot+dmUrelsq;
end
madd=mVsqtot./Uinit.^2;
where Uinit is the input flow, a vector with 1 m/s in the X direction, rho is hardcoded to 1000 kg/m^3 (water) and mVsqtot is an accumulator for the sum of dm*Vrel^2.
Unfortunately this does not give me correct result of 2.86 kg for a cylinder with radius of 30.165 mm, length 1 m, but rather close to half of that, 1.4361 kg. What am I doing wrong? Is it correct that the viscosity set is neglected with potentialFoam, but what density is used?
potentialFoam能求附加质量吗
最近做了wagner的楔形物如水的模拟,然后做了球鼻艏入水的二维模拟,需要求一下附加质量力,看了一些文章,也在openfoam里找了下,不知道用哪个求解器。我就看了下potentialFoam那个cylinder算例,那个能求吗。
各位大佬们,本人现在研究方向是slamming,能推荐一些书籍看看么,并且最近才接触openfoam,感觉要学的好多!