这里已经有关于向量场的了。现在我想多了解一点。
如何使每个箭头具有相同的固定长度,并按颜色定义值的大小?
还是不可能把流线画成格鲁图呢?如果可能的话,我该怎么做?
现在我有这个,需要升级它。
set term pngcairo
set title 'Navier-Stokes Equation'
set terminal png size 1280,720
set output 'vec.png'
plot 'vec' u 1:2:($3/$5):($4/$5) w vec t 'Vector Field'
更新
多亏了
对于如何编写代码来实现这个概念,我有点困惑。所以,如果有人知道如何解决我的问题,请帮帮我。
我有一个want循环来运行一个模拟,在这个过程中,我想看到一个质量的流体在t_domain给出的时间内的扩散。时间t是以dt的步数递增的,这两个变量都是类型为double的变量,但是我只想在Del_t时间间隔内拍摄扩散过程的快照。我的意思是,如果我运行1000秒的模拟,时间t增加0.01秒,我想每10秒拍摄一次扩散过程的快照。
现在,我无法弄清楚,在文件的第9步中,打印密度的代码部分。有人能给我建议吗?
while(t<t_domain)
{
//runs a function to
我只想逐行从文件中读取,并将每一行附加到数组中。我以前做过,没有任何问题,我不知道为什么这一次它不工作!打印函数只显示文件的第一行!!这是我的代码:
keyword_array = []
with open('local directory\\C0577785c.txt') as my_keywordfile:
for keyword in my_keywordfile.readline().strip("[]").strip("'").split(","):
keyword_ar
我想要生成一个2D纹理,首先生成一些2D频域系数,然后使用ifft生成纹理。为了计算法线映射,我想取这个二维函数的导数。我只想知道以下内容是否正确:
in 1D I would just do a loop for (k = 0; k < N, k++) coefficients[k] *= i * k; where i is sqrt(-1)
in 2D I should do this? for (m=0;m<N;m++) for (n=0;n<N;n++) coefficients[m,n] *= (i*m)*(i*n) -> -m*n;
如果存在,我需要确定2d Bezier曲线上的拐点(曲率变化的点),参数化为t,0 <= t <= 1。我最初的方法是沿着曲线取样,评估二阶导数,找出导数符号变化的点。
2DVector curvature1, curvature2;
for (double t = 0, t <= 1.0; t += STEP) {
curvature1 = bezier.CurvatureAt(t);
curvature2 = bezier.CurvatureAt(t + (STEP/2.0 >= 1.0 ? 0 : t + STEP/2.0));
if (i
我正在研究一个耦合的DEM-CFD模型。我使用Fipy作为它的CFD部分,我的几何形状是圆柱体。如何将Gmsh中生成的网格导入到Fipy?
我需要生成一个圆柱几何体,并通过在Fipy中生成3D三角形网格来对其进行离散化。我试图在Fipy的帮助下找到一个命令来生成柱面几何中的3D三角形网格,但我找不到任何与3D相关的东西。我在Gmsh软件中生成了网格,该软件基于Fipy的帮助而被Fipy接受。我现在有网格文件,但是我不能把它导入到Fipy中。
每当我试图将它导入到模型中时,我都会遇到这个错误Gmsh version must be >= 2.0。然而,我的Gmsh版本是4。
我也使用了Gm
我正在使用fipy来解决一个计算域中具有孤立区域的扩散问题。请参见示意图。下图:
在fipy.FaceVariable的帮助下,在@Daniel Wheeler的帮助下,该问题在Fipy下进行了建模,并使用了fipy.FaceVariable定义的可变系数。
然而,计算速度无法满足我的命令,这比使用有限差分法的cython代码要慢得多。如果我想加速fipy计算,我该怎么做?以下是我的演示代码:
from pylab import *
from numpy import *
import fipy
from scipy.spatial import Delaunay
from fipy.vari
我已经生成了一个gmsh网格,我通过Gmsh2D导入到了Gmsh2D中。我想把表面的表面固定下来,但是我不知道该怎么做。
在我在文档中找到的FiPy示例中,建议给出某些行的名称,以便稍后处理它们。在我将网格导入fipy之后,如何做到这一点?
// note: if you do not use any labels, all Cells will be included.
Physical Surface("Outer") = {1};
Physical Surface("Middle") = {2};
Physical Surface("Inner&
我是第一次接触fipy,所以如果我问一些显而易见的问题,请原谅我的无知。但是我无法运行已经存在的(并且在其他机器上工作的)脚本,获取在安装过程中可能遗漏了某些东西的EnvironmentError: Gmsh version must be >= 2.0.?Fipy手册对gmsh的实现有点含糊。它只提供了下载文件的链接,可以通过它生成网格,但它没有说明gmsh应该安装在哪里(所以我的gmsh-3.0.3-Windows是默认的,也就是Program Files (x86) )。它是否应该安装在某个特定的目录中,也许与fipy相同?(我真的为一个愚蠢的问题道歉,但如果你知道该问什么问题,你
我运行了以下代码来模拟二维网格中圆柱体周围的流动:
from fipy import CellVariable, FaceVariable, Grid2D, DiffusionTerm, ImplicitSourceTerm, PeriodicGrid2DTopBottom, DistanceVariable, Viewer
from fipy.tools import numerix
L = 1.0
N = 50
dL = L / N
viscosity = 1
U = 1.
#0.8 for pressure and 0.5 for velocity are typical relaxa
通过FiPy.grid3D("")导入.msh2网格类型,在PFC软件中进行网格化有什么问题?该网格从openFOAM中导出,并转化为FiPy的可读.msh。嵌入式FiPy版本为3.1,Python版本为2.7。
IPython 3.1.0 -- An enhanced Interactive Python.
? -> Introduction and overview of IPython's features.
%quickref -> Quick reference.
help -> Python's own h