|November 22, 2006, 09:03||
Hi, people! I wrote my own
Join Date: Mar 2009
Location: Nizhnij Novgorod, Russia
Posts: 52Rep Power: 10
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?
Thank in advance.
|November 22, 2006, 09:49||
I don't know much about using
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 316Rep Power: 16
I don't know much about using finite volume solvers for coupled EM, but I have written a solver (not in OpenFOAM) to solve electric fields coupled to reacting flow. I don't think that you can count on the segregated solution approach to always be stable, at least not without some underrelaxation. The reason is simple. You have (at least) six fields...the D-vector and B-vector fields that all couple strongly to each other. D is based on the curl of B. B is based on the curl of D. That coupling is all external. So the linear system of an equation for D or B is only going to have a diagonal and right-hand side.
d(D0)/dt = function(B0,B1,B2)
d(B0)/dt = function(D0,D1,D2)
So there is no linear system to solve. You are stuck either doing explicit timestepping or an iterative scheme based only on successive substitution. Both are potentially unstable depending on your choice of timestep and/or underrelaxation factors. As a rule, unstructured problems are usually more difficult to solve and less stable than structured problems. Vector components aligned with grid faces, no skew terms...all nice things that vanish with the grid structure. Moving to unstructured meshes may simply have moved your solver out of the basin of convergence you were enjoying before. You can try underrelaxing the solution more or reducing the timestep.
A better solution to this stability problem is to solve the fully coupled system. I.e., solve for all components of D and B from a block structured system. However, that capability does not currently exist in OpenFOAM.
|Thread||Thread Starter||Forum||Replies||Last Post|
|solving problem !!! need help !!!||bne||FLUENT||0||June 19, 2007 07:36|
|Problem with solving for unstructured mesh||evgenii||OpenFOAM Running, Solving & CFD||5||November 27, 2006 05:37|
|Problem with solving for unstructured mesh||evgenii||OpenFOAM Running, Solving & CFD||0||November 22, 2006 09:18|
|Problem with solving for unstructured mesh||evgenii||OpenFOAM Running, Solving & CFD||0||November 22, 2006 09:13|
|problem in solving "wave generation" problem||san||FLUENT||2||April 3, 2006 23:37|