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.

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

Question to Rw in the strain-only parametrization #118

Open
opinner opened this issue Feb 18, 2024 · 0 comments
Open

Question to Rw in the strain-only parametrization #118

opinner opened this issue Feb 18, 2024 · 0 comments

Comments

@opinner
Copy link

opinner commented Feb 18, 2024

In the function shearstrain(), one finds the code snippet

# Strain only parameterization
# Use assumed shear/strain ratio of 3
Rw = 3
eps_st[iwin] = eps_strain(eps0, Nm, N0, Sst, Sstgm, Rw, f)
krho_st[iwin] = diffusivity(eps_st[iwin], Nm, Gam=Gam0)

Inside eps_strain(), the function aspect_ratio_correction_st(Rw) is called, which evaluates to $1$ for Rw = 3, making the whole function superfluous. But observations point to possible higher values of $R_\omega = 7\pm3$ (Kunze et al., 2006 [p.1560]).
This value changes aspect_ratio_correction_st(Rw = 7) = 2.69. Currently, I manually multiply the end result for the dissipation rate with this factor, to change my $R_\omega$ value.

Would it be possible to make Rw_strain = 3 an optional argument of the shearstrain() function? Or does that only clutter the argument list, as a workaround is possible?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant