/[ascend]/trunk/models/johnpye/fprops/solve_ph.c
ViewVC logotype

Log of /trunk/models/johnpye/fprops/solve_ph.c

Parent Directory Parent Directory | Revision Log Revision Log


Links to HEAD: (view) (download) (as text) (annotate)
Sticky Revision:

Revision 2664 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Fri Jan 18 06:02:15 2013 UTC (10 years, 8 months ago) by jpye
File length: 9687 byte(s)
Diff to previous 2662
Trying to debug fprops_triple_point for Toluene with pengrob correlation. Something strange is happening with fratio.

Revision 2662 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Fri Jan 18 00:34:51 2013 UTC (10 years, 8 months ago) by jpye
File length: 9689 byte(s)
Diff to previous 2661
Added some C tests for sat and solve_ph routines.
Added parameter 'source' to fprops_fluid, to allow data with a particular source to be specified (optionally). This allows forced testing explicitly of the RPP data even when Helmholtz data is also available.
Tweaked some convergence parameters in sat.c and solve_ph.c to get convergence across broader range of parameters for more fluids (more testing required).

Revision 2661 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Jan 17 00:15:12 2013 UTC (10 years, 8 months ago) by jpye
File length: 9291 byte(s)
Diff to previous 2660
Update GPL header, debugging issue with (p,h) for Water.

Revision 2660 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Jan 16 05:57:03 2013 UTC (10 years, 8 months ago) by jpye
File length: 9359 byte(s)
Diff to previous 2658
Working on updating rankine_fprops and associated models to work with new fprops2 code.
Some issue discovered with (p,h) for water (added python/solve_ph1.py to check it).
Next cunit tests to drive some ASCEND models embedding FPROPS (note use of slvreq for this).

Revision 2658 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Jan 16 02:37:20 2013 UTC (10 years, 8 months ago) by jpye
File length: 9144 byte(s)
Diff to previous 2654
Working on problem with (p,h) solution with new FPROPS code.

Revision 2654 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Jan 16 02:04:07 2013 UTC (10 years, 8 months ago) by jpye
File length: 9142 byte(s)
Diff to previous 2651
Merged (copied) fprops2 branch to trunk. This copy overwrites the trunk version of fprops,
and history is preserved via the fprops2 sequence. I hope that's OK...

Revision 2651 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Dec 13 07:29:48 2012 UTC (10 years, 9 months ago) by jpye
File length: 7996 byte(s)
Diff to previous 2308
Fixing GPL header, removing postal address (rpmlint incorrect-fsf-address)

Revision 2308 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Sep 2 03:56:14 2010 UTC (13 years ago) by jpye
File length: 8062 byte(s)
Diff to previous 2290
Turn off FPE in solve_ph.c (otherwise causes problems on Ubuntu 10.04 with rankine_co2 model).

Revision 2290 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 15 10:11:54 2010 UTC (13 years, 1 month ago) by jpye
File length: 8060 byte(s)
Diff to previous 2288
Trying to catch errors associated with pressure below triple point for CO2.

Revision 2288 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 15 06:16:13 2010 UTC (13 years, 1 month ago) by jpye
File length: 8031 byte(s)
Diff to previous 2287
Fixed some more FPE issues with pt->l going off the end of the array!
satcvgc working everywhere now except exactly at sat_p for p=p_t.

Revision 2287 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 15 01:19:21 2010 UTC (13 years, 1 month ago) by jpye
File length: 7845 byte(s)
Diff to previous 2271
Trying to track down an error that's started occurring with CO2 saturation curves.

Revision 2271 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 8 12:17:00 2010 UTC (13 years, 1 month ago) by jpye
File length: 7661 byte(s)
Diff to previous 2270
Removed p_t, p_c from HelmholtzData, added routines that calculate those values directly.
'satcvgc.py' now runs without errors for 100% of data points.

Revision 2270 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 8 09:07:43 2010 UTC (13 years, 1 month ago) by jpye
File length: 7611 byte(s)
Diff to previous 2269
Seem to have fixed the solve_ph problems, but new problems in fprops_sat_p have arisen.

Revision 2269 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Sun Aug 8 05:40:41 2010 UTC (13 years, 1 month ago) by jpye
File length: 7439 byte(s)
Diff to previous 2268
Avoid stepping outside saturation region in fprops_sat_hf.

Revision 2268 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Fri Aug 6 07:31:35 2010 UTC (13 years, 1 month ago) by jpye
File length: 7394 byte(s)
Diff to previous 2267
Improved first guess has eliminated many of the errors, still problems though.

Revision 2267 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Fri Aug 6 00:54:48 2010 UTC (13 years, 1 month ago) by jpye
File length: 6590 byte(s)
Diff to previous 2266
Adding better starting guess for T,rho when solving (p,h) in compressed subcritical liquid region.

Revision 2266 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Aug 5 13:53:01 2010 UTC (13 years, 1 month ago) by jpye
File length: 6182 byte(s)
Diff to previous 2262
The (p,h) function is converging almost everywhere. The home-made Newton iteration wasn't starting well for supercritical states
when the starting guess had rho = rho_c, so changing slightly away from fixed almost all the problems.

Revision 2262 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Aug 5 08:34:43 2010 UTC (13 years, 1 month ago) by jpye
File length: 5436 byte(s)
Diff to previous 2261
Deleted sat2.c and sat2.h, that code is not being used anywhere in ASCEND/FPROPS now. REFPROP algorithm replaced by Akasaka approach.
Moving sat3.c code into sat.h.

Revision 2261 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Aug 5 07:41:31 2010 UTC (13 years, 1 month ago) by jpye
File length: 5437 byte(s)
Diff to previous 2260
First attempt at Akasaka algorithm. But something's badly broken in the test suite, need to go find what.

Revision 2260 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Aug 5 06:28:01 2010 UTC (13 years, 1 month ago) by jpye
File length: 5336 byte(s)
Diff to previous 2259
Error messages to detect temperatures below critical point.

Revision 2259 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Thu Aug 5 03:37:14 2010 UTC (13 years, 1 month ago) by jpye
File length: 5096 byte(s)
Diff to previous 2257
Working on convergence of (p,h) over wide range.

Revision 2257 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Aug 4 09:45:52 2010 UTC (13 years, 1 month ago) by jpye
File length: 4758 byte(s)
Diff to previous 2256
Converging in supercritical region now.
Fixed some nan errors in helmholtz, needs cleaning up.

Revision 2256 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Aug 4 08:33:48 2010 UTC (13 years, 1 month ago) by jpye
File length: 4327 byte(s)
Diff to previous 2255
Converging... for one little case at least.

Revision 2255 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Aug 4 04:43:18 2010 UTC (13 years, 1 month ago) by jpye
File length: 3805 byte(s)
Diff to previous 2254
Adding bespoke (p,h) iteration code, but looks like some errors with the partial derivatives.

Revision 2254 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Wed Aug 4 02:56:52 2010 UTC (13 years, 1 month ago) by jpye
File length: 2088 byte(s)
Diff to previous 2252
For some reason the fdfsolver is not varying the guess values...

Revision 2252 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Tue Aug 3 14:30:53 2010 UTC (13 years, 1 month ago) by jpye
File length: 2061 byte(s)
Diff to previous 2251
Fixed solve_ph for saturation region.
Still not implemented 'nonsolver' for non-saturation cases.

Revision 2251 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Tue Aug 3 14:12:37 2010 UTC (13 years, 1 month ago) by jpye
File length: 2043 byte(s)
Diff to previous 2248
Added binding for solve_ph, but code is not yet completed.

Revision 2248 - (view) (download) (as text) (annotate) - [select for diffs]
Modified Tue Aug 3 08:48:59 2010 UTC (13 years, 1 month ago) by jpye
File length: 1941 byte(s)
Diff to previous 2247
More work on adding (p,h) solver in FPROPS.

Revision 2247 - (view) (download) (as text) (annotate) - [select for diffs]
Added Tue Aug 3 07:51:40 2010 UTC (13 years, 1 month ago) by jpye
File length: 1669 byte(s)
Some work towards (p,h) solver and derivatives function for FPROPS.

This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, enter a numeric revision.

  Diffs between and
  Type of Diff should be a

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22