# Hello Foamers.....

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

 July 12, 2011, 01:15 Hello Foamers: need mhdFoam help #1 Senior Member   Tushar Chourushi Join Date: Jul 2009 Location: IIT-Indore, India Posts: 319 Blog Entries: 1 Rep Power: 10 I am working on mhdFoam from past few weeks, my problem comes when I do change the transportProperties: since I want to check for the air as a conducting fluid so I am putting the Properties of Air values in it.. and studying the effect of variation on boundary layer when it flows over the flat plate... In my case velocity is quite small of 4 m/s... when I run my case it runs for few iteration and then courant no. keeps on increasing after few iteration and even the variables U, P, B, and pB stop converging after some iteration.. The boundary condition in my case are: 1) velocity inlet and pressure outlet.. 2) Far-field at the surrounding and no-slip at walls... Even i tried with the other BC's: 1) pressure inlet and velocity outlet... 2) Far-field at the surrounding and no-slip at walls... But the solution didn't converge after some iteration and shows the following errors: #### solver stops after some run ####### Courant Number mean: 4.08481e+77 max: 3.35924e+82 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 505337, No Iterations 1001 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.1354e+07, No Iterations 1001 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.30113e+06, No Iterations 1001 DICPCG: Solving for p, Initial residual = 1, Final residual = 6.32227, No Iterations 1001 time step continuity errors : sum local = 1.50302e+103, global = -3.0307e+97, cumulative = -3.0307e+97 DICPCG: Solving for p, Initial residual = 0.999369, Final residual = 4.20609, No Iterations 1001 time step continuity errors : sum local = 2.56876e+109, global = 2.06169e+102, cumulative = 2.06166e+102 #0 Foam::error:rintStack(Foam::Ostream&) in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so" #4 Foam::fvMatrix::solve(Foam::dictionary const&) in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libfiniteVolume.so" #5 in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/plasmairFoam" #6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6" #7 in "/home/iitgn/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/plasmairFoam" Floating point exception I would like to know, Has anyone tried this case earlier? Last edited by Tushar@cfd; July 12, 2011 at 01:44.

 July 12, 2011, 02:01 #2 Senior Member   Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 14 Does your case converge without the MHD-part?

 July 12, 2011, 02:59 Hi bernhard, #3 Senior Member   Tushar Chourushi Join Date: Jul 2009 Location: IIT-Indore, India Posts: 319 Blog Entries: 1 Rep Power: 10 I got your point.. You are right the equations won't be sufficient for this case as it won't reach convergence.. I need to modify as per my problem case... Thank you..

 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 21kalee OpenFOAM Running, Solving & CFD 7 December 19, 2012 18:18 ternik OpenFOAM 2 October 26, 2010 11:35 hk318i OpenFOAM 16 June 9, 2010 08:30 Miroo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 0 November 20, 2009 11:46 21kalee OpenFOAM Running, Solving & CFD 0 December 29, 2007 08:36

All times are GMT -4. The time now is 16:56.