我想把图片中的红色框的一部分输出,然后在paraview中做后处理,为了方便,我建立了一个同样的参数(黄色框),想用类似湍流模型里边的方法试一下,但是写不出来
我还尝试了另一种方法:在createfield.h里面定义了一个参数,这样会写出一个文件,但是写出的文件不会更新。
其中m是自定义的一个参数
benqing
帖子
-
-
@ct 没有,这个不该出现负数的吧
-
@ct 后来我直接用的k作为模化的湍动能,用UPrime2Mean的结果计算瞬时的resloved湍动能
-
@benqing 解决了,只算一个就可以,算流量比较方便
Q=A×V,A=(Ly-h)×Lz=(3h-h)×4.5h=9h×h,V=1,Q=9×0.028×0.028×1=0.007056,所以我的结果跟@心里的小漩涡 的结果不太一样,后续看看计算结果 -
@心里的小漩涡 小伙伴还在论坛吗?最近有兴趣用fluent跑跑这个case,但是这里面对周期性边界上mass flow rate和pressure Gradient设定不是很理解,这个数值是如何算出来的?
-
@李东岳 在 雷诺应力与URANS 中说:
我一直在找相关的sci,你这面有文章见过类似的操作么
URANS的话,没有见到过,这也只是我自己的想法。我用过PANS模型,有一些文献是将两者相加的,其实也没有理论的介绍,就是后处理的方法。
我找找看有没有URANS相关的文献! -
@李东岳 两者相加吧,我认为若URANS可以较好的解析流场的话,uPrime2Mean会占更大比重,但是不能忽略R。理论方面的话,可以参考PANS模型。
-
这个问题我解决了,宏C_WALL_DIST在高版本里面确实不能用了,从CFDonline上找到了一个可以用的代码,代码如下:
#include "udf.h" #include "prox.h" static cxboolean wall_dist_set = FALSE;//cxboolean布尔运算,与或非,对错 DEFINE_ON_DEMAND(set_wall_dist_udm0) { #if !RP_HOST//表示非host节点的话,执行后面的代码,fluent并行的话,会启动一个 host节点和指定个数的node节点,host节点主要处理GUI和结果收集之类的任务node节点负责分区计算任务 Domain *domain; Thread *t; cell_t c; if (! wall_dist_set) { domain = Get_Domain(1);//ROOT_DOMAIN_ID? Alloc_Storage_Vars(domain, SV_RTMP_0, SV_NULL);//Alloc_Storage_Vars 用于更改分配 Calc_Cell_Wall_Distance_New(domain, SV_RTMP_0); thread_loop_c(t,domain) { begin_c_loop(c,t) { C_UDMI(c,t,0) = C_TMP0(c,t); } end_c_loop(c,t) } wall_dist_set = TRUE; } #endif /* !RP_HOST */ } DEFINE_ON_DEMAND(reset_udm0) { #if !RP_HOST Domain *domain; Thread *t; cell_t c; domain = Get_Domain(1); thread_loop_c(t,domain) { begin_c_loop(c,t) { C_UDMI(c,t,0) = 0.0; } end_c_loop(c,t) } #endif /* !RP_HOST */ wall_dist_set = FALSE; }
说实话这个代码我不能完全看懂,但是主要思想就是用DEFINE_ON_DEMAND定义两个命令reset_udm0和set_wall_dist_udm0,前者将壁面到最近壁面的距离reset为0,后者可以求出想要的结果,经过测试两个命令应该都是必须的,另外我放的代码和CFDonline上有一点不同,domain = Get_Domain(1);这一行在原代码里面是domain = Get_Domain(ROOT_DOMAIN_ID);注意修改,下图是我测试的一个水翼网格的结果,结果是对的,
再啰嗦几句代码使用方法,这个代码直接复制到一个.c文件中,在fluent中添加memory的数量为2,编译.c文件并load,在calculate之前,点Execute on Demand,并先后运行reset_udm0和set_wall_dist_udm0,然后就可以calculate了
-
@bestucan 谢谢,这个问题解决了,需要设置UDS的个数,不论有没有标量方程
-
你们的问题解决了吗,我遇到了同样的问题
-
想在fluent中用UDF实现,非线性雷诺应力,能从CFDonline上找到一个相关的代码,CFDonline链接,我在fluent中能够编译成功,但是在开始计算的时候出现问题,报错:
============================================================================== Node 0: Process 13836: Received signal SIGSEGV. ============================================================================== MPI Application rank 0 exited before MPI_Finalize() with status 2 The fl process could not be started.
测试发现是C_UDSI这部分的问题,在我的理解中UDS是用于建立标量方程的,这个里面采用UDSI计算了UU等雷诺应力分量,有考虑在UDS对话框中添加UDS,但是对这些量并没有输运方程,不知道应该如何处理,希望有人帮忙解决,多谢!
代码贴在下面/* turbulence model: nonlinear k-epsilon turbulence model by Ehrhard (1999) */ /* UDF (User Defined Function) Routine for FLUENT */ /* see also the paper: */ /* Implementation and comparison of different turbulence models for three dimensional wall jets with Fluent */ /* FLUENT CFD Forum 2005 Bad Nauheim, Germany */ /* Heschl Ch*., Sanz W.**, Klanatsky P.* */ /* * FachhochschulstudiengçÏge Burgenland GmbH, SteinamangerstraŽße 21, A-7423 Pinkafeld */ /* ** TU Graz, Institut of Thermal Turbomachinery and Machine Dynamics, Inffeldgasse 25, A-8010 Graz */ /* Contact: christian.heschl@fh-burgenland.at, wolfgang.sanz@tugraz.at, peter.klanatsky@fh-burgenland.at */ /* Die Autoren ?ernhemen keine Verantwortung oder Haftung f? den Inhalt der unten angef?rten UDF-Routinen */ #include "udf.h" #include "math.h" /* Turbulence model constants */ const real C_1=-0.2; const real C_2=0.4; const real C_5=0.0; const real CE_1=1.44; const real CE_2=1.92; const real SIG_K=1.0; const real SIG_E=1.3; const real C_MU_ke=0.09; const real kappa=0.42; const real E_log=9.793; /* User-defined scalars */ enum { UU, VV, WW, UV, UW, VW }; DEFINE_SOURCE(u_source, c, t, dS, eqn) { real source; dS[eqn]=0.0; source= - C_R(c,t) * ( C_UDSI_G(c,t,UU)[0] + C_UDSI_G(c,t,UV)[1] + C_UDSI_G(c,t,UW)[2] );//C_R(c,t) = rho density return source; } DEFINE_SOURCE(v_source, c, t, dS, eqn) { real source; dS[eqn]= 0.0; source= - C_R(c,t) * ( C_UDSI_G(c,t,UV)[0] + C_UDSI_G(c,t,VV)[1] + C_UDSI_G(c,t,VW)[2] ); return source; } DEFINE_SOURCE(w_source, c, t, dS, eqn) { real source; dS[eqn]= 0.0; source= - C_R(c,t) * ( C_UDSI_G(c,t,UW)[0] + C_UDSI_G(c,t,VW)[1] + C_UDSI_G(c,t,WW)[2] ); return source; } DEFINE_TURBULENT_VISCOSITY(user_mu_t,c,t) { real S,W,c_mu,T,mu_t; real S11, S12, S13, S21, S22, S23, S31, S32, S33; real W11, W12, W13, W21, W22, W23, W31, W32, W33; S11 = 0.5*( C_DUDX(c,t)+C_DUDX(c,t) ); S12 = 0.5*( C_DUDY(c,t)+C_DVDX(c,t) ); S13 = 0.5*( C_DUDZ(c,t)+C_DWDX(c,t) ); S21 = 0.5*( C_DVDX(c,t)+C_DUDY(c,t) ); S22 = 0.5*( C_DVDY(c,t)+C_DVDY(c,t) ); S23 = 0.5*( C_DVDZ(c,t)+C_DWDY(c,t) ); S31 = 0.5*( C_DWDX(c,t)+C_DUDZ(c,t) ); S32 = 0.5*( C_DWDY(c,t)+C_DVDZ(c,t) ); S33 = 0.5*( C_DWDZ(c,t)+C_DWDZ(c,t) ); W11 = 0.5*( C_DUDX(c,t)-C_DUDX(c,t) ); W12 = 0.5*( C_DUDY(c,t)-C_DVDX(c,t) ); W13 = 0.5*( C_DUDZ(c,t)-C_DWDX(c,t) ); W21 = 0.5*( C_DVDX(c,t)-C_DUDY(c,t) ); W22 = 0.5*( C_DVDY(c,t)-C_DVDY(c,t) ); W23 = 0.5*( C_DVDZ(c,t)-C_DWDY(c,t) ); W31 = 0.5*( C_DWDX(c,t)-C_DUDZ(c,t) ); W32 = 0.5*( C_DWDY(c,t)-C_DVDZ(c,t) ); W33 = 0.5*( C_DWDZ(c,t)-C_DWDZ(c,t) ); T=C_K(c,t)/C_D(c,t);//C_K(c,t) = turbulent kinitic energy(TKE); C_D(c,t) = epsilon S=T*sqrt( 2*(S11*S11+S12*S12+S13*S13+S21*S21+S22*S22+S23*S23 +S31*S31+S32*S32+S33*S33) ); W=T*sqrt( 2*(W11*W11+W12*W12+W13*W13+W21*W21+W22*W22+W23*W23 +W31*W31+W32*W32+W33*W33) ); c_mu=MIN( 1./(0.9*pow(S,1.4)+0.4*pow(W,1.4)+3.5) , 0.15 ); mu_t=c_mu*T*C_K(c,t); return mu_t; } DEFINE_ADJUST(rsm_adjust,domain) { Thread *t; cell_t c; real T; real mu_t, c_mu; real S11, S12, S13, S21, S22, S23, S31, S32, S33; real W11, W12, W13, W21, W22, W23, W31, W32, W33; real S, W, Sij, Wij; real C_3, C_4, C_6, C_7; real P11, P22, P33, Pkk; real tau_w, u_tauw, u_star, y_star, Gk, u_mag; /* Set the turbulent viscosity */ thread_loop_c(t,domain) if (FLUID_THREAD_P(t)) { begin_c_loop(c,t) { S11 = 0.5*( C_DUDX(c,t)+C_DUDX(c,t) ); S12 = 0.5*( C_DUDY(c,t)+C_DVDX(c,t) ); S13 = 0.5*( C_DUDZ(c,t)+C_DWDX(c,t) ); S21 = 0.5*( C_DVDX(c,t)+C_DUDY(c,t) ); S22 = 0.5*( C_DVDY(c,t)+C_DVDY(c,t) ); S23 = 0.5*( C_DVDZ(c,t)+C_DWDY(c,t) ); S31 = 0.5*( C_DWDX(c,t)+C_DUDZ(c,t) ); S32 = 0.5*( C_DWDY(c,t)+C_DVDZ(c,t) ); S33 = 0.5*( C_DWDZ(c,t)+C_DWDZ(c,t) ); W11 = 0.5*( C_DUDX(c,t)-C_DUDX(c,t) ); W12 = 0.5*( C_DUDY(c,t)-C_DVDX(c,t) ); W13 = 0.5*( C_DUDZ(c,t)-C_DWDX(c,t) ); W21 = 0.5*( C_DVDX(c,t)-C_DUDY(c,t) ); W22 = 0.5*( C_DVDY(c,t)-C_DVDY(c,t) ); W23 = 0.5*( C_DVDZ(c,t)-C_DWDY(c,t) ); W31 = 0.5*( C_DWDX(c,t)-C_DUDZ(c,t) ); W32 = 0.5*( C_DWDY(c,t)-C_DVDZ(c,t) ); W33 = 0.5*( C_DWDZ(c,t)-C_DWDZ(c,t) ); T=C_K(c,t)/C_D(c,t); S=T*sqrt( 2*(S11*S11+S12*S12+S13*S13+S21*S21+S22*S22+S23*S23 +S31*S31+S32*S32+S33*S33) ); W=T*sqrt( 2*(W11*W11+W12*W12+W13*W13+W21*W21+W22*W22+W23*W23 +W31*W31+W32*W32+W33*W33) ); Sij=S11*S11+S12*S12+S13*S13+S21*S21+S22*S22+S23*S23+S31*S31+S32*S32+S33*S33; Wij=W11*W11+W12*W12+W13*W13+W21*W21+W22*W22+W23*W23+W31*W31+W32*W32+W33*W33; c_mu=MIN( 1./(0.9*pow(S,1.4)+0.4*pow(W,1.4)+3.5) , 0.15 ); mu_t=c_mu*T*C_K(c,t); C_3=2.0-exp(-SQR(S-W)); C_4=-32.0*SQR(c_mu); C_6=-16.0*SQR(c_mu); C_7=16.0*SQR(c_mu); C_UDSI(c,t,UU)= C_1*mu_t*T* (S11*S11+S12*S12+S13*S13-1./3.*Sij) + 2.*C_2*mu_t*T*(W12*S12+W13*S13) + C_3*mu_t*T*(W12*W12+W13*W13-1./3.*Wij) - 2.*C_4*mu_t*T*T * (W12*(S11*S12+S12*S22+S13*S23)+W13*(S11*S13+S12*S23+S13*S33)) + C_6*mu_t*T*T* (S11*Sij) + C_7*mu_t*T*T* (S11*Wij);//User-Defined Scalar(UDS) cell variables C_UDSI(c,t,VV)= C_1*mu_t*T* (S12*S12+S22*S22+S23*S23-1./3.*Sij) + 2.*C_2*mu_t*T*(W23*S23-W12*S12) + C_3*mu_t*T*(W12*W12+W23*W23-1./3.*Wij) + 2.*C_4*mu_t*T*T * (W12*(S11*S12+S12*S22+S13*S23)-W23*(S12*S13+S22*S23+S23*S33)) + C_6*mu_t*T*T* (S22*Sij) + C_7*mu_t*T*T* (S22*Wij); C_UDSI(c,t,WW)= C_1*mu_t*T* (S13*S13+S23*S23+S33*S33-1./3.*Sij) - 2.*C_2*mu_t*T*(W13*S13+W23*S23) + C_3*mu_t*T*(W13*W13+W23*W23-1./3.*Wij) + 2.*C_4*mu_t*T*T * (W13*(S11*S13+S12*S23+S13*S33)+W23*(S12*S13+S22*S23+S23*S33)) + C_6*mu_t*T*T* (S33*Sij) + C_7*mu_t*T*T* (S33*Wij); C_UDSI(c,t,UV)= C_1*mu_t*T* (S11*S12+S12*S22+S13*S23) + C_2*mu_t*T*(W12*(S22-S11)+W13*S23+W23*S13) + C_3*mu_t*T*(W13*W23) + C_4*mu_t*T*T* (W12*(S11*S11+S13*S13-S22*S22-S23*S23)-W13*(S12*S13+S22*S23+S23*S33)-W23*(S11*S13+S12*S23+S13*S33)) + C_6*mu_t*T*T* (S12*Sij) + C_7*mu_t*T*T* (S12*Wij); C_UDSI(c,t,UW)= C_1*mu_t*T* (S11*S13+S12*S23+S13*S33) + C_2*mu_t*T*(W13*(S33-S11)+W12*S23-W23*S12) -C_3*mu_t*T*(W12*W23) + C_4*mu_t*T*T* (W13*(S11*S11+S12*S12-S23*S23-S33*S33)-W12*(S11*S13+S22*S23+S23*S33)+W23*(S11*S12+S12*S22 +S13*S23)) + C_6*mu_t*T*T* (S13*Sij) + C_7*mu_t*T*T* (S13*Wij); C_UDSI(c,t,VW)= C_1*mu_t*T* (S12*S13+S22*S23+S23*S33) + C_2*mu_t*T*(W23*(S33-S22)-W12*S13-W13*S12) + C_3*mu_t*T*(W12*W13) + C_4*mu_t*T*T* (W23*(S12*S12+S22*S22-S13*S13-S33*S33)+W12*(S11*S13+S12*S23+S13*S33)+W13*(S11*S12+S12*S22+S13*S23)) + C_6*mu_t*T*T* (S23*Sij) + C_7*mu_t*T*T* (S23*Wij); C_UDMI(c,t,0)= -2.*mu_t*S11 + 2./3.*C_K(c,t) + C_UDSI(c,t,UU); C_UDMI(c,t,1)= -2.*mu_t*S22 + 2./3.*C_K(c,t) + C_UDSI(c,t,VV); C_UDMI(c,t,2)= -2.*mu_t*S33 + 2./3.*C_K(c,t) + C_UDSI(c,t,WW); C_UDMI(c,t,3)= -2.*mu_t*S12 + C_UDSI(c,t,UV); C_UDMI(c,t,4)= -2.*mu_t*S13 + C_UDSI(c,t,UW); C_UDMI(c,t,5)= -2.*mu_t*S23 + C_UDSI(c,t,VW); } end_c_loop(c,t) } }
-
@yfleung 这个具体没写过,找找定义模仿现有程序试试吧
-
@yfleung postChannel这个程序很好改,你可以自定义添加量,比如k和epsilon等等;求utau的话,需要壁面剪切应力(tau_wall),u_tau=sqrt(tau_wall/密度)--其实在of中不需要除以密度,计算壁面剪切应力用wallShearStress,但是看代码的话,wallShearStress采用的速度是瞬时速度,其实这个工具适合计算定常,如果非定常的话,我修改过程序采用的UMean;这个时候得到的是面上的数据,你自己再做合适的后处理就行了
-
@yfleung 对于槽道流有一个非常好用的工具,叫postChannel,在channel395那个算例能用,这个工具可以做空间(槽道流中应该是流向和展向)的平均,比较容易用,具体你可以自己下看下,如果有问题可以再问。
-
@fangyuanaza 没有遇到过,你去网上找找原因吧,可以看看是不是版本原因,或者这个字典里其他的位置的错误,我用的2.3.0,感觉还是设定的问题
-
@fangyuanaza 我当时用face这个类型,你试试吧
-
@Mania 还有另外一种方法是安装两个版本的openfoam,然后后处理的时候用低版本的操作
-
@Mania 抱歉昨天忙了一天没有看到消息,没回复。
他的版本是比较低的版本,大概是2.3.0左右的,所以高版本下会出问题,你可以找下那个连接下有没有高版本的
或者另外一个帖子下有其他的链接,你可以试一下 -
@Mania 不知道版本有没有影响,在2017版本里面公式内的变量需要用大括号{},你可以看下tecplot help,下面这张图是17版本的help的提示
所以你换这个写的话,应该是{V34}=...,你可以试一下。
另外,其实你这个做三个面的平均的话,完全可以用matlab,非常简单,几行代码的事 -
一年多了还是没有解决这个算例,不知道有没有人做过相类似的,重新描述下这个算例,如图所示一个后台阶流动,但是是在图中左下(0,0,0)位置为原点绕Z轴旋转,主要流动参数,
雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$,其中$U_b$是体积平均速度;
旋转数$U_b=\frac{\Omega(H-h)}{U_b}$,$\Omega$是角速度;
台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$,这个与槽道流DNS结果是完全一致的,也就是说台阶前是完全发展的槽道流。目前存在以下问题:
1、求解器选择,SRF类还是需要自己修改求解器?
2、合理的进口条件,不旋转的相对好给,但是旋转的case,速度进口边界条件和初始值不知道如何给合理。
3、无量纲化的流向速度为什么会小于1,如下图所示。以我的经验,槽道流里面最大速度比体积平均速度肯定是大的,也就是$\frac{U_{max}}{U_b}$是大于1的,但是为什么在旋转流动里面就会于1呢?
参考文献
1、DNS原始文献:Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation
2、图片来自:Development and Application of Rotation and Curvature Correction,是篇博士论文 -
@liujunCFD 湍流层流粘度比这种说法来自哪,没有听过。湍流粘度比是一个估计值,指的是湍流粘度比流体粘度 $\mu_t/\mu $,一般取10就可以,这里有个解释可以看下
-
@liujunCFD CFDonline计算工具,可以自己对照一下
https://www.cfd-online.com/Tools/turbulence.php -
@chszkc 在我的理解里是不行的,设定了cyclic边界后就是周期性边界条件了,在ployMesh里面的boundary里面,inlet/outlet就必须设定为cyclic,而zerogradient适用于wall和patch,不能设定为cyclic,两者矛盾,应该计算不了。
-
@Stan同学 是的,这个跟paraview是不一样的,paraview有internalMesh这个选项,tecplot是没有的
你做振动的话,使用FSI吗? -
@Stan同学 不太清楚你算的什么算例,tecplot显示出来的应该是只有边界,要处理内部场的话可以尝试做个slice,或者处理一个Q准则等内部流场的涡结构,显示一下看看有没有结果
-
@Stan同学 把第一个zone--Step 1737显示试一下
-
@小考拉 我用的第二个
遍历网格单元的代码不知道你指的是什么意思
EXE = $(FOAM_USER_APPBIN)/spanwiseAvg这个就是编译一个程序,可以直接在终端调用,确实跟求解器的编译方法是一样的,其实求解器也是一个程序,类似于win里面的.exe的应用程序,只不过在这里需要通过终端调用 -
想计算LES下的耗散率,搜到cfdonline的一个程序CreateTurbulenceFieldsLES,是这个createTurbulenceFieldsLES2.zip ,这个是cfdonline链接,里面对于耗散率的计算非常简单,贴下代码
volScalarField epsilon ( IOobject ( "epsilon", runTime.timeName(), mesh,IOobject::NO_READ, IOobject::NO_WRITE ), turbulence->epsilon() );
对于这个不是很理解,箭头应该是指针吧,但是我不能很好的找到计算公式,貌似是采用的层流下的计算方式,因为在这个程序另一个文件createFields.H下有这么一处代码:
singlePhaseTransportModel laminarTransport(U, phi); autoPtr<incompressible::turbulenceModel> turbulence ( incompressible::turbulenceModel::New(U, phi, laminarTransport) );
然后我找了下层流laminar.C下计算的方法
tmp<volScalarField> laminar::epsilon() const { return tmp<volScalarField> ( new volScalarField ( IOobject ( "epsilon", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::NO_WRITE ), mesh_, dimensionedScalar ( "epsilon", sqr(U_.dimensions())/dimTime, 0.0 ) ) ); }
两个具体问题,我找的laminar.C这个是不是计算epsilon的公式,第二个问题层流下计算耗散率epsilon的公式为什么是这样的,按照定义应该是下图这样的,Eij就是指耗散率,应该是脉动速度的乘积,希望得到解答,谢谢!
-
@shepherd 其实我也不太懂,感觉是后处理的问题,我算过那个beamInCrossFlow的算例,结果也是只有流体域的网格变化,固体域的后处理不知道怎么弄
-
@shepherd 你截图没截全,现在提示的是warning,你看下是不是有error的提示,目录下应该有allrun的,运行不了吗?
-
@shepherd 感谢回复,希望后面有什么问题可以继续交流!
-
@bestucan 谢谢,感谢前辈详细的指导,我先好好研究研究下PPT
-
@炎冰灬灬 好的,谢谢
-
@炎冰灬灬 你好,我对fsifoam下的固体求解器有些疑问,发了一个帖子,希望能给予一些指导,在这就不重复问题了,这个是问题链接,先谢谢了
https://www.cfd-china.com/topic/3931/fsifoam中的固体求解器 -
在调研OpenFOAM中流固耦合的模块,看到固体求解器这部分卡住了,
从fsi的package可以看到,现在有两种固体求解器分别是unsIncrTotalLagrangianSolid和unsTotalLagrangianSolid,对于我一个固体方程小白来说,看代码还有点困难,找了一些资料,都是直接讲固体方程式是什么样的,完全没有这两个求解器的介绍,我甚至不知道这个求解去的名字应该如何理解,uns-total-lagrangian--**全部拉格朗日?
还有就是目前的关于结构方程的介绍与现在ansys中的完全不同,是方法不一样还是我没完全理解,有没有人给指条明路。 -
比如,我想计算k=0.5*(uu+vv+ww),其实这个是很简单的公式,写个小程序的话也不难,但是感觉程序写多了会很乱,我想能不能直接用paraview计算,然后输出成openfoam格式的场文件?
-
已经做过测试,tecplot可以直接导入openfoam未合并的processor的数据
-
@小考拉 paraview能不能做不清楚,但是可以写程序实现,下面是cfdonline上的程序,我用过,可以实现,你可以研究一下,源程序好像只有速度压力这些东西,alpha.water的话需要自己添加进去
https://www.cfd-online.com/Forums/openfoam-post-processing/106966-spanwise-averaging-three-dimensional-field.html -
@xiapfCFD
应该是这个位置,你试一下 -
@东岳 如果不考虑合并的话,处理数据,包括提取线、面上的速度,压力场什么的都可以在decomposedPar中实现,只是做云图这些处理,paraview用不太习惯,不知道有没有可能不做合并处理,可以用tecplot这些处理结果呢
-
@benqing 试了一个流场nut,第二种方法也不行,提示下面错误
而且工作量巨大 -
@东岳 所有时间都不行。
已经有人提交bug反馈了,结论是有两种方法,
一是swak4Foam引起的,注释掉就能用,但是我的算例并未用这个库,所以这个解决不了;另外一种是说AMI分到两个processor上引起的,解决方法是,先临时将AMI边界改成其他边界,合并完成之后再修改回来;下面是链接
我的算例文件比较大,用了40个核来计算,如果要改的话,有点太麻烦了,所以先上这里来问一下,看有没有人能够解决。 -
我用的是2.3.0,算例时双流道的泵,有周期面,有动静交界面,动静交界面采用AMI处理,计算可以进行,但是在重构流场(reconstructPar)的时候出错,错误信息如下,有人遇到过这个问题吗,CFDonline上有很多人遇到这个问题,但是没有解决方案
-
@一二 都是用matlab做的
-
emm,画出来了,打扰了,下次一定先自己好好研究在来提问
-
@bestucan 你好,又遇到一个图。这样的图,画散点图好弄,怎么加上线和面呢?
-
@bestucan 感谢回复和指导
虽然不太一样,但是基本满足我的期望 -
我主用matlab,用过origin,这两个软件能做吗?
给个方向我可以自己研究。 -
准确的应该是wallShearStress,计算的是壁面剪切应力tau_wall
求助:OpenFOAM空化流动计算,空化源项如何输出
LES中关于Resolved TKE和Modelled TKE求解比例的问题
LES中关于Resolved TKE和Modelled TKE求解比例的问题
PK! OpenFOAM VS Fluent计算效率
PK! OpenFOAM VS Fluent计算效率
雷诺应力与URANS
雷诺应力与URANS
UDF求计算域中某网格到壁面的最近距离(C_WALL_DIST宏)
UDF,非线性雷诺应力,USDI?
UDF求计算域中某网格到壁面的最近距离(C_WALL_DIST宏)
UDF,非线性雷诺应力,USDI?
请问各位大佬,如何对数据在指定平面进行空间平均?
请问各位大佬,如何对数据在指定平面进行空间平均?
请问各位大佬,如何对数据在指定平面进行空间平均?
openfoam后处理按照网格信息提取数据
openfoam后处理按照网格信息提取数据
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
SRF求解器的使用
关于k-omega计算的问题
关于k-omega计算的问题
关于cyclic边界条件
Tecplot无法显示内部流体
Tecplot无法显示内部流体
Tecplot无法显示内部流体
paraView-求展向平均问题
关于CreateTurbulenceFieldsLES疑问
关于流固耦合求解器fsiFoam的问题
关于流固耦合求解器fsiFoam的问题
fsiFoam中的固体求解器
fsiFoam中的固体求解器
FsiFoam并行问题
FsiFoam并行问题
fsiFoam中的固体求解器
paraview可以直接将calculator的结果输出成场文件吗
有AMI边界的算例,reconstructPar时出错
paraView-求展向平均问题
tecplot作图求助
有AMI边界的算例,reconstructPar时出错
有AMI边界的算例,reconstructPar时出错
有AMI边界的算例,reconstructPar时出错
有AMI边界的算例,reconstructPar时出错
这种图怎么做
这种图怎么做
这种图怎么做
这种图怎么做
这种图怎么做
输出壁面摩擦速度