JPoptimization_github/Ra_Reitz_load_data.m at master · unghee/JPoptimization_github · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
clear all
close all

%%
mechanism={'MFC'};

fuel_sim={'afteroptimize_3'};
fuel_name={'JetA'};
date = {'05_03_2017'};
equi =1;

pressure = [40];
for k=1:length(pressure)
    pressure_text{k}=['P',num2str(pressure(k)),'atm'];
end
num_cases=25;


for i=1:length(mechanism)
    for j=1:length(fuel_sim)

%         phi=1 simulation results
        for k=1:length(pressure)
        location=['C:\Users\unghee\Dropbox\post_process','\',mechanism{i},'\',fuel_name{1},'_',num2str(pressure(k)),'atm','_','phi',num2str(equi),'_',date{1},'\',fuel_sim{j}];

        sim.(mechanism{i}).(fuel_sim{j}).(pressure_text{k})=read_ignition_delay(location,num_cases);

        end


    end
end

directory=['C:\Users\unghee\Dropbox\post_process','\',mechanism{1},'\',fuel_name{1},'_',num2str(pressure(1)),'atm','_','phi',num2str(equi),'_',date{1}];

cd(directory);
cd(fuel_sim{1});
filename = ['simulation_result','_',fuel_sim{1},'_',num2str(pressure(1)),'atm','_','phi','_',num2str(equi),'.mat'];
save(filename,'sim')