Can we get the mass balance from the model output U,V,W?

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
cathyyangfeng
Posts: 24
Joined: Fri Aug 29, 2008 4:13 am
Location: Virginia Institute of Marine Sciences

Can we get the mass balance from the model output U,V,W?

#1 Unread post by cathyyangfeng »

I recently confused about mass conservation problem in ROMS.

I use the model output U,V,W on each grid, and calculated Pu/Px + Pv/Py + Pw/Pz problems on designed model grid. I could not get a completely zero on each grid. Can we get an idealized mass conservation from this way in calculation?

linzhenhua
Posts: 64
Joined: Mon Oct 17, 2005 2:02 am
Location: Institute of Oceanology,Chinese Academy of Sciences

Re: Can we get the mass balance from the model output U,V,W?

#2 Unread post by linzhenhua »

cathyyangfeng wrote:I recently confused about mass conservation problem in ROMS.

I use the model output U,V,W on each grid, and calculated Pu/Px + Pv/Py + Pw/Pz problems on designed model grid. I could not get a completely zero on each grid. Can we get an idealized mass conservation from this way in calculation?

I think the total volume should be conserved because the vertical velocity is calculated according to volume conservation rule(Pu/Px + Pv/Py + Pw/Pz=0). The non-zero value is possibly because you use different discretization scheme than that used in ROMS.

cathyyangfeng
Posts: 24
Joined: Fri Aug 29, 2008 4:13 am
Location: Virginia Institute of Marine Sciences

Re: Can we get the mass balance from the model output U,V,W?

#3 Unread post by cathyyangfeng »

Yes, I agree with what you said. I did not use the central finite difference as the model did.

I believe the model is volume conserved. But if we select a small region and use Int(velocity x dArea). why usually we could not get a totally zero?

Int(velocity x dArea) = Int(u*dydz) + Int(v*dxdz) + Int(w*dydz) at each boundary place.

If we interested in the mass flux goes in from different boundaries, what's going to be the right way for calculation?

linzhenhua
Posts: 64
Joined: Mon Oct 17, 2005 2:02 am
Location: Institute of Oceanology,Chinese Academy of Sciences

Re: Can we get the mass balance from the model output U,V,W?

#4 Unread post by linzhenhua »

The volume for a specific region is not conserved because there is a free surface.

Post Reply