Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
@SnoWolf2018 @玉玺
多谢!!
@李东岳 :140: ,效果很好。还能把点调大点,调成正方形,然后就可以没有空隙了。pointtype, pointsize这些
@李东岳 感觉好像能行了:140: 感谢李老师!
@QuietCheng 结果怎么样了?
代码中的描述的作用是为了使求解稳定,在原有格式的基础上减去了 fvm::Sp(fvc::surfaceIntegrate(faceFlux), vf) ,而计算收敛的时候该项为0,所以理论上应该不会影响精度。
工科末流985硕士毕业,想在武汉找CFD相关的工作,家人们有没有推荐的公司。
@王孟飞 这个问题我也遇到了,不知道怎么解决
@紫色星光 在设置节点数的小窗口,往下拉,有个copy parameter,不勾选。如果不勾选还会联动,那就是切块的问题了。
谢谢各位,解决了:mihu:
大家好! 我想创建两个等腰直角三角形的Y型网格划分,2d不能直接创建Y-BLOCK,只查到了先O型划分再移动节点,间接建立Y Grid,但是这样就会影响我相邻块的关联,有什么办法能实现其划分吗?如下图,红色区域就是两个直角等腰三角形,其相邻的块都是矩形块,我用O型划分了一个,但是周围的块的关联已经很奇怪了。在此先谢过各位前辈了!
cfd2fd2a-2a3a-457e-90da-8343185b0c87-image.png
@chengan-wang 或许你可以搞一个局部细加密。
之前有点问题,最终如下 \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}
@东岳 谢谢老师,我已经找到网格生成规律了,从左至右,由下至上地规律依次生词。满足误之前需求工作,reNumberMesh有需要我再去琢磨吧,再谢!
@lyc 你好,我很少遇到。即便遇到都是重启就能解决的那种。
不知道你遇见过没,第一个值太小比如-15次方,基本不以收敛,或者在多孔介质壁面上云图不连续,?
这个是不是用到了基于物理量+动网格技术。。。
@东岳 万分感谢!
@东岳 后来就只用snappyHexMesh画的网格做了一些简单的模拟……:zoule:
@CYW 谢谢。在interFoam里加入过能量方程,但是在加物理量比如普康特数啥的出现了点问题。。。我对照着compressibleInterFoam研究研究,,有问题再向你请教。不好意思没帮到你什么忙
东岳大神、各路大神好: 最近在模拟一个上下振动混合器中的单相三维流场,如图。 考虑到计算时间和成本问题,想采用moving frame of reference的方法来代替动网格进行振动的模拟。目前设置inlet速度为振动速度,壁面和出口都设置的是wall,速度也是振动速度,而板和轴的速度为0.但得到的流场似乎不太对,想请教下边界条件的设置存在哪些问题?WeChat Image_20190929112939.jpg