Variability and mass-loss are common phenomena in massive OB-type stars. It is argued that they are caused by violent strange mode instabilities identified in corresponding stellar models. We present a systematic linear stability analysis with respect to radial perturbations of massive OB-type stars with solar chemical composition and masses between 23 to 100 M$_{\odot}$. For selected unstable stellar models, we perform non-linear simulations of the evolution of the instabilities into the non-linear regime. Finite amplitude pulsations with periods in the range between hours and 100 days are found to be the final result of the instabilities. The pulsations are associated with a mean acoustic luminosity which can be the origin of a pulsationally driven wind. Corresponding mass-loss rates lie in the range between 10$^{-9}$ and 10$^{-4}$ M$_{\odot}$ yr$^{-1}$ and may thus affect the evolution of massive stars.

A linear stability analysis of models for evolved primordial stars with masses between 150 and 250 M$_{\odot}$ is presented. Strange mode instabilities with growth rates in the dynamical range are identified for stellar models with effective temperatures below log T$_{\rm{eff}}$ = 4.5. For selected models the final fate of the instabilities is determined by numerical simulation of their evolution into the non-linear regime. As a result, the instabilities lead to finite amplitude pulsations. Associated with them are acoustic energy fluxes capable of driving stellar winds with mass-loss rates in the range between 7.7 $\times$ 10$^{-7}$ and 3.5 $\times$ 10$^{-4}$ M$_{\odot}$ yr$^{-1}$.

Numerical simulations of the evolution of strange-mode instabilities into the non-linear regime have been performed for a wide range of stellar parameters for Wolf-Rayet stars. It has been shown that the Wolf-Rayet models reach radial velocities which amount to up to 30% of their escape velocity. The acoustic luminosities suggest a connection to the observed mass loss. Most of the models show a jump in the mean effective temperature after reaching the non-linear regime. This jump is related to the run of the opacity.

We analyse the presence of nonradial oscillations in Cepheids, a problem which has not been theoretically revised since the work of Dziembowsky (1977) and Osaki (1977). Our analysis is motivated by a work of Moskalik et al. (2004) which reports the detec tion of low amplitude periodicities in a few Cepheids of the large Magellanic cloud. These newly discovered periodicities were interpreted as nonradial modes.} {Based on linear nonadiabatic stability analysis, our goal is to reanalyse the presence and stability of nonradial modes, taking into account improvement in the main input phys ics required for the modelling of Cepheids.} {We compare the results obtained from two different numerical methods used to solve the set of differential equations: a matrix method and the Ricatti method.} {We show the limitation of the matrix method to find low order p-modes ($l<6$), because of their dual character in evolved stars such as Cepheids. For higher order p-modes, we find an excellent agreement between the two methods.} {No nonradial instability is found below $l=5$, whereas many unstable nonradial modes exist for higher orders. We also find that nonradial modes remain unstable, even at hotter effective temperatures than the blue edge of the Cepheid instability strip, where no radial pulsations are expected.

Strange mode instabilities in the envelopes of massive stars lead to shock waves, which can oscillate on a much shorter timescale than that associated with the primary instability. The phenomenon is studied by direct numerical simulation using a, with respect to time, implicit Lagrangian scheme, which allows for the variation by several orders of magnitude of the dependent variables. The timestep for the simulation of the system is reduced appreciably by the shock oscillations and prevents its long term study. A procedure based on domain decomposition is proposed to surmount the difficulty of vastly different timescales in various regions of the stellar envelope and thus to enable the desired long term simulations. Criteria for domain decomposition are derived and the proper treatment of the resulting inner boundaries is discussed. Tests of the approach are presented and its viability is demonstrated by application to a model for the star P Cygni. In this investigation primarily the feasibility of domain decomposition for the problem considered is studied. We intend to use the results as the basis of an extension to two dimensional simulations.

The evolution of strange mode instabilities into the non linear regime has been followed by numerical simulation for an envelope model of a massive star having solar chemical composition, M=50M_sun, T_eff=10^4K and L=1.17*10^6 L_sun. Contrary to previously studied models, for these parameters shocks are captured in the H-ionisation zone and perform rapid oscillations within the latter. A linear stability analysis is performed to verify that this behaviour is physical. The origin of an instability discovered in this way is identified by construction of an analytical model. As a result, the stratification turns out to be essential for instability. The difference to common stratification instabilities, e.g., convective instabilities, is discussed.