Opened 16 years ago
Last modified 15 years ago
#273 new enhancement
Implement Kinematic Viscosity
Reported by: | ole | Owned by: | steve |
---|---|---|---|
Priority: | normal | Milestone: | ANUGA enhancements |
Component: | Functionality and features | Version: | |
Severity: | normal | Keywords: | |
Cc: |
Description
Ted and Rudy to provide formula for KV forcing term. Tests will show eddy formations.
Change History (4)
comment:1 Changed 16 years ago by ole
comment:2 Changed 16 years ago by ole
New idea: Using the existing architecture we could compute the second order derivatives in two iterations by taking the gradients of gradients. In other words store the first order derivatives when computed and then use the same algorithm to compute the second order derivatives. Information beyond immediate neighbours would therefore be implicit.
comment:3 Changed 16 years ago by ole
In phone conversation with Stephen Roberts early September 2008 he suggested that we used the second order derivative already implemented in the smoothing algorithm in ANUGA. In other words we could use the Finite-Element solver to approximate the viscosity term for use with the Finite-volume scheme.
Moreover, the majority of time spent in the smoothing algorithm has to do with the least squares process, so the second order derivative on its own is likely to be rather fast.
comment:4 Changed 15 years ago by ole
- Milestone set to ANUGA enhancements
- Owner changed from ole to steve
Stephen Roberts reckons it'll be a bit involved as one needs to approximate the second order derivative and this will involve going beyond immediate neighbours. Obviously, this will take some effort to get right especially near the boundaries.