#include "udf.h"
#include "unsteady.h"
#include "profile.h"
#include <stdio.h>
#include <stdlib.h>

#define MAX_POINTS 4000
#define MAX_TIMESTEPS 997
#define FILENAME "velocity_data.csv"

typedef struct {
double x, y, z;
double u[MAX_TIMESTEPS];
} PointData;

void read_data(PointData* points) {
FILE* file = fopen(FILENAME, "r");
if (file == NULL) {
printf("Error: Could not open file.\n");
exit(1);
}

for (int i = 0; i < MAX_POINTS; i++) { fscanf(file, "%lf,%lf,%lf", &points[i].x, &points[i].y, &points[i].z); for (int j = 0; j < MAX_TIMESTEPS; j++) { fscanf(file, ",%lf", &points[i].u[j]); } } fclose(file);

}

DEFINE_PROFILE(unsteady_velocity, thread, position)
{
real x[ND_ND];
real t = CURRENT_TIME;
int time_step = N_TIME;
PointData points[MAX_POINTS];
face_t f;

read_data(points); begin_f_loop(f, thread) { F_CENTROID(x, f, thread); for (int i = 0; i < MAX_POINTS; i++) { if (x[0] == points[i].x && x[1] == points[i].y && x[2] == points[i].z) { F_PROFILE(f, thread, position) = points[i].u[time_step]; break; } } } end_f_loop(f, thread)

}

1710316733826.jpg

使用滑移网格fluent报错：Note: zone-surface: cannot create surface from sliding interface zone.

x3d
用这个导出，导入到blender

不知道当时的问题是否解决了，我去年做的一部分工作里面也需要相关的处理。当时采用的Matlab对导出的数据进行处理，不过在数据处理之前需要将 x,y当中一些极小的值变成0，比如10的负十几次方的值，因为这些fluent输出的值可能会对处理产生影响，然后导入程序就可以变成矩阵形式了，程序大概如下。供以后有需要的朋友们参考。

data1 = xlsread('D:\User.....');

% 获取X、Y坐标和密度值的向量
X = data1(:, 1);
Y = data1(:, 2);
density = data1(:, 3);

% 确定X和Y的唯一值，并按升序排列
unique_x = unique(X);
unique_y = unique(Y);

% 创建行列编号映射表
x_mapping = containers.Map(unique_x, 1:length(unique_x));
y_mapping = containers.Map(unique_y, 1:length(unique_y));

% 初始化矩阵
num_rows = length(unique_y);
num_cols = length(unique_x);
densityMatrix1 = zeros(num_rows, num_cols);

% 将密度值填充到矩阵中
for i = 1:length(X)
row = y_mapping(Y(i));
col = x_mapping(X(i));
densityMatrix1(row, col) = density(i);
end

a80fed28-45ef-4507-a08d-13abf19c7ccc-image.png 这个udf具体该怎么编写和使用？

loop的是面网格上的所有网格。你的需求应该能实现，最简单的办法是，每个面网格中心去读一遍文档，把散点距离自己最近的那个点的速度赋给所在网格。但是会有其他问题，散点疏密是不是和面网格疏密比较匹配，也许一个面网格配多个散点或者反之。

想设置一个壁面接触角随时间增加减小，用公式描述
QMF9EZM9NI7B{$AS17ASF{M.jpg
但是模拟过程中不到0.01s就已经完全铺开了
8d0dfcc4-1ada-4bec-8a31-d4000f59d818-0~FKZ}KBX0WEF3X7@2%@NTD.png
接触角为定值的时候没有问题，请问接触角随时间增加减小的设置哪里出了问题，该如何设置

用Fluent大涡模拟计算充分发展的槽道湍流，流向和展向为周期性边界条件，上下为壁面，如何设置流体进出口呢？看了之前的一些帖子，需要进口设定一定压力梯度吗？如何定义压力梯度啊

替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程，现在模型验证一直没弄好。首先是编写了一个简单的UDF，采用热平衡模型，导入fluent后提醒我定义比热的那一块有错误，求各位前辈帮我看一下；其次，我添加了Boussinesq项后，优化了一点点液相云图中熔融线一直平行的问题，但是与论文中的还是有差距，且我的熔融线上下厚度保持一致（实际上应该下面厚很多），UDF中还需要添加什么项来解决这个问题呢？以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图，求各位帮我解惑，不甚感激。2.png 3.pn替代文字g 4.png

我想建立一个水槽，水流从水槽一边冲过去，然后又从另一边冲过来，两边的入口和出口边界要随时间变化，有没有大神指导一下，怎么实现？谢谢！

请参考无痛苦ns方程笔记的这一节：R、sigma、prime2Mean、以及雷诺应力相关量

不清楚商软那面怎么回事，本身我也没太关注这个，我查了一下，OpenFOAM本身这面的streamline看起来是你说的patiline：

This function object generates streamline data by sampling a set of user-specified fields along a particle track, transported by a user-specified velocity field.

后处理软件paraview的streamline是如何定义的还没找到。

关闭超线程，选择msmpi并行方式，计算速度明显改善。但感觉还有很大的提升空间，目前先这样将就着用了

我模拟的对象是反应堆内的稳态流动传热，由于堆芯内组件数量庞大，精细建模复杂，因此采用多孔介质模型对堆芯区域进行简化。

为了验证这种等效的合理性和误差，我先对单个组件栅元进行多孔介质简化，并将多孔介质模拟结果与未简化的模型结果对比验证多孔介质模型的合理性。

单个组件栅元的截面图如下：
图片1.png

其二维轴对称模型如下图：
f510f5de-0d5e-4389-b577-ad2d3defe29c-image.png

模拟结果
(1)未采用多孔介质模型
我用了二维轴对称模型来模拟，设置固体流体内热源，入口质量流量和出口压力边界条件，流体和固体分别加载体积热源P1和P2（W/m3），流体和固体的导热系数均为常数，得到的稳态温度场和压损分布如下：
bca065f1-7cc4-4385-8113-d7e0acd6d23e-image.png
流体出口温度700℃，固体由于带内热源温度高于流体，固体最高温度为799℃，固体平均温度为756℃，组件压损为2400Pa

(2)采用多孔介质模型模拟
由于固体带内热源，流固存在较大温差，采用非热平衡模型
根据栅元几何结构指定：孔隙率γ=0.1 ，界面换热系数hfs=5000W/M2K（与未采用多孔介质模型算得的换热系数一致），界面面积密度Afs=6.71m-1

通过流速-压力实验关系确定粘性损失项系数和惯性损失项系数，只考虑x方向的流动损失，y方向系数放大100倍

别给多孔介质流体和多孔介质固体设置了内热源，分别为P1×孔隙率和P2×（1-孔隙率），初始化后检查内热源加载总量无误；流体固体的热导率输入值和上个模型一致

得到的模拟结果：
af1efa10-a95f-4e49-a69c-63b4b4357ca0-image.png
采用多孔介质非热平衡模型模拟得到的流体温度是正确的，压损也是正确的，但固体温度分布明显不对，固体域温度热点仅为725℃，平均温度为711℃，均远小于实际模型算得的固体温度值。

我自己的猜测：
多孔介质模型将单组件变为流固均匀介质后，不能模拟固体区域沿导热路径带来的固体温升？刚好我这个栅元模型中固体由于导热产生的温升非常高，如果是这样，说明多孔介质非热平衡模型不适用于我所研究的问题，因为简化后会显著低估固体温度热点。

请问我的理解是正确的么

参数波动跟监测点位置有关，放在分离区跟放在直管中间是不一样的，另外，感觉稳态计算是监测面上的通量吧，通量守恒就可以认为收敛了吧。

@黄昏晓为谁破 好的感谢前辈回复

请问我想使用dpm模型模拟液滴的凝固，应该选什么类型得到颗粒呢？我看几种类型的粒子没有考虑凝固过程，求教！！！

@benf 不客气不客气，欢迎引用我们组的论文~这个模型有外国组已经验证过了，所以如果你有兴趣编应该可以原样复现

我觉得应该搞清楚液滴在真空状态下的变化，然后针对性的设置环境条件

格式不对噻