Page 356 - Process Modelling and Simulation With Finite Element Methods

Electrokinetic Flow 343

...
% Solve dynamic problem
fem.sol=femtime(fem, ...
'tlist', 0:0.1000000000000000l:tau, ...
' atol ,
{'u',0.01,'v',0.01,'p',Inf,'Y',0.001,'phi',0.001,'lmu',0.001,
.
.
.
'lrnv',0.001},
'rtol', 0.01, ...
'jacobian','equ', ...
.
'mass', 'full', . .
'ode', 'fldaspk', ...
1 odeopt ,
I
[I
struct('InitialStep',{ },'MaxOrder',{5},'MaxStepi,{[]}, ...
...
'maxiter',{6},'estrat',{O},'complex',{O}),
'out', 'sol', ...
'stop', 'on', . .
.
'init', init, ...
'report', 'on', ...
'timeind','auto', ...
' index2v' , [I , . . .
'indexlv' [I,. . .
,
'consistent', 'bweuler', ...
'krylov', 'direct', ...
'context','local', ...
.
'sd', 'Off', . .
'nullfun', 'flnullorth', ...
'blocksize',5000, ...
...
'solcomp',{~p','u','v','Y','phi','lmu','lmv~},
'linsolver', 'matlab', ...
'uscale', 'auto') ;
% Save current fern structure for restart purposes
pppppapppppppppp~epppWB~~p~ppppp~ppp~pppppp~~p~pp~~~ppppp~ppp~~pppp
~
~
~
o
0
o
o
o
~
.
~
~
0
0
~
0
~
0
0
.
0
0
0
0
0
o
~
0
o
o
o
.
b
o
o
~
~
~
o
~
.
o
.
o
o
~
0
~
~
~
0
0
~
~
~
D
~
0
0
0
~
0
0
SPP
000
f emO=f ern;
feml=fem;
for k=l:2*cycles-l
swap=phia;
phia=phib;
phib=swap;
...
% Solve nonlinear problem
fem.sol=femnlin(fem,. . .
.
'out', 'sol', . .
'stop', 'on', ...
'init', init, ...
'report', 'on', ...
'context', 'local', ...
.
'sd', 'Off', . .
'nullfun', 'flnullorth', ...
'blocksize',5000, ...
~solcomp',(~p~,'u','v','phi','lmu','lmv~),..
...
'linsolver','matlab',
.
'bsteps' 0, . .
,
'ntol', 9.9999999999999995e-007, ...
'hnlin' 'off I,. . .
,