Issue 
A&A
Volume 646, February 2021



Article Number  L5  
Number of page(s)  7  
Section  Letters to the Editor  
DOI  https://doi.org/10.1051/00046361/202040259  
Published online  05 February 2021 
Letter to the Editor
Surface effects and turbulent pressure
Assessing the GasΓ_{1} and ReducedΓ_{1} empirical models
^{1}
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 place Jules Janssen, 92195 Meudon, France
email: kevin.belkacem@obspm.fr
^{2}
Univ. of Appl. Sciences Technikum Wien, Dept. Applied Mathematics and Physics, Höchstädtplatz 6, 1200 Wien, Austria
^{3}
WolfgangPauliInstitute c/o Faculty of Mathematics, University of Vienna, OskarMorgensternPlatz 1, 1090 Wien, Austria
^{4}
MaxPlanckInstitut für Sonnensystemforschung, JustusvonLiebig Weg 3, 37077 Göttingen, Germany
Received:
29
December
2020
Accepted:
14
January
2021
The application of the full potential of stellar seismology is made difficult by the improper modelling of the uppermost layers of solarlike stars and their influence on the modelled frequencies. Our knowledge of these socalled ‘surface effects’ has improved thanks to the use of 3D hydrodynamical simulations, however, the calculation of eigenfrequencies relies on empirical models for the description of the Lagrangian perturbation of turbulent pressure, namely: the reducedΓ_{1} model (RGM) and the gasΓ_{1} model (GGM). Starting from the fully compressible turbulence equations, we derived both the GGM and RGM models by using a closure to model the flux of turbulent kinetic energy. We find that both models originate from two terms: the source of turbulent pressure due to compression produced by the oscillations and the divergence of the flux of turbulent pressure. We also demonstrate that they are both compatible with the adiabatic approximation and, additionally, that they imply a number of questionable assumptions, mainly with respect to mode physics. Among other hypotheses, it is necessary to neglect the Lagrangian perturbation of the dissipation of turbulent kinetic energy into heat and the Lagrangian perturbation of buoyancy work.
Key words: waves / convection / Sun: oscillations
© K. Belkacem et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Systematic differences between observed and modelled eigenfrequencies is a longstanding problem in stellar seismology. For the Sun and solarlike stars, this has been identified and found to be associated with an improper modelling of their uppermost layers (e.g. Brown 1984; Dziembowski et al. 1988; ChristensenDalsgaard & Thompson 1997). To circumvent this problem, wellchosen combinations of frequencies (e.g. Roxburgh & Vorontsov 2003) or adhoc corrections (Kjeldsen et al. 2008; Ball & Gizon 2014; Sonoi et al. 2015) are commonly used. Nevertheless, to be able to exploit all the information contained in the observed frequencies, the physics underlying the surface effects must be understood.
To this end, Rosenthal et al. (1999) used 3D hydrodynamical simulations, which allowed them to account for the mean turbulent pressure (as well as convective backwarming; see Trampedach et al. 2017) in the equilibrium structure. Their work was followed by Piau et al. (2014), Magic & Weiss (2016), Houdek et al. (2017), Jørgensen & Weiss (2019), Schou & Birch (2020) for the Sun and Sonoi et al. (2015, 2017), Ball et al. (2016), Trampedach et al. (2017), Jørgensen et al. (2017, 2018, 2019, 2021), Manchon et al. (2018), Jørgensen & Angelou (2019), Houdek et al. (2019), Mosumgaard et al. (2020) for solarlike stars. A drawback of this approach is the need to compute adiabatic eigenfrequencies using empirical models to describe the Lagrangian perturbation of turbulent pressure.
For this purpose, Rosenthal et al. (1999) introduced two empirical models. The reducedΓ_{1} model (RGM), which is based on the argument that the perturbation of turbulent pressure is negligible compared to the perturbation of gas pressure, and the gasΓ_{1} model (GGM), which is based on a consideration that the perturbations of gas pressure and turbulent pressure are in phase. The authors favoured the GGM as it permits them to better reproduce the observed solar frequencies. Most of the abovementioned works then applied the GGM, except in a few cases (e.g. Jørgensen & Weiss 2019), where the RGM was guided by the nonadiabatic calculation from Houdek et al. (2017). However, the latter was recently challenged by Schou & Birch (2020), using the eigenfunctions as inferred directly from 3D numerical simulations. Therefore, the issue of surface effects remains and it demands extensive theoretical investigation to gain insights into the physics of the problem.
In this work, we aim at deriving and then assessing the theoretical validity of the GGM and RGM empirical models.
2. Equation governing the turbulent pressure
When including the mean turbulent pressure in the adiabatic oscillation equations, the Lagrangian perturbation of turbulent pressure is to be prescribed. This is the role of the GGM and RGM models, which are extensively described in Appendix A. To derive them, a first step is to express the equation for turbulent pressure.
2.1. Averaged equation for turbulent pressure
We start with the equation for the averaged Reynolds stresses (e.g. Canuto 1997; Gatski & Bonnet 2013):
where the overbar denotes the Reynolds average and the tilde denotes the Favre average (see Appendix A for a definition), u_{i} is the i component of the velocity field, P_{g} is the gas pressure (it includes contributions from the radiative field as well as body forces), ρ is the density, τ_{jk} is the viscous stress tensor, and are the Reynolds stresses, where is the ith component of the velocity fluctuation around its Favre average (). We also employ the notations, ∂_{t} ≡ ∂/∂t and ∂_{i} ≡ ∂/∂x_{i}, Einstein’s notation for repeated indices, and the pseudoLagrangian derivative is defined by .
We then consider the rr component of Eq. (1) (where r is the radial coordinate) and identify the Reynolds average with the horizontal average. Thus, we obtain:
where is the turbulent pressure, , , and . Except for notational differences, Eq. (2) is strictly equivalent to Eq. (16b) in Canuto (1997). The first term of the righthand side of Eq. (2) is a source of turbulent pressure due to the compression produced by radial oscillations. The second term corresponds to the pressure work (or the buoyancy work), which is a source of turbulent pressure in convective regions. The three following terms in brackets correspond to transport terms. Finally, the last two terms are dissipative terms.
Several assumptions are now needed. The molecular diffusion flux () is neglected compared to the other fluxes appearing in Eq. (2). This is justified as we are considering fully developed turbulence with high Reynolds numbers (e.g. Nordlund et al. 2009). The last term of Eq. (2) is assumed to be proportional to the rate of dissipation of turbulent kinetic energy into heat, that is, . This standard assumption is made possible because dissipation by molecular forces occurs at almost isotropic small scales (e.g. Pope 2000). Finally, the pressuredilatation term () is neglected because it scales as the square of the turbulent Mach number (Sarkar 1992) (the maximum is about 0.3 in the Sun). By adopting those approximations, Eq. (2) is reduced to:
where is the dissipation rate of turbulent kinetic energy into heat.
2.2. Modelling the transport of turbulent pressure
The flux of turbulent pressure, in Eq. (3), is difficult to model in the uppermost layers of solarlike stars for which the downgradient approximation fails (see Canuto 2009; Kupka & Muthsam 2017). To overcome this issue, we adopt the closure initially proposed by Canuto (1992) (see also Canuto 2011) and supported by 3D direct numerical simulations (Kupka & Muthsam 2007). It reads:
where is the flux of turbulent kinetic energy, is the flux of turbulent dissipation, is the specific turbulent kinetic energy, and c is a parameter to be specified.
From Eq. (4), we then obtain:
where ω ≡ ϵ/k is the turbulence frequency and
The parameter α, which can be understood as the degree of anisotropy of the flux of turbulent kinetic energy, will be supposed to be known and obtained directly from the solar 3D numerical simulation (see Sect. 3.2).
Now, using Eq. (5), we can rewrite Eq. (3) as:
where one still needs a prescription for the flux of turbulent dissipation. To do so, we use the equation governing ϵ. A standard procedure consists in constructing a parametric ϵequation as is done in twoequation models (e.g. k − ϵ models). Following this approach, we can write (e.g. Pope 2000; Wilcox 2006; Gatski & Bonnet 2013):
The first term in the righthand side of Eq. (8) is related to the production of turbulent kinetic energy by compression. The second term represents for the effect of bulk compressions and expansions onto ϵ (Coleman & Mansour 1991). The third term is related to production of turbulent kinetic energy by the buoyancy work, and the last term is related to both the viscous destruction and production due to vortex stretching (Gatski & Bonnet 2013). The coefficients , and are discussed in the following section.
3. Recovering the gasΓ_{1} and reducedΓ_{1} models
Here, we derive an expression for the Lagrangian perturbation of turbulent pressure and make a number of assumptions and discuss their validity.
3.1. Perturbation of turbulent pressure
Using Eq. (7) together with Eq. (8), we obtain the desired expression for the equation governing turbulent pressure:
where we used the averaged continuity equation (Eq. (A.5)) and Φ ≡ r_{ii}/P_{t} is the anisotropy factor.
To go further, we assume the following hypothesis: (H1) The Lagrangian perturbation of the pressurevelocity fluctuations is neglected, that is, . (H2) The turbulence frequency (ω) is supposed to vary on a length scale that is much larger than the length scale of F_{iir}. Accordingly, the second term of the righthand side of Eq. (5) can be neglected and, thus, the last term of Eq. (9) vanishes. (H3) The Lagrangian perturbation of α is neglected, that is, δα = 0. (H4) The Lagrangian perturbation of the dissipation of turbulent kinetic energy into heat (δϵ) is neglected. (H5) The perturbation of the buoyancy work, , is neglected. (H6) The Lagrangian perturbation of density is assumed to be real.
Perturbing Eq. (9) and applying H1 to H6, we get:
where δP_{t} and δρ are the pseudoLagrangian perturbations of turbulent pressure and density, respectively. As an introduction to the next section, we can already mention that Eq. (10) allows us to recover the RGM and GGM models.
Tracing back the origin of Eq. (10) shows that it results from two terms in Eq. (2): the source of turbulent pressure due to compression produced by the oscillations and the divergence of the flux of turbulent pressure. It is also important to mention that points H1 to H6 are fully consistent with the assumptions made to obtain the equation for the perturbation of gas pressure (δP_{g}, see A.9 and Appendix B). Moreover, given H1, H3, and H6 and further neglecting the perturbation of the radiative flux, we find that H5 is equivalent to adopting the adiabatic limit (see Appendix C). Hence, both equations for δP_{t} (Eq. (10)) and δP_{g} (Eq. (A.9)) are compatible with the adiabatic approximation. Conversely, the adiabatic approximation on its own is not sufficient to derive Eqs. (10) and (A.9).
The first assumption, H1, is essentially equivalent to neglecting the perturbation of the convective flux () because is proportional to the convective flux as shown by Canuto (1997) (see his Eq. (35c)). Such an assumption is, however, not strictly valid because the perturbation of the convective flux does not vanish even in the adiabatic limit (e.g. Sonoi et al. 2017). H2 is difficult to properly assess because we need to determine ϵ and the simplest possible way to do this is to consider ϵ ∝ k^{3/2}/H_{p} (e.g. Pope 2000), where the dissipation lengthscale is assumed to scale as the pressure scaleheight (H_{p}). Using, the solar 3D simulation described in Belkacem et al. (2019), we may readily find that H2 is valid near the superadiabatic peak but not in the quasiadiabatic regions. For H3, it is equivalent to assuming that perturbations of horizontal and vertical contributions of turbulent kinetic energy adjust instantaneously to each other and are thus in phase. A look at the thirdorder equation on fluxes (see Canuto 1997) demonstrates that the situation is much more complex because many terms are capable to introduce some redistribution and, in turn, phase shifts, when perturbed. Concerning H4 and H5, it consists of neglecting the Lagrangian perturbation of both the turbulent dissipation and the buoyancy work. For mode damping, these two contributions exactly compensate the contribution of turbulent pressure in the limit of a vanishing flux of turbulent kinetic energy and with Γ_{3} − 1 = 2/3, where Γ_{3} ≡ (∂lnT/∂lnρ)_{s} (e.g. Ledoux & Walraven 1958; Grigahcène et al. 2005). More recently, Belkacem et al. (2019) demonstrated using the normal modes of a 3D solar hydrodynamical simulation that the perturbation of both terms plays an essential role for the mode damping rates. Hence, as modal surface effects also partly rely on the phase differences between the perturbations of density and turbulent pressure, the impacts of those contributions to surface effects ought to be assessed definitively. Finally, H6 assumes that adding the mean turbulent pressure to the hydrostatic equilibrium only introduces a negligible phase shift to the perturbation of density in the adiabatic limit. These remarks lead to questions regarding the validity of both the RGM and GGM models because, even in the adiabatic limit, they introduce oversimplifying hypotheses regarding the properties of turbulent convection and mode physics.
3.2. Application to solar pmode frequencies
Now we go on to investigate how Eq. (10) permits us to recover the GGM and RGM models. One prerequisite is the specification of the coefficients, and . A commonly accepted value for the first coefficient is (e.g. Pope 2000; Wilcox 2006; Gatski & Bonnet 2013). For , we adopt the model of Coleman & Mansour (1991) which gives:
where n ≃ 0.75 is the exponent of the viscosity law on temperature. Indeed, Eq. (11) is aimed at accounting for the effect of bulk compressions and expansions onto ϵ due to the dependence of the viscosity on temperature (Coleman & Mansour 1991). We note that Eq. (11) was obtained using the rapid distorsion theory (see Hunt & Carruthers 1990; Cambon et al. 1993, for reviews) and under the assumption of adiabatic compression. In the quasiadiabatic regions of the solar convection zone, those approximations are relatively applicable because the modal period is much shorter than both the typical turnover time scale and the thermal time scale. However, in the superadiabatic layers, those time scales become of comparable magnitude and the validity of Eq. (11) becomes questionable. Therefore, the adopted value of used in this work must be considered a guideline, rather than a firm value.
In this way, we computed the frequency differences (for radial modes) between the observed frequencies, taken from Broomhall et al. (2009) and Davies et al. (2014), and theoretical frequencies computed with a classical shooting method. For the latter, we used the RGM, GGM, and the model developed in this work using several values of the closure coefficient, c (see Eq. (4)). Theoretical frequencies have been obtained by integrating Eqs. (A.7)–(A.9), complemented by either Eq. (A.12) for the RGM, Eq. (A.13) for the GGM, and Eq. (10) for the model proposed in this article. The equilibrium model had been obtained by patching a CESTAM model together with a solar ANTARES 3D simulation (see Belkacem et al. 2019, for details) using the same methodology as described in Sonoi et al. (2015). The 3D model is however not exactly solar because the effective temperature is 5750 ± 18 K and the chemical composition is (X, Y, Z) = (0.7373, 0.2427, 0.0200) with Grevesse et al. (1993) mixture. To allow for a comparison with observed solar frequencies, the 3D model is patched with a CESTAM model with a helium abundance of 0.2485 and the resulting frequencies are rescaled to match the standard solar frequency Hz (G is the gravitational constant, M_{⊙} the solar mass, and R_{⊙} the solar radius). Such a procedure is sufficient for our purposes, even if a frequency shift remains at low frequencies, because we are interested in differential effects. From this patched model, all the equilibrium quantities have been inferred.
Figure 1 shows that both the frequency differences obtained using the GGM and the RGM models can be recovered by adjusting the value of the closure coefficient, c, introduced in Eq. (4). First, we considered what guidance for c we might take from existing 3D numerical simulations. A value of c ≈ 0.6 was derived from four direct numerical simulations of fully compressible convection, two of which were discussed in Kupka & Muthsam (2007). Those simulations were done with the ASCIC3 code (Muthsam et al. 1995), which is used for such work because it does not introduce, directly or indirectly, subgrid scale viscosities: dissipation takes place by way of explicit, physical viscosity and a very small contribution from time integration only. The data were collected for a total of 12 different model configurations (five of which have been discussed in detail in Muthsam et al. 1995, 1999). They covered a range in Prandtl number Pr from 0.1 to 1 at a ‘zone Rayleigh number’ of 10^{5} to 10^{6}. For a horizontal domain width eight times the size of this zone depth and with three to four granules found along each horizontal direction, this implies a ‘granulation diameterbased Rayleigh number’ that is ∼16 to 50 times larger or an Ra in the range of 10^{6} to 5 × 10^{7}. This yields a product of Ra and Pr in the range of 10^{5} to 5 × 10^{7}. That (squared) ratio of the thermal diffusion time scale to the buoyancy time scale agrees with results for the upper part of the solar convection zone (Kupka et al. 2020) despite the fact that the convection zones are much more shallow than in the solar case, where additionally Pr ≪ 10^{−6}. Kupka & Muthsam (2007) found c ≈ 0.6 in all cases where sufficient numerical resolution had been ensured, also for cases not shown therein. Confirming those results for much lower values of Pr and a higher ratio of total flux to radiative flux would be useful. Recalling the discussion in Belkacem et al. (2019) on computing the dissipation rate from realistic solar simulations, we have to point out here that large eddy simulations, whether they use hyperviscosity, a Smagorinskytype subgrid scale model, or a Riemann solver, are not the best tools for computing or ϵ; those quantities depend on viscosity related processes which peak near small multiples of the grid scale. This could heavily bias computations of those quantities by the numerical method used. We thus consider conclusions based on low Pr direct numerical simulations for the physical range of Ra ⋅ Pr of interest as the safer way to estimate c. In comparison, the ratio α from Eq. (5) can be safely estimated from solar granulation simulations, as the contributions to this quantity peak at length scales resolved in those simulations.
Fig. 1.
Frequency difference for radial solar modes between the observed frequencies, as inferred by Broomhall et al. (2009) and Davies et al. (2014), and the modelled frequencies computed using the GGM model, the RGM model, and the model of the present work (see text for details). 
4. Conclusions
By using the Reynolds and Favreaveraged, fully compressible NavierStokes equations, we demonstrate that it is possible to develop a model which recovers both the RGM and GGM empirical models. Interestingly enough, this is based on a relation that has been shown to originate from a compensation between the source of turbulent pressure due to compression produced by the oscillations and the divergence of the flux of turbulent pressure. We show the RGM and GGM models are compatible with the adiabatic approximation but this also implies more drastic and unrealistic physical assumptions regarding turbulence and mode physics.
A comparison with solar frequencies shows that while recovering the RGM and GGM, the results are sensitive to the closure coefficients (c in Eq. (4) but also appearing in Eq. (10)). To consolidate the value of these parameters, only direct 3D numerical simulations prove useful. Unfortunately, this is still out of reach for the Sun due to the limits on our current numerical capacities. Hence, extrapolations from more accessible parameter ranges remain necessary.
It is thus difficult to draw a conclusion on which of the two models is more appropriate in application. What complicates the matter further, given the abovementioned assumptions which are needed to recover these empirical models, we can safely conclude none of them are firmly physically grounded, even in the adiabatic limit. However, we may still quantify the hypotheses (H1 to H6) on which both RGM and GGM models rely, and more precisely, their individual effect on mode frequencies. This will be made possible by either using a realistic treatment of turbulent convection based on a timedependent and nonadiabatic treatment or by using normal modes of direct 3D numerical simulations. For the latter, dedicated longduration simulations (in order to have a sufficient statistic and to resolve the normal modes) with a large spatial extension (to have a sufficient number of normal modes) need to be computed and this must be done by resolving all spatial scales to obtain an accurate estimate of turbulent dissipation. For the former, a qualitative leap forward is needed because current 1D formalisms based on the mixinglength theory all have their shortcomings, including issues related to free parameters and questionable physical assumptions (see Houdek & Dupret 2015, for details), which prevent them assuring that we properly grasp the physics of the problem.
Acknowledgments
F. Kupka is grateful to the Austrian Science Fund FWF for support through projects P29172N and P33140N and support from European Research Council (ERC) Synergy Grant WHOLESUN #810218.
References
 Ball, W. H., & Gizon, L. 2014, A&A, 568, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ball, W. H., Beeck, B., Cameron, R. H., & Gizon, L. 2016, A&A, 592, A159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belkacem, K., Kupka, F., Samadi, R., & GrimmStrele, H. 2019, A&A, 625, A20 [EDP Sciences] [Google Scholar]
 Broomhall, A. M., Chaplin, W. J., Davies, G. R., et al. 2009, MNRAS, 396, L100 [Google Scholar]
 Brown, T. M. 1984, Science, 226, 687 [Google Scholar]
 Cambon, C., Coleman, G. N., & Mansour, N. N. 1993, J. Fluid Mech., 257, 641 [Google Scholar]
 Canuto, V. M. 1992, ApJ, 392, 218 [Google Scholar]
 Canuto, V. M. 1997, ApJ, 482, 827 [Google Scholar]
 Canuto, V. M. 2009, in Turbulence in Astrophysical and Geophysical Flows, eds. W. Hillebrandt, & F. Kupka, 756, 107 [Google Scholar]
 Canuto, V. M. 2011, A&A, 528, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 ChristensenDalsgaard, J., & Thompson, M. J. 1997, MNRAS, 284, 527 [NASA ADS] [Google Scholar]
 Coleman, G. N., & Mansour, N. N. 1991, Phys. Fluids A, 3, 2255 [Google Scholar]
 Davies, G. R., Broomhall, A. M., Chaplin, W. J., Elsworth, Y., & Hale, S. J. 2014, MNRAS, 439, 2025 [Google Scholar]
 Dziembowski, W. A., Paterno, L., & Ventura, R. 1988, A&A, 200, 213 [NASA ADS] [Google Scholar]
 Favre, A. 1969, SIAM, Philadelphia, 231 [Google Scholar]
 Gatski, T., & Bonnet, J. 2013, Compressibility, Turbulence and High Speed Flow (Elsevier Science & Technology Books) [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse, 15 [Google Scholar]
 Grigahcène, A., Dupret, M.A., Gabriel, M., Garrido, R., & Scuflaire, R. 2005, A&A, 434, 1055 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Houdek, G., & Dupret, M.A. 2015, Liv. Rev. Sol. Phys., 12, 8 [Google Scholar]
 Houdek, G., Trampedach, R., Aarslev, M. J., & ChristensenDalsgaard, J. 2017, MNRAS, 464, L124 [Google Scholar]
 Houdek, G., Lund, M. N., Trampedach, R., et al. 2019, MNRAS, 487, 595 [Google Scholar]
 Hunt, J. C. R., & Carruthers, D. J. 1990, J. Fluid Mech., 212, 497 [Google Scholar]
 Jørgensen, A. C. S., & Angelou, G. C. 2019, MNRAS, 490, 2890 [Google Scholar]
 Jørgensen, A. C. S., & Weiss, A. 2019, MNRAS, 488, 3463 [Google Scholar]
 Jørgensen, A. C. S., Weiss, A., Mosumgaard, J. R., Silva Aguirre, V., & Sahlholdt, C. L. 2017, MNRAS, 472, 3264 [Google Scholar]
 Jørgensen, A. C. S., Mosumgaard, J. R., Weiss, A., Silva Aguirre, V., & ChristensenDalsgaard, J. 2018, MNRAS, 481, L35 [Google Scholar]
 Jørgensen, A. C. S., Weiss, A., Angelou, G., & Silva Aguirre, V. 2019, MNRAS, 484, 5551 [Google Scholar]
 Jørgensen, A. C. S., Montalbán, J., Angelou, G. C., et al. 2021, MNRAS, 500, 4277 [Google Scholar]
 Kjeldsen, H., Bedding, T. R., & ChristensenDalsgaard, J. 2008, ApJ, 683, L175 [Google Scholar]
 Kupka, F., & Muthsam, H. J. 2007, in Convection in Astrophysics, eds. F. Kupka, I. Roxburgh, & K. L. Chan, IAU Symp., 239, 86 [Google Scholar]
 Kupka, F., & Muthsam, H. J. 2017, Liv. Rev. Comput. Astrophys., 3, 1 [Google Scholar]
 Kupka, F., Fabbian, D., Zaussinger, F., & Gizon, L. 2020, Phys. Rev. Lett., submitted [Google Scholar]
 Ledoux, P., & Walraven, T. 1958, Handbuch der Physik, 51, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Magic, Z., & Weiss, A. 2016, A&A, 592, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Manchon, L., Belkacem, K., Samadi, R., et al. 2018, A&A, 620, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosumgaard, J. R., Jørgensen, A. C. S., Weiss, A., Silva Aguirre, V., & ChristensenDalsgaard, J. 2020, MNRAS, 491, 1160 [Google Scholar]
 Muthsam, H. J., Goeb, W., Kupka, F., Liebich, W., & Zoechling, J. 1995, A&A, 293, 127 [Google Scholar]
 Muthsam, H. J., Göb, W., Kupka, F., & Liebich, W. 1999, New Astron., 4, 405 [Google Scholar]
 Nordlund, Å., & Stein, R. F. 2001, ApJ, 546, 576 [Google Scholar]
 Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2 [Google Scholar]
 Piau, L., Collet, R., Stein, R. F., et al. 2014, MNRAS, 437, 164 [Google Scholar]
 Pope, S. B. 2000, Turbulent Flows (Cambridge University Press) [Google Scholar]
 Rosenthal, C. S., ChristensenDalsgaard, J., Houdek, G., et al. 1995, in Helioseismology, ESA Spec. Publ., 376, 459 [Google Scholar]
 Rosenthal, C. S., ChristensenDalsgaard, J., Nordlund, Å., Stein, R. F., & Trampedach, R. 1999, A&A, 351, 689 [NASA ADS] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sarkar, S. 1992, Phys. Fluids A, 4, 2674 [Google Scholar]
 Schou, J., & Birch, A. C. 2020, A&A, 638, A51 [CrossRef] [EDP Sciences] [Google Scholar]
 Sonoi, T., Samadi, R., Belkacem, K., et al. 2015, A&A, 583, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sonoi, T., Belkacem, K., Dupret, M. A., et al. 2017, A&A, 600, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Trampedach, R., Aarslev, M. J., Houdek, G., et al. 2017, MNRAS, 466, L43 [Google Scholar]
 Wilcox, D. 2006, Turbulence Modeling for CFD (DCW Industries), 1 [Google Scholar]
Appendix A: The gasΓ_{1} (GGM) and reducedΓ_{1} (RGM) empirical models
When turbulent pressure is included in the mean stratification for computing the classical adiabatic oscillations, it is necessary to determine the perturbation of turbulent pressure and, subsequently, the perturbation of total pressure (e.g. Rosenthal et al. 1999). To illustrate this scenario, let us begin by considering the mass and momentum conservation equations, which read:
where P_{g} is the gas pressure, g_{i} is the ith component of the gravitational acceleration, τ_{jk} is the viscous stress tensor, ρ is the density, and u_{i} the ith component of the velocity field. We also employ the notation ∂_{t} ≡ ∂/∂t and ∂_{i} ≡ ∂/∂x_{i} as well as Einstein’s notation for repeated indices.
Equations (A.1) and (A.2) are averaged using both a classical Reynolds average and a densityweighted average (also commonly named to as Favre average, see Favre 1969). For a quantity X, the Reynolds average is defined as:
and the Favre average is defined, for a quantity Y, by
As we consider a compressible flow, using a Reynolds average for density and gas pressure and a Favre average for the velocity field greatly simplifies the averaged equations (e.g. Canuto 1997; Nordlund & Stein 2001; Belkacem et al. 2019). Applying this procedure for Eqs. (A.1) and (A.2) gives
where are the Reynolds stresses (where is the kth component of the velocity fluctuation around its Favre average), the pseudoLagrangian derivative is defined by . To derive Eqs. (A.5) and (A.6), gravity fluctuations have been neglected and it has been assumed that viscosity does not affect the mean momentum equation and thus the largescale flow (). It is important to notice that this approximation does not mean that viscous dissipation is neglected, because it appears in the equations governing the turbulent quantities (Canuto 1997).
We now identify the Reynolds average with the horizontal average so that in Eq. (A.6) only the rr component of the Reynolds stress remains. Turbulent pressure is thus defined by . In addition, we split the mean quantities such that (where is the timeaverage of and δX is the pseudoLagrangian perturbation corresponding to the radial oscillations). Therefore, from Eqs. (A.5) and (A.6), we get the desired oscillation equations in the pseudoLagrangian frame:
where σ is the angular frequency, ξ_{r} is the radial component of the eigendisplacement () and δP_{tot} = δP_{g} + δP_{t} is the Lagrangian perturbation of total pressure. In addition, in the adiabatic limit, we also have the thermodynamic relation:
which has been derived by using a number of approximations that are explicitly stated in Appendix B. Therefore, except for the boundary conditions, the system is not closed because we must specify δP_{t} or equivalently δP_{tot}. We note that in the following the subscript ‘0’ for denoting the temporally and horizontally averaged quantities will be dropped for ease of notation.
To solve Eqs. (A.7) and (A.8), we need to specify a prescription for the perturbation of total pressure. It is thus necessary to express the perturbation of turbulent pressure with the perturbation of density, so that:
where 𝒜 is to be determined. Therefore, using Eq. (A.10), we formally write:
Two empirical models have previously been introduced by Rosenthal et al. (1999). The first is the RGM, which consists of the argument that δP_{t} = 0 or 𝒜 = 0, so that Eq. (A.11) becomes:
is called the reduced Γ_{1}. This approximation was initially introduced by Rosenthal et al. (1995) based on the observation that some nonlocal mixinglength theory shows that density and gas pressure perturbations are almost in phase quadrature with the perturbation of turbulent pressure. The authors, therefore, considered that the real part of δP_{t} can be neglected. The second model is the GGM and has been introduced by Rosenthal et al. (1999). It consists of arguing the opposite, namely, that the perturbation of gas pressure and turbulent pressure are in phase. Hence, using 𝒜 = Γ_{1} together with Eqs. (A.9), (A.10), and (A.11), we obtain:
As can be seen with Eqs. (A.12) and (A.13), both the RGM and GGM present the major advantage of being easily implemented for computing adiabatic oscillation while including the effect of turbulent pressure, provided that the mean turbulent pressure is prescribed. For GGM, it is sufficient to replace the gas pressure by the total pressure in the classical adiabatic oscillation equations while for the RGM, it is also necessary to replace Γ_{1} by the reduced given by Eq. (A.12).
Appendix B: Relation between perturbation of gas pressure and density
We follow the derivation proposed by Rosenthal et al. (1999), which starts with the equation governing specific entropy. It reads:
where D/Dt ≡ ∂/∂t + u_{i}∂_{i}, ρ is the density, T is the temperature, s is the specific entropy, is the kcomponent of the radiative flux, and τ_{ik} is the viscous stress tensor.
Equation (B.1) can be recast in terms of gas pressure and density by using the thermodynamic identity:
where P_{g} is the gas pressure, (Γ_{3} − 1)≡(∂lnT/∂lnρ)_{s}. Then, using Eq. (B.2) together with the conservation of mass, after averaging, Eq. (B.1) finally becomes:
which is strictly equivalent to Eq. (11) of Rosenthal et al. (1999). Following the same authors, it is assumed that timevarying parts (i.e. the Lagrangian perturbation) of the last four terms of Eq. (B.3) are zero. Consequently, the perturbation of Eq. (B.3) permits us to recover Eq. (A.9), that is
which is the classical relation for adiabatic oscillations. We note that we assumed , which is quite an accurate approximation, as verified by our solar numerical 3D simulation.
The assumptions made to derive Eq. (B.4) from Eq. (B.3), which consists of assuming that the Lagrangian perturbation of the last four terms of Eq. (B.3) are zero, require further discussion. Indeed, as recognized by Rosenthal et al. (1999), those approximations are quite radical. It is nevertheless useful to go a step further and to explain the underlying physical assumptions. To that end, let us recast Eq. (B.3) in the following form:
where, for sake of simplicity and without loss of meaning, we assumed that thermodynamic quantities are timeindependent. To recover Eq. (B.4) from Eq. (B.5), the Lagrangian perturbations of the last six terms must be neglected. More precisely:

The perturbation of the buoyancy work (), which also appears as a source of turbulent kinetic energy and thus a sink of thermal energy, is considered to be null.

The perturbation of the divergence of the terms and are set to zero. For the former, considering a perfect gas, it is proportional to the convective (enthalpy) flux because of the relation:
where ℛ is the ideal gas constant. Concerning , it is also essentially proportional to the convective flux as shown by Canuto (1997) using a polytropic relation. Therefore, we can conclude that neglecting the perturbation of those two terms is equivalent to neglecting the perturbation of the convective heat flux.

The perturbation of the pressurestrain rate, , is neglected. This can be justified as it scales as the squared turbulent Mach number of the rate of dissipation of turbulent kinetic energy into heat (Sarkar 1992). Hence, it can be neglected because we are considering turbulent flow at relatively low turbulent Mach numbers.

The perturbations of the radiative flux as well as the dissipation rate of turbulent energy into heat (i.e. ) are also considered to be negligible.
Appendix C: Averaged equation for specific entropy
Furthermore, we now look towards inferring the averaged entropy equation. To do so, rather than starting directly with the entropy equation, we consider the equation governing the enthalpy. It reads:
where D/Dt ≡ ∂/∂t + u_{i}∂_{i} and is the kcomponent of the radiative flux. Now, averaging Eq. (C.1) leads to:
where is the kcomponent of the convective flux. Then using the thermodynamic relation, , valid at leading order, and identifying the Reynolds average to the horizontal average gives:
Neglecting the pressurestrain rate as for the equation governing turbulent pressure and adopting the same notations, we finally obtain:
Perturbing Eq. (C.4) thus permits us to derive an expression for the perturbation of the buoyancy work. It reads:
Using the same approximations as described in Sect. 3.1 to derive Eq. (10), immediately leads to δs = 0.
All Figures
Fig. 1.
Frequency difference for radial solar modes between the observed frequencies, as inferred by Broomhall et al. (2009) and Davies et al. (2014), and the modelled frequencies computed using the GGM model, the RGM model, and the model of the present work (see text for details). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.