本文先容OpenFOAM中samping操作基础。
functionObject不错在求解流程中输出指定的物理量信息。当求解计较结束后,不错选拔使用Sampling获得特定位置物理信息,此时主要使用postProcess门径来斥逐。
欺骗postProcess进行采样,需要界说采样字典文献,这些文献位于system旅途下,文献不错淘气定名。当采样计较结束后,算例旅途中将会新建一个名为postProcessing的新文献夹,在该文献夹内部包含与字典文献同名的子目次,其中存储了采样获得的恶果数据。
底下贯串子际算例来进行形容。
如下图所示,算例文献system文献夹中包含两个字典文献probeDict及sampleDict,其中界说了取样操作所需的参数。
图片
1、probesDict字典文献
该文献中界说了需要在计较域中进行取样的点。本算例取样点的位置如下图所示。
图片
probesDict文献本质如下所示。
FoamFile{ version 2.0; format ascii; class dictionary; object probesDict;}// * * * * * * * * * * * * * * * * * //// 指定类型为probestype probes; // 指定需要索取的物理场fields( p U); // 取样点的坐标列表probeLocations( (1.0 0 0) (1.25 0 0) (1.5 0 0) (1.75 0 0) (2.0 0 0) (2.0 -0.25 0) (2.0 -0.5 0) (2.0 0 .25 0) (2.0 0 .5 0));
文献准备结束后不错使用高唱:
postProcess -func probeDict -time 2000
这里指定参数-time 2000,则只输出2000 s期间的数据,若不指定此参数, 武汉市科佳机电进出口有限责任公司则会每一期间的数据。
开动结束后, 首页-汉盛奋皮具有限公司算例文献夹中的postProcessing旅途下会多出一个名为probesDict的子文献夹,湖北碧山机械股份有限公司其中包含了2000 s期间监测点位置的U与p的值。
图片
关于本算例,2000/U文献中的本质如下所示,列举了所检测的9个点位置的速率重量:
# Probe 0 (1 0 0)# Probe 1 (1.25 0 0)# Probe 2 (1.5 0 0)# Probe 3 (1.75 0 0)# Probe 4 (2 0 0)# Probe 5 (2 -0.25 0)# Probe 6 (2 -0.5 0)# Probe 7 (2 0.25 0)# Probe 8 (2 0.5 0)# Time 0 1 2 3 4 5 6 7 8 2000 (1.0960031 -0.61045141 0) (1.0003352 -0.35057338 1.1984886e-20) (0.99511675 -0.22368641 -1.4884685e-20) (0.99587925 -0.17217799 6.8574687e-21) (0.99434334 -0.13173595 0) (0.85681371 -0.10854681 3.6035834e-20) (0.95135192 -0.1125737 2.353505e-20) (1.0221783 -0.12625997 -1.5616405e-20) (1.041875 -0.12290482 -6.6065831e-20)
文献2000/p本质如下所示,其列举了9个定位置的压力值:
# Probe 0 (1 0 0)# Probe 1 (1.25 0 0)# Probe 2 (1.5 0 0)# Probe 3 (1.75 0 0)# Probe 4 (2 0 0)# Probe 5 (2 -0.25 0)# Probe 6 (2 -0.5 0)# Probe 7 (2 0.25 0)# Probe 8 (2 0.5 0)# Time 0 1 2 3 4 5 6 7 8 2000 -0.60129439 -0.061042557 -0.020796443 -0.011252459 -0.0028472253 0.020868043 0.041151249 -0.030775763 -0.050775396
2、sampleDict字典文献
该文献界说了一条线,建筑设计之后对该线上的物理量进行索取。
FoamFile{ version 2.0; format ascii; class dictionary; object sampleDict;}// * * * * * * * * * * * * * * * //// 指定类型为sets,频繁不错为points及linestype sets;libs ("libsampling.so"); // 插值阵势,不错是cellPoint、cellPointFace、// cellPointWallModified或cellinterpolationScheme cellPoint; // 文献存储神志setFormat raw;surfaceFormat raw; // 指定待取样的物理场// 这里指定了速率与壁面剪切应力// 本算例中并不行索取到壁面剪切应力// 物理量不存在也没相干系fields( U wallShearStress); // 界说靠拢sets( profile0 { // 指定几何类型为线 type lineCellFace; axis distance; // 两个点细目一条线 start ( 0.75150 0.04767 0 ); end ( 0.76168 0.14715 0 ); } );/////// 若不明晰重要字选项,不错使用banana步调
实践以下高唱索取数据:
postProcess -func sampleDict –time 2000
高唱实践结束后会在postProcessing文献夹中增多子文献夹sampleDict,并在其中存储采样数据。
图片
不错使用gnuplot绘图速率重量沿直线长度的离别:
set xlabel 'length(m)'set ylabel 'velocity(m/s)'set key outsideplot 'profile0_U.xy' u 1:2 w l lw 2 lc 7 t 'u_x','' u 1:3 w l lw 2 t 'u_y','' u 1:4 w l lw 2 t 'u_z'
如下图所示。
图片
关于采选不同的数据集,其存储的文献体式有所不同。
如点集采样,关于标量场,其文献存储体式为:
#POINT_COORDINATES (X Y Z) SCALAR_VALUE0 0 0.05 13.3109950 0 0.1 1 9.263215...
若为矢量场建筑设计,则其存储体式为:
#POINT_COORDINATES (X Y Z) VECTOR_COMPONENTS (X Y Z)0 0 0.05 0 0 13.3109950 0 0.1 0 0 19.263215...
若采样集为线,则存储体式为:
#AXIS_COORDINATE SCALAR_VALUE0 13.3109950.0015 1 9.263215...
若为矢量场,则其存储体式为:
#AXIS_COORDINATE VECTOR_COMPONENTS (X Y Z)0 0 0 13.3109950.0015 0 0 19.263215...
若采样的为名义,则存储体式
#POINT_COORDINATES (X Y Z) SCALAR_VALUE0 0 0.05 13.3109950 0 0.1 19.263215...
若为矢量场,则其存储体式为:
首页-微俊哲咖啡有限公司#POINT_COORDINATES (X Y Z) VECTOR_COMPONENTS (X Y Z)0 0 0.05 0 0 13.3109950 0 0.1 0 0 19.263215...本站仅提供存储管事,通盘本质均由用户发布,如发现存害或侵权本质,请点击举报。