正在加载视频...

视频加载失败

返回首页

skye crosshair| #VALORANT code: 0;s;1;P;c;8;u;7FFF0AFF;h;0;d;1;b;1;m;1;0t;6;0l;1;0o;1;0a;1;0f;0;1o;2;1a;0.61;1m;0;1e;0.34

anime_daily's profile picture

Anime Daily

171,992 subscribers

192,132 次观看 • 1 年前 •via X (Twitter)

游戏#VALORANT
LIVE

Anya Rossi• Live Now

Private livecam show

0 条评论

暂无评论

原始帖子的评论将显示在这里

相关视频

Sakura crosshair| #VALORANT code: 0;p;0;P;c;8;u;FF8D89FF;o;0.157;d;1;b;1;a;0.484;m;1;0t;4;0l;2;0o;2;0a;1;0e;0.39;1o;3;1a;1;1m;0;1e;0.177
0:11

Sakura crosshair| #VALORANT code: 0;p;0;P;c;8;u;FF8D89FF;o;0.157;d;1;b;1;a;0.484;m;1;0t;4;0l;2;0o;2;0a;1;0e;0.39;1o;3;1a;1;1m;0;1e;0.177

Anime Daily

983,781 次观看 • 1 年前

Piplup crosshair | #VALORANT code: 0;s;1;P;c;4;o;1;d;1;z;6;a;0;f;0;m;1;0t;10;0l;10;0v;0;0g;1;0o;10;0a;0;0f;0;1t;6;1l;0;1v;3;1g;1;1o;0;1a;1;1s;0.053;1e;0.206
0:15

Piplup crosshair | #VALORANT code: 0;s;1;P;c;4;o;1;d;1;z;6;a;0;f;0;m;1;0t;10;0l;10;0v;0;0g;1;0o;10;0a;0;0f;0;1t;6;1l;0;1v;3;1g;1;1o;0;1a;1;1s;0.053;1e;0.206

Anime Daily

268,702 次观看 • 1 年前

このクロスヘアバフってるから使ってみて!! 0;c;1;s;1;P;u;0000FFFF;o;0.2;f;0;m;1;0t;1;0l;3;0v;2;0g;1;0o;2;0a;1;0f;0;1b;0;S;d;0 #VALORANT #おすすめのクロスヘア教えて
0:10

このクロスヘアバフってるから使ってみて!! 0;c;1;s;1;P;u;0000FFFF;o;0.2;f;0;m;1;0t;1;0l;3;0v;2;0g;1;0o;2;0a;1;0f;0;1b;0;S;d;0 #VALORANT #おすすめのクロスヘア教えて

ぽるぴ

368,999 次观看 • 3 年前

6MANS GAUNTLET COMPLETE 🏆 3-1 Rank H 🏆 3-0 Rank E 🏆 3-0 Rank D 🏆 3-1 Rank C 🏆 3-1 Rank B 🏆 3-1 Rank B+ ☠️ 0-3 Rank A 🏆 3-1 Rank X 🏆 3-2 Rank S
0:10

6MANS GAUNTLET COMPLETE 🏆 3-1 Rank H 🏆 3-0 Rank E 🏆 3-0 Rank D 🏆 3-1 Rank C 🏆 3-1 Rank B 🏆 3-1 Rank B+ ☠️ 0-3 Rank A 🏆 3-1 Rank X 🏆 3-2 Rank S

GGM1 Coconut

85,633 次观看 • 6 个月前

Hail, hail Mother Mahalakshmi! H a i l M o t h e r L a k s h m i, b y w h o s e g r a c e h a p p i n e s s, p r o s p e r i t y a n d u n b r o k e n g o o d f o r t u n e s h o w e r u p o n t h e h o m e a n d c o u r t y a r d O Mother Lakshmi,seated on🙌 Grok
0:15

Hail, hail Mother Mahalakshmi! H a i l M o t h e r L a k s h m i, b y w h o s e g r a c e h a p p i n e s s, p r o s p e r i t y a n d u n b r o k e n g o o d f o r t u n e s h o w e r u p o n t h e h o m e a n d c o u r t y a r d O Mother Lakshmi,seated on🙌 Grok

🌺Sai🌺

52,635 次观看 • 8 天前

Hail, hail Mother Mahalakshmi! H a i l M o t h e r L a k s h m i, b y w h o s e g r a c e h a p p i n e s s, p r o s p e r i t y a n d u n b r o k e n g o o d f o r t u n e s h o w e r u p o n t h e h o m e a n d c o u r t y a r d O Grok
0:10

Hail, hail Mother Mahalakshmi! H a i l M o t h e r L a k s h m i, b y w h o s e g r a c e h a p p i n e s s, p r o s p e r i t y a n d u n b r o k e n g o o d f o r t u n e s h o w e r u p o n t h e h o m e a n d c o u r t y a r d O Grok

🌺Sai🌺

24,789 次观看 • 22 天前

░░▒▓█▓▒░▒█▓▒░▒ 2 0 2 6. 9. 6 1 2 : 0 0 C O M I N G S O O N ░▒▓█▓▒░▒▓▒░▒█▓
0:13

░░▒▓█▓▒░▒█▓▒░▒ 2 0 2 6. 9. 6 1 2 : 0 0 C O M I N G S O O N ░▒▓█▓▒░▒▓▒░▒█▓

アイドルマスター公式

752,074 次观看 • 13 天前

for(float i,z,d,f;i++<1e2;o+=vec4(4,6,8.+z,0)/f-min(dFdx(z)*r.y+z,0.)/exp(d*d/.1)){vec3 p=z*(FC.rgb*2.-r.xyy)/r.y,c=p;p.z+=8.;c.z*=3.;for(f=1.;f++<9.;c+=sin(c.yzx*f+z+t*.5)/f);z+=min(f=.1+abs(.2*c.y+abs(p.y+.8)),d=max(length(p)-3.,.9-length(p-vec3(-1,1,3))))/7.;}o=tanh(o/2e3);
0:12

for(float i,z,d,f;i++<1e2;o+=vec4(4,6,8.+z,0)/f-min(dFdx(z)*r.y+z,0.)/exp(d*d/.1)){vec3 p=z*(FC.rgb*2.-r.xyy)/r.y,c=p;p.z+=8.;c.z*=3.;for(f=1.;f++<9.;c+=sin(c.yzx*f+z+t*.5)/f);z+=min(f=.1+abs(.2*c.y+abs(p.y+.8)),d=max(length(p)-3.,.9-length(p-vec3(-1,1,3))))/7.;}o=tanh(o/2e3);

Xor

166,930 次观看 • 1 年前

A R o s s C h e e k Goal Puts The Knappers 1-0 Up
0:16

A R o s s C h e e k Goal Puts The Knappers 1-0 Up

Knaphill FC

13,764 次观看 • 2 年前

E K D A N T B H A L S U R A T,M U S H A K D A M A N R A T🚩🙏 JAY SHRI GANESHAY NAMAH 🙏 J A Y S H R I G A N E S H A Y N A M A H 🙏 MANGAL MURTI MORAYA ✨ M A N G A L M U R T I M O R A Y A ✨ GANPATI BAPPA MORAYA ❤️ G A N P A T I B A P P A M O R A Y A ❤️ C H A N D R A S U R Y A K O C H H O D I D E E P A K J A L A T 🙌 Grok
0:15

E K D A N T B H A L S U R A T,M U S H A K D A M A N R A T🚩🙏 JAY SHRI GANESHAY NAMAH 🙏 J A Y S H R I G A N E S H A Y N A M A H 🙏 MANGAL MURTI MORAYA ✨ M A N G A L M U R T I M O R A Y A ✨ GANPATI BAPPA MORAYA ❤️ G A N P A T I B A P P A M O R A Y A ❤️ C H A N D R A S U R Y A K O C H H O D I D E E P A K J A L A T 🙌 Grok

राधा

56,393 次观看 • 1 个月前

Orchard vec3 p,v=normalize(FC.rgb*2.-r.xyx),c=v/v.y;c.z+=.5*t;for(float z,i,b,g,m;i++<5e1;z+=.8*max(b=length((p.y-m)/1e2/(abs(sin(c.xz/.1))-.05/v.y)),min(4.-m,g=length(sin(p.xz)+1.-.1*(1.+sin(p.y-p.zx*.5))*m))-b),o.rgb+=(.7-v)/(g+b))p=z*v+1.,p.z-=t,m=abs(++p.y);o=tanh(o/5e2);
0:12

Orchard vec3 p,v=normalize(FC.rgb*2.-r.xyx),c=v/v.y;c.z+=.5*t;for(float z,i,b,g,m;i++<5e1;z+=.8*max(b=length((p.y-m)/1e2/(abs(sin(c.xz/.1))-.05/v.y)),min(4.-m,g=length(sin(p.xz)+1.-.1*(1.+sin(p.y-p.zx*.5))*m))-b),o.rgb+=(.7-v)/(g+b))p=z*v+1.,p.z-=t,m=abs(++p.y);o=tanh(o/5e2);

Xor

12,216 次观看 • 7 个月前

for(float i=0.,z=0.,d=0.,s=0.;i+++z<1e2;){vec3 q=z*normalize(vec3(FC.xy*2.-r,r.y));q.z+=t*4.;q.xy*=rotate2D(q.z*.1);for(s=1.;s<32.;s/=.5)q+=(abs(cos(q.yzx*s))-.68)/s;z+=d=.005+abs((length(q.yx)-6.)*q.x)/8.;o.rgb+=hsv(.3+sin(q.x)*.2,.5,.2)*exp(-z*.1)/d;}o=tanh(o/2e2);#つぶやきGLSL
0:15

for(float i=0.,z=0.,d=0.,s=0.;i+++z<1e2;){vec3 q=z*normalize(vec3(FC.xy*2.-r,r.y));q.z+=t*4.;q.xy*=rotate2D(q.z*.1);for(s=1.;s<32.;s/=.5)q+=(abs(cos(q.yzx*s))-.68)/s;z+=d=.005+abs((length(q.yx)-6.)*q.x)/8.;o.rgb+=hsv(.3+sin(q.x)*.2,.5,.2)*exp(-z*.1)/d;}o=tanh(o/2e2);#つぶやきGLSL

Yohei Nishitsuji

41,292 次观看 • 13 天前

"Waveform" for(float i,d,z,f;i++<9e1;o+=(cos(z*.5+t+vec4(0,2,4,3))+1.3)/d/z){vec3 p=z*normalize(FC.rgb*2.-r.xyy)+1.;f=max(-p.y,0.);p.y+=f+f;for(d=1.;d<3e1;d+=d)p.y+=cos(p*d+2.*t*cos(d)+z).x/d;z+=d=(.1*f+abs(p.y-1.)/(++f*f)+max(d=p.z+3.,-d*.1))/8.;}o=tanh(o/9e2);
0:20

"Waveform" for(float i,d,z,f;i++<9e1;o+=(cos(z*.5+t+vec4(0,2,4,3))+1.3)/d/z){vec3 p=z*normalize(FC.rgb*2.-r.xyy)+1.;f=max(-p.y,0.);p.y+=f+f;for(d=1.;d<3e1;d+=d)p.y+=cos(p*d+2.*t*cos(d)+z).x/d;z+=d=(.1*f+abs(p.y-1.)/(++f*f)+max(d=p.z+3.,-d*.1))/8.;}o=tanh(o/9e2);

Xor

46,883 次观看 • 1 年前

#Dodgers Bobby Miller gets through his first rehab inning. 2 strike outs, 1 ground out, one hit batter. 23 P / 13 S / 1 H / 0 ER / 0 BB / 2 K / 98mph FB / 85 mph changeup / 80 mph curveball. Control a bit shacky, but good first inning.
0:14

#Dodgers Bobby Miller gets through his first rehab inning. 2 strike outs, 1 ground out, one hit batter. 23 P / 13 S / 1 H / 0 ER / 0 BB / 2 K / 98mph FB / 85 mph changeup / 80 mph curveball. Control a bit shacky, but good first inning.

Dodgers_After_Duty

56,967 次观看 • 1 个月前

Following many two-dimensional turbulence animations, it is time to share computer code! Implementation of a vorticity-streamfunction method with a pseudo-spectral discretization + third-order Runge-Kutta for time integration (our CFD class semester project). The MATLAB code: M = 256; % number of points N = M; Lx = 2*pi; Ly = 2*pi; nu = 5e-4; % kinematic viscosity Sc = 0.7; % Schmidt number beta = 0; % meridional gradient of Coriolis parameter ar = 0.02; %random number amplitude b = 1; % mean scalar gradient CFLmax = 0.8; tend = 200; % end time x=linspace(0,Lx,M+1); x(end)=[]; dx = Lx/M; kx=[0:M/2 -M/2+1:-1]*2*pi/Lx; y=linspace(0,Ly,N+1); y(end)=[]; ky=[0:N/2 -N/2+1:-1]*2*pi/Ly; dy = Ly/N; time = 0; index_kmax = ceil(M/3); kmax = kx(index_kmax); filter = ones(M,N); filter(index_kmax+1:2*index_kmax+3,index_kmax+1:2*index_kmax+3)=0; rng(64); [u, v, omega, psi, ddx, ddy, idel2, kk, k2]=deal(zeros(M,N)); for j=1:N ddx(:,j)=1i*kx; end for i=1:M ddy(i,:)=1i*ky; end for i=1:M for j=1:N idel2(i,j)=-kx(i)^2-ky(j)^2; end end idel2=1./idel2; idel2(1,1)=0; for i=1:M for j=1:N kk(i,j)=kx(i)^2+ky(j)^2; k2(i,j)=kx(i)^2+ky(j)^2; if kk(i,j) >= 6^2 && kk(i,j) <= 7^2 % forcing kk(i,j) = -kk(i,j); end if kk(i,j) <= 2^2 kk(i,j) = 8*kk(i,j); % large-scale dissipation end end end for i=1:M for j=1:N u(i,j) = cos(2*x(i))*sin(2*y(j))+ar*rand; v(i,j) = -sin(2*x(i))*cos(2*y(j))+ar*rand; end end uhat = fft2(u); vhat = fft2(v); omegahat = ddx.*vhat - ddy.*uhat; % make vorticity phi = rand(size(u)); phihat = fft2(phi); ncid = netcdf.create(' 'CLOBBER'); dimid_x = netcdf.defDim(ncid, 'x', M); dimid_y = netcdf.defDim(ncid, 'y', N); dimid_time = netcdf.defDim(ncid, 'time', netcdf.getConstant('NC_UNLIMITED')); varid_x = netcdf.defVar(ncid, 'x', 'NC_FLOAT', [dimid_x]); varid_y = netcdf.defVar(ncid, 'y', 'NC_FLOAT', [dimid_y]); varid_u = netcdf.defVar(ncid, 'u', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_v = netcdf.defVar(ncid, 'v', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_omega = netcdf.defVar(ncid, 'vorticity', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_phi = netcdf.defVar(ncid, 'scalar', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_dissipation = netcdf.defVar(ncid, 'dissipation', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_time = netcdf.defVar(ncid, 'time', 'NC_FLOAT', [dimid_time]); netcdf.endDef(ncid); netcdf.putVar(ncid, varid_time, 0, 1, time); netcdf.putVar(ncid, varid_x, 0, M, x); netcdf.putVar(ncid, varid_y, 0, N, y); netcdf.putVar(ncid, varid_phi, [0, 0, 0], [M, N, 1], phi); netcdf.putVar(ncid, varid_u, [0, 0, 0], [M, N, 1], u); netcdf.putVar(ncid, varid_v, [0, 0, 0], [M, N, 1], v); netcdf.putVar(ncid, varid_omega, [0, 0, 0], [M, N, 1], omega); netcdf.putVar(ncid, varid_dissipation, [0, 0, 0], [M, N, 1], 0*omega); dt = 0.5*min([dx dy]); nstep = 1; while time < tend psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); facto = exp(-nu*8/15*dt*kk); factp = exp(-nu/Sc*8/15*dt*k2); r0o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r0p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = facto.*(omegahat + dt*8/15*r0o); % update omega phihat = factp.*(phihat + dt*8/15*r0p); % update phi %%%% Substep 2 psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); r1o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r1p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = omegahat + dt*(-17/60*facto.*r0o + 5/12*r1o); phihat = phihat + dt*(-17/60*factp.*r0p + 5/12*r1p); facto = exp(-nu*(-17/60+5/12)*dt*kk); factp = exp(-nu/Sc*(-17/60+5/12)*dt*k2); omegahat = omegahat.*facto; phihat = phihat.*factp; %%%% Substep 3 psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; % max(max(abs(real(ifft2(1i*ddx.*uhat+1i*ddy.*vhat))))) % divergence u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); r2o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r2p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = omegahat + dt*(-5/12*facto.*r1o + 3/4*r2o); phihat = phihat + dt*(-5/12*factp.*r1p + 3/4*r2p); facto = exp(-nu*(-5/12+3/4)*dt*kk); factp = exp(-nu/Sc*(-5/12+3/4)*dt*kk); omegahat = omegahat.*facto; phihat = phihat.*factp; phihat = filter.*phihat; omegahat = filter.*omegahat; time = time + dt; nstep = nstep + 1; CFL = max(max(abs(u)))/dx*dt+max(max(abs(v)))/dy*dt; if mod(nstep,20)==0 phi = real(ifft2(phihat)); omega = real(ifft2(omegahat)); dissipation = 2*nu*(real(ifft2(ddx.*uhat)).^2 + real(ifft2(ddy.*uhat)).^2 + real(ifft2(ddx.*vhat)).^2 + real(ifft2(ddy.*vhat)).^2); eta = (nu^3/mean(dissipation,'all'))^0.25; subplot(221); pcolor(x,y,omega'); title('Vorticity'); shading flat; axis equal tight; colorbar; drawnow subplot(222); pcolor(x,y,phi'); title('Scalar'); shading flat; axis equal tight; colorbar; drawnow subplot(223); pcolor(x,y,dissipation'); title('Dissipation'); shading flat; axis equal tight; colorbar; drawnow subplot(224); pcolor(x,y,u); title('u-velocity'); shading flat; axis equal tight; colorbar; drawnow fprintf(1,'step = %d time = %g dt = %g CFL = %g kmax*eta = %g %g\n', nstep, time, dt, CFL, eta.*kmax, eta.*kmax/sqrt(Sc)); [~, dim_time_len] = netcdf.inqDim(ncid,dimid_time); netcdf.putVar(ncid, varid_time, [dim_time_len], [1], time); netcdf.putVar(ncid, varid_phi, [0, 0, dim_time_len], [M, N, 1], phi); netcdf.putVar(ncid, varid_u, [0, 0, dim_time_len], [M, N, 1], u); netcdf.putVar(ncid, varid_v, [0, 0, dim_time_len], [M, N, 1], v); netcdf.putVar(ncid, varid_omega, [0, 0, dim_time_len], [M, N, 1], omega); netcdf.putVar(ncid, varid_dissipation, [0, 0, dim_time_len], [M, N, 1], dissipation); netcdf.sync(ncid); end dt = CFLmax/CFL*dt; %0.0005; end netcdf.close(ncid);
0:10

Following many two-dimensional turbulence animations, it is time to share computer code! Implementation of a vorticity-streamfunction method with a pseudo-spectral discretization + third-order Runge-Kutta for time integration (our CFD class semester project). The MATLAB code: M = 256; % number of points N = M; Lx = 2*pi; Ly = 2*pi; nu = 5e-4; % kinematic viscosity Sc = 0.7; % Schmidt number beta = 0; % meridional gradient of Coriolis parameter ar = 0.02; %random number amplitude b = 1; % mean scalar gradient CFLmax = 0.8; tend = 200; % end time x=linspace(0,Lx,M+1); x(end)=[]; dx = Lx/M; kx=[0:M/2 -M/2+1:-1]*2*pi/Lx; y=linspace(0,Ly,N+1); y(end)=[]; ky=[0:N/2 -N/2+1:-1]*2*pi/Ly; dy = Ly/N; time = 0; index_kmax = ceil(M/3); kmax = kx(index_kmax); filter = ones(M,N); filter(index_kmax+1:2*index_kmax+3,index_kmax+1:2*index_kmax+3)=0; rng(64); [u, v, omega, psi, ddx, ddy, idel2, kk, k2]=deal(zeros(M,N)); for j=1:N ddx(:,j)=1i*kx; end for i=1:M ddy(i,:)=1i*ky; end for i=1:M for j=1:N idel2(i,j)=-kx(i)^2-ky(j)^2; end end idel2=1./idel2; idel2(1,1)=0; for i=1:M for j=1:N kk(i,j)=kx(i)^2+ky(j)^2; k2(i,j)=kx(i)^2+ky(j)^2; if kk(i,j) >= 6^2 && kk(i,j) <= 7^2 % forcing kk(i,j) = -kk(i,j); end if kk(i,j) <= 2^2 kk(i,j) = 8*kk(i,j); % large-scale dissipation end end end for i=1:M for j=1:N u(i,j) = cos(2*x(i))*sin(2*y(j))+ar*rand; v(i,j) = -sin(2*x(i))*cos(2*y(j))+ar*rand; end end uhat = fft2(u); vhat = fft2(v); omegahat = ddx.*vhat - ddy.*uhat; % make vorticity phi = rand(size(u)); phihat = fft2(phi); ncid = netcdf.create(' 'CLOBBER'); dimid_x = netcdf.defDim(ncid, 'x', M); dimid_y = netcdf.defDim(ncid, 'y', N); dimid_time = netcdf.defDim(ncid, 'time', netcdf.getConstant('NC_UNLIMITED')); varid_x = netcdf.defVar(ncid, 'x', 'NC_FLOAT', [dimid_x]); varid_y = netcdf.defVar(ncid, 'y', 'NC_FLOAT', [dimid_y]); varid_u = netcdf.defVar(ncid, 'u', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_v = netcdf.defVar(ncid, 'v', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_omega = netcdf.defVar(ncid, 'vorticity', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_phi = netcdf.defVar(ncid, 'scalar', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_dissipation = netcdf.defVar(ncid, 'dissipation', 'NC_FLOAT', [dimid_x, dimid_y, dimid_time]); varid_time = netcdf.defVar(ncid, 'time', 'NC_FLOAT', [dimid_time]); netcdf.endDef(ncid); netcdf.putVar(ncid, varid_time, 0, 1, time); netcdf.putVar(ncid, varid_x, 0, M, x); netcdf.putVar(ncid, varid_y, 0, N, y); netcdf.putVar(ncid, varid_phi, [0, 0, 0], [M, N, 1], phi); netcdf.putVar(ncid, varid_u, [0, 0, 0], [M, N, 1], u); netcdf.putVar(ncid, varid_v, [0, 0, 0], [M, N, 1], v); netcdf.putVar(ncid, varid_omega, [0, 0, 0], [M, N, 1], omega); netcdf.putVar(ncid, varid_dissipation, [0, 0, 0], [M, N, 1], 0*omega); dt = 0.5*min([dx dy]); nstep = 1; while time < tend psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); facto = exp(-nu*8/15*dt*kk); factp = exp(-nu/Sc*8/15*dt*k2); r0o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r0p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = facto.*(omegahat + dt*8/15*r0o); % update omega phihat = factp.*(phihat + dt*8/15*r0p); % update phi %%%% Substep 2 psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); r1o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r1p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = omegahat + dt*(-17/60*facto.*r0o + 5/12*r1o); phihat = phihat + dt*(-17/60*factp.*r0p + 5/12*r1p); facto = exp(-nu*(-17/60+5/12)*dt*kk); factp = exp(-nu/Sc*(-17/60+5/12)*dt*k2); omegahat = omegahat.*facto; phihat = phihat.*factp; %%%% Substep 3 psihat = -idel2.*omegahat; uhat = ddy.*psihat; vhat = -ddx.*psihat; % max(max(abs(real(ifft2(1i*ddx.*uhat+1i*ddy.*vhat))))) % divergence u = real(ifft2(uhat)); v = real(ifft2(vhat)); omegadx = real(ifft2(ddx.*omegahat)); omegady = real(ifft2(ddy.*omegahat)); r2o = -fft2(u.*omegadx+v.*omegady)+beta*vhat; r2p = -fft2(u.*real(ifft2(ddx.*phihat))+v.*real(ifft2(ddy.*phihat)))+b*vhat; omegahat = omegahat + dt*(-5/12*facto.*r1o + 3/4*r2o); phihat = phihat + dt*(-5/12*factp.*r1p + 3/4*r2p); facto = exp(-nu*(-5/12+3/4)*dt*kk); factp = exp(-nu/Sc*(-5/12+3/4)*dt*kk); omegahat = omegahat.*facto; phihat = phihat.*factp; phihat = filter.*phihat; omegahat = filter.*omegahat; time = time + dt; nstep = nstep + 1; CFL = max(max(abs(u)))/dx*dt+max(max(abs(v)))/dy*dt; if mod(nstep,20)==0 phi = real(ifft2(phihat)); omega = real(ifft2(omegahat)); dissipation = 2*nu*(real(ifft2(ddx.*uhat)).^2 + real(ifft2(ddy.*uhat)).^2 + real(ifft2(ddx.*vhat)).^2 + real(ifft2(ddy.*vhat)).^2); eta = (nu^3/mean(dissipation,'all'))^0.25; subplot(221); pcolor(x,y,omega'); title('Vorticity'); shading flat; axis equal tight; colorbar; drawnow subplot(222); pcolor(x,y,phi'); title('Scalar'); shading flat; axis equal tight; colorbar; drawnow subplot(223); pcolor(x,y,dissipation'); title('Dissipation'); shading flat; axis equal tight; colorbar; drawnow subplot(224); pcolor(x,y,u); title('u-velocity'); shading flat; axis equal tight; colorbar; drawnow fprintf(1,'step = %d time = %g dt = %g CFL = %g kmax*eta = %g %g\n', nstep, time, dt, CFL, eta.*kmax, eta.*kmax/sqrt(Sc)); [~, dim_time_len] = netcdf.inqDim(ncid,dimid_time); netcdf.putVar(ncid, varid_time, [dim_time_len], [1], time); netcdf.putVar(ncid, varid_phi, [0, 0, dim_time_len], [M, N, 1], phi); netcdf.putVar(ncid, varid_u, [0, 0, dim_time_len], [M, N, 1], u); netcdf.putVar(ncid, varid_v, [0, 0, dim_time_len], [M, N, 1], v); netcdf.putVar(ncid, varid_omega, [0, 0, dim_time_len], [M, N, 1], omega); netcdf.putVar(ncid, varid_dissipation, [0, 0, dim_time_len], [M, N, 1], dissipation); netcdf.sync(ncid); end dt = CFLmax/CFL*dt; %0.0005; end netcdf.close(ncid);

Computational Fluid Dynamics group

226,726 次观看 • 3 年前

a=(x,y,k=x/12-12.5,d=abs(cos(k/2)+sin(y/2)*cos(k/5)/2+1))=&gt;point((q=x/4+99+(d+0)*k*(.5+cos(d*2-t+y/8)))*.8*cos(c=y/69-t/8+d/6)+200,q/2*cos(c+d/3+y/69)+200) t=0,draw=$=&gt;{t||createCanvas(w=400,w);background(0).stroke(w,46);for(t+=PI/30,i=4e4;i--;)a(i%300,i/300)}#つぶやきProcessing
0:16

a=(x,y,k=x/12-12.5,d=abs(cos(k/2)+sin(y/2)*cos(k/5)/2+1))=&gt;point((q=x/4+99+(d+0)*k*(.5+cos(d*2-t+y/8)))*.8*cos(c=y/69-t/8+d/6)+200,q/2*cos(c+d/3+y/69)+200) t=0,draw=$=&gt;{t||createCanvas(w=400,w);background(0).stroke(w,46);for(t+=PI/30,i=4e4;i--;)a(i%300,i/300)}#つぶやきProcessing

ア

18,708 次观看 • 1 年前

19:51 | 🔴🔥#URGENTE #Providencia Se declara INCENDIO EN HOSPITAL CALVO MACKENNA. Se despacha B-1, M-5, Q-7, RX-14, H-4, S-1, B-21, B-22, BT-18 en apoyo a B-13, B-3, B-14, Q-8, MX-13 (Ricardo Hernández) Red de Emergencia ONG Carlos González Baeza.
0:15

19:51 | 🔴🔥#URGENTE #Providencia Se declara INCENDIO EN HOSPITAL CALVO MACKENNA. Se despacha B-1, M-5, Q-7, RX-14, H-4, S-1, B-21, B-22, BT-18 en apoyo a B-13, B-3, B-14, Q-8, MX-13 (Ricardo Hernández) Red de Emergencia ONG Carlos González Baeza.

Alejandro Alegria 🌴

59,184 次观看 • 2 年前