skip to main content

Title: Limits on Mode Coherence Due to a Non-static Convection Zone
The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of DAV convection zones is very sensitive to temperature (i.e., MCZ∝T−90eff) leading to the possibility that even "small amplitude" pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is slightly out of phase with that required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. We compute the size of this effect for specific examples and discuss the results in the context of Kepler and K2 observations.
; ;
Award ID(s):
Publication Date:
Journal Name:
21st European Workshop on White Dwarfs
Sponsoring Org:
National Science Foundation
More Like this
  1. We consider high-order strongly nonlinear long wave models expanded in a single small parameter measuring the ratio of the water depth to the characteristic wavelength. By examining its dispersion relation, the high-order system for the bottom velocity is found stable to all disturbances at any order of approximation. On the other hand, systems for other velocities can be unstable and even ill-posed, as signified by the unbounded maximum growth. Under the steady assumption, a new third-order solitary wave solution of the Euler equations is obtained using the high-order strongly nonlinear system and is expanded in an amplitude parameter, which is different from that used in weakly nonlinear theory. The third-order solution is shown to well describe various physical quantities induced by a finite-amplitude solitary wave, including the wave profile, horizontal velocity profile, particle velocity at the crest and bottom pressure. For numerical computations, the first- and second-order strongly nonlinear systems for the bottom velocity are considered. It is shown that finite difference schemes are unstable due to truncation errors introduced in approximating high-order spatial derivatives and, therefore, a more accurate spatial discretization scheme is necessary. Using a pseudo-spectral method based on finite Fourier series combined with an iterative scheme formore »the inversion of a non-local operator, the strongly nonlinear systems are solved numerically for the propagation of a single solitary wave and the head-on collision of two counter-propagating solitary waves of finite amplitudes, and the results are compared with previous laboratory measurements.« less
  2. We consider a strongly nonlinear long wave model for large amplitude internal waves in a three-layer flow between two rigid boundaries. The model extends the two-layer Miyata–Choi–Camassa (MCC) model (Miyata, Proceedings of the IUTAM Symposium on Nonlinear Water Waves , eds. H. Horikawa & H. Maruo, 1988, pp. 399–406; Choi & Camassa, J. Fluid Mech. , vol. 396, 1999, pp. 1–36) and is able to describe the propagation of long internal waves of both the first and second baroclinic modes. Solitary-wave solutions of the model are shown to be governed by a Hamiltonian system with two degrees of freedom. Emphasis is given to the solitary waves of the second baroclinic mode (mode 2) and their strongly nonlinear characteristics that fail to be captured by weakly nonlinear models. In certain asymptotic limits relevant to oceanic applications and previous laboratory experiments, it is shown that large amplitude mode-2 waves with single-hump profiles can be described by the solitary-wave solutions of the MCC model, originally developed for mode-1 waves in a two-layer system. In other cases, however, e.g. when the density stratification is weak and the density transition layer is thin, the richness of the dynamical system with two degrees of freedom becomes apparent and new classesmore »of mode-2 solitary-wave solutions of large amplitudes, characterized by multi-humped wave profiles, can be found. In contrast with the classical solitary-wave solutions described by the MCC equation, such multi-humped solutions cannot exist for a continuum set of wave speeds for a given layer configuration. Our analytical predictions based on asymptotic theory are then corroborated by a numerical study of the original Hamiltonian system.« less
  3. For Rayleigh-Bénard-Poiseuille flows, thermal stratification resulting from a wall-normal temperature gradient together with an opposing gravitational field can lead to buoyancy-driven instability. Moreover, for sufficiently large Reynolds numbers, viscosity-driven instability can occur. Two higher-order-accurate methods based on the full and linearized Navier-Stokes equations were developed for investigating the temporal stability of such flows. The new methods employ a spectral discretization in the homogeneous directions. In the wall-normal direction, the convective and viscous terms are discretized with fifth-order-accurate biased and fourth-order-accurate central compact finite differences. A fourth-order-accurate explicit Runge-Kutta method is employed for time integration. To validate the methods, the primary instability was investigated for different combinations of the Reynolds and Rayleigh number. The results from these primary stability investigations are consistent with linear stability theory results from the literature with respect to both the onset of the instability and the dependence of the temporal growth rate on the wave angle. For the cases with buoyancy-driven instability, strong linear growth is observed for a broad range of spanwise wavenumbers. The largest growth rates are obtained for a wave angle of 90deg. For the cases with viscosity-driven instability, the linear growth rates are lower and the first mode to experience nonlinear growthmore »is a higher harmonic with half the wavelength of the fundamental.« less
  4. Context. Before reaching their quiescent terminal white-dwarf cooling branch, some low-mass helium-core white dwarf stellar models experience a number of nuclear flashes which greatly reduce their hydrogen envelopes. Just before the occurrence of each flash, stable hydrogen burning may be able to drive global pulsations that could be relevant in shedding some light on the internal structure of these stars through asteroseismology, similarly to what occurs with other classes of pulsating white dwarfs. Aims. We present a pulsational stability analysis applied to low-mass helium-core stars on their early white-dwarf cooling branches going through CNO flashes in order to study the possibility that the ε mechanism is able to excite gravity-mode pulsations. We assess the ranges of unstable periods and the corresponding instability domain in the log g  −  T eff plane. Methods. We carried out a nonadiabatic pulsation analysis for low-mass helium-core white-dwarf models with stellar masses between 0.2025 and 0.3630  M ⊙ going through CNO flashes during their early cooling phases. Results. We found that the ε mechanism due to stable hydrogen burning can excite low-order ( ℓ  = 1, 2) gravity modes with periods between ∼80 and 500 s for stars with 0.2025 ≲  M ⋆ / M ⊙  ≲ 0.3630 locatedmore »in an extended region of the log g  −  T eff diagram, with effective temperature and surface gravity in the ranges 15 000 ≲  T eff  ≲ 38 000 K and 5.8 ≲ log g  ≲ 7.1, respectively. For the sequences that experience multiple CNO flashes, we found that with every consecutive flash, the region of instability becomes wider and the modes are more strongly excited. The magnitudes of the rate of period change for these modes are in the range of ∼10 −10 –10 −11  [s/s]. Conclusions. Since the timescales required for these modes to reach amplitudes large enough to be observable are shorter than their corresponding evolutionary timescales, the detection of pulsations in these stars is feasible. Given the current problems in distinguishing some stars that populate the same region of the log g  −  T eff plane, the eventual detection of short-period pulsations may help in the classification of such stars. Furthermore, if a low-mass white dwarf star were found to pulsate with low-order gravity modes in this region of instability, it would confirm our result that such pulsations can be driven by the ε mechanism. In addition, confirming a rapid rate of period change in these pulsations would support the idea that these stars actually experience CNO flashes, as has been predicted by evolutionary calculations.« less
  5. We examine the linear stability of a homogeneous gas–solid suspension of small Stokes number particles, with a moderate mass loading, subject to a simple shear flow. The modulation of the gravitational force exerted on the suspension, due to preferential concentration of particles in regions of low vorticity, in response to an imposed velocity perturbation, can lead to an algebraic instability. Since the fastest growing modes have wavelengths small compared with the characteristic length scale ( $U_{g}/{\it\Gamma}$ ) and oscillate with frequencies large compared with ${\it\Gamma}$ , $U_{g}$ being the settling velocity and ${\it\Gamma}$ the shear rate, we apply the WKB method, a multiple scale technique. This analysis reveals the existence of a number density mode which travels due to the settling of the particles and a momentum mode which travels due to the cross-streamline momentum transport caused by settling. These modes are coupled at a turning point which occurs when the wavevector is nearly horizontal and the most amplified perturbations are those in which a momentum wave upstream of the turning point creates a downstream number density wave. The particle number density perturbations reach a finite, but large amplitude that persists after the wave becomes aligned with the velocity gradient.more »The growth of the amplitude of particle concentration and fluid velocity disturbances is characterised as a function of the wavenumber and Reynolds number ( $\mathit{Re}=U_{g}^{2}/{\it\Gamma}{\it\nu}$ ) using both asymptotic theory and a numerical solution of the linearised equations.« less