#SOURCE YPOS 2014 HOW TO#

init_vertices ( nvert ) i_vert = 0 for x in xpos : for y in ypos : editor. open ( mesh, 2, 2 ) # Add the vertices (nodes) editor. linspace ( 0, height, ny ) # Create the mesh and open for editing mesh = df. linspace ( 0, length, nx ) ypos = numpy. ''' # The number of mesh entities nvert = nx * ny ncell = 2 * ( nx - 1 ) * ( ny - 1 ) # Positions of the vertices xpos = numpy. Should give exactly the same results as using the built in dolfin.RectangleMesh() class So, for posterity, here is a reimplementation of dolfin.RectangleMesh in a much more cumbersome (and flexible) way: import numpy import dolfin as df def create_mesh ( length, height, nx, ny, show = False ): '''

#SOURCE YPOS 2014 HOW TO#

While playing with implementing a simple solver for fluid flow in a tank by use of FEniCS to solve the Laplace equation I needed to create my own mesh in code.įEniCS is quite well documented, but I had to look at the source code for some of the mesh conversion routines to find out how to build a mesh from scratch. In tutorial but I have to change it to laplacian(y)=y0.I may end up doing a large project using FEniCS which is a collection of tools for automated, efficient solution of differential equations using the finite element method. In a rectangular domain that y0 is a velocity distribution that is obtained by solving the other problem in this domain and I should import it to this solution. Also M_PI is pi number.Ĭould you plz give me procedure to solve this equation: Maybe that's not the case in your version. Also, I added the 'Foam::' because otherwise I cannot compile it.

I use zero boudary condition in a 1 by 1 square. The second change was to have a source term which isĬompatible with my boundary conditions. You may also be able to do it in a better way. The first change is to be able to read and write your source term and see if you evaluate it right or wrong. S = Foam::sin(M_PI*xPos) + Foam::sin(M_PI*yPos)) My T file is giving values for inlet condition but not for the outlet, lowerwall, etc. << ' ClockTime = ' << runTime.elapsedClockTime() << ' s' Info<< 'ExecutionTime = ' << runTime.elapsedCpuTime() << ' s' Info<< 'Time = ' << runTime.timeName() << nl << endl įor (int nonOrth=0 nonOrth<=nNonOrthCorr nonOrth++) Info<< ' Calculating temperature distribution ' << endl įor (runTime++ !runTime.end() runTime++) LowerwallT = (2.0 -faceCentresLowerwall.y())*(1.0 - faceCentresLowerwall.y()) Label lowerwallPatchID = mesh.boundaryMesh().findPatchID('lowerwall') įvPatchScalarField& lowerwallT = T.boundaryField() Ĭonst vectorField& faceCentresLowerwall = mesh.Cf().boundaryField() UpperwallT = faceCentresUpperwall.y()*(4.0 - faceCentresUpperwall.y()) Label upperwallPatchID = mesh.boundaryMesh().findPatchID('upperwall') įvPatchScalarField& upperwallT = T.boundaryField() Ĭonst vectorField& faceCentresUpperwall = mesh.Cf().boundaryField() OutletT = 5.0 - faceCentresOutlet.y()*(2.0 - faceCentresOutlet.y()) Label outletPatchID = mesh.boundaryMesh().findPatchID('outlet') įvPatchScalarField& outletT = T.boundaryField() Ĭonst vectorField& faceCentresOutlet = mesh.Cf().boundaryField() InletT = faceCentresInlet.y()*(1.0 - faceCentresInlet.y()) Label inletPatchID = mesh.boundaryMesh().findPatchID('inlet') įvPatchScalarField& inletT = T.boundaryField() Ĭonst vectorField& faceCentresInlet = mesh.Cf().boundaryField() (I know this is a crude way of doing it). I however incorporated the values as a list and tried it out and it is working fine. Could you point out where exactly I am doing a mistake. But this time I got the non uniform values for inlet alone and repeated the same procedure for outlet, lowerwall and upperwall, but didn't see any change.

Solve poisson equation just add a source term - CFD Online Discussion Forums

2022.11