Overleg:Shallow water equations

Uit EurosWiki
Versie door Wgodlieb (overleg | bijdragen) op 22 dec 2008 om 18:10
(wijz) ← Oudere versie | Huidige versie (wijz) | Nieuwere versie → (wijz)
Naar navigatie springen Naar zoeken springen

Ik probeer de vergelijkingen er ook nog even netje bij te maken, maar dat is nog best lastig. Willem 22 Dec 2008 17:10 (CET)

WX,YtDXDYX=ux,YWx,Y*DYxux+1,YWx+1,YDYx+1+vX,y*WX,y*DXvX,y+1*WX,y+1*DX


ux,YtDXDYxWxY=impulsnormaal+impulstangentieel+zwaartekracht+kromming+coriolis+getijkracht+frictie

impulsnormaal=Wx1,Yux1,y2DYx1Wx+1,Yux+1,y2DYx+12

impulstangentieel=Wx,yvx,y2DXWx+1,yux+1,y2DX

dwdt[x][y]=u[x][y]*uW[x][y]*dyo[x]u[x+1][y]*uW[x+1][y]*dyo[x+1]+(v[x][y]*vW[x][y]v[x][y+1]*vW[x][y+1])*DX;W[x][y]+=dwdt[x][y]*DT/DX/dym[x];dudt[x][y]=(uW[x1][y]*u[x1][y]*u[x1][y]*dyo[x1]uW[x+1][y]*u[x+1][y]*u[x+1][y]*dyo[x+1])/2.0+//impulstransport.5*g*(H[x1][y]H[x][y])*dyo[x]*uW[x][y]+//potentieleenergie(LOW[x][y]*(v[x][y]+v[x1][y])*(u[x][y]+u[x][y1])ROW[x][y]*(v[x][y+1]+v[x1][y+1])*(u[x][y+1]+u[x][y]))/4.0*DX+//tangentieelimpulstransportu[x][y]/R*dyo[x]*DX*vopu[x][y]*sino[x]+//curvature2*omega*vopu[x][y]*sino[x]*uW[x][y]*DX*dyo[x]+//corioliskrachtuGM[x][y]*uW[x][y]*DX*dyo[x]+//maangetijkrachtu[x][y]*0.01*DX*dyo[x];//frictiedvdt[x][y]=(vW[x][y1]*v[x][y1]*v[x][y1]vW[x][y+1]*v[x][y+1]*v[x][y+1])/2.0*DX+.5*g*(H[x][y1]H[x][y])*vW[x][y]*DX+(LOW[x][y]*(u[x][y]+u[x][y1])*(v[x][y]+v[x1][y])*dyo[x]LBW[x][y]*(u[x+1][y]+u[x+1][y1])*(v[x+1][y]+v[x][y])*dyo[x+1])/4.0+uopv[x][y]/R*dym[x]*DX*uopv[x][y]*sinm[x]+2*omega*uopv[x][y]*sinm[x]*vW[x][y]*DX*dym[x]+vGM[x][y]*vW[x][y]*DX*dym[x]+v[x][y]*0.01*DX*dym[x];dudt[x][y]/=uW[x][y]*DX*dyo[x];dvdt[x][y]/=vW[x][y]*DX*dym[x];