% Reverse_test.m
%
% Model exported on Sep 20 2011, 13:28 by COMSOL 4.2.0.187.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');

model.modelPath('D:\Eigene Dateien\Simulations\COMSOL\3-Influence Grid\Simulations\Test');

model.modelNode.create('mod1');

model.geom.create('geom1', 2);

model.mesh.create('mesh1', 'geom1');

model.physics.create('spf', 'LaminarFlow', 'geom1');

model.study.create('std1');
model.study('std1').feature.create('stat', 'Stationary');

model.geom('geom1').feature.create('sq1', 'Square');
model.geom('geom1').feature('sq1').set('type', 'solid');
model.geom('geom1').feature('sq1').set('base', 'corner');
model.geom('geom1').feature('sq1').set('pos', {'-0.65' '-0.4'});
model.geom('geom1').feature('sq1').set('l', '0.75');
model.geom('geom1').run('sq1');
model.geom('geom1').run;

model.material.create('mat1');
model.material('mat1').name('Helium');
model.material('mat1').set('family', 'air');
model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.66');
model.material('mat1').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]');
model.material('mat1').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]');
model.material('mat1').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]');
model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat1').propertyGroup('def').func('eta').set('funcname', 'eta');
model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('eta').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'144.0' '900.0' '3.7251756E-6+6.83450863E-8*T^1-5.07299333E-11*T^2+2.36688744E-14*T^3'});
model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat1').propertyGroup('def').func('Cp').set('funcname', 'Cp');
model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('Cp').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'144.0' '900.0' '5200.0'});
model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat1').propertyGroup('def').func('rho').set('funcname', 'rho');
model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.004003/8.314/T');
model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.004003/8.314/T,pA)'; 'T' 'd(pA*0.004003/8.314/T,T)'});
model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat1').propertyGroup('def').func('k').set('funcname', 'k');
model.material('mat1').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('k').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('k').set('pieces', {'144.0' '900.0' '0.0294900023+5.07655059E-4*T^1-4.22501605E-7*T^2+2.1209438E-10*T^3'});
model.material('mat1').propertyGroup('def').addInput('temperature');
model.material('mat1').propertyGroup('def').addInput('pressure');
model.material('mat1').set('family', 'air');

model.mesh('mesh1').feature.create('ftri1', 'FreeTri');
model.mesh('mesh1').run;

model.result.create('pg1', 'PlotGroup2D');
model.result('pg1').run;
model.result('pg1').feature.create('mesh1', 'Mesh');
model.result('pg1').run;
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature('surf1').set('descr', ['Geschwindigkeit, Betrag Netzqualit' native2unicode(hex2dec('00e4'), 'Cp1252') 't']);
model.result('pg1').feature('surf1').set('expr', 'qual');
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;

model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('s1', 'Stationary');
model.sol('sol1').feature('s1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 25);
model.sol('sol1').feature('s1').feature.create('d1', 'Direct');
model.sol('sol1').feature('s1').feature('d1').set('linsolver', 'pardiso');
model.sol('sol1').feature('s1').feature('fc1').set('linsolver', 'd1');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 25);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.result.create('pg2', 2);
model.result('pg2').set('data', 'dset1');
model.result('pg2').feature.create('surf1', 'Surface');
model.result('pg2').feature('surf1').set('expr', {'spf.U'});
model.result('pg2').set('frametype', 'spatial');
model.result('pg2').name('Geschwindigkeit (spf)');
model.result.create('pg3', 2);
model.result('pg3').set('data', 'dset1');
model.result('pg3').feature.create('con', 'Contour');
model.result('pg3').feature('con').set('expr', {'p'});
model.result('pg3').set('frametype', 'spatial');
model.result('pg3').name('Druck (spf)');
model.result('pg3').feature('con').set('number', 40);

model.sol('sol1').runAll;

model.result('pg2').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').set('data', 'dset1');
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature('surf1').set('colortablerev', 'on');

mphplot(model,'pg1')