Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

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

Open
utterances-bot opened this issue Mar 14, 2023 · 6 comments
Open

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

utterances-bot opened this issue Mar 14, 2023 · 6 comments

Comments

@utterances-bot
Copy link

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 utterances
Copy link

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 utterances
Copy link
Owner

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
Copy link

parkermac commented Mar 15, 2023 via email

@ahornerdevine utterances
Copy link

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 utterances
Copy link
Owner

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
Copy link

parkermac commented Mar 16, 2023 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
4 participants
@parkermac @ajleeson @utterances-bot @ahornerdevine