CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem with solving for unstructured mesh (http://www.cfd-online.com/Forums/openfoam-solving/59922-problem-solving-unstructured-mesh.html)

 evgenii November 22, 2006 08:53

Hi, people! I wrote my own

Hi, people!

I wrote my own solver to solve electromagnetic problems. This program solves set of curl Maxwell equations. and now I am testing it solving problems for resonators. I used it for quite complicated geometry as cross-like volume, it works fine even for quite coarse mesh.

Then I decided to test it for cylindrical geometry, but I failed. I got really unstructured mesh and for this mesh solution exponentionally grows to infinity. Everything that I changed is structured ortogonal mesh to unstructured non-ortogonal mesh.

Could anyone give me advice or any assumption how this problem can be solved?

Efimenko Evgeny

 dmoroian November 22, 2006 11:34

If by "really unstructured" yo

If by "really unstructured" you ment stretched, then this could be the problem.
Try running checkMesh first to see the quality of your mesh.

Dragos

 evgenii November 22, 2006 15:53

Thank you for advise! But I am

Thank you for advise! But I am not common with this utility. Can you tell me something about following output?

Checking geometry...
Boundary openness in x-direction = 4.78574e-20
Boundary openness in y-direction = -3.1552e-20
Boundary openness in z-direction = 2.21605e-18
Boundary closed (OK).
Max cell openness = 6.35275e-22 Max aspect ratio = 1.0484. All cells OK.

Minumum face area = 3.08352e-07. Maximum face area = 1.81143e-06. Face area magnitudes OK.

Min volume = 3.22811e-10. Max volume = 1.72756e-09. Total volume = 5.02447e-05. Cell volumes OK.

Mesh non-orthogonality Max: 19.1576 average: 3.7965
Non-orthogonality check OK.

Face pyramids OK.

Max skewness = 24.4986 percent. Face skewness OK.

Minumum edge length = 0.000486087. Maximum edge length = 0.00151374.

All angles in faces are convex or less than 10 degrees concave.

Face flatness (1 = flat, 0 = butterfly) : average = 0.999998 min = 0.998494
All faces are flat in that the ratio between projected and actual area is > 0.8

Geometry check done.

Regards,
Efimenko Evgeny

 dmoroian November 23, 2006 03:44

Well, it basically says that y

Well, it basically says that your mesh is good. That means your choice of discretisation is probably not adequate.
I am not familiar with Maxwell equations (the curl form). Are they parabolic or hyperbolic, linear or nonlinear?
State your discretisation schemes for the different terms. Also, what solvers are set for different variables? (system/fvSchemes, system/fvSolution)

Dragos

 olwi November 24, 2006 21:36

A simple discretisation that w

A simple discretisation that works fine on cartesian grids will fail badly on non-structured grids. I have worked only with liquid-metal MHD, which is easier than the full Maxwell's. Already then you need to make sure you conserve a lot of things. You may for example have to compute certain terms on the faces, to assure that face fluxes are conservative. Then cell centre values can be reconstructed from face fluxes (there's a reconstruct method for that). It's not trivial, and I think you would probably have to skim the literature for algorithms.

If your new geometry is cylinder symmetric, however, you can still make an orthogonal grid in cylindrical coordinates... Maybe your simple discretization might be enough then.

/Ola

 evgenii November 27, 2006 05:37

Hi, all! Thanks greatly for

Hi, all!