include convective driving in gyre

[anwesha108]

Hi,

I am working with the driving of g-modes in DAVs.. now, the way that nad works is by 'frozen convection' where FROZEN_PESNELL_4 gives the driving but it is ignoring the lagrangian perturbation of the convective flux.
while scheme1
conv_term = lambda*c_rad*(3._WP + dc_rad)/(c_1*alpha_omg*omega_c**2)
completely ignores the convective flux,
scheme4
conv_term = lambda*(c_lum*(3._WP + dc_lum) - (c_lum - c_rad))/(c_1*alpha_omg*omega_c**2)
just ignores the lagrangian variation of the convective flux.

I was wondering, if there was any way to include the lagrangian variation of the convective flux to the equations. I tried by modifying the second equation in gyre_nad_eqns.fpp, to

conv_term = lambda*c_lum*(3._WP + dc_lum)/(c_1*alpha_omg*omega_c**2)
but that is not producing any different result than the original one (scheme4)..

can I get some help on how to go about this one so that I can include that term in calculation.
Thank you!

[rhtownsend]

Hi,

I am working with the driving of g-modes in DAVs.. now, the way that nad works is by 'frozen convection' where FROZEN_PESNELL_4 gives the driving but it is ignoring the lagrangian perturbation of the convective flux.
while scheme1
conv_term = lambda*c_rad*(3._WP + dc_rad)/(c_1*alpha_omg*omega_c**2)
completely ignores the convective flux,
scheme4
conv_term = lambda*(c_lum*(3._WP + dc_lum) - (c_lum - c_rad))/(c_1*alpha_omg*omega_c**2)
just ignores the lagrangian variation of the convective flux.

I was wondering, if there was any way to include the lagrangian variation of the convective flux to the equations. I tried by modifying the second equation in gyre_nad_eqns.fpp, to

conv_term = lambda*c_lum*(3._WP + dc_lum)/(c_1*alpha_omg*omega_c**2)
but that is not producing any different result than the original one (scheme4)..

can I get some help on how to go about this one so that I can include that term in calculation.
Thank you!

Hi Anu --

As a start to answering your questions, it might be useful to be more specific with terminology. For instance, "FROZEN_PESNELL_4" neglects the the Lagrangian perturbation to the convective luminoisity. This is *not* the same as neglecting the perturbation to the convective flux. You can see this from equation 21.4 of Unno et al (1989):

delta L = 4 pi r^2 (delta F_r + 2 xi_r/r F_r)

This applies to the total flux/luminosity, *and* their respective radiative and convective components. So, setting delta L_conv = 0 means that

delta F_r,conv = - 2 xi_r/r F_r,conv

i.e., the Lagrangian perturbation to the convective flux is not zero!

cheers,

Rich