search

ajleeson / research_blog

Weekly updates of research progress
0 stars 0 forks source link

research_blog/2023/03/14/lwsrc-idealized-run #25

Open utterances-bot opened 1 year ago

utterances-bot commented 1 year ago

Looking at LwSrc in Idealized Runs | Weekly Research Update

First order of business: Happy Pi Day! Does a tart count as a pie? This week I tested out a couple of different combinations of mixing coefficients and advectoin schemes, but still nothing seems to improve conditions at Birch Bay WWTP. I then focused my attention on the idealized estuary. I re-ran the model for one day with point sources discharging at a rate of 1 m3/s (rather than 1000 m3/s). I also disabled LtracerSrc so that the inflowing water enters at the same temperature and salinity as the water already in the grid cell. The results are not as extreme as what we observe at point sources in LiveOcean, however there are still interesting vertical velocity patterns in the water column. More details below. Mixing Coeff

https://ajleeson.github.io/research_blog/2023/03/14/lwsrc-idealized-run.html

parkermac commented 1 year ago

Aurora, small note: I think in the code (where it adds the tracer) turning off LTracerSrc only has an effect if you are not using MPDATA. PM

ajleeson commented 1 year ago

Hi Parker,

Thanks for checking in! I ran a quick test, and it seems that LTracerSrc = F works even when using MPDATA. Below is a salinity section across a freshwater point source with LTracerSrc = F: ltracersrcF

And here is one for LTracerSrc = T: ltracersrcT

MPDATA was used in both cases.

Given these results, I think the idealized estuary results in the blog post shouldn't be impacted by buoyancy effects.

Aurora

parkermac commented 1 year ago

Got it. Nice checking.

On Wed, Mar 15, 2023 at 10:11 AM ajleeson @.***> wrote:

Hi Parker,

Thanks for checking in! I ran a quick test, and it seems that LTracerSrc = F works even when using MPDATA. Below is a salinity section across a freshwater point source with LTracerSrc = F: [image: ltracersrcF] https://user-images.githubusercontent.com/15829099/225387168-b0e665d2-aabf-4b22-a021-793e6f71cedd.png

And here is one for LTracerSrc = T: [image: ltracersrcT] https://user-images.githubusercontent.com/15829099/225387176-a3169bd6-5465-48c3-bf51-24eeb3cf217f.png

MPDATA was used in both cases.

Given these results, I think the idealized estuary results in the blog post shouldn't be impacted by buoyancy effects.

Aurora

— Reply to this email directly, view it on GitHub https://github.com/ajleeson/research_blog/issues/25#issuecomment-1470428459, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC5WRXZUGIB5HYIUE2X2NETW4HZ3HANCNFSM6AAAAAAV25AL4U . You are receiving this because you commented.Message ID: @.***>

--

Parker MacCready

Leo Maddox Endowed Professor in Oceanography

Research Professor, School of Oceanography I live and work on the land and waters of the Squaxin Island Tribe. Email: @.** URL: faculty.washington.edu/pmacc LiveOcean Daily Forecasts http://faculty.washington.edu/pmacc/LO/LiveOcean.html Cell: (360) 359-1936, Office: (206) 685-9588 pronouns: he/him*

ahornerdevine commented 1 year ago

First, we had quiche in celebration of Pi day (and an apple pie), so I definitely think tart counts. Second, staring at the video in Fig 6 it looks like the 'relative' v-velocity at the source is always negative over the whole water column. Can you try subtracting the tidal velocity (eg the velocity profile at a location away from the source) from the velocity profile at the source? This makes it look like the addition of this very small amount of fluid is significantly influencing the horizontal momentum, which it shouldn't. Is there any way that the initial horizontal velocity of the added fluid is wrong? If you require that the initial velocity of the added fluid matches that of the ambient fluid (which is what I'd do), then you are effectively adding a tiny bit of momentum.

ajleeson commented 1 year ago

Hi Alex,

Thanks for the input! It is very strange that the horizontal velocity at the source is different than elsewhere. I'll subtract the tidal velocity and see what it looks like. Like you mentioned, maybe this will give us more clues about how the sources are implemented in ROMS.

Aurora

parkermac commented 1 year ago

Note that the velocity from a mooring extraction is just pulled in by nearest-neighbor, so it is not interpolated to the rho-grid.

In contrast, the section making code (if you are using my new version) does interpolate everything.

PM

On Thu, Mar 16, 2023 at 9:24 AM ajleeson @.***> wrote:

Hi Alex,

Thanks for the input! It is very strange that the horizontal velocity at the source is different than elsewhere. I'll subtract the tidal velocity and see what it looks like. Like you mentioned, maybe this will give us more clues about how the sources are implemented in ROMS.

Aurora

— Reply to this email directly, view it on GitHub https://github.com/ajleeson/research_blog/issues/25#issuecomment-1472296941, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC5WRX7GRFXCMZQFKQCO57TW4M5CXANCNFSM6AAAAAAV25AL4U . You are receiving this because you commented.Message ID: @.***>

--

Parker MacCready

Leo Maddox Endowed Professor in Oceanography

Research Professor, School of Oceanography I live and work on the land and waters of the Squaxin Island Tribe. Email: @.** URL: faculty.washington.edu/pmacc LiveOcean Daily Forecasts http://faculty.washington.edu/pmacc/LO/LiveOcean.html Cell: (360) 359-1936, Office: (206) 685-9588 pronouns: he/him*