[haze_1986]

Hi all, just would like to verify certain parts of the code:
http://www.tfd.chalmers.se/~hani/kur...vistReport.pdf

From the tutorial it was mentioned that h is the mean surface height and h0 is the deviation from the mean surface height? Should it be the other way round? In the code h is used in most of the momentum equations and theoretically this should be the mean plus the deviation.

fvVectorMatrix hUEqn
(
fvm::ddt(hU)
+ fvm::div(phiv, hU)
);

On the other hand, in the mometum predictor, h+h0 is used in the code. Theoretically, we need to use the deviation from the mean inside the grad, which based on the tutorial definition should be h0.

solve(hUEqn + (F ^ hU) == -magg*h*fvc::grad(h + h0));

Is there something I am missing here?

[loic_d]

Hi all,
I faced this problem last month.
This tutorial contains a lot of errors, h is the water height and h0 the bed height.

[RoderickZ]

Originally Posted by loic_d

Hi all,
I faced this problem last month.
This tutorial contains a lot of errors, h is the water height and h0 the bed height.

With water height (h) you refer to the elevation of the free surface relative to the geoid and h_0 is the bathymetry ?
or h the height from topography to free surface?

(https://users.oden.utexas.edu/~arbog.../dawson_v2.pdf)

[randolph]

In OpenFOAM, your coordinate origin should be defined somewhere below the bottom floor. h0 is your bottom elevation (vertical distance from the bottom to your coordinate origin). h is your water column depth (from your water surface to the top of your bottom floor). And the htotal is a dummy; as defined as htotal = h0+h, which is your total water surface elevation (water surface to your origin).

Further information can be found at: A high-order arbitrarily unstructured finite-volume model of the global atmosphere: Tests solving the shallow-water equations.


Thanks,
Rdf

[TomS]

You will find an example of a shallowFoam application and details of verification and set-up here:

https://www.mts-cfd.com/ballater

Tom