<div dir="ltr">Hi Paul,<div><br></div><div>Yes, the Orlanski open BCs are suitable for variable phase velocity, but it works perfectly only in theory. In reality, a number of numerical studies demonstrated some problems with the Orlanski scheme. For more details, you may see, e.g., 'Revisiting open boundary conditions from the point of view of characteristic variables' by E. Blayo *, L. Debreu, 2005, etc.</div><div><br></div><div>From my own experience with MITgcm applied for a real object (with realistic forcings, etc), the Orlanski scheme does not work very well, and the Dirichlet+sponge condition is the best option available in the model, even better than the Stevens OBCs implemented in MITgcm. It' s just my experience, maybe in other cases it will be the opposite.</div><div><br></div><div>---</div><div>Kind regards, </div><div>Stanislav Martyanov</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">пн, 18 апр. 2022 г. в 09:28, Burns, Paul <<a href="mailto:P.Burns2@exeter.ac.uk">P.Burns2@exeter.ac.uk</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Jody,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thanks for your message. </div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
My understanding is that the Orlanski scheme is suitable for cases where a spectrum of waves is excited. In Orlanski (1976) it shows that the scheme uses a variable phase velocity (varying in time, location along boundary, and with each prognostic variable)
to try and account for different waves reaching the boundary (shown in equation (2.2)). Orlanski also performed tests, computing the evolution of a collapsing bubble in stratified flow, which excited many different internal waves, with the results showing
only small distortion using the boundary condition (2.2). In the same paper it is also said that the sponge condition is not perfect and may create small distortions. It seems that neither boundary condition is perfect. </div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
My understanding is that the parameter Cmax that the users sets in MITgcm sets only the maximum phase speed that the Orlanski scheme accounts for, but that the scheme still computes variable phase speeds below this limit.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I've performed some initial tests using the Orlanski scheme and I'm not seeing obvious reflections at the boundary, with flow structures apparently leaving the domain without big distortions. But I wish to confirm things further than visual inspection, which
isn't sufficient alone.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Cheers,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Paul</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="gmail-m_-1911444968858732682Signature">
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div style="margin:0px;font-size:15px;color:rgb(32,31,30);text-align:start;background-color:rgb(255,255,255)">
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="color:rgb(12,100,192)"><u> </u></span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">Dr Paul Burns MInstP FRMetS</span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">Research Fellow</span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">Geophysical and Astrophysical Fluid Dynamics Group</span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">Department of Mathematics</span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">Harrison Building</span></p>
<p style="font-size:11pt;font-family:Calibri,sans-serif;margin:0px"><span style="margin:0px;color:rgb(31,78,121)">University of Exeter</span></p>
</div>
<br>
</div>
</div>
</div>
</div>
<div id="gmail-m_-1911444968858732682appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-1911444968858732682divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" target="_blank">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Jody Klymak <<a href="mailto:jklymak@uvic.ca" target="_blank">jklymak@uvic.ca</a>><br>
<b>Sent:</b> 17 April 2022 17:43<br>
<b>To:</b> <a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a> <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Subject:</b> Re: [MITgcm-support] Orlanski radiation condition parameters</font>
<div> </div>
</div>
<div style="overflow-wrap: break-word;">
<div style="background-color:rgb(255,235,156);width:100%;border-style:solid;border-color:rgb(156,101,0);border-width:1pt;padding:2pt;font-size:10pt;line-height:12pt;font-family:Calibri;color:black;text-align:left">
<span style="color:rgb(156,101,0);font-weight:bold">CAUTION:</span> This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe.</div>
<br>
<div><br>
<div><br>
Hi Paul,</div>
<div><br>
</div>
<div>Others should correct me if I am mistaken, but I don’t think Orlanski will work particularly well for this case because you do not have one dominant wave speed in your problem. Rather, you might try a gentle sponge. </div>
<div><br>
</div>
<div>Cheers, Jody<br>
<br>
</div>
<div><br>
<blockquote type="cite">
<div>On Apr 12, 2022, at 11:58 AM, Burns, Paul <<a href="mailto:P.Burns2@exeter.ac.uk" target="_blank">P.Burns2@exeter.ac.uk</a>> wrote:</div>
<br>
<div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Hi,</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
This is a bit basic, but I would like to check my implementation of my lateral boundary conditions. I am running a 2D simulation and prescribe (Dirichlet conditions) a Garrett-Munk spectrum of waves on my left boundary, which then propagate into my domain.
I would like to allow the motions in the interior solution to radiate away through my right boundary, so I've asked the model to use the Orlanski condition on the right boundary. </div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
In Orlanski (1976) it states that the (variable) phase velocity C must be less than the maximum numerical velocity, which for my setup is dx/dt = 25/1 = 25 m/s. So in data.obcs I've set Cmax = 20 m/s. This is quite different to the value used in the 'Gravity
Plume On a Continental Slope' tutorial, which has a similar domain size to mine, but this is because of my much smaller timestep due to my fine vertical grid required for my features of interest. I realise I am using large aspect ratio cells, but this is
to keep the simulation feasible with the compute resources I have. I will try and test sensitivity of my results to horizontal resolution.</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Can someone confirm how you fix the remaining parameter of the Orlanski scheme? T<span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">hat is, </span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<div style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgba(128,128,128,0.5)">
<span style="font-variant-ligatures:no-common-ligatures"><span> </span>cVelTimeScale</span></div>
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Many thanks,</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Paul</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<div id="gmail-m_-1911444968858732682x_Signature">
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<div style="margin:0px;font-size:15px;color:rgb(32,31,30);text-align:start;background-color:rgb(255,255,255)">
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="color:rgb(12,100,192)"><u> </u></span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">Dr Paul Burns MInstP FRMetS</span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">Research Fellow</span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">Geophysical and Astrophysical Fluid Dynamics Group</span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">Department of Mathematics</span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">Harrison Building</span></div>
<div style="margin:0px;font-size:11pt;font-family:Calibri,sans-serif">
<span style="margin:0px;color:rgb(31,78,121)">University of Exeter</span></div>
</div>
<br>
</div>
</div>
</div>
</div>
<span style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">_______________________________________________</span><br style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<span style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">MITgcm-support
mailing list</span><br style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<a href="mailto:MITgcm-support@mitgcm.org" style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">MITgcm-support@mitgcm.org</a><br style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<a href="https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=04%7C01%7CP.Burns2%40exeter.ac.uk%7C361fb4c681ad47b6905908da20918ac9%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C637858106991748874%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=s8yjYgBQE1fdMLE562bef7JM2U4FDu6JXJizyQilCPk%3D&reserved=0" style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>