[Sethu]

Dear all,

Im trying to plot min/max of pressure (p) in my domain using pyFoamPlotwatcher.py along with massflow rate. Im able to plot for massflow rate, but min(p) values are not getting plotted.

My log file looks like this,

//log_file_sample_start//
...
ExecutionTime = 46.73 s ClockTime = 51 s

MassFlows: outlet = 0.000254 inlet = -0.000254
fieldMinMax minmaxdomain output:
min(p) = -9.27757 at position (0.0549444 -0.0108001 -1.09092e-21)
max(p) = 13.3839 at position (0.207022 0.0251462 -1.63172e-21)
min(U) = (9.83972 -1.23842 0) at position (-0.0198086 0.000156695 -2.90727e-21)
max(U) = (10.0882 0.0457922 5.69625e-23) at position (0.000794772 0.00125782 -3.55445e-21)

//log_file_sample_end//


My customRegexp file is written like this,

massFlow
{
accumulation first;
enabled yes;
expr "MassFlows: outlet = (%f%) inlet = (%f%)";
name Custom01_massFlow;
theTitle "Massflow rate";
titles
(
outlet
inlet
);
type regular;
}

min(p)
{
expr "min(p) = (%f%) at position \(%f% %f% %f%\)";
name Custom02_min(p);
theTitle "min_p";
titles
(
min_p
);
type regular;
}


As I mentioned earlier, massflow values are plotted but min(p) value is not getting plotted. Moreover, i didnt receive any error.

Pls give me your suggestions.


Sethu

[arieljeds]

Hi, I'm having trouble using the customRegexp and swak4Foam function objects (following from one of the OFW presentation examples). I'm getting the following error when I try to run a case with some swak4Foam expressions in controlDict and a customRegexp file:

PyFoam WARNING on line 144 of file /usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Execution/ParallelExecution.py : which can not find a match for waveFoam . Hoping for the best 
Traceback (most recent call last):
  File "/usr/local/bin/pyFoamPlotRunner.py", line 5, in <module>
    pkg_resources.run_script('PyFoam==0.6.5', 'pyFoamPlotRunner.py')
  File "/usr/lib/python2.7/dist-packages/pkg_resources.py", line 499, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/usr/lib/python2.7/dist-packages/pkg_resources.py", line 1235, in run_script
    execfile(script_filename, namespace, namespace)
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/EGG-INFO/scripts/pyFoamPlotRunner.py", line 5, in <module>
    PlotRunner()
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Applications/PlotRunner.py", line 60, in __init__
    **kwargs)
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Applications/PyFoamApplication.py", line 424, in __init__
    result=self.run()
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Applications/PlotRunner.py", line 135, in run
    jobId=self.opts.jobId)
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Execution/GnuplotRunner.py", line 171, in __init__
    plottingImplementation=plottingImplementation)
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Execution/GnuplotRunner.py", line 74, in __init__
    plottingImplementation=plottingImplementation)
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/Execution/AnalyzedCommon.py", line 278, in createPlots
    endTime=custom.end))
  File "/usr/local/lib/python2.7/dist-packages/PyFoam-0.6.5-py2.7.egg/PyFoam/LogAnalysis/RegExpLineAnalyzer.py", line 71, in __init__
    self.exp=re.compile(self.strExp,reFlags)
  File "/usr/lib/python2.7/re.py", line 190, in compile
    return _compile(pattern, flags)
  File "/usr/lib/python2.7/re.py", line 242, in _compile
    raise error, v # invalid expression
sre_constants.error: unbalanced parenthesis
Exception in thread Thread-3 (most likely raised during interpreter shutdown):
Traceback (most recent call last):
  File "/usr/lib/python2.7/threading.py", line 551, in __bootstrap_inner
  File "/usr/lib/python2.7/threading.py", line 757, in run
  File "/usr/lib/python2.7/threading.py", line 403, in wait
  File "/usr/lib/python2.7/threading.py", line 258, in wait
<type 'exceptions.TypeError'>: 'NoneType' object is not callable

My controlDict is as follows:

  application		waveFoam;
    
    startFrom		latestTime; //startTime; 
     
    startTime		0;
    
    stopAt		endTime;
    
    endTime		300;
    
    deltaT		0.05;
    
    writeControl	adjustableRunTime;
    
    writeInterval	0.5;
    
    purgeWrite		0;
    
    writeFormat		ascii;
    
    writePrecision	6;
    
    writeCompression	off; 
    
    timeFormat		general;
    
    timePrecision	6;
    
    runTimeModifiable	true;
    
    adjustTimeStep	yes; 
    
    maxCo		0.9; // increase the Co as much as possible 
    
    maxAlphaCo		0.9;
    
    maxDeltaT		1;
  
  // **************************************************************************** //
  
   
   libs
   (
   	"libOpenFOAM.so"
   	"libsimpleSwakFunctionObjects.so"
   	"libswakFunctionObjects.so"
   	"libgroovyBC.so"
   );
   
   functions
   {

        forces
     {
	type forces;
	functionObjectLibs ("libforces.so");
	patches (cylinder); 
	rhoName rho;
	rhoInf 1028;
	CofR (0 0 0);
	outputControl timeStep;//outputTime;// 
	outputInterval 100;
	log true;
      }
      
      
      velocityStatistics
      {
      
     	type swakExpression;
     	verbose true; 
     	outputControlMode timeStep; 
     	outputInterval	   1;
     	valueType	internalField;
     	expression	"mag(U)";
     	accumulations (
     		weightedQuantile0.1
     		weightedAverage
     		weightedQuantile0.9
     		max
     		);
      }
      
      minStreamFunctionPos 
      {
	 type 	swakExpression; 
	 verbose true; 
	 outputControlMode	timeStep; 
	 outputInterval		1; 
	 expression "streamF";
	 valueType   internalField; 
	 accumulations (
	    	 min
	 	);
	 expressions "minPosition(streamF)";
      }
     	

   }

And finally my customRegexp file is:

velocity {
   	 theTitle  "Velocity statistics";
    	 expr "Expression velocityStatistics :  weightedQuantile0.1=(.+) weightedAverage=(.+) 
   	       weightedQuantile0.9=(.+) max=(.+)";
   	titles (
    		"10 %"
    		average
    		"90 %"
    		maximum
   	 );
  	 ylabel "[m/s]";
}

I can't see an unbalanced paranthesis... can anyone help me on this?

Thanks in advance

[RomainBou]

hello dear friends,

I am trying to plot the output of a faceSource functionObject that computes the total massflow at inlet.

Here is my log for one timestep:

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.0116454, Final residual = 5.6463e-08, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.00353414, Final residual = 3.34245e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00429777, Final residual = 8.15157e-09, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0390386, Final residual = 0.000158842, No Iterations 4
time step continuity errors : sum local = 4.00815e-08, global = 3.11859e-13, cumulative = 1.62584e-11
DICPCG:  Solving for PotE, Initial residual = 7.94541e-07, Final residual = 7.94541e-07, No Iterations 0
ExecutionTime = 5.58 s  ClockTime = 6 s

faceSource debit_sortant output:
    sum(inlet) of phi = -0.00164543

I am using pyFoam 0.6.5 and here is my customRegexp file :

massflowplot
{
  enabled   yes;
  expr     "    sum(inlet) of phi = (%f%)";
  name Custom_flow_plot;
  theTitle "Custom_flow_plot";
  titles
   (
    massflow
   );
  type regular;
}

Then i run

pyFoamPlotRunner.py mySolver

and it says "Reading regular expressions from /home/romain/OpenFOAM/romain-3.0.1/run/Cuve2_variations/Cuve2Cas1_variation_B/customRegexp".
But it doesn't plot and I don't understand why. Any idea ? Do I have to write something else somewere else ? Is the regexp correct ? thanks a lot for any help !

[Mars409]

Hi All,

How do you get pyFoamPlotWatcher to plot the added particle count and the escape count from solver.log?

The typical time step output for particles look like below:


Time = 0.875

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
Cloud: kinematicCloud
Current number of parcels = 5045
Current mass in system = 0.00647515
Linear momentum = (-0.002015723 -0.001078755 -0.0004498539)
|Linear momentum| = 0.002330069
Linear kinetic energy = 0.01457947
injector_1:
number of parcels added = 6546
mass introduced = 0.008395572
Parcel fate (number, mass) : patch cyclone_inletX
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch cyclone_outlet
- escape = 1501, 0.001920422
- stick = 0, 0


I wrote a Python script to print to a CSV file and plot it using Gnuplot:

#!/usr/bin/python

import re
import sys
from subprocess import call

if len(sys.argv) > 1 :
        TEXTFILE = sys.argv[1]
else : 
        TEXTFILE = 'solver.log'

if len(sys.argv) > 2 :
    CSVFILE = sys.argv[2]
else :
    CSVFILE = 'timeAddsEscapes.csv'
    
f1=open(CSVFILE, 'w+')

STATE0 = 'TIME'
STATE00 = 'SEARCH_ADDED'
STATE1 = 'SEARCH_OUTLET'
STATE2 = 'SEARCH_ESCAPE'

RE0 = 'Time = (\d+.\d+)'
RE00 = r'.*parcels added\s*=\s*(\d+)'
RE1 = r'\s*Parcel fate.*\:.*outlet$'
RE2 = r'\s*-\s*escape\s*=\s*(\d+)'



state = STATE0

print >>f1, "Time","Added", "Escapes"
with open(TEXTFILE, 'r') as f:
	for line in f:
		if state == STATE0 :
			pat = RE0
		elif state == STATE00 :
			pat = RE00
		elif state == STATE1 :
			pat = RE1
		else :
			pat = RE2
		match = re.match(pat, line)
		if not match :
			continue
		#if state == STATE0 or state == STATE2 :
		#print line
		if state == STATE0 :
			time = float(match.group(1))
			state = STATE00
		elif state == STATE00 :
			added = int(match.group(1))
			state = STATE1
		elif state == STATE1 :
			state = STATE2
		else :
			escapes = int(match.group(1))
			state = STATE0
			if added > 0 : 
				# print "Time = ", time, " added = ", added, " escapes = ",escapes
				print >>f1, time, added, escapes
f1.close()

call(["plotEscapes.sh", CSVFILE])

Here's a shell script, called by the name 'plotEscapes.sh' as in the above Python script, to plot the CSV values using Gnuplot:

#!/bin/bash
if (( $# != 1 )); then
    echo "Illegal number of parameters"
    cat <<EOT 
Usage:
    plotEscapes.sh <CSV filename>
EOT
exit
fi


gnuplot -persist << EOT
	set term X11 size 700, 900
	set title 'Particle escapes & Particles added vs Time'
	set xlabel 'Time/s'
	set ylabel 'Particle count'
	set grid
	set multiplot layout 2,1 columnsfirst
	plot for [i=2:3] '$1' using 1:i lt 1 with lines title columnhead at end 


	set logscale y
	plot for [i=2:3] '$1' using 1:i lt 2 with lines title columnhead at end
	unset multiplot
EOT