How to do Interface between Aspen Hysys and Matlab and perform Hybrid optimization

Hi

Is there anyone who can me help me please. I am doing LNG plant optimization. To perform the optimization I first did the simulation on Hysys and interface HYSYS with Matlab. The model is still now in steady-state condition. When I am using the compressor outlet pressure as the decision variables, The simulation is working and it can minimize the energy consumption.

But, when I am using the refrigerant mass flow rate or outlet pressure of throttle valves as decision variables, the simulation get's unconverged and the simulation stopped.

In addition to this, when I am using refrigerant mass fraction as decision variables and run the optimization it shows the following error "Invoke Error, Dispatch Exception: No such interface supported"

If I the convert the simulation into dynamic mode, the problem will be solved ?

whenever, i am trying to convert the model into dynamic it's showing that notification in this picture.

[Mario Malic]

Hi Ahmed,

I can't write to you privately, so I am commenting here.

Unfortunately, I can't help with your issue, since I do not use Aspen Hysis. If you want other people to help you (who know mentioned software) with your issue, you have to put specify the problem in detail.

[Rasel Ahmed]

Thank you @mario. Actually, I already solve this issue. Currently, i am facing another problem. Please check the following link. Can you please help me about it ?

https://www.mathworks.com/matlabcentral/answers/595414-how-can-i-use-the-refrigerant-mass-fraction-mole-fraction-value-for-optimization

[Mario Malic]

Well, I believe that's related to that software itself, not MATLAB, can only instruct you to try to set the value that MATLAB sets, manually into software and figure out why doesn't it work, or is it even supposed to work.

[Rasel Ahmed]

The problem is with Matlab. I need to keep the summation of 4 mass fraction value as 1. I have to write it as a constraint. the summation can't be more than or less than one. In each iteration matlab will find the optimum combination of the mass fraction, so the summation will be one. I stated it as following but it's not working:

M = Spreadsheet_name.Cell('B8').CellValue + Spreadsheet_name.Cell('B9').CellValue... 
+ Spreadsheet_name.Cell('B10').CellValue + Spreadsheet_name.Cell('B11').CellValue 
if M < 1, or M >1 
    M =1;
end 

Can you please make the correction in this code.

[Mario Malic]

Why would you even calculate this, if you want to set it to 1?

M = 1;

"But, when I am using the refrigerant mass flow rate or outlet pressure of throttle valves as decision variables, the simulation get's unconverged and the simulation stopped." - maybe you are setting values in your simulation that make your simulation unsolvable.

"Invoke Error, Dispatch Exception: No such interface supported" - maybe you are setting property value that cannot be changed.

I'd suggest that you read the Aspen documentation first. Do the change you want in the Aspen software, not in MATLAB, and see if it works.

[Rasel Ahmed]

The simulation is ok. I already checked everything.

it's not working in this way. I think I have to build some correlation among the mass fraction and convert them to code so that I can change the mass fractions by keeping the summation of mass fraction as 1.

[Mario Malic]

The rest of the code that is missing, is important then. You are saying that M has to be 1 or sum of those 4 cell values must be 1. You have to affect those cell values, so that their summation equals to 1. You have to check your constraints in optimization.

[Rasel Ahmed]

Hi Mario,

I think I couldn't explain the problem properly, here I do:

I have interfaced HYSYS and MATLAB using a spreadsheet.

Now, Hysys calculates values of variables x1, x2, x3, x4 and send to MATLAB, which optimizes them and send back the values to Hysys.

The sum of these four variables should be equal to one during optimization, that is x1+x2+x3+x4=1, that is if value of x1 is reduced by MATLAB then the values of remaining variables should be increased to make the sum equals 1.

I am not figuring out a way to correctly write in MATLAB. How should I formulate it?

[Mario Malic]

I am not sure how and if that's possible. Work through examples in your chosen algorithm, that'll give you an idea I am not sure if it's possible to.

Try with this (not sure if it will/should work), see fmincon(or ga) and input arguments Aeq = [1 1 1 1] ; beq = 1.

I won't reply to this thread anymore, as question is not specific..

[Walter Roberson]

if M < 1, or M >1 
  M =1;
end

Wrong way of doing or

if M < 1 || M > 1
  M = 1;
end

However, there are only two ways that the or could be false.

If M is 1 already then the test would be false, leaving M as 1. In which case M = 1 unconditionally would be the same. Unless, that is, you use the test to take other actions such as normalizing values, in which case M ~= 1 would be a better test.

The other way the test could fail is if M is nan. Comparison to nan is always false.

[Rasel Ahmed]

thanks walter for your comment. I will apply this logic in matlab

[Ting Liang]

Hi, Rasel Ahmed,

Excuse me.

Can you share me some documents or references about how to do optimization by connecting MATLAB and HYSYS? My email is [e-mail].

Thanks.

Best regards

Ting Liang

[Rasel Ahmed]

Hi liang,

I am sharing the codes to interface hysys with matlab, and providing another link of the details code. Have you tried to use mixed refrigerant mass fraction / mass flow as decision variables ?. I am facing some problem that I couldn't change the mixed refrigerant mass fraction values.

MyObject=actxserver('Hysys.Application');
global MySimCase;
FileNamePath = 'C3MR7v6';
MySimCase = MyObject.SimulationCases.Open([cd,strcat('\',...
    FileNamePath,'.hsc')]);
MySimCase.Visible = true;
MyOperations = MySimCase.Flowsheet.Operations; % Getting data from our simcase operations
MyFlowSheet = MySimCase.get('flowsheet'); %accessing to flowsheet
MyMaterialStreams = MyFlowSheet.get('MaterialStreams'); %accessing to simulated material streams
MyEnergyStreams = MyFlowSheet.get('EnergyStreams'); % Energy streams

https://www.mathworks.com/matlabcentral/answers/569845-hello-i-d-like-to-ask-you-something-concerning-optimization-communicating-with-aspen-hysys

[Ting Liang]

Hi, Rasel Ahmed,

Thanks a lot for your sharing.

Sorry, I am not working on mixed refrigerant, but on air, sorry that I cannot help with your problem now, but hope we can communicate with and help each other more later.

Thanks.

Best regards

Ting Liang

Answer 1

[bilal kazmi]

if you provide any help in this please share some file as iam also in search of specific file

[Rasel Ahmed]

Hi kazmi,

I am sharing the codes to interface hysys with matlab, and providing another link of the details code. Have you tried to use mixed refrigerant mass fraction / mass flow as decision variables ?. I am facing some problem that I couldn't change the mixed refrigerant mass fraction values.

MyObject=actxserver('Hysys.Application');

global MySimCase;

FileNamePath = 'C3MR7v6';

MySimCase = MyObject.SimulationCases.Open([cd,strcat('\',...

FileNamePath,'.hsc')]);

MySimCase.Visible = true;

MyOperations = MySimCase.Flowsheet.Operations; % Getting data from our simcase operations

MyFlowSheet = MySimCase.get('flowsheet'); %accessing to flowsheet

MyMaterialStreams = MyFlowSheet.get('MaterialStreams'); %accessing to simulated material streams

MyEnergyStreams = MyFlowSheet.get('EnergyStreams'); % Energy streams

https://www.mathworks.com/matlabcentral/answers/569845-hello-i-d-like-to-ask-you-something-concerning-optimization-communicating-with-aspen-hysys