nm_ldd_init_tttt_test
nm命令
nm
命令是一个用于查看二进制文件、库文件、可执行文件中符号信息的命令行工具。它通常用于分析程序的二进制文件，以获取关于程序结构的信息。
以OpenFOAM中的一个.o文件为例，OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/applications/solvers/DNS/dnsFoam/dnsFoam.o
用nm命令可以查看其中的符号，nm -C dnsFoam.o。
1 | 0000000000000000 V vtable for Foam::DimensionedField<double, Foam::volMesh> |
以下是 nm 命令的一些常用输出格式：
U：未定义符号，表示该符号在当前目标文件中被引用，但并未在当前目标文件中定义，需要在链接时从其他地方解析。T：文本段（代码段）符号，表示对应的符号在代码段中，通常是可执行程序的一部分。D：数据段符号，表示对应的符号在数据段中，通常是全局变量。B：BSS段符号，表示对应的符号在BSS段中，通常是未初始化的全局变量。W：弱符号，表示对应的符号是一个可以被重定位或被其他模块覆盖的符号。S：特殊符号，通常用于标记一些特殊的符号，比如函数的开始和结束。- 地址：符号的内存地址，可以是十六进制数。
- 符号名：函数名、变量名等。
TODO: 结合C++的type_trait，typeid等函数来解释符号表。引入函数签名等。
一个示例演示： 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
int main() {
int x = 10;
double y = 20.5;
// 使用typeid获取变量的类型信息
const std::type_info& typeInfoX = typeid(x);
const std::type_info& typeInfoY = typeid(y);
// 输出类型信息的名称
std::cout << "Type of x: " << typeInfoX.name() << std::endl;
std::cout << "Type of y: " << typeInfoY.name() << std::endl;
return 0;
}
这种获取类型信息的方式允许程序在运行时获取对象的实际类型，也就是RTTI(Run-Time Type Information)
编译后，运行结果： 1
2Type of x: i
Type of y: d
nm -C a.out查看其符号表可以发现，
1 | 0000000000003d80 V typeinfo for double@CXXABI_1.3 |
NOTE：实际在运用typeid的时候，一般会将其转换为人可读的（用cxxabi.h库来实现）。比如下面：
Ask gpt:
cxxabi.h是 C++ 标准库中的一个头文件，提供了用于操作 C++ ABI（Application Binary Interface，应用程序二进制接口）的一些函数。C++ ABI 定义了在不同编译器和平台下，C++ 编译器如何生成二进制代码以及如何在不同模块之间进行交互。这包括函数名的重整、异常处理、虚函数表等方面的规范。
cxxabi.h中的函数可以用于对符号（比如函数名）进行反解析，将符号的字符串表示形式转换为其内部的实际表示形式，或者反之。其中最常用的函数可能就是
abi::__cxa_demangle()，它用于将符号的字符串表示解析成人类可读的形式。举例来说，如果你有一个 C++ 函数名的字符串（例如 "_ZN5MyClass3fooEv"），你可以使用
__cxa_demangle()将其解析成对应的函数名（例如 "MyClass::foo()"）。总的来说，
cxxabi.h提供了一些工具，使得在运行时可以对 C++ 符号进行更灵活的处理，这对于某些需要动态获取或者修改函数信息的场景可能会非常有用。
所以，修改后的代码：
1 |
|
得到的结果： 1
2Type of x: int
Type of y: double
注意：不用的编译器，如gnu，clang，MSVC的typeid行为可能不用，详细可见： parallel101/course: 高性能并行编程与优化 - 课件 (github.com)
最近用nm干了啥
最近用nm来查看这个库中有没有某个符号（报错出现了undefined reference）。最后发现是代码的CMake找错了库
= =!
U：未定义符号，表示该符号在当前目标文件中被引用，但并未在当前目标文件中定义，需要在链接时从其他地方解析。
更新版：OpenFOAM在vscode的代码跳转设置
以前写过一些（理解还不太够），主要对法1的具体操作进行更新，
代码跳转的本质：找到这个文件
- 找到所有OpenFOAM中（编译后）的
lnInclude所在路径
find . -name "lnInclude" > tttt.findlnInclude
将OpenFOAM所有的lnInclude都找到。将所有的lnInclude链接都加在includePath中（仅仅加上Make/options下链接到的库还不够，代码中还是有不能跳转的地方）。
VScode: ctrl + shift + P， 搜索C++。
同理，代码不用编译都可以跳转，在windows下一样的。只要找到合适的lnInclude路径。
linux下OpenFOAM7的
lnInclude。/home/tttt/OpenFOAM前面替换一下即可。一个例子：
#include <typeinfo>是标准库中的头文件，要跳转这个，可以加上C++编译器的includePath，方法如下：
gcc -E -v - 查看这个命令打印信息中的这些：
1 | /usr/lib/gcc/x86_64-linux-gnu/11/include |
将这些粘贴到includePatch中（注意修改前缀）
1 | /home/tttt/OpenFOAM/OpenFOAM-7/src/regionModels/thermalBaffleModels/lnInclude/** |
Windows下（没编没有lnInclde，链到上一个路径也可以，注意windows和linux的文件分隔符），这个是没有编译，在windows中直接解压从GitHub下载的OpenFOAM包。
1 | D:\02_DevProject\00_OpenFOAM\OpenFOAM-7-master\src\regionModels\thermalBaffleModels\** |
设置后的效果：
类似的可以对Eigen进行同种操作
- Eigen代码跳转
设置： 1
/Users/tongyanjun/tttt_file/zz_lib_tttt/eigen-3.4.0/**
测试代码： 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
//这里eigen库可以成功跳转了
using Eigen::MatrixXd;
int main()
{
MatrixXd m(2,2);
m(0,0) = 3;
m(1,0) = 2.5;
m(0,1) = -1;
m(1,1) = m(1,0) + m(0,1);
std::cout << m << std::endl;
}
用doxygen生成类似OpenFOAM的继承关系图和接口手册
生成用的指令
doxygen -g Doxygen.config： 会在当前路径下生成一个Doxygen.configDoxygen.config文件里面有个选项HAVE_DOT推荐打开（不打开就是OpenFOAM那种类图，成员变量和函数较多的时候就不打开了），打开了里面有详细的函数。
更多设置就要研究这个Doxygen.config啦
doxygen Doxigen.config：当前路径生成了html- 如果在linux下
firefox ./html/index.html打开，windows直接网页打开，就可以看到了
生成后的（代码中按照一定的注释规范，doxygen会自动识别）：
这个图就和OpenFOAM的API Guide一样啦
The OpenFOAM Source Code Guide | OpenFOAM v7 | The OpenFOAM Foundation
以前刚开始学OpenFOAM，既不会C++也不会Doxygen，看API Guide简直是痛苦-0 -。 现在自己写个类试一试用doxygen生成一下，对阅读OpenFOAM的API Guide也有帮助（最近在写一个类，感觉类的设计简直是个规则类怪谈= =!，等哪天整理整理发出来乐呵一下）
通过API Guide来看OpenFOAM
以一个OpenFOAM求解器常用的类来看，createFields.H中的fvMesh。
1 | // createFields.H |
fvMesh 这个类继承了PrimitiveMesh
和polyMesh。
继续看PrimitiveMesh
和polyMesh可以发现在顶层的PrimitiveMesh，有许多函数定义为virtual xxx xxxx = 0;定义成了纯虚的接口，意味着实现应该去他的子类去找。当然，子类也可以用父类已经实现了的函数，所以有时在OpenFOAM找具体在哪儿实现的要套好几层。
生成doxygen手册测试用代码
1 | /** |
安装zsh以及oh my zsh
oh my zsh可以在WSL中安装，用vscode和windows自带的终端打开就是oh my zsh的界面！- 安装
oh my zsh之前先安装zsh
CuSO4Gem/ohmyzsh-offline: install oh my zsh offline (github.com)
学习在fluent和OpenFOAM中设置压力边界条件
update:
[CFD] When and Why do I need Operating Pressure, Temperature and Density? - YouTube
为什么需要Operating Pressue
补充关于压力出入口边界设置的资料：
ANSYS FLUENT 12.0 User's Guide - 7.3.3 Pressure Inlet Boundary Conditions (enea.it)
！！[CFD] Pressure-Inlet Boundary Conditions - YouTube
讲的超好！里面有对不同速度工况下（不可压，可压，其中可压又分为亚音速和超音速）压力入口的设置
一个视频理解Ansys Fluent里的压力设置及含义（表压or静压）_哔哩哔哩_bilibili
- Operating Condition
- Operating Pressure[Pa] +Gauge Pressure(Relative Pressure, Static Pressure) \(P_{absolute} = P_{operating}+P_{gauge}\) The gauge pressure is what we enter in the Boundary conditions panels. Set the gauge pressure is used to avoid problems with round-off errors which occur when the dynamic pressure differences in a fluid are small compared to the absolute pressure level
- Reference pressure location, is where the absolute static pressure is known
Reference pressure in OpenFOAM -- CFD Online Discussion Forums (cfd-online.com)
1 | In others the reference is set only if required. |
这个算例里面有句话值得注意：over-specify！(setting p and U at the same side)
ASK GPT:
In OpenFOAM, when setting up a boundary condition with a constant mass flow rate velocity, you usually need to also specify an appropriate pressure boundary condition to ensure the problem is well-posed and physically accurate. The choice of pressure boundary condition depends on the specific problem and flow conditions. Here are a few options you can consider:
Inlet/Outlet Pressure: You can use the
inletOutletboundary condition for pressure at the inlet or outlet. This condition allows the pressure to be determined based on the velocity field and the specified mass flow rate, ensuring that the flow is consistent with the desired mass flow rate.1
2
3
4
5
6inletOutlet
{
type inletOutlet;
inletValue uniform (0 0 0); // Specify the velocity components
value $internalField; // To set pressure equal to internalField
}Fixed Value or Fixed Gradient: Depending on the flow problem, you might be able to set a fixed value or fixed gradient for the pressure at the boundary. However, ensure that the boundary conditions do not lead to unphysical results.
1
2
3
4
5fixedValue
{
type fixedValue;
value $internalField; // Specify the desired pressure value
}1
2
3
4
5fixedGradient
{
type fixedGradient;
gradient uniform (0 0 0); // Specify a gradient if needed
}Zero Gradient: For some cases, you might use a zero gradient condition for pressure at the boundary. This assumes that the pressure gradient is zero, which might be suitable if the boundary is far from the region of interest or if you have a long enough domain.
1
2
3
4zeroGradient
{
type zeroGradient;
}Pressure Inlet: If you have a known pressure value at the inlet, you can directly specify that value. However, this might affect the mass flow rate, so ensure consistency with your velocity boundary condition.
1
2
3
4
5
6pressureInlet
{
type pressureInlet;
p uniform 101325; // Specify the pressure value
U uniform (0 0 0); // Specify the velocity components
}
Remember to select a pressure boundary condition that is consistent with your desired mass flow rate velocity condition and that reflects the physics of your problem. Always validate your choice through simulations and comparisons with experimental or analytical results if available.
c_cpp_properties的设置
ubuntu 下 GCC/G++ 的 include 搜索路径查看与设置_ubuntu g++ clude_许野平的博客-CSDN博客
echo 'main(){}'|gcc -E -v -
1 | #include <...> search starts here: |
OpenFOAM和fluent的SIMPLE计算设置对比【还没写完】
预备知识
- (Dynamic viscosity)动力粘度：\(\mu\) \(，这个\)\(\mu\)与牛顿内摩擦公式有关：\(\tau = \mu du/ dy\)， 单位为\(Pa\cdot s\)或\(N\cdot s/m^2\)
- (Kinematic viscosity)运动粘度：\(\nu = \mu / \rho\) ，单位为\(m^2/s\)
在OpenFOAM中，一般设置运动粘度： 1
nu [ 0 2 -1 0 0 0 0 ] 0.01;
theory.pdf (dyfluid.com) 对于simpleFoam来说，OpenFOAM求解的方程是这样的：
$$ \[\begin{align} \nabla\cdot\bf{U}=0 \\ \nabla \cdot (\mathbf{U}\mathbf{U})-\nabla \cdot(\nu \nabla \mathbf{U})=-\nabla p \end{align}\] $$ 对于不可压流体时，OpenFOAM对上述方程整个除了密度\(\rho\) 。所以会在边界条件的p文件内看到p的量纲是\(m^2/s^2\)：
1 | dimensions [0 2 -2 0 0 0 0]; |
压力\(p\) 原量纲为：\(kg / m\cdot s^2\) ，除以密度\(kg/m^3\)
\(k-\epsilon\) 模型
对于求解采用\(k-\epsilon\) 模型时，湍流粘度\(\nu_{t}\) 与湍动能\(k\) 和湍流耗散率\(\epsilon\) 结合，有公式： \[ \nu_t = C_{\mu} \frac{k^2}{\epsilon} \] 标准 k-ε 模型： 常数： \[ \begin{align} C_{mu} = 0.0 \\ C_1 = 1.44 \\ C_2 = 1.92 \\ α_k = 1 \\ α_ε = 0.76923 \end{align} \]
Turbulence dissipation rate -- CFD-Wiki, the free CFD reference (cfd-online.com) 湍流耗散率单位：\(m^2/s^3\)
关于边界flowRateInletVelocity
flowRateInletVelocity
value uniform (0 0 0) -- CFD Online Discussion Forums
(cfd-online.com) OpenFOAM:
API Guide: flowRateInletVelocityFvPatchVectorField Class Reference -
Example of the boundary condition specification for a volumetric
flow rate: 1
2
3
4
5
6
7
8<patchName>
{
type flowRateInletVelocity;
volumetricFlowRate 0.2;
extrapolateProfile yes;
value uniform (0 0 0);
}
example of the boundary condition specification for a mass flow rate:
1
2
3
4
5
6
7
8
9<patchName>
{
type flowRateInletVelocity;
massFlowRate 0.2;
extrapolateProfile yes;
rho rho;
rhoInlet 1.0;
value uniform (0 0 0);
}检测入口和出口处的流量相同的 Difference between specified and calculated mass flow rate -- CFD Online Discussion Forums (cfd-online.com)
边界条件的\(k-\epsilon\) 模型该怎么取
calculation of k, epsilon and omega -- CFD Online Discussion Forums (cfd-online.com)
Turbulence free-stream boundary conditions -- CFD-Wiki, the free CFD reference (cfd-online.com)
- 看Henry等人写的关于这部分的内容
Notes on CFD: General Principles - 7.3 Inlet turbulence
- 看这里的估算方法
测试数学公式
LaTex: $$ \[\begin{align} \int_{V_p} \nabla \cdot (\rho \mathbf{U} \phi)dV & = \sum_f {\mathbf{S} \cdot (\rho \mathbf{U} \phi )_f } \notag \\ & = \sum_f \mathbf{S} \cdot (\rho \mathbf{U} )_f \phi_f \notag \\ & = \sum_f{F \phi_f} \label{Convection term-1} \end{align}\] $$
Markdown:
\[ \int_{V_p} \nabla \cdot (\rho \mathbf{U} \phi)dV & = \sum_f {\mathbf{S} \cdot (\rho \mathbf{U} \phi )_f } \notag \\ & = \sum_f \mathbf{S} \cdot (\rho \mathbf{U} )_f \phi_f \notag \\ & = \sum_f{F \phi_f} \]
\[ \phi_c = \phi_f+\psi \]
测试结果：LaTeX公式可以过，markdown多行的解析还有问题