fipy

2熱度

2回答

我是fipy的新手，所以如果我問一些顯而易見的問題，請原諒我的無知。但是我無法運行一個已經存在的（並且在其他機器上運行）腳本，得到EnvironmentError: Gmsh version must be >= 2.0.我可能在安裝過程中遺漏了一些東西？關於gmsh的實施，Fipy手冊有點模糊。它只提供下載文件的鏈接，並且網格可以由它生成，但是它沒有說應該安裝WHERE gmsh（所以我的默

0熱度

1回答

如何提高Fipy效率

我正在使用fipy解決計算域中孤立區域的擴散問題。參見原理圖。如下： where there is no flux between isolated BC, and flux exists at periodic BC 這個問題在Fipy模型下，在@Daniel Wheeler的幫助下通過fipy.FaceVariable定義varing coeff。但是，計算速度無法滿足我的命令，這比使用有限

2熱度

1回答

翻譯2- d的PDE和條件成準確Fipy碼

enter image description here 圖像包括控制方程，初始和邊界condition.It descripe板和流體之間的熱傳遞的問題。我不知道如何使用fipy編碼包含var的二維問題和邊界條件。這是我的嘗試。 from fipy import * import numpy as np #constant Pe=2400 le_L=1/20000 L_l=2000

0熱度

1回答

一般邊界條件

我很感謝設置一般邊界條件的幫助，-grad(y) + g(y) = 0其中g是未知的某些函數y。這裏有一個簡單的例子1D，我不能去上班： N=3 h=1./(float(N)-1.) mesh = Grid1D(nx=N, dx=h) c=CellVariable(mesh=mesh,value=0.5) ## Dirichlet boundary conditions #c.co

2熱度

1回答

在Python中解決PDE源代碼

我正在使用FiPy解決受生物學啓發的問題。本質上，我想表示一個2D平面，在不同的點我有源和匯。源以固定的速率發射襯底（不同的源可以具有不同的速率）並且阱以固定的速率消耗襯底（不同的阱可以具有不同的速率）。我的代碼： import numpy.matlib from fipy import CellVariable, Grid2D, Viewer, TransientTerm, Diffusio

1熱度

1回答

FiPy - 域延伸/幀增長

我想在深度（h（t））隨時間變化的域上求解簡單的擴散方程（dT/dt = K d2T/dx2）。因此得到的等式是： dT/dt = K/h^2 d2T/dx2 + z/h dh/dt dT/dz 其中z現在是固定的0-> 1域。新術語是幀平流，我試圖包括它，但我正在努力與空間依賴係數。當我把它給了對流項之外： mesh.cellCenters[0]*PowerLawConvectionTe

0熱度

1回答

使用FiPy和Mayavi的解決3D

擴散方程我感興趣的解決， \frac{\delta \phi}{\delta t} - D \nabla^2 \phi - \alpha \phi - \gamma \phi = 0 以下是工作，但我有幾個問題：是否有可能通過FiPy提高性能？儘管計算時間很長，但我覺得這裏的nx, ny, nz箱子非常小。我不明白爲什麼陣列 X, Y, and Z是如此之大。請注意，在第一幀中，我們放大

2熱度

3回答

您如何在Fipy中指定Neumann（法向面的固定通量）邊界條件？

我如何明確地設置fipy網格中的邊界面的流量法線爲一個特定值，而不會限制面內的流量分量？（1）通量的固定分量垂直於邊界面，或（2）作爲面的通量的完整規定。默認的fipy條件是前者（值= 0），但顯式方法（faceGrad.constrain）是後者。通過將以下代碼添加到fipy diffusion.mesh20x20示例的結尾並注意不同的臉部漸變結果可以理解此問題。 facesNeumann =

1熱度

1回答

存儲術語數組，而不是fipy對象

我是新來fipy，所以我道歉，如果這是一個愚蠢的問題（和this似乎並沒有幫助我）。但有存儲在人類可讀（或Python可讀）形式fipy對象的方式，除了在上面的問題提出？這隻適用於單元格變量。如果我想做一些更看中/自定義繪製比什麼是默認fipy觀衆，我該怎麼辦呢？採取例如一個簡單的一維擴散： from fipy import * # USER-DEFINED PARAMETERS nx

1熱度

1回答

FiPy內部進行邊界條件

我是一個新人，FiPy，我解決了泊松方程的三維體積內的潛力。它適用於表面邊界條件，但現在我需要在裏面放置一個導體。這將是一個恆定的潛在表面，我意識到你不能將potential.constrain用於內部表面。該文檔建議使用ImplicitSourceTerm以及定義表面的掩膜，但不清楚如何使用它來限制電位恆定，或者等同於約束電場垂直於表面。這可能嗎？感謝您的任何幫助。