Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (geom_testpardiso) Browse

fennecs/matlab/dispespicFluid.m
50e38dd6641fgeom_testpardiso

Last Change
Raw File

dispespicFluid.m
View Options

function f=dispespicFluid(M)
fieldstep=1;
f=uifigure('Name','Fluid data');
mf=uipanel(f,'Position',[5 50 f.Position(3)-10 f.Position(4)-55]);
mf.AutoResizeChildren='off';
m=uipanel(f,'Position',[5 5 f.Position(3)-10 40]);
sgtitle(mf,sprintf('t=%0.5e s',M.t2d(fieldstep)))
sld = uislider(m,'Position',[10 30 0.6*m.Position(3) 3]);
sld.Value=fieldstep;
sld.Limits=[1 length(M.t2d)];
edt = uieditfield(m,'numeric','Limits',[1 length(M.t2d)],'Value',1);
edt.Position=[sld.Position(1)+sld.Position(3)+25 5 40 20];
edt.RoundFractionalValues='on';
Printbt=uibutton(m,'Position',[edt.Position(1)+edt.Position(3)+10 5 40 20],'Text', 'Save');
%Playbt=uibutton(m,'Position',[Printbt.Position(1)+Printbt.Position(3)+10 5 40 20],'Text', 'Play/Pause');
sld.ValueChangingFcn={@updatefigdata,edt,mf};
edt.ValueChangedFcn={@updatefigdata,sld,mf};
Printbt.ButtonPushedFcn={@plotGridButtonPushed};
[R,Z]=meshgrid(M.rgrid,M.zgrid);
Rinv=1./R;
Rinv(:,1)=0;
PlotEspic2dfluiddata(mf,M,fieldstep);
function PlotEspic2dfluiddata(fig,M,fieldstep)
%PlotEspic2dgriddata Plot the 2d data of espic2d at time step fieldstep
sgtitle(fig,sprintf('t=%0.5e s',M.t2d(fieldstep)))
ax1=subplot(2,2,1,'Parent',fig);
surface(ax1,M.zgrid,M.rgrid,M.N(:,:,fieldstep),'edgecolor','none');
xlim(ax1,[M.zgrid(1) M.zgrid(end)])
ylim(ax1,[M.rgrid(1) M.rgrid(end)])
xlabel(ax1,'z [m]')
ylabel(ax1,'r [m]')
title(ax1,'Position')
c = colorbar(ax1);
c.Label.String= 'n[m^{-3}]';
%c.Limits=[0 max(M.N(:))];
%caxis(ax1,[0 max(M.N(:))]);
view(ax1,2)
%set(ax1,'colorscale','log')
ax2=subplot(2,2,2,'Parent',fig);
surface(ax2,M.zgrid,M.rgrid,M.fluidUR(:,:,fieldstep)./M.vlight,'edgecolor','none');
%plot(ax2,M.zgrid,data.pot(:,5))
xlabel(ax2,'z [m]')
ylabel(ax2,'r [m]')
colormap(ax2,'jet')
c = colorbar(ax2);
%c.Limits=[min(M.fluidUR(:,:,:)) max(M.fluidUR(:,:,:))];
c.Label.String= 'U_r [m/s]';
title(ax2,'Radial velocity')
grid(ax2, 'on')
%caxis(ax2,[min(M.fluidUR(:)) max(M.fluidUR(:))])
view(ax2,2)
ax3=subplot(2,2,3,'Parent',fig);
surface(ax3,M.zgrid,M.rgrid,abs(M.fluidUTHET(:,:,fieldstep).*Rinv'),'edgecolor','none')
%mean(M.fluidUTHET(:,:,end).*Rinv')
xlabel(ax3,'z [m]')
ylabel(ax3,'r [m]')
colormap(ax3,'jet')
c = colorbar(ax3);
c.Label.String= '\omega_\theta [1/s]';
%c.Limits=[min(M.fluidUTHET(:)) max(M.fluidUTHET(:))];
%caxis(ax3,[min(M.fluidUTHET(:)) max(M.fluidUTHET(:))])
title(ax3,'Azimuthal velocity')
set(ax3,'colorscale','log')
grid(ax3, 'on')
view(ax3,2)
ax4=subplot(2,2,4,'Parent',fig);
surface(ax4,M.zgrid,M.rgrid,M.fluidUZ(:,:,fieldstep)./M.vlight,'edgecolor','none')
xlabel(ax4,'z [m]')
ylabel(ax4,'r [m]')
colormap(ax4,'jet')
c = colorbar(ax4);
%c.Limits=[min(M.fluidUZ(:)) max(M.fluidUZ(:))];
c.Label.String= 'U_z [m/s]';
title(ax4,'Axial velocity')
grid(ax4, 'on')
%caxis(ax4,[min(M.fluidUZ(:)) max(M.fluidUZ(:))])
view(ax4,2)
linkaxes([ax1 ax2 ax3 ax4],'xy')
end
function plotGridButtonPushed(btn,ax)
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
f=figure();
PlotEspic2dgriddata(f,M,sld.Value);
f.PaperOrientation='landscape';
[~, name, ~] = fileparts(M.file);
print(f,sprintf('%sfluid%d',name,sld.Value),'-dpdf','-fillpage')
end
function updatefigdata(control, event, Othercontrol, fig)
if strcmp(event.EventName,'ValueChanged')
fieldstep=floor(control.Value);
control.Value=fieldstep;
else
fieldstep=floor(event.Value);
end
Othercontrol.Value=fieldstep;
sgtitle(fig,sprintf('t=%0.5e s',double((fieldstep-1)*M.it1)*M.dt))
%% update Position histogram
ax1=fig.Children(end);
ax1.Children(1).CData=M.N(:,:,fieldstep);
ax1.Children(1).ZData=M.N(:,:,fieldstep);
view(ax1,2)
%% update Radial velocity
fig.Children(end-2).Children(1).CData=M.fluidUR(:,:,fieldstep)./M.vlight;
fig.Children(end-2).Children(1).ZData=M.fluidUR(:,:,fieldstep)./M.vlight;
%% update Azimuthal velocity
fig.Children(end-4).Children(1).CData=abs(M.fluidUTHET(:,:,fieldstep).*Rinv');
fig.Children(end-4).Children(1).ZData=abs(M.fluidUTHET(:,:,fieldstep).*Rinv');
%% update Axial velocity
fig.Children(end-6).Children(1).CData=M.fluidUZ(:,:,fieldstep)./M.vlight;
fig.Children(end-6).Children(1).ZData=M.fluidUZ(:,:,fieldstep)./M.vlight;
drawnow limitrate
end
end
c4science · Help