function out = model
%
% helm2143.m
%
% Model exported on Dec 15 2011, 21:43 by COMSOL 4.2.0.187.

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


model = ModelUtil.create('Model');

model.modelPath('D:\Program Files\comsol_matlab');

model.name('helm2043.mph');

model.comments(['Helmholtz Coil\n\nA Helmholtz coil is a parallel pair of identical circular coils spaced one radius apart and wound so that the current flows through both coils in the same direction. This winding results in a uniform magnetic field between the coils with the primary component parallel to the axes of the two coils.']);

model.param.set('J0', '1[A/m^2]', 'Current density in coil');

%mod8.param.set('J0', '2[A/m^2]', 'Current density in coil');
%mod8.sol('sol1').run;
%mod8.result.export('data1').run;
%mod8来自在mat command窗口中输入的mod8=model

model.modelNode.create('mod1');

model.geom.create('geom1', 3);
model.geom('geom1').feature.create('wp1', 'WorkPlane');
model.geom('geom1').feature('wp1').geom.feature.create('sq1', 'Square');
model.geom('geom1').feature('wp1').geom.feature.create('sq2', 'Square');
model.geom('geom1').feature.create('rev1', 'Revolve');
model.geom('geom1').feature.create('sph1', 'Sphere');
model.geom('geom1').feature('wp1').geom.feature('sq1').set('pos', {'-0.4' '0.2'});
model.geom('geom1').feature('wp1').geom.feature('sq1').set('base', 'center');
model.geom('geom1').feature('wp1').geom.feature('sq1').set('size', '0.05');
model.geom('geom1').feature('wp1').geom.feature('sq2').set('pos', {'-0.4' '-0.2'});
model.geom('geom1').feature('wp1').geom.feature('sq2').set('base', 'center');
model.geom('geom1').feature('wp1').geom.feature('sq2').set('size', '0.05');
model.geom('geom1').feature('rev1').selection('input').set({'wp1.sq1' 'wp1.sq2'});
model.geom('geom1').run;

model.selection.create('sel1', 'Explicit');
model.selection('sel1').geom('geom1', 3, 2, {'exterior'});
model.selection('sel1').set([2 3]);
model.selection('sel1').name('Coil');

model.material.create('mat1');
model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat1').propertyGroup('def').func.create('cs', 'Analytic');

model.physics.create('mf', 'InductionCurrents', 'geom1');
model.physics('mf').feature.create('ecd1', 'ExternalCurrentDensity', 3);
model.physics('mf').feature('ecd1').selection.set([2 3]);

model.mesh.create('mesh1', 'geom1');
model.mesh('mesh1').feature.create('ftet1', 'FreeTet');
model.mesh('mesh1').feature('ftet1').feature.create('size1', 'Size');
model.mesh('mesh1').feature('ftet1').feature('size1').selection.geom('geom1', 3);
model.mesh('mesh1').feature('ftet1').feature('size1').selection.set([2 3]);

model.view('view1').set('renderwireframe', true);
model.view('view2').axis.set('xmin', '-1.1382290124893188');
model.view('view2').axis.set('xmax', '1.1382290124893188');

model.material('mat1').name('Air');
model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'});
model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'});
model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T');
model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'});
model.material('mat1').propertyGroup('def').func('rho').set('plotargs', {'pA' '' ''; 'T' '' ''});
model.material('mat1').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'});
model.material('mat1').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)');
model.material('mat1').propertyGroup('def').func('cs').set('args', {'T'});
model.material('mat1').propertyGroup('def').func('cs').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'});
model.material('mat1').propertyGroup('def').func('cs').set('plotargs', {'T' '' ''});
model.material('mat1').propertyGroup('def').set('relpermeability', {'1' '0' '0' '0' '1' '0' '0' '0' '1'});
model.material('mat1').propertyGroup('def').set('relpermittivity', {'1' '0' '0' '0' '1' '0' '0' '0' '1'});
model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.4');
model.material('mat1').propertyGroup('def').set('electricconductivity', {'0[S/m]' '0' '0' '0' '0[S/m]' '0' '0' '0' '0[S/m]'});
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)]' '0' '0' '0' 'k(T[1/K])[W/(m*K)]' '0' '0' '0' 'k(T[1/K])[W/(m*K)]'});
model.material('mat1').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]');
model.material('mat1').propertyGroup('def').addInput('temperature');
model.material('mat1').propertyGroup('def').addInput('pressure');

model.physics('mf').feature('ecd1').set('Je', {'-J0*z/sqrt(x^2+z^2)'; '0'; 'J0*x/sqrt(x^2+z^2)'});

model.mesh('mesh1').feature('size').set('hauto', 6);
model.mesh('mesh1').feature('ftet1').feature('size1').set('hmax', '0.05');
model.mesh('mesh1').feature('ftet1').feature('size1').set('hminactive', false);
model.mesh('mesh1').feature('ftet1').feature('size1').set('hcurveactive', false);
model.mesh('mesh1').feature('ftet1').feature('size1').set('hnarrowactive', false);
model.mesh('mesh1').feature('ftet1').feature('size1').set('hgradactive', false);
model.mesh('mesh1').run;

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

model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').attach('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
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.create('i1', 'Iterative');
model.sol('sol1').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature.create('sg1', 'SORGauge');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('sg1', 'SORGauge');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature.create('kp1', 'KrylovPreconditioner');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature('kp1').feature.create('sg1', 'SORGauge');
model.sol('sol1').feature('s1').feature.remove('fcDef');

model.result.dataset('dset1').selection.named('sel1');
model.result.dataset('dset1').selection.geom('geom1', 3, 2, {'exterior'});
model.result.dataset('dset1').selection.set([2 3]);
model.result.dataset.create('ps1', 'ParSurface');
model.result.create('pg1', 'PlotGroup3D');
model.result('pg1').feature.create('slc1', 'Slice');
model.result('pg1').feature.create('arwv1', 'ArrowVolume');
model.result('pg1').feature.create('surf1', 'Surface');
model.result.export.create('data1', 'Data');

model.sol('sol1').feature('st1').name('Compile Equations: Stationary');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('i1').set('linsolver', 'fgmres');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature('kp1').set('prefun', 'bicgstab');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature('kp1').set('iterm', 'iter');
model.sol('sol1').runAll;

model.result.dataset('dset1').selection.named('sel1');
model.result.dataset('dset1').selection.geom('geom1', 3, 2, {'exterior'});
model.result.dataset('dset1').selection.set([2 3]);
model.result.dataset('ps1').set('parmax1', '0.2');
model.result.dataset('ps1').set('parmax2', '0.2');
model.result.dataset('ps1').set('expry', 'sqrt(s1*s1-s2*s2)');
model.result.dataset('ps1').set('exprz', 's2');
model.result('pg1').name('Magnetic Flux Density (mf)');
model.result('pg1').set('title', 'Slice: Magnetic flux density norm (T) Arrow: Magnetic field Surface: 1 ');
model.result('pg1').set('titleactive', false);
model.result('pg1').feature('slc1').set('quickplane', 'xy');
model.result('pg1').feature('slc1').set('quickznumber', '1');
model.result('pg1').feature('arwv1').set('expr', {'mf.Hx' 'mf.Hy' 'mf.Hz'});
model.result('pg1').feature('arwv1').set('descr', 'Magnetic field');
model.result('pg1').feature('arwv1').set('xnumber', '24');
model.result('pg1').feature('arwv1').set('ynumber', '10');
model.result('pg1').feature('arwv1').set('znumber', '1');
model.result('pg1').feature('arwv1').set('arrowscale', '0.5');
model.result('pg1').feature('surf1').set('expr', '1');
model.result('pg1').feature('surf1').set('unit', '');
model.result('pg1').feature('surf1').set('descr', '1');
model.result('pg1').feature('surf1').set('coloring', 'uniform');
model.result('pg1').feature('surf1').set('color', 'white');
model.result.export('data1').set('data', 'ps1');
model.result.export('data1').set('expr', {'mf.normB'});
model.result.export('data1').set('unit', {'T'});
model.result.export('data1').set('descr', {'Magnetic flux density norm'});
model.result.export('data1').set('filename', 'D:\Program Files\comsol_matlab\out1923.txt');
model.result.export('data1').run;

out = model;