CFD中文网

请问OpenFOAM有关于波浪与柔性结构相互作用的求解器吗？

下图是我基于不同的密度场(xh)进行速度场(U)和伴随速度场(Ua)的计算结果，xh有些小差异，但是伴随速度场计算出来完全不一样，右图计算出错了。
a712471d-a2cf-4ca5-8ca9-88d11620cffc-image.png

速度场计算是基于SA模型进行，在UEqn方程中增加了源项：

tmp<fvVectorMatrix> tUEqn ( fvm::div(phi, U) + turbulence->divDevReff(U) + fvm::Sp(alpha, U) == fvOptions(U) ); fvVectorMatrix& UEqn = tUEqn.ref(); UEqn.relax();

其中alpha是由密度场xh插值进行的，alphaMax=2.5e5，qu=0.01。

alpha=alphaMax*qu*(1-xh)/(qu+xh);

对于伴随速度场的计算，和UEqn类似，如下所示：

volVectorField adjointTransposeConvectiona((fvc::grad(U) & Ua)); tmp<fvVectorMatrix> tUaEqn ( fvm::div(-phi, Ua) + adjointTransposeConvectiona + turbulence->divDevReff(Ua) + fvm::Sp(alpha, Ua) == fvOptions(Ua) ); fvVectorMatrix& UaEqn = tUaEqn.ref(); UaEqn.relax();

请大家帮忙看看，会有哪些可能性导致代码计算出现异常。

@李东岳 在 SRFSimpleFOAM求解出来负压力 中说：

simplefoam求解器重要的是压差，而不是真实的值。因此负压是完全可能的。

那请问李老师，按理说流体经过这样一个通道，肯定会产生压损，我算过静止的沿程也是压力降低，可是这个压力居然沿程增高，明显不符合物理常识了。。。。8cbd17a4-0a42-4f6a-a369-3300fb3e7dca-image.png

@bestucan 多谢老师回复。这台工作站两个cpu，跨cpu线程计算的，不知道跟上传下载流量有关没，还是因为我这边是远程网络连接导致的。

我在指定跨cpu的线程计算时，如8和16线程，cpu很大概率上会有线程闲置或者无法满负荷运行，极少的情况下会100%运行，该情况下效率比较高；如果指定的是单cpu的线程时都是满负荷运行，但是在核数较多时效率反而较低，比如说32核。

如果是内存带宽受限的话，会体现在那个数据上呢，多谢老师解惑

我参考cantera官网示例中的diffusion_flame_batch.py，设置一批应变率变化时的一维对冲火焰解。通过反复调用前一计算结果作为初始解，来逐渐增加或降低应变率，从而得到计算结果。当我按照代码设置应变率上升时（即和源代码中一样设置strain_factor = 1.25），边界结果无明显问题。但当我设置应变率下降时（将strain_factor改为0.8），在燃料侧的grid=0处，其结果中的组分组成随应变率下降而逐渐偏离边界条件。我的燃料设置为Y_O2:0.2，Y_NH3:0.1417，Y_N2:0.6583。取应变率下降循环第50次的结果而言，grid=0处的O2仅为0.15左右，NH3也偏差很大。这个可能是什么原因引起的？

示例中的源代码部分如下：

# PART 3: STRAIN RATE LOOP # Compute counterflow diffusion flames at increasing strain rates at 1 bar # The strain rate is assumed to increase by 25% in each step until the flame is # extinguished strain_factor = 1.25 # Exponents for the initial solution variation with changes in strain rate # Taken from Fiala and Sattelmayer (2014) exp_d_a = - 1. / 2. exp_u_a = 1. / 2. exp_V_a = 1. exp_lam_a = 2. exp_mdot_a = 1. / 2. # Restore initial solution file_name, entry = names("initial-solution") f.restore(file_name, name=entry) # Counter to identify the loop n = 0 # Do the strain rate loop while np.max(f.T) > temperature_limit_extinction: n += 1 print('strain rate iteration', n) # Create an initial guess based on the previous solution # Update grid f.flame.grid *= strain_factor ** exp_d_a normalized_grid = f.grid / (f.grid[-1] - f.grid[0]) # Update mass fluxes f.fuel_inlet.mdot *= strain_factor ** exp_mdot_a f.oxidizer_inlet.mdot *= strain_factor ** exp_mdot_a # Update velocities f.set_profile('velocity', normalized_grid, f.velocity * strain_factor ** exp_u_a) f.set_profile('spread_rate', normalized_grid, f.spread_rate * strain_factor ** exp_V_a) # Update pressure curvature f.set_profile('lambda', normalized_grid, f.L * strain_factor ** exp_lam_a) try: # Try solving the flame f.solve(loglevel=0)

:146: :146:

foamyQuadMesh是之前测试性的代码，后来给删了。你编译时候强制略过好了。或者自动停止之后，手动继续编译别的

@random_ran 大佬好久没来了。

你应该类似输出U、position的方式来输出曳力或者曳力系数。需要改一些代码。但是我最近备课。没法更详细解释了。不过论坛里好像最近有几个做拉格朗日的挺活跃。比如这位大佬可能知道。 @youhaoyu

最近debug不了算例了，我要回国讲课了，备课比较紧张.. 看看有没有别的大佬给你康康吧

请教一下各位大佬，我在做双自由度涡激振动的时候，设置了如下的dynamicMeshDict文件，openFoam动网格程序是如何识别水平方向和竖直方向的弹簧的呢？还是说我这样的配置程序根本无法识别，但我这样的设置程序能正常计算，而且两个方向都有运动，不过模拟结果确实差蛮多。我需要怎么去调整呢？
d55f0aa65c8502fd9d43f11f8aea289.png

@youhaoyu
2ae0aafe-3a6e-4664-b101-52d186f51c46-image.png 这行代码我没对应上:xinlei:

感觉和网格有关，可以试试优化一下网格

感谢大佬百忙之中抽空解惑。对于壁面剪切力矢量，$\tau\cdot\mathbf{n}$在壁面上投影即为该壁面的剪切力，这点我是认同的。那其他方向的分量的意义是什么？是否为壁面所受到的压应力？还是说该分量没有意义，只是因为shearstress在近壁面网格体心处的值插值到壁面处得到的假值？

好的谢谢老师

各位大佬好，在使用OpenFOAM中的粒子信息输出过程中，有这么一条信息
1e907062-a87a-4c95-b9fc-297bafcb26e3-image.png
上图中标白色部分的Parcel fate system(number, mass)是什么意思呢？是逃离了系统吗？

不好讲故事是比较困境了，哎..

@Elibathe康 这位美女2018年之后就没来过了 :143:

算例是twoSpheresGlowinskiMPI

似乎解决了这个问题，把文本中的单引号改为双引号后，现在切换版本运行目前没啥问题:135:

@benqing 在 SRFSimpleFoam里边界条件的设置？ 中说：

你好，我最近准备使用SRF求解器，但是应用过程中一直存在问题，请问明白前述的问题了吗？
或者，有没有关于SRF求解器的资料可以分享一下？谢谢！

您好，我最近也困于SRF求解器，请问您现在怎么样了，有找到相关资料可以分享一下吗？谢谢