Welcome, Guest
Username: Password: Remember me

TOPIC: Sharp salt-wedge interface

Sharp salt-wedge interface 1 month 1 week ago #45640

  • giuliaperego
  • giuliaperego's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 22
  • Thank you received: 1
Hello.
I am trying to model a salt wedge in a river estuary.
I would like to have a sharper interface between high salinity and freshwater in my results. While here I have a high mixing and several kilometers of intermediate salinity (see picture).
What can I change ?
I attach my .cas file.
Thank you very much
Giulia
Attachments:
The administrator has disabled public write access.

Sharp salt-wedge interface 3 weeks 1 day ago #45663

  • giuliaperego
  • giuliaperego's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 22
  • Thank you received: 1
Hello,
I would be very grateful if anybody could give me some advice
The administrator has disabled public write access.

Sharp salt-wedge interface 3 weeks 1 day ago #45665

  • Sokolov
  • Sokolov's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 38
  • Thank you received: 13
Hello,

There are too many things to play with :)

Check your boundary condition:

PRESCRIBED ELEVATIONS = 0;0
PRESCRIBED FLOWRATES = 0;3.0

Where water should go on the first boundary? To start with you can try to set river discharge on one end and set free sea level on both end (check you .cli file, and set code 4 for the level)

Good luck
The administrator has disabled public write access.

Sharp salt-wedge interface 2 weeks 4 days ago #45683

  • pham
  • pham's Avatar
  • OFFLINE
  • Administrator
  • Posts: 1559
  • Thank you received: 602
Hello Giulia,

If I were you, I would read the recommendation appendix (E) and would try a few tips (maybe not suitable to your model, but let's try):
- IMPLICITATION FOR DEPTH = 0.55 (default), 0.6 or a value up to 0.9
- IMPLICITATION FOR VELOCITIES = 0.55 (default), 0.6 or a value up to 0.9

I would also change the vertical turbulence model rather than constant viscosity along the vertical:
e.g. VERTICAL TURBULENCE MODEL = 3 (k-epsilon). You can also try = 2 (mixing length) and try different mixing length models (perhaps with damping function). I do not know if they could suit your model, but you can also try = 5 (Spalart-Allmaras) or = 7 (k-omega). You can also try what is available in the GOTM library.

Hope this helps,

Chi-Tuan
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.