使用Mesh2D构造函数重新创建fipy Grid2D
Mesh2D是fipy库中的一个类，用于创建二维网格。它可以通过Mesh2D构造函数来创建。
构造函数的语法如下：
Mesh2D(nx, ny, dx=1.0, dy=1.0)参数说明：
nx：网格在x方向上的单元格数ny：网格在y方向上的单元格数dx：单元格在x方向上的尺寸，默认为1.0dy：单元格在y方向上的尺寸，默认为1.0
使用Mesh2D构造函数可以重新创建fipy的Grid2D对象，Grid2D是fipy库中用于表示二维网格的类。
以下是使用Mesh2D构造函数重新创建fipy Grid2D的示例代码：
from fipy import Grid2D, Mesh2D
nx = 10
ny = 10
dx = 0.1
dy = 0.1
mesh = Mesh2D(nx, ny, dx, dy)
grid = Grid2D(mesh)在上述示例中，我们创建了一个10x10的网格，每个单元格的尺寸为0.1x0.1。然后使用Mesh2D对象创建了一个Grid2D对象。
Mesh2D和Grid2D在计算领域中常用于有限差分方法（Finite Difference Method）的数值计算，用于解决偏微分方程等问题。它们可以用于模拟物理现象、流体力学、热传导等领域的计算。
腾讯云提供了云计算相关的产品和服务，其中与网格计算相关的产品是腾讯云弹性容器实例（Elastic Container Instance，ECI），它提供了一种快速部署容器化应用的方式，可以用于网格计算等场景。您可以通过访问腾讯云弹性容器实例的官方文档了解更多信息：腾讯云弹性容器实例
相关·内容
我在fipy中使用了上面的代码，但是颜色条的范围总是不变(总是显示-1,1)。我的版本是Python 3.6.7和matlablib 3.1.1。
#create a viewer to see the results
if __name__ == '__main__':
viewer = Viewer(vars=psi, datamin=-0.01, datamax=0)
viewer.plot()
我尝试了以下方法:1、我更改了"C:\Users\Lenovo\fipy\examples\diffusion“中的代码，例如，更
浏览 9提问于2019-10-20得票数 0
回答已采纳
我对用fipy解微分方程很感兴趣。
当我使用Grid2D时，以下代码工作正常。
from fipy import *
mesh = Grid2D(nx=3, ny=3)
#mesh = Grid3D(nx=3, ny=3, nz=3)
phi = CellVariable(name='solution variable', mesh=mesh, value=0.)
phi.constrain(0, mesh.facesLeft)
phi.constrain(10, mesh.facesRight)
coeff = CellVariable(mesh=mesh, value=1
浏览 0提问于2018-03-01得票数 0
from fipy import CellVariable, Grid2D, DiffusionTerm
nx = 20 # grid size on coordinate axes
dx = 1 # grid spacing
ny = nx; dy = dx; L = dx*nx
mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)
phi = CellVariable(name = "phi", mesh = mesh, value = 0.5)
eq = (0. == DiffusionTerm( coeff
浏览 13提问于2021-02-27得票数 2
我正在尝试使用python中的fipy包来解决pde。我为同样的代码编写的代码已经附上。
!pip install fipy
from fipy import *
mesh= Grid2D(nx=0.001,dx=100,ny=0.0005,dy=100)
phil=CellVariable(name='Sol variable',mesh=mesh)
phil.constrain(0,mesh.facesBottom)
phil.constrain(1,mesh.facesTop)
n=1.7*10**(-6)*((0.026*numerix.exp(phil/0.026)+
浏览 1提问于2020-05-03得票数 1
我有下面的偏微分方程
我想知道如何用Fipy python来表示源术语。我试过以下几种方法
from fipy import *
nx = 50
ny = 1
dx = dy = 0.025 # grid spacing
L = dx * nx
mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)
phi = CellVariable(name="solution variable", mesh=mesh, value=0.)
convCoeff = ((10.,), (10.,))
Gamma =
浏览 5提问于2019-09-18得票数 0
回答已采纳
1回答
我想知道如何在Fipy中将狄拉克增量函数表示为源项。我想解下面的方程式 ? 我已经尝试了以下代码 from fipy import *
nx = 50
ny = 1
dx = dy = 0.025 # grid spacing
L = dx * nx
mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)
phi = CellVariable(name="solution variable", mesh=mesh, value=0.)
Gamma=1
delta=1 # I want knowing how to make this rig
浏览 17提问于2019-09-21得票数 1
回答已采纳
2回答
二维网格类在C++中的实现
、、、、
我想要实现一个2D网格类，这是内存效率，因为它将用于蒙特卡罗模拟。然而，在这个阶段，这实际上与2D数组类相同。不用说，我是C++的新手。我写了以下代码。
#define MESH2D_H_INCLUDED
class Mesh2D
{
private:
double* mesh;
size_t rows;
size_t columns;
public:
/* constructor */
Mesh2D(size_t Nx, size_t Ny){
rows = Nx;
columns = Ny;
浏览 2提问于2020-07-03得票数 2
回答已采纳
1回答
我正在使用Fipy解决一个包含两个参数或常量的PDE系统，所以我想知道是否也可以在Fipy中估计这些参数，或者其他哪些库更适合这样做。
注意:我知道scipy有一些功能(用于MLE的optimize.minimize)，但我不确定是否足够将它们应用于Fipy的代码。
更新:对于下面的PDE系统，我想估计两个未知参数："Beta“和"m”
在Fipy中求解此PDE的函数如下所示：
import scipy as sci
import fipy as fipy
import numpy as np
from fipy import *
# Grid
nx = 100
n
浏览 1提问于2017-03-01得票数 0
回答已采纳
1回答
我正在尝试使用FiPy重现Matlab示例。但是，我无法重现Matlab的结果(使用Fipy时我是新手)。我认为问题可能出在函数或边界条件声明中。此外，对于初始条件，我没有找到关于t的变量的setValue属性，所以我用非常笨拙的循环结构覆盖了它。下面是我的测试代码：
from fipy import *
from fipy import Grid2D, CellVariable,TransientTerm, DiffusionTerm
from fipy.tools import numerix
import matplotlib.pyplot as plt
import numpy as
浏览 7提问于2021-09-21得票数 0
1回答
有一个数独解算器方法，可以检查行/列位置的解决方案，并在测试数字是否不正确后进行回溯。它先打印No Solution，然后打印Abort Trap Error: 6。我正在使用一个应该返回TRUE的测试用例。Row和Col最初以0的形式传入参数
可能的问题是什么？我想这可能是我递增row和col值的方式。它读入单行文件并转换为2D网格
下面是测试输入: 4.....8.5.3..........7......2.....6.....8.4......1.......6.3.7.5..2.....1.4......
#include <stdio.h>
#include <st
浏览 15提问于2018-02-23得票数 0
回答已采纳
考虑一下这个小pybind11包装器+测试：
setup.py
from pybind11.setup_helpers import Pybind11Extension
from pybind11.setup_helpers import build_ext
from setuptools import setup
setup(
name="wrapper",
ext_modules=[Pybind11Extension("wrapper", ["wrapper.cpp"])],
cmdclass={"build
浏览 4提问于2021-05-08得票数 11
1回答
我用FiPy来解决一个受生物学启发的问题。
本质上，我想代表一个二维平面，在不同的点，我有源和汇。源以固定的速率发射基板(不同的源可以有不同的速率)，汇以固定的速率消耗基板(不同的汇可以有不同的速率)。我的代码：
import numpy.matlib
from fipy import CellVariable, Grid2D, Viewer, TransientTerm, DiffusionTerm, ImplicitSourceTerm, ExplicitDiffusionTerm
from fipy.tools import numerix
from time import *
nx
浏览 2提问于2016-04-03得票数 2
回答已采纳
我正在尝试实现下面的示例 https://nbviewer.jupyter.org/github/barbagroup/CFDPython/blob/master/lessons/14_步骤_11.ipynb，但是我面临着一些问题。我认为主要的问题是，我在边界条件以及定义方程中的项时遇到了一些问题。 PDE包括： ? 初始条件和边界条件是： ? 我知道我的变量是速度矢量(u，v)和压力(p)的两个分量。按照示例并使用FiPy，我对PDE进行了如下编码： import numpy as np
import pandas as pd
import matplotlib.pypl
浏览 112提问于2021-02-17得票数 1
回答已采纳
如果我使用以下命令创建一个fipy 2D网格： mesh = Grid2D(nx=3, ny=3, dx=1., dy=1.) 我得到了一个合理的矩形单元网格。但是，如果我取这个网格的值并使用它们来构造一个普通的Mesh2D mesh2d = Mesh2D(mesh.vertexCoords, mesh.faceVertexIDs, mesh.cellFaceIDs) 我得到了一个三角形的网格，其中一些三角形似乎是无效的。如果我这样做： model = CellVariable(mesh=mesh2d, value=1.) 并使用matplotlib查看器查看，一些三角形是未填充的，我认为这
浏览 16提问于2019-04-24得票数 0
回答已采纳
1回答
我是个新手。我正在尝试使用fipy解决以下一组耦合的pdes (变量是p，n和psi) :-
qDn∇2n − qun∇.(n∇ψ) = q(n-10**11)/10**(-6), in Domain 1
qDp∇2p + qup∇.(p∇ψ) = -q(p-10**21)/10**(-6), in Domain 1
e∇2ψ = −(p − n- N) in Domain 1
e∇2ψ = 0 in Domain 2, where Domain 1 is the region defined by y<=h and Domain 2 is the region defined by
浏览 4提问于2020-05-22得票数 0
这个问题跟在这个()后面。
弗雷德里科建议我使用OOP，并给出了我最初问题的解决方案。现在，我想使用子模块(cf4.中子模块的优点)。
这是我的密码：
module grids
implicit none
type, abstract, public :: grid
contains
procedure (grid_init), deferred :: init
end type grid
type, public, extends(grid) :: grid2D
real, dimension(:,:), allocata
浏览 2提问于2022-05-30得票数 1
我正试图从一个地区获得价值。我知道有一种方法可以通过引入坐标来调用这些值：
var([(x1, x2, x3) #xvalues
(y1, y2, y3)]) #yvalues
但是，例如，如果我以前在var.setValue (接受区域)中使用'where‘attr来设置值，那么就会对它们有多少个点以及它们的坐标感到困惑。
var.setValue(20., where=(Y <= y1) & ((X <= x1) & (X > x2)))
是否有任何方法可以像在var.setValue中那样“查询”一个区域(where=‘.’)来获取这些
浏览 2提问于2021-03-19得票数 0
回答已采纳
我想要解一个含源项的流体流动问题。但是，我无法按照我所希望的那样，在最后用源项写出方程。是否有方法将FaceVariable转换为cellVariable。
我收到了这个错误
TypeError:系数不能是FaceVariable。
代码如下
from fipy import CellVariable, FaceVariable, Grid2D, DiffusionTerm, Viewer, ExponentialConvectionTerm
from fipy.tools import numerix
from fipy.meshes.gmshMesh import Gmsh2D
i
浏览 1提问于2018-08-17得票数 0
1回答
我能用FiPy或.模块？对于半径为“R”的圆柱体，相关的网格在其上创建O‘网格扫描样。从表面上看，“R”大于最外层的网格人脸中心，因此两者之间存在差异(图1.)而且"R“很难被解决。
我需要提到的所有外部网格面中心围绕缸壁(而不是顶部和底部的墙壁)。网格系统由"FiPy Gmsh“导入(FiPy版本= 3.1和Python2.7)。FiPy 'solver.mesh.getFaceCenters()‘为网格人脸中心获取一个简单的数组。我认为外六角形网格层与圆筒壁有关的外面可以通过以下代码来处理：
reduce(numpy.logical_and, (FORMU
浏览 1提问于2021-03-07得票数 0
回答已采纳
在我看来，FiPys对Gmsh的支持取决于所安装的FiPy版本，这听起来很合理。我在哪里可以找到关于哪个版本的gmsh在哪个版本的fipy中仍然受支持的信息，这样我就不必弄清楚哪些命令有效，哪些命令只会引发错误。 摘录FiPy文档： FiPy works with arbitrary polygonal meshes generated by Gmsh. FiPy provides two wrappers classes (Gmsh2D and Gmsh3D) enabling Gmsh to be used directly from python. The classes can be
浏览 14提问于2020-01-13得票数 0
回答已采纳
1回答
安装的最佳方法是什么(在ubuntu16.04机器上)？
我曾经尝试过conda create --name FIPY --channel conda-forge python=3 fipy numpy scipy matplotlib pysparse gmsh，但结果却产生了很多这样的冲突：
Collecting package metadata (current_repodata.json): done
Solving environment: failed with repodata from current_repodata.json, will retry with next re
浏览 9提问于2022-09-27得票数 0
回答已采纳
我已经尝试了几个小时来安装FiPy，我已经安装了Pip和许多其他的东西来让它工作。Pip成功地安装了许多我需要的东西，但是我不能让它在PySparse或FiPy上工作。为什么我尝试安装PySparse，我得到一个错误：
$ pip install pysparse
Downloading/unpacking pysparse
Could not find a version that satisfies the requirement pysparse (from versions: 1.1.1-dev, 1.2-dev, 1.2-dev202, 1.2-dev203, 1.2-dev21
浏览 2提问于2014-07-07得票数 2
2回答
我正在尝试创建一个网格，据我所知，这个问题来自于在它自己内部使用一个模板类的指针，这是合法的，直到我尝试使用它，当编译器抱怨的时候。我正在寻找一种方法来使用指针的一个类的模板在其自身使用使用的指针，并做我稍后会做的事情。我用g++版本5编译，我使用的编译命令是g++ *.cpp -o main -std=c++11，我得到的错误将跟随在代码片段之后。 struct Vector2D
{
Vector2D( ) { }
Vector2D( int x , int y ): x( x ) , y( y ) { } ;
int x , y ;
} ;
tem
浏览 102提问于2019-01-28得票数 2
回答已采纳
2回答
FiPy不工作
、、
在Mac (OSX10.10.2和Python2.7.6)上重新安装FiPy之后，我尝试了一些文档中的示例，但是它们没有工作。因此，我使用内置的测试例程(从python中键入"fipy.test()“)测试fipy.test，但是经过一长串好的doctest之后，我得到了这个日志，有两个失败：。
我也尝试了一个新的LTS Ubuntu安装，但我也遇到了同样的问题。有什么解决办法吗？在我看来，这看起来像是FiPy模块的问题，而不是平台的问题。
浏览 5提问于2015-04-03得票数 2
回答已采纳
3回答
最近，我通过FiPy将FiPy包安装到了我的Macbook上，其中包含了所有的依赖项。我在Python中将FiPy和NumPy作为包调用没有问题。
现在我已经开始工作了，我想看看这些例子。但是，我无法在我的计算机中找到“基本目录”或FiPy目录。
我怎样才能找到基本目录？如果我已经通过Macports安装了所有这些，那么我甚至有基本目录吗？
作为注释，我正在使用Python27。
求你了救命！谢谢。
浏览 4提问于2014-09-10得票数 1
回答已采纳
1回答
我正在使用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
浏览 1提问于2016-12-01得票数 0
1回答
根据FiPy网站的说法，你需要“至少一个求解器”，它们作为一个可能的“求解器”。我肯定已经安装了SciPy --我可以用import scipy导入它。我还没有正式安装带有setup.py或类似的东西的fipy，我只是把代码放在一个名为FiPy的文件夹中。然后，当我进入Python shell并输入import fipy时，我得到：
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "fipy/__init__.py", line
浏览 0提问于2018-05-10得票数 1
我对fipy并不熟悉，所以如果这是个愚蠢的问题，我很抱歉(而且似乎帮不了我)。但是，除了上述问题中的建议之外，是否有一种方法可以以人类可读的(或python可读的)形式存储fipy对象？这仅适用于单元格变量。如果我想做一些比默认的fipy查看器更花哨的/定制的绘图，我怎么做呢？
举一个简单的一维扩散为例：
from fipy import *
# USER-DEFINED PARAMETERS
nx = 100
dx = 0.1
D = 1.0
bound1 = 30
bound2 = 70
# PREPARED FOR SOLUTION
mesh = Grid1D(nx=nx, dx=dx
浏览 4提问于2017-09-06得票数 2
回答已采纳
我试图在Python3.6上安装FiPy。遇到了麻烦。我遵循了NIST网站的安装说明：。
但是，在最后一步pip install fipy中，我得到了一个我不明白的错误：
(MYFIPYENV) >pip安装fipy收集fipy缓存条目反序列化失败，条目忽略Cache条目反序列化失败，条目忽略
下载FiPy-3.1.3.tar.gz
命令python setup.py egg_info: Traceback (最近一次调用)：File ""，第1行，除了ImportError，e：^ SyntaxError:无效语法
命令"python setup.py egg
浏览 0提问于2018-02-28得票数 0
回答已采纳
1回答
为了了解FiPy是如何工作的，我想解决具有固定端点的：
w''''(x) = q(x,t), w(0) = w(1) = 0, w'(0) = w'(1) = 0.
为了简单起见，让q(x,t) = sin(x)。
如何在FiPy中定义和解决它？如何针对方程中唯一的自变量指定源项sin(x)？
from fipy import CellVariable, Grid1D, DiffusionTerm, ExplicitDiffusionTerm
from fipy.tools import numerix
nx = 50
dx = 1/n
浏览 0提问于2019-05-08得票数 1
回答已采纳
1回答
管状三维扫描啮合
、、
我想sweep_like三维啮合一个空心圆柱，在Fipy中有任何网格模块吗？如何通过Fipy或python中的其他库实现这一目标？
浏览 2提问于2020-04-23得票数 0
回答已采纳
是否有可能连接两个(或更多)网格？示例：
import fipy
meshA = fipy.Grid2D(nx=5, ny=3)
meshB = fipy.Grid2D(nx=5, ny=6)
mesh = fipy.vstack(meshA, meshB)
1D或2D网格也是如此。
中间问题:是否可以定义偏移量？
浏览 1提问于2015-09-21得票数 1
回答已采纳
1回答
我使用fipy来建模线性化的，它本质上是
我假设我可以将f(x)建模为一个边界条件。如果epsilon(x)是一个常量，那么fipy可以处理这个问题：
phi = CellVariable(mesh)
dielectric_solvent = 80.0
dielectric_inner = 4.0
LHS = (DiffusionTerm(coeff = dielectric_solvent))
RHS = phi
eq = LHS == RHS
dr = np.linalg.norm(mesh.faceCenters, axis=0)
mask = (dr<.5)
浏览 1提问于2014-04-09得票数 5
回答已采纳
在Windows10上，我使用当前的FiPy版本3.4.2.1运行Python3.7。
如果我运行circle.py，我会得到：
Traceback (most recent call last):
File "circle.py", line 182, in <module>
exec(fipy.tests.doctestPlus._getScript())
File "<string>", line 39, in <module>
File "C:\python37\lib\site-pac
浏览 14提问于2021-01-05得票数 0
我正在尝试将网格文件(.msh)从GMSH4.7.0导入到FiPy==3.4.2.1 (Windows10)(Python2.7.9)中。该网格是使用开源有限元软件(Cast3M)创建的，然后将其作为.MED文件导出到Gmsh。Gmsh使用相同数量的元素和节点完美地查看了我的网格，没有任何问题。然后，我继续从Gmsh中将网格导出为.msh文件。(文件-导出-msh文件-保存)。我在做：
mesh=fp.Gmsh3D('cylinder.msh') #After importing fipy
所以我得到了这个错误：
Traceback (most recent call last
浏览 95提问于2020-11-13得票数 0
回答已采纳
2回答
Lambda函数作为参数
、、、、
我第一次和兰巴斯合作。我应该编写一个函数walk()，它以lambda函数作为参数。
在标题中，我声明函数为：
template<class T>
void walk(T operation) const;
我们应该在.inl中定义函数，我已经这样做了：
template<class T>
void Sea::Grid2D<T>::walk(T operation) const{
for(auto a : Sea::Grid2D<T>::grid){
operation(a);
}
}
我的问题就在这一点
浏览 5提问于2021-06-14得票数 1
回答已采纳
我正在尝试用FiPy实现一个三元合金的相场凝固模型。我已经看过FiPy网站上提供的大多数阶段领域示例，我的模型类似于。
浓度的演化方程如下所示：
对于C_1和C_2 (C_3是溶剂)，应该解决这个问题。将浓度方程与相变方程耦合，得到D_i = D_i(phi)，h= h(phi)。
在所有三个项中，(C_i)的求解变量都存在非线性依赖关系，这使得我不确定如何在FiPy中定义它们。第一项(red)是具有非线性系数的扩散项，这应该很好，但我应该如何定义反扩散项和相位变换项？
我试着将它们定义为具有非线性系数的扩散项和对流项，但没有成功。因此，我希望得到一些关于我如何定
浏览 0提问于2016-10-05得票数 2
OpenFOAM部分：
我在OpenFOAM中建立了一个空心圆柱形的网格。在该模型中，流动将从外圆柱体壁向内筒壁流动；右旋定律用于啮合点如下：
我曾尝试过用OpenFOAM编译的paraFoam绘制人脸法线，但我做不到。可以用字形法线绘制边界网格面(位于顶部、底部、内壁和外圆柱壁上)，但不能绘制内部网格面。
Q OpenFOAM:是否有任何函数代码或方法在paraFoam中绘制它们？
Mesh文件已经导出以在python中绘制，如上图所示。由此可见，与其他旋转面法线相比，起始端重叠面上的网格面具有相反的方向。
Q OpenFOAM:这些面的法线方向会导致错误的结果吗？或者结果被控制，只靠右手
浏览 13提问于2021-04-09得票数 0
通过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
浏览 13提问于2021-02-28得票数 0
我正在使用FiPy解决一个包含积分项的PDE系统。但是，在FiPy中不允许使用集成术语。有人能帮我弄一下这个系统吗？System of PDE
浏览 22提问于2021-04-21得票数 0
2回答
我想在Python3.3 (Windows7或Linux/ubuntu-mint)上安装Fipy，可以吗？我需要使用Fipy来改进我的python3.3代码。我必须把它翻译成python2.7吗？有谁知道吗？有什么建议吗？
浏览 0提问于2012-12-19得票数 3
考虑以下方案：
[<Struct>]
type Grid2D<'T> =
val RowLength : int
val Data : 'T[]
new(rowLength, data) = { RowLength = rowLength; Data = data }
member this.Item
with get(rowIndex, columnIndex) =
this.Data.[rowIndex * this.RowLength + columnIndex]
浏览 5提问于2015-07-04得票数 7
回答已采纳
1回答
我有一个由18个线性(假设)PDE组成的系统。然而，它们都是非常笨拙的PDE。每次我尝试使用FiPy来解决耦合PDE系统时，都会出现一个错误。我自己解决了其中的大部分问题，但我似乎无法通过这一次。
首先，您将需要一个方程式列表，这些方程将在程序后面添加：
import time
from fipy import Variable, FaceVariable, CellVariable, Grid1D, Grid2D, ExplicitDiffusionTerm, TransientTerm, DiffusionTerm, Viewer,PowerLawConvectionTerm, Impli
浏览 1提问于2019-04-26得票数 1
回答已采纳
1回答
from fipy import *
nx = 50
dx = 1.
mesh = Grid1D(nx=nx, dx=dx)
phi = CellVariable(name="solution variable",
mesh=mesh,
value=0.)
D = 1.
valueLeft = 1
valueRight = 0
phi.constrain(valueRight, mesh.facesRight)
phi.constrain(valueLeft, mesh.facesLeft)
浏览 0提问于2015-04-09得票数 1
我试图在python中使用FiPy来解决pde问题。它显示的是RuntimeError: Factor正好是奇异的。这不是在这些方程中的情况，而是在我计算的大多数方程中。
(1+αv)uv/(1+h(1+αv)u) +D∇^2 u
(1+αv)uv/(1+h(1+αv)u) -v+D∇^2v
from fipy import *
nx=ny=100
dx=dy=0.25
L=dx*nx
dt=0.01
sigma=10.0
h=0.1
C=0.8
alp=0.57 #alpha
mesh =Grid2D(dx=dx,dy=dy,nx=nx,ny=ny)
u=CellVariable(nam
浏览 14提问于2020-07-31得票数 0
回答已采纳
1回答
我希望使用FiPy解决扩散方程，并阅读了他们的一些文档，但似乎找不到任何与编写扩散项相关的内容，该扩散项包括作为自变量(即空间)的函数的附加项。我找到的最接近的东西是在FAQ上，他们建议将额外的术语重写为ConvectionTerm。然而，我相信这只适用于附加项是解变量的函数而不是自变量的情况。例如，我正在尝试求解一个具有以下扩散项的一维扩散方程(其中导数是w.r.t.到自变量x，y是解变量)： D * sin(x) * Div_x {sin(x) * Grad_x {y}} 我觉得这是一个非常简单的表达式，但是我找不到如何用FiPy表示法来表达它。任何帮助都将不胜感激！ 精确代码： fro
浏览 12提问于2019-09-15得票数 2
回答已采纳
我想在下面的图片中对柱状骨折进行网格处理：enter image description here 网格中的单元具有6个面的六面体形式，8个nodes.This网格将在fipy中用于求解扩散方程。(径向流动)有人能帮助回答我关于使用gmsh代码的网格的问题吗？因为这是与fipy兼容的唯一方法。我很感谢你的帮助。非常感谢。
浏览 70提问于2020-11-10得票数 0
我正在研究二维网格上微生物生物量(b1)分布的生物学模型。从生物量中产生一种蛋白质(p1)。生物质能在网格上扩散，而蛋白质则不扩散。只有当产生一定量的蛋白质(p > p_lim)时，生物量才会扩散。我试图通过使用一个虚拟单元变量z与扩散系数相乘来实现这一点，并且只在p> p_lim的单元格中将其设置为0到1。
条件很好，当在单元中达到p的临界量时，z被设置为1，然后发生扩散。然而，扩散仍然不能与我想要的速度工作，因为在计算扩散时，使用的是面变量，而不是单元本身的值。Z的面始终是带有z=1的单元格的平均值，而它的相邻的单元是z=0的均衡器。然而，我希望扩散以原来的速度工作，即使相邻的
浏览 2提问于2020-11-22得票数 0
回答已采纳
3回答
在我的中，我允许用户创建/指定一个首选的空间解析方法和一个首选的解析方法。
下面是一个空间分区的例子。
struct SpatialBase
{
virtual void a() = 0;
virtual void b() = 0;
...
};
struct Grid2D : public SpatialBase { ... override a() and b() ... };
struct QuadTree : public SpatialBase { ... override a() and b() ... };
struct World
{
std
浏览 3提问于2013-08-22得票数 0
回答已采纳
如何显式地将fipy网格中边界面法向的通量设置为特定值，而不限制该面内的通量组件？
Neumann边界条件可以被指定为：(1)垂直于边界面的通量的固定分量，或(2)作为该面上通量的完整规范。默认的fipy条件是前者(值= 0)，但显式方法(faceGrad.constrain)是后者。这个问题可以通过将以下代码添加到fipy 示例的末尾并注意不同的面部渐变结果来理解。
facesNeumann = mesh.exteriorFaces & ~facesTopLeft & ~facesBottomRight
print 'grad(phi) with default Ne
浏览 48提问于2016-09-18得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
