Difference between revisions of "Conservation"
| (One intermediate revision by the same user not shown) | |||
| Line 1: | Line 1: | ||
Volume conservation is not enforced in SCHISM but the error is generally small. Mass conservation is enforced with upwind/TVD transport. | Volume conservation is not enforced in SCHISM but the error is generally small. Mass conservation is enforced with upwind/TVD transport. | ||
| − | It's important to use proper method in checking the volume conservation. In particular, you should make use of the flag consv_check in param.in, instead of calculating flows using the outputs yourself. This is because the velocity outputs have interpolated from side centers to nodes. Strictly speaking, checking volume conservation in SCHISM entails a calculation of semi-implicit fluxes to be consistent with the original formulation | + | It's important to use proper method in checking the volume conservation. In particular, you should make use of the flag consv_check in param.in, instead of calculating flows using the outputs yourself. This is because the velocity outputs have interpolated from side centers to nodes. Strictly speaking, checking volume conservation in SCHISM entails a calculation of semi-implicit fluxes to be consistent with the original formulation (see transport solver in the online manual). |
Another thing to be aware of is the boundary condition. For example, the imposed river flow is conserved downstream, but if you also impose the elevation b.c., this is not the case anymore. The elev. b.c. as an essential b.c. over-rules the flow b.c. (a Neumann b.c.). | Another thing to be aware of is the boundary condition. For example, the imposed river flow is conserved downstream, but if you also impose the elevation b.c., this is not the case anymore. The elev. b.c. as an essential b.c. over-rules the flow b.c. (a Neumann b.c.). | ||
| + | |||
| + | [https://ccrm.vims.edu/yinglong/SVN_large_files/Mass-conservation-in-SCHISM-and-test-cases.pdf Here] is a document that describes mass conservation with settling velocity. | ||
Latest revision as of 15:51, 18 November 2025
Volume conservation is not enforced in SCHISM but the error is generally small. Mass conservation is enforced with upwind/TVD transport.
It's important to use proper method in checking the volume conservation. In particular, you should make use of the flag consv_check in param.in, instead of calculating flows using the outputs yourself. This is because the velocity outputs have interpolated from side centers to nodes. Strictly speaking, checking volume conservation in SCHISM entails a calculation of semi-implicit fluxes to be consistent with the original formulation (see transport solver in the online manual).
Another thing to be aware of is the boundary condition. For example, the imposed river flow is conserved downstream, but if you also impose the elevation b.c., this is not the case anymore. The elev. b.c. as an essential b.c. over-rules the flow b.c. (a Neumann b.c.).
Here is a document that describes mass conservation with settling velocity.
