標簽[fipy] - 堆棧內存溢出

[英]Python Fipy explicit terms with derivatives and variable diffusion constant

我正在嘗試求解四個耦合 PDE（兩個標量密度場和一個矢量場 - 但我以分量方式求解）。其中第一個方程中的第二項由下式給出現在，在 Fipy 代碼中（\rho 和矢量分量分別用 m、px 和 py 表示）我將它們寫下來： c 方程中第二項的 Fipy 實現是否正確？代碼沒有出現任何錯誤，但我 ...

如何使用 FiPy 模擬矩形環上的熱擴散？

[英]How to simulate a heat diffusion on a rectangular ring with FiPy?

我剛開始使用 FiPy 求解 PDE 並在矩形環狀銅體上進行熱擴散實驗。有時這是一個 plot 的模擬結果。我使用Grid2D()作為網格，使用CellVariable.constrain()指定邊界條件。綠點是外表面的中心，其中T = 273.15 + 25 ( K )，藍點是內表面的中 ...

使用 anaconda 安裝 fipy 的最佳方法是什么？

[英]What is the best way to install fipy with anaconda?

使用 anaconda3 安裝 fipy的最佳方法是什么（在 ubuntu 16.04 機器中）？我試過conda create --name FIPY --channel conda-forge python=3 fipy numpy scipy matplotlib pysparse gms ...

有效地表示 FiPy 中的“零級”可飽和水槽且沒有負濃度

[英]Representing "zero-order" saturable sinks in FiPy efficiently and without negative concentration

我正在嘗試使用 FiPy 來解決具有小圓形水槽（粒子濃度 = n ）的簡單 2D 粒子擴散場景的瞬態行為。我希望水槽是可飽和的，即水槽內應該有最大可能的損失率。同時，濃度不得以 go 為負。接收器的精確行為應該是：當進入水槽的粒子總速率<=最大速率時，水槽吸收了它里面的所有粒子。這就 ...

python 中的 2d fipy PDE 求解器中的 stackoverflow 錯誤

[英]stackoverflow error in 2d fipy PDE solver in python

我正在嘗試解決二維中的三個耦合 PDE 和兩個 ODE（耦合）的系統。該問題試圖解決粒子在介質中移動的情況。介質由場速度分量 vx、vy 和密度 m 給出，並且有一個粒子在該介質中移動，position Rx, Ry。它運行了一段時間，但隨后拋出錯誤：“致命的 Python 錯誤：無法從堆棧溢 ...

Fipy 錯誤：''因子完全是奇異的''，當應用 Neumann 邊界條件時

[英]Fipy error:’’The Factor is exactly singular’’, when applying Neumann boundary conditions

我們正在嘗試解決 Fipy 中的一維耦合連續性-泊松問題。當應用 Dirichlet 條件時，它給出了正確的結果，但是當我們將邊界條件更改為更接近我們的問題的 Neumann 條件時，它給出了“因子完全奇異”錯誤。非常感謝任何幫助。代碼如下（0<x<2.5）： ...

在 Fipy 中求解的線性系統

[英]the linear system being solved in Fipy

我是否通過緩存 Matrix 和 RHS 並運行求解/掃描 1 次迭代來獲得線性系統？另外，是否可以手動設置正在使用的矩陣？謝謝 ...

如何在gmsh上從圓心到邊緣進行擴散？ [Python]

[英]How to make a diffusion from the center of the circle to the edge on gmsh? [python]

我正在嘗試使用 gmsh 工具在 python 上模擬營養物質在腫瘤中的擴散作為空間和時間的函數。所以我需要初始濃度位於中心，然后擴散性將其向外擴散。這是完整的代碼，參數還不正確，但我已經在嘗試獲得從內到外的擴散率，而不是相反。這是我的代碼：非常感謝！ ...

通量邊界條件一維擴散

[英]flux boundary condition 1D diffusion

我目前正在學習如何使用 fipy，並且在嘗試運行下面的代碼時遇到了問題。問題是從右邊界流出的靜止一維擴散。我在其他帖子中讀到了不適用於 fipy 中的這個問題。我不斷收到錯誤消息： RuntimeError: Factor 是完全奇異的但是對於兩個 Neumann 邊界條件，這 ...

一維平流-擴散系統中不連續的fipy問題

[英]fipy problem with discontinuity in an 1D advection-diffusion system

我正在嘗試對涉及可變平流速度的一維平流擴散問題進行建模。它涉及一塊向下流動的冰塊。我遇到了在速度變化的網格點處發生非物理不連續的問題（附圖）。該代碼適用於恆定速度的情況。任何幫助將不勝感激頂部繪制平流系數作為深度的函數中間繪制： cie 板的溫度分布作為深度的函數。兩條邊的邊界條件分別為 ...

無法在 FIPY 中重現球形網格的結果

[英]Unable to reproduce results for spherical grid in FIPY

我正在嘗試重現這篇文章的結果（考慮到 jeguyer 在回答中提供的修復）。但是當我嘗試運行完全相同的代碼時，我收到以下警告：生成的 X_ca 數組全部由 -inf 組成。有什么辦法可以解決嗎？先感謝您。 ...

在 FiPy 中解決人口平衡 PDE 時遇到問題

[英]Having trouble solving Population balance PDE in FiPy

我是使用 FiPy 的初學者，也是一般使用有限體積法的初學者。如果我的問題很幼稚，請原諒我我正在嘗試求解具有如下邊界條件的PDE形式的人口平衡方程Boundary Conditions x 的取值范圍是 0.5e-6 到 1000e-6，我想使用均勻網格。這是我目前擁有的代碼，但無法正常工作。 ...

在 FiPy package 中掃描時選擇 dt 會改變解的時域？

[英]Chosen dt in sweep in FiPy package changes time domain of solution?

我目前正在嘗試解決具有 18 個 PDE 和 ODE 的大型 PDE 系統（這就是我使用 FiPy 而不是 scipy 中的 odeint 的原因）。其中一個等式非常簡單：dROS/dt = 0.06，我實現如下：所需的時間范圍是從 o 到 3650，所以我用這種方式求解，使用掃描（我在這里只 ...

FiPy 求解偶偏微分方程導致 IndexError

[英]FiPy Solving Couple PDEs Results in IndexError

我是 FiPy 的新手，我正在嘗試解決一對耦合 PDE，描述了一種材料對裝在圓柱體中的吸附劑的吸附。材料流經圓柱體，其濃度處於 2 個相中的任何一個：流動相或固定相。我將流動相的網格設為沿圓柱體長軸的一維網格。我將固定相的網格作為圓柱體橫截面積的一維圓柱體網格。兩相是耦合的，我試圖求解流動 ...

如何用 FiPy 表示這些非標准 pdes？

[英]How to represent these non-standard pdes with FiPy?

我正在嘗試使用 FiPy 求解一組偏微分方程，但不確定如何表示方程：請參閱此處的偏微分方程。給我帶來麻煩的術語是紅色矩形和藍色矩形中的術語。紅色矩形是對流項乘以 function，取決於 x、T1 和 T2；藍色矩形中的一個是瞬態項乘以另一個 function，具體取決於 x、T1 和 T2 ...

使用 FiPy 進行反向擴散

[英]Reverse diffusion with FiPy

在我的鳥類遷徙中，我想代表鳥類在非洲度過一段時間后返回歐洲。為了實現這一點，我需要能夠“逆轉”鳥類的擴散。我試圖為擴散設置一個負系數，但它當然不起作用。所以我想知道是否有可能用 FiPy 實現這一點。 ...

FiPy 中隨時間演化的平流系數，平流 - 擴散方程 (2D)

[英]Evolutive advection coefficient with time in FiPy, advection - diffusion equation (2D)

我正在嘗試開發一個代表鸛從法國遷移到中非的 model。我選擇了 go 和 FiPy。我想將時間設置為完整的遷移時間：首先，鸛在法國隨機擴散，但是當時間達到極限值時，開始遷移並出現平流系數。我首先嘗試了擴散系數，它起作用了：結果如下：前五張快照：然后擴散系數變化：但是當我嘗試相同 ...

如何使用 FiPy 在同一側求解具有 Dirichlet 和 Neumann 邊界條件的 1D Poisson 方程

[英]How to solve 1D Poisson equation with Dirichlet and Neumann boundary conditions on the same side using FiPy

顯然，即使 Dirichlet 和 Newmann 邊界條件在同一側，具有恆定源的一維泊松方程也有唯一解。但是我想不出在 FiPy 中解決這個問題的方法。請指教。我試過了：這不是在facesLeft gradValue 但是，它可以在facesRight gradValue ...

運行 stokesCavity.py 示例時出現查看器錯誤

[英]viewer error when running the stokesCavity.py example

使用 python-3.9 和 python-fipy-3.4.2.1-1 運行 Manjaro stable。剛開始接觸FiPy，最終對寫單相流和兩相流代碼感興趣。自然地，我嘗試使用以下命令運行examples/flow/stokesCavity.py （從所有第一個文本中刪除）： pyth ...

FIPY 求解 PDE 和 ODE

[英]FIPY Solving a PDE coupled with an ODE

我有一組耦合方程，其中主偏微分方程（時間和位置 z 的函數）給出為：第二個方程是這樣的：其中k_m = f(q)和q^* = f(c) 。正如你所看到的，第二個方程是一個常微分方程（ q不直接依賴於空間）。我發現很難編寫代碼來耦合這兩個方程。到目前為止，對於我忽略第二個方程並采 ...