CFD中文网

好的好的

我发现这个欧拉场到拉格朗日场，拉格朗日场到欧拉场的映射，有各种不同的方法..

Sorry, Chinese input missing....
0_1476947519913_oie_yyxHcXChIqBJ.png

有个问题想问问大伙，如题

@cloud1990 在 VOF求解自由水面问题，求解空气区域吗？ 中说：

Image

Left: Correct jet shape predicted by TruVOF technique used in FLOW-3D. Right: Incorrect jet shape predicted by pseudo-VOF technique used by other CFD codes.

Image
Left: FLOW-3D‘s TruVOF technique predicts jet impingement on wall and some outflow. Right: Pseudo-VOF methods don’t predict realistic jetting of fluid on side walls.

从结果上看，flow3D 所谓的TruVOF是有问题的，air velocity 和　air resistance = 0?　这个可能比较适合粗网格计算，作为一个tricky的存在．但理论上，这个TruVOF是无法解决wind wave的问题．

想做重叠网格的入水仿真，然后前景网格导入进背景网格后，显示云图时出现弹体之外的网格缺失了。是什么情况，这是负体积么d05b8f7a-a478-4841-8d82-e89bd4f376fb-image.png 290f4ec5-7312-4d2f-bc4a-9233470b86dd-image.png

简单的点线面操作撤销是没问题的，检查拓扑或者是操作包含的数量太多ICEM就没办法了，解决这个问题的最好办法就是你在做包含多步操作的问题时候先另存为一下是最保险的，问了一下学长他也是用这个方法的

@李东岳 好的，这里是一片净土

@李东岳 我这试了下system下创建fieldAverage文件并复制上述代码，删掉functions和fieldAverage1可以计算完成后处理，controDict不需要挂载，完整的文件内容如下

/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 9 \\/ M anipulation | ------------------------------------------------------------------------------- Description Calculates and writes the time averages of given list of fields. \*---------------------------------------------------------------------------*/ type fieldAverage; libs ( "libfieldFunctionObjects.so" ); writeControl writeTime; timeStart 0; timeEnd 100; resetOnOutput false; fields ( U { mean on; prime2Mean on; base time; } );

之前有点问题，最终如下
\par $C_p$采用4次多项式分段拟合
\begin{equation}\label{equ:NS_Cp}
C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
\end{equation}
\par静焓
\begin{equation}\label{equ:NS_H}
H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
\end{equation}
\par熵
\begin{equation}\label{equ:NS_S}
S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
\end{equation}
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equNSUs}
U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
\end{equation}
\begin{equation}\label{equNSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equNSTsta}
T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
\end{equation}
$\bullet$求解静压(已知总温、总压和静温)
\par由p等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equNS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation}

这个是哪国人的名字哟

是否为容差 matchTolerance 太大导致的衔接不上~

你好，想问一下你升阻力和文献对的上吗，我算的误差有12%

@Samuel-Tu 可行蓝色可以设置为internal面，顺便统计一下风量

@李东岳 @lsprxd @赵一铭 对于标准壁面函数我是这样理解的：在求解流场的时候还是用原来的k方程和epsilon方程来求解，只不过是将在壁面第一层网格（有可能是网格面上的值，有可能是网格中心的值）上的信息用壁面函数求解的结果来修正，而且壁面函数其实只求解壁面第一层网格的信息，第二层网格之后就是输运方程求解的结果了。不知我这样理解是否正确？

0_1478308139674_QQ图片20161105090838.png
公式18分行了还是显示不下。
大家都怎么处理的？用的是latex

0_1470732751048_QQ图片20160809165158.png

谢谢各位，解决了:mihu:

dos2unix命令是用来将DOS格式的文本文件转换成UNIX格式下的文本文件，这两个系统下的文本文件的差异在于换行标志，DOS下是用\r\n断行，UNIX下采用\n作为断行。将windows下面写好的一个文本文件转移到linux下，无法直接运行，打开会发现出现很多^M的标志，可以采用dos2unix命令进行清除。dos2unix命令起且仅起这个作用。因此，对于windows下生成的mesh文件，如果转移到linux下能够正常运行跑通算例，则完全不需要dos2unix命令。
同样，checkMesh也是，如果算例都能跑，则不需要这一步骤hhhh:papa:

@youmengtian
你用的什么显示器 :cheeky:

这个啊，先排错，
看看是“只有20可以”，还是“数字小了才可以”。
比如试试19，21，试试200。

错误定的越准，描述才越准，也越容易搜到前人的经验