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

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.

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
Note: See TracTickets for help on using tickets.