Clear Filters
Clear Filters

Error using System-Level Refrigeration Cycle (2P)

16 views (last 30 days)
while shifting from r134a to ammonia i get this strange error message :"The values in the table grid vector are not strictly monotonic".
the same cycle run normally with r600a or with r1234yf. can you guys help me solving this issue? I'd really like to test the ammonia behaviour.

Answers (2)

Yifeng Tang
Yifeng Tang on 30 Jun 2023
Hi Carlo,
We believe the error has to do with the fact that NH3 has a non-monotonic supercritical line. It needs to be fixed in the future.
For now, the easiest workaround is to use the Two-Phase Fluid Properties (2P) block instead, and get NH3 properties from REFPROP or Coolprop for a pressure range below the supercritical line. You can use the twoPhaseFluidTables command to generate the tables needed. You will need either REFPROP or Coolprop on your computer though.
I tried this:
NH3props = twoPhaseFluidTables([10,2000],[0.01,10],25,25,60,'Ammonia','C:\Program Files (x86)\REFPROP\')
gcb % select the Two-Phase Fluid Properties (2P) block when you do this
And it seems to work OK.
carlo ferro
carlo ferro on 8 Jul 2023
Thank you for your reply. In fact we have followed this approach and now it seems to be working. Now I should control the 2P cooling block in another way: limit the maximum mechanincal power avaiable to a certain value and varying the mass flow of the gas according to that. do you think it is possible to achieve this goal?
Yifeng Tang
Yifeng Tang on 11 Jul 2023
That I am not sure.
I would imagine it's possible, since you can use the pwr output from the block as a feedback signal and adjust the S input (load fraction). You may need to add a PS transfer function to provide some time delay.
On the other side, the gas mass flow on the condenser side can be used to control the refrigerant pressure / sat. temperature in the condenser. It probably also has an effect on other parts of the cycle, but I'm not sure.

Sign in to comment.

carlo ferro
carlo ferro on 3 Jul 2023
Hi Yifeng,
thank you for your reply. We are trying with Coolprop due to the fact that is opensource, and it seems to be an excel addon. Do you think that copying the data of nh3 from the excel table into the interface mask would work or do we have to replace all the function script?
carlo ferro
carlo ferro on 3 Jul 2023
i have another question:
what do you mean in line 2 with
"gcb % select the Two-Phase Fluid Properties (2P) block when you do this"
how can I paste the Two-Phase Fluid Properties (2P) block in the command line 2 ?
Yifeng Tang
Yifeng Tang on 3 Jul 2023
gcb will give you the path name of the currently selected block in Simulink. So you go to your model, select the Two-Phase Fluid Properties (2P) block by left click on it with your mouse, switch to MATLAB, run gcb command, or assign the output of gcb to a variable.
I know, it's confusing for the first time, lol. Same thing as this short example:

Sign in to comment.


Find more on Two-Phase Fluid Library in Help Center and File Exchange




Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!