Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the computation of eigenfrequencies for equilibrium models including turbulent pressure
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The space-borne missions have provided a wealth of highly accurate data. However, our inability to properly model the upper-most region of solar-like stars prevents us from making the best of these observations. This problem is called surface effect and a key ingredient to solve it is turbulent pressure for the computation of both the equilibrium models and the oscillations. While 3D hydrodynamic simulations help to include properly the turbulent pressure in the equilibrium models, the way this surface effect is included in the computation of stellar oscillations is still subject to uncertainties. We aim at determining how to properly include the effect of turbulent pressure and its Lagrangian perturbation in the adiabatic computation of the oscillations. We also discuss the validity of the gas-gamma model (GGM) and reduced gamma model (RGM) approximations, which have been used to compute adiabatic oscillations of equilibrium models including turbulent pressure. We use a patched model of the Sun with an inner part constructed by a 1D stellar evolution code (CESTAM) and an outer part by the 3D hydrodynamical code (CO$^5$BOLD). Then, the adiabatic oscillations are computed using the ADIPLS code for the GGM and RGM and with the MAD code imposing the adiabatic condition on an existing time-dependent convection (TDC) formalism. We show that the computation of the oscillations using the TDC formalism in the adiabatic limit improves significantly the agreement with the observed frequencies compared to the GGM and RGM. Of the components of the turbulent pressure perturbation, the perturbation of the density and advection term is found to contribute most to the frequency shift. We propose a formalism to evaluate the frequency shift due to the inclusion of the term with the turbulent pressure perturbation in the variational principle in order to extrapolate our result to other stars.
rate research
Read More
We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the structure models. The eigenfrequencies obtained are compared to those of 1-D models relying on the usual phenomenologies of convection and to observations of the MDI instrument aboard SoHO. We also investigate how the magnetic activity could change the eigenfrequencies and the solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects a 1.2 kG horizontal field. All models and observed eigenfrequencies are fairly close below 3 mHz. Above 3 mHz the eigenfrequencies of the phenomenological convection models are above the observed eigenfrequencies. The frequencies of the models based on the 3-D simulations are slightly below the observed frequencies. Their maximum deviation is ~ 3 mu Hz at 3 mHz but drops below 1 mu Hz at 4 mHz. Replacing the hydrodynamic by the magnetohydrodynamic simulation increases the eigenfrequencies. The shift is negligible below 2.2 mHz and then increases linearly with frequency to reach ~ 1.7 mu Hz at 4 mHz. The impact of the simulated activity is a 14 milliarcsecond shrinking of the solar layers near the optical depth unity.
The use of the full potential of stellar seismology is made difficult by the improper modeling of the upper-most layers of solar-like stars and their influence on the modeled frequencies. Our knowledge on these emph{surface effects} has improved thanks to the use of 3D hydrodynamical simulations but the calculation of eigenfrequencies relies on empirical models for the description of the Lagrangian perturbation of turbulent pressure: the reduced-$Gamma_1$ model (RGM) and the gas-$Gamma_1$ model (GGM). Starting from the fully compressible turbulence equations, we derive both the GGM and RGM models using a closure to model the flux of turbulent kinetic energy. It is found 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. It is also demonstrated that they are both compatible with the adiabatic approximation but also imply a number of questionable assumptions mainly regarding mode physics. Among others hypothesis, one has to neglect the Lagrangian perturbation of the dissipation of turbulent kinetic energy into heat and the Lagrangian perturbation of buoyancy work.
The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger.
We have constructed models for the sun at three stages of its evolution: a zero-age main sequence model, the present sun, and a subgiant model. For each model, the turbulent pressure and turbulent kinetic energy were calculated from 3-d radiative hydrodynamical simulations (described in the poster by Robinson et al.), and inserted into the 1-d stellar models. We note that in these simulations, the turbulent pressure is not a free parameter, but can be computed from the resulting velocity field. We show the calculated p-mode frequencies for the model of the present sun, with and without turbulent pressure, and compare them to the observed solar frequencies. When the turbulent pressure is included in the models, the calculated frequencies are brought closer to the observed frequencies in the sun by up to two $mu Hz$, strictly from structural effects. The effect of including turbulent pressure on p-mode frequencies is also shown for the zero-age main sequence model. Our models also suggest that the importance of turbulent pressure increases as the star evolves into the subgiant region. We discuss the importance of also including realistic turbulence as well as radiation in the non-adiabatic calculation of oscillation frequencies.
We present models of turbulent mixing at the boundaries between hot (T~10^{6-7} K) and warm material (T~10^4 K) in the interstellar medium, using a three-dimensional magnetohydrodynamical code, with radiative cooling. The source of turbulence in our simulations is a Kelvin-Helmholtz instability, produced by shear between the two media. We found, that because the growth rate of the large scale modes in the instability is rather slow, it takes a significant amount of time (~1 Myr) for turbulence to produce effective mixing. We find that the total column densities of the highly ionized species (C IV, N V, and O VI) per interface (assuming ionization equilibrium) are similar to previous steady-state non-equilibrium ionization models, but grow slowly from log N ~10^{11} to a few 10^{12} cm^{-2} as the interface evolves. However, the column density ratios can differ significantly from previous estimates, with an order of magnitude variation in N(C IV)/N(O VI) as the mixing develops.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
