fipy

    2熱度

    1回答

    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回答

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

    2熱度

    2回答

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

    1熱度

    1回答

    我想在深度(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

    2熱度

    3回答

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

    1熱度

    1回答

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

    0熱度

    1回答

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

    2熱度

    1回答

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

    1熱度

    1回答

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

    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