# Problem with solving for unstructured mesh

 Register Blogs Members List Search Today's Posts Mark Forums Read

 November 22, 2006, 09:03 Hi, people! I wrote my own #1 Member   Efimenko Evgenii Join Date: Mar 2009 Location: Nizhnij Novgorod, Russia Posts: 52 Rep Power: 10 Sponsored Links 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? Thank in advance. Efimenko Evgeny

 November 22, 2006, 09:49 I don't know much about using #2 Senior Member   Michael Prinkey Join Date: Mar 2009 Location: Pittsburgh PA Posts: 334 Rep 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.

 November 24, 2006, 10:17 Which form of Maxwell's equati #3 Member   E. David Huckaby Join Date: Mar 2009 Posts: 57 Rep Power: 10 Which form of Maxwell's equations have you implemented ?

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post bne FLUENT 0 June 19, 2007 07:36 evgenii OpenFOAM Running, Solving & CFD 5 November 27, 2006 05:37 evgenii OpenFOAM Running, Solving & CFD 0 November 22, 2006 09:18 evgenii OpenFOAM Running, Solving & CFD 0 November 22, 2006 09:13 san FLUENT 2 April 3, 2006 23:37