파일:Airflow-Obstructed-Duct.png
From Wikipedia, the free encyclopedia
미리 보기 크기: 800 × 571 픽셀 다른 해상도: 320 × 229 픽셀 | 640 × 457 픽셀 | 1,024 × 731 픽셀 | 1,270 × 907 픽셀
원본 파일 (1,270 × 907 픽셀, 파일 크기: 85 KB, MIME 종류: image/png)
파일 설명
이 그림은 벡터 그래픽 버전(SVG)이 있습니다. 래스터(비트맵) 그림 대신 벡터 그래픽 그림을 사용하는 것이 좋습니다.
File:Airflow-Obstructed-Duct.png → File:N S Laminar.svg
|
설명Airflow-Obstructed-Duct.png |
A simulation using the navier-stokes differential equations of the aiflow into a duct at 0.003 m/s (laminar flow). The duct has a small obstruction in the centre that is parallel with the duct walls. The observed spike is mainly due to numerical limitations. This script, which i originally wrote for scilab, but ported to matlab (porting is really really easy, mainly convert comments % -> // and change the fprintf and input statements) Matlab was used to generate the image. %Matlab script to solve a laminar flow %in a duct problem %Constants inVel = 0.003; % Inlet Velocity (m/s) fluidVisc = 1e-5; % Fluid's Viscoisity (Pa.s) fluidDen = 1.3; %Fluid's Density (kg/m^3) MAX_RESID = 1e-5; %uhh. residual units, yeah... deltaTime = 1.5; %seconds? %Kinematic Viscosity fluidKinVisc = fluidVisc/fluidDen; %Problem dimensions ductLen=5; %m ductWidth=1; %m %grid resolution gridPerLen = 50; % m^(-1) gridDelta = 1/gridPerLen; XVec = 0:gridDelta:ductLen-gridDelta; YVec = 0:gridDelta:ductWidth-gridDelta; %Solution grid counts gridXSize = ductLen*gridPerLen; gridYSize = ductWidth*gridPerLen; %Lay grid out with Y increasing down rows %x decreasing down cols %so subscripting becomes (y,x) (sorry) velX= zeros(gridYSize,gridXSize); velY= zeros(gridYSize,gridXSize); newVelX= zeros(gridYSize,gridXSize); newVelY= zeros(gridYSize,gridXSize); %Set initial condition for i =2:gridXSize-1 for j =2:gridYSize-1 velY(j,i)=0; velX(j,i)=inVel; end end %Set boundary condition on inlet for i=2:gridYSize-1 velX(i,1)=inVel; end disp(velY(2:gridYSize-1,1)); %Arbitrarily set residual to prevent %early loop termination resid=1+MAX_RESID; simTime=0; while(deltaTime) count=0; while(resid > MAX_RESID && count < 1e2) count = count +1; for i=2:gridXSize-1 for j=2:gridYSize-1 newVelX(j,i) = velX(j,i) + deltaTime*( fluidKinVisc / (gridDelta.^2) * ... (velX(j,i+1) + velX(j+1,i) - 4*velX(j,i) + velX(j-1,i) + ... velX(j,i-1)) - 1/(2*gridDelta) *( velX(j,i) *(velX(j,i+1) - ... velX(j,i-1)) + velY(j,i)*( velX(j+1,i) - velX(j,i+1)))); newVelY(j,i) = velY(j,i) + deltaTime*( fluidKinVisc / (gridDelta.^2) * ... (velY(j,i+1) + velY(j+1,i) - 4*velY(j,i) + velY(j-1,i) + ... velY(j,i-1)) - 1/(2*gridDelta) *( velY(j,i) *(velY(j,i+1) - ... velY(j,i-1)) + velY(j,i)*( velY(j+1,i) - velY(j,i+1)))); end end %Copy the data into the front for i=2:gridXSize - 1 for j = 2:gridYSize-1 velX(j,i) = newVelX(j,i); velY(j,i) = newVelY(j,i); end end %Set free boundary condition on inlet (dv_x/dx) = dv_y/dx = 0 for i=1:gridYSize velX(i,gridXSize)=velX(i,gridXSize-1); velY(i,gridXSize)=velY(i,gridXSize-1); end %y velocity generating vent for i=floor(2/6*gridXSize):floor(4/6*gridXSize) velX(floor(gridYSize/2),i) = 0; velY(floor(gridYSize/2),i-1) = 0; end %calculate residual for %conservation of mass resid=0; for i=2:gridXSize-1 for j=2:gridYSize-1 %mass continuity equation using central difference %approx to differential resid = resid + (velX(j,i+ 1)+velY(j+1,i) - ... (velX(j,i-1) + velX(j-1,i)))^2; end end resid = resid/(4*(gridDelta.^2))*1/(gridXSize*gridYSize); fprintf('Time %5.3f \t log10Resid : %5.3f\n',simTime,log10(resid)); simTime = simTime + deltaTime; end mesh(XVec,YVec,velX) deltaTime = input('\nnew delta time:'); end %Plot the results mesh(XVec,YVec,velX) |
날짜 | 2007년 2월 24일 (원본 올리기 일시) |
출처 | en.wikipedia에서 공용으로 옮겨왔습니다. |
저자 | 영어 위키백과의 User A1 |
라이선스
Public domainPublic domainfalsefalse |
이 작품은 저작자인 영어 위키백과의 User A1에 의해 퍼블릭 도메인으로 공개된 작품입니다. 이 공개 선언은 전 세계적으로 유효합니다. 만약 저작권의 포기가 법률적으로 가능하지 않은 경우, User A1은 이 작품을 법적으로 허용되는 한도 내에서 누구나 자유롭게 어떤 목적으로도 제한 없이 사용할 수 있도록 허용합니다.Public domainPublic domainfalsefalse |
기존 올리기 기록
The original description page was here. All following user names refer to en.wikipedia.
- 2007-02-24 05:45 User A1 1270×907×8 (86796 bytes) A simulation using the navier-stokes differential equations of the aiflow into a duct at 0.003 m/s (laminar flow). The duct has a small obstruction in the centre that is paralell with the duct walls. The observed spike is mainly due to numerical limitatio
이 파일에 묘사된 항목
다음을 묘사함
24 2 2007
image/png
86,796 바이트
907 화소
1,270 화소
파일 역사
날짜/시간 링크를 클릭하면 해당 시간의 파일을 볼 수 있습니다.
날짜/시간 | 섬네일 | 크기 | 사용자 | 설명 | |
---|---|---|---|---|---|
현재 | 2007년 5월 2일 (수) 00:52 | 1,270 × 907 (85 KB) | Smeira | {{Information |Description=A simulation using the navier-stokes differential equations of the aiflow into a duct at 0.003 m/s (laminar flow). The duct has a small obstruction in the centre that is paralell with the duct walls. The observed spike is mainly |
이 파일을 사용하는 문서
이 파일을 사용하고 있는 모든 위키의 문서 목록
다음 위키에서 이 파일을 사용하고 있습니다:
- anp.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ar.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ba.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- bg.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- bn.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ca.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ckb.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- cs.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- de.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- en.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- en.wikiquote.org에서 이 파일을 사용하고 있는 문서 목록
- es.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- fa.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- he.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- hif.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- hi.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- hr.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- hy.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- id.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- jv.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ko.wikiversity.org에서 이 파일을 사용하고 있는 문서 목록
- map-bms.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- ms.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- mwl.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- pt.wikipedia.org에서 이 파일을 사용하고 있는 문서 목록
- Isaac Newton
- Equação diferencial
- Equações de Navier-Stokes
- Equação diferencial linear
- Equação diferencial de Bernoulli
- Equação diferencial de d'Alembert
- Decaimento exponencial
- Equação de Laplace
- Equação diferencial parcial
- Equação de Poisson
- Equação do calor
- Lema de Grönwall
- Teorema de Picard-Lindelöf
- Método de Runge-Kutta
- Equação de Mason-Weaver
- Equação do pêndulo
- Equação de onda
- Método Multigrid
- Teoria de Sturm-Liouville
이 파일의 더 많은 사용 내역을 봅니다.