CFD中文网

@李东岳 最后还是用了后处理的方法进行计算，又要做实验又要做模拟有点赶，用的是老板之前研究的方法几何平均粒径、几何标准差与浓度的关系，算了浓度再来得到颗粒数关系。不过布朗运动要加进去我感觉得用群体平衡方法。

@卡洛 请问你是用的哪种方法？

之前有点问题，最终如下
\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}

@bestucan 明白了，感谢指点 :）

是做巷道通风仿真么，带条缝的附壁风筒？

@李东岳
已经看到了，非常感谢您:xiexie:

启动的时候在启动命令后面加上后缀试试

-driver x11

@李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
这篇文章对三种重构方法做了对比。

你是不是想把曲面展开为一个平面？https://blog.csdn.net/mrbaolong/article/details/106949340

@wwzhao

简单的感受了一下两种赋值方式，效率确实很大提高。

#include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07

这个软件使用的太少了。好像论坛也没什么人可以讨论 :134:

@HaoC51547 不需要 原来是多少就是多少

@wsy11 谢谢您~我现在已经可以成功构建网格啦^^

@李东岳 我就记得我之前搞过，后来坑太深了，就没搞了。

这些涉及到图形处理的各个层，又是平台不同，又是商业公司夹私货，那些个组件的名字本来就乱七八糟，缩写规则还不一样。想短时间了解就更乱。

大概就是，可以软件渲染（OSMesa），可以硬件渲染（EGL）。
其实都是基于Open GL，就这三个英语词空格隔开一起搜都一大堆讲更多术语关系的文章。

paraview 默认用软件渲染，想要硬件渲染把 EGL 编译选项打开就行了，这就有个忘了打开报错的：
https://discourse.paraview.org/t/build-paraview-for-nvidia-gpu/8309
但这个是单GPU的，不知道多GPU的怎么搞，

我的GPU是当显卡用的，一开paraview就调用GPU，当然是为“图形界面渲染"的,不是为“paraview渲染”的，这俩还都是opengl下面的，那些名词太乱了。

我试着编了一下 pvserver能多线程跑在GPU上，
Screenshot from 2022-06-11 15-52-46.png
就是这个方法：
https://github.com/Kitware/ParaView/blob/master/Documentation/dev/build.md

当然编译之前已经安装好nvidia自带的驱动了。

nvidia 早些年自己出过 paraview 的插件：
https://www.nvidia.com/en-us/data-center/index-paraview-plugin/
不过是针对 paraview 5 的

solidworks导入icem时候有些地方确实有点问题，我现在先用sw把大框架画出来，自己在icem里面再编辑一下小框架，最后再切分一下，有点麻烦。你现在已经解决把，有啥好办法吗

@李东岳 谢谢李老师

他这个聊天记录我没太懂。什么是拖布结构？

好专业好专业，看的我蒙蒙的

:146: :146: :146:

@李东岳 好嘞，谢谢老师解惑

请问大家，有个ICEM19.4导出的.blk文件，在新版本（2022）中导入显示版本不匹配，但是.tin文件就可以，有没有什么好办法可以查看.blk文件？