Velocity in diferent layers

MOHID Water forum. Questions and discussion around MOHID Water
gdreguero@hotmail.com
MOHID Beginner
Posts: 7
Joined: Thu Jan 01, 1970 00:00

Velocity in diferent layers

Postby gdreguero@hotmail.com » Thu Jan 28, 2010 14:31

Hi,
I´m trying to simulate the currents from an South Atlantic area. As first step was considered, as forcing of the system, velocities in X and Y with time-varying amplitudes (IN_TIDES archive). The domain is divided in 20 layers, but the velocity results are the same for any of them. Do you have any idea of what I am doing wrong?

Thanks a lot.

administrator@.
Site Admin
Posts: 349
Joined: Thu Jan 01, 1970 00:00
Location: Instituto Superior Tecnico - Lisbon - Portugal
Contact:

Postby administrator@. » Mon Feb 15, 2010 09:05

Hi,

Can you post some more information about this. What options have you selected? Can you post the Hydrodynamic, Turbulence and tide file?

Best regards,
Luis

gdreguero@hotmail.com
MOHID Beginner
Posts: 7
Joined: Thu Jan 01, 1970 00:00

Postby gdreguero@hotmail.com » Thu Feb 25, 2010 17:06

Hi Luis,
then I attach the contents of the files you asked.

Hydrodynamic_1.dat (the rest of keywords are zero)

BIHARMONIC_COEF : 1.00E+09
DISCRETIZATION : 2
DT_OUTPUT_TIME : 600
EVOLUTION : Solve_Equations
HMIN_ADVECTION : 0.5
HMIN_CHEZY : 0.1
HORIZONTALADVECTION : 1
HORIZONTALDIFFUSION : 1
IMPLICIT_HORADVECTION : 0
IMPLICIT_VERTADVECTION : 1
IMPLICIT_VERTDIFFUSION : 1
INITIAL_ELEVATION : 0
INITIAL_ELEVATION_VALUE : 0
INITIAL_VEL_U : 0
INITIAL_VEL_V : 0
LOCAL_DENSITY : 0
LOCAL_SOLUTION : 1
RESIDUAL : 0
SLIPPING_CONDITION : 1
TIDE : 1
TIDEPOTENTIAL : 0
TIME_SERIE : 1
TIME_SERIE_LOCATION : G:\\MOHID-Preliminares\\GeneralData\\TimeSerieLocation_02feb.dat
UP_CENTER : 1
UPSTREAM : Upwind
VELNORMALBOUNDARY : 1
VELTANGENTIALBOUNDARY : 1
VERTICALADVECTION : 1
VERTICALDIFFUSION : 1
VMIN_CHEZY : 0.1
VOLUMEVARIATION : 1
OUTPUT_TIME : 0 600


Turbulence_1.dat
VISCOSITY_V : 0.0010
VISCOSITY_H : 10.0000


I made several runs changing values for RADIATION, LOCAL_SOLUTION, NULL_BOUND_HORADV, VELNORMALBOUNDARY, VELTANGENTIALBOUNDARY, VISCOSITY_H and VISCOSITY_V.

Tide Components (The original archive have 3 more points like is shown bellow, with diferent locations inside grid.)

<begingauge>
NAME : velocidad 1
LONGITUDE: -63 12 0
LATITUDE : -49 36 0
GRID_I : 99
GRID_J : 2
!METRIC_X : 485621.9
!METRIC_Y : 4509678
REF_LEVEL : 0
TIME_REF : 0
EVOLUTION : Time Serie
EVOLUTION_VEL : Time Serie
TIME_SERIE_FILE : G:\\MOHID-Preliminares\\GeneralData\\VelConst-Minutes.dat
LEVEL_COLUMN : 2
VELU_COLUMN : 3
VELV_COLUMN : 4
<endgauge>


TimeSerie: VelConst-Minutes.dat

TIME_UNITS : MINUTES
SERIE_INITIAL_DATA : 2002 10 1 0 0 0

<BeginTimeSerie>
0 1 1 0
30 1 1 0
60 1 1 0
90 1 1 0
120 1 1 0
...
<EndTimeSerie>

administrator@.
Site Admin
Posts: 349
Joined: Thu Jan 01, 1970 00:00
Location: Instituto Superior Tecnico - Lisbon - Portugal
Contact:

Postby administrator@. » Wed Mar 03, 2010 18:00

Hi,

It is a bit tricky to do this. What you're doing is imposing velocities in the X direction with 1 m/s along the entire open boundary, and setting the water level to zero. This is a completely unrealistic situation and I don't think the model will handle it well.
Is your objective to impose a velocity field along the boundary, with different velocities at different grid cells and layers?

Luis

gdreguero@hotmail.com
MOHID Beginner
Posts: 7
Joined: Thu Jan 01, 1970 00:00

Postby gdreguero@hotmail.com » Fri Mar 05, 2010 14:19

Hi Luis, thanks for your respond. Yes, that is exactly what I intent to do. what is the way to do it?

Thanks again. Paula.

administrator@.
Site Admin
Posts: 349
Joined: Thu Jan 01, 1970 00:00
Location: Instituto Superior Tecnico - Lisbon - Portugal
Contact:

Postby administrator@. » Wed Mar 31, 2010 13:16

Hi Paula,

Sorry for the late answer. You may find some more info here:

http://www.mohid.com/forum/viewtopic.php?t=1124&highlight=

Hope this helps,
Luis


Return to “MOHID Water”

Who is online

Users browsing this forum: No registered users and 6 guests