Thatcher - Harleman time lagThatcher - Harleman time laghttps://oss.deltares.nl/c/message_boards/find_thread?p_l_id=1806765&threadId=14509252020-06-03T03:42:46Z2020-06-03T03:42:46ZRE: Thatcher - Harleman time lagKonstantinos Matsoukishttps://oss.deltares.nl/c/message_boards/find_message?p_l_id=1806765&messageId=14779842018-03-07T10:01:22Z2018-03-07T10:01:22ZDear Erik,<br /><br />Thank you for your help and suggestions. Indeed I am using a sigma model.<br />I have managed to surpass the issue by setting the horizontal diffusion equal to zero at the last two grid lines at the open boundaries.<br />In this way, salinity is not constrained anymore by the boundaries. <br /><br />Thanks,<br />KostasKonstantinos Matsoukis2018-03-07T10:01:22ZRE: Thatcher - Harleman time lagErik de Goedehttps://oss.deltares.nl/c/message_boards/find_message?p_l_id=1806765&messageId=14721062018-03-01T15:25:17Z2018-03-01T15:25:17ZDear Kostas,<br /><br />From your figure it is clear that the problem occurs at outflow. However, at outflow the TH boundary condition does not play a role. (NB. This only has some effect at inflow.) I guess that the problem is due to horizontal diffusion at theopen boundaries. Do you apply a sigma-model? (NB. In a Z-model the horizontal diffusion at open boundaries is neglected. However, in a sigma model this isn't the case.) <br />Assuming that you are using a sigma model, you can reduce the horizontal diffusion near open boundaries by applying spatially varying horizontal diffusion values. As a simple test you might do a simulation with the horizontal diffusion set to zero. <br /><br />Wondering whether the explanation above is valid or not and with kind regards,<br /><br />Erik de Goede, DeltaresErik de Goede2018-03-01T15:25:17ZRE: Thatcher - Harleman time lagKonstantinos Matsoukishttps://oss.deltares.nl/c/message_boards/find_message?p_l_id=1806765&messageId=14712042018-02-28T17:53:22Z2018-02-28T17:53:22ZDear Erik,<br /><br />Thank you very much for your reply.It is very explanatory.<br />I did a simulation without TH and the results look good (attached an example of salinity contours and flow vectors I get without TH).<br />However, it seems that the plume is forced in this case by the prescribed boundary conditions that I give for salinity equal to 30.<br />My flow boundaries are: Riemann at the north and south, water elevation equal to zero at the offshore. <br />I was looking for an appropriate way to free the transport at these boundaries and allow the contours to surpass if possible the boundaries.<br />Could this be done by implementing the TH condition and what are the restrictions in this case please? Could the TH be applied for example in a water level boundary?<br />Is there any other method to free the boundaries for salinity?<br /><br />Thank you,<br />KostasKonstantinos Matsoukis2018-02-28T17:53:22ZRE: Thatcher - Harleman time lagErik de Goedehttps://oss.deltares.nl/c/message_boards/find_message?p_l_id=1806765&messageId=14710562018-02-28T14:58:11Z2018-02-28T14:58:11ZDear Kostas,<br />The Thatcher Harleman boundary condition is meant to suppress unphysical currents near open boundaries. At first, try a simulation without TH boundary conditions. If unphysical currents occur, try positive TH boundary values; for example of 60 (minutes). At a last alternative try a negative TH value, which corresponds to a so-called Neumann boundary condition (d c / d x = 0).<br /><br />In general, the TH boundary conditions have a marginal effect on the model results. If all options for the TH boundary condition give unrealistic results, then other options come into play (e.g., locally adaptation of the depth or the horizontal diffusivity). This is of later interest, At first, try zero or positive values for the TH boundary conditions. Success!<br /><br />Erik de Goede, DeltaresErik de Goede2018-02-28T14:58:11ZThatcher - Harleman time lagKonstantinos Matsoukishttps://oss.deltares.nl/c/message_boards/find_message?p_l_id=1806765&messageId=14509242018-02-07T15:20:48Z2018-02-07T13:54:33ZHi,<br /><br />I am trying to model salinity flow in a rectilinear grid with DELFT3D-FLOW. The model is forced only by a river discharge coming from the western boundary. <br />I want my north,south and offshore boundaries to be open so in my first test case, I prescribed the Riemann invariant (taken zero) at the north and south boundaries and water level equal to zero at the offshore boundary. For salinity calculations, at first a uniform value equal to 30psu is set as initial conditions and 30 at the other boundaries too (north,south,offshore). This simple test seemed to work fine at first with the flow vectors extracting from the channel to the outer field and developing a logical pattern.<br />But the transport boundaries are required to be of a free exit as well. Following the instructions given by the manual, this is done by activating the Thatcher-Harleman boundary condition and introducing a negative value of the required time of residence in the relationship as this ensures salinity transport to be calculated exclusively from the model and not influenced from the outside (ie free exit). In this case though, this transport condition results in a flow parallel to the north and south boundaries instead of being an exit flow, while a reverse flow occurs at the offshore boundary caused obviously by inappropriate density gradients. <br />In other words, applying a free exit condition for the salinity equation at the boundaries makes the flow to become parallel to the boundary or reversed back into the model area instead of getting through the boundary unobstructed. Is there a reason for this to happen ? <br />Any suggestions of how this could be achieved, that is, to have both flow and salinity transport free exit conditions at the open boundaries and the proper physical results as well? <br />Are there any other conditions that could be applied with or without the Thatcher-Harlemann option to overcome this issue?<br /> <br /><br />Thank you,<br />KostasKonstantinos Matsoukis2018-02-07T13:54:33Z