CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Convergence check on forces

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2011, 17:43
Default Convergence check on forces
  #1
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 16
andrea.pasquali is on a distinguished road
Hi,
I'm trying to add the convergence check on forces function following the residuals convergence check in simpleFoam.
For example, for the forces pressure, I added those lines:
Quote:
scalar eqnResidualFF = 1, maxResidualFF = 0;

eqnResidualFF = solve
(
fm.first().first() += sum(fN)
).initialResidual();

maxResidualFF = max(eqnResidualFF, maxResidualFF);
After compilation I got the error:
Quote:
error: invalid use of void expression
What I'm wronging? I think it does not like the expression:
Quote:
+=
Could anyone help me to go ahead?

Thanks

Andrea
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   March 16, 2011, 06:24
Default
  #2
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 16
andrea.pasquali is on a distinguished road
Hi,
I solved my problems and I did the new lib to check the error on the forces / moments.
At now I calculate the increment between it i and it i-1:

err = || (force)[i] - (forces)[i-1] ||

I attached the plots for increment, residuals and force.

I'd like to stop the calculation when the force becomes more or less constant (2000-3000 its) without waiting for residuals convergence check (because it is possible my Ux Uy Uz residuals don't decrease under the convergence declared in fvSolution).

I hoped to do that checking the increment but as you can see it is very dispersive.
Is it possible to do what I want to do?

Thanks in advance for any help and suggestion

Andrea
Attached Images
File Type: png increment.png (19.0 KB, 20 views)
File Type: png residulas.png (8.3 KB, 20 views)
File Type: png force.png (5.5 KB, 18 views)
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Forces in V1.6 terrybarnaby OpenFOAM Post-Processing 72 September 2, 2015 17:49
Forces calulated through pressure LVDH OpenFOAM Post-Processing 2 February 26, 2010 04:15
Calculate forces without hydrostatic pressure geir_oye FLUENT 4 November 12, 2009 10:12
changes to forces in 1.6 linnemann OpenFOAM Running, Solving & CFD 0 July 30, 2009 09:49
convergence problem with SIMPLER NURAY KAYAKOL Main CFD Forum 1 February 24, 1999 14:43


All times are GMT -4. The time now is 23:30.