On the basis of consistent stellar evolution and nonlinear stellar pulsation calculations the Cepheid V1033 Cyg is shown to be the post--main sequence star at the first crossing of the instability strip during gravitational contraction of the helium core. The observed light variability of V1033 Cyg is due to radial oscillations in the fundamental mode. The best agreement (within one percent) between recent observations and the theoretical estimate of the period change rate was obtained for the evolutionary sequence with stellar mass $M=6.3M_\odot$ and helium and heavier element fractional abundances $Y=0.28$ and $Z=0.022$, respectively. The age of the star, the luminosity, the radius, the effective temperature and the surface gravity are $t_\mathrm{ev}=5.84\times 10^7$ yr, $L=2009L_\odot$, $R=45.6R_\odot$, $T_\mathrm{eff}=5726 \mathrm{K}$, $\log g=1.92$.

The evolutionary calculations for population I stars with masses on the main sequence 5 M_\odot <= M_0 <= 6.1 M_\odot and initial fractional abundances of helium Y_0=0.28 and heavier elements Z_0=0.02 were carried out to the stage of central helium exhaustion. Selected core helium--burning models were used as initial conditions for solution of the equations of hydrodynamics and time--dependent convection describing radial pulsations of Cepheids. In the Hertzsprung--Russel diagram the evolutionary tracks are shown to cross the red edge of the instability strip for M_0 > 5.1M_\odot. The grid of hydrodynamic models of core helium--burning Cepheids on the stages of the second and the third crossings of the instability strip was computed. The pulsation period \Pi and the rate of period change \dot\Pi were determined as a continuous function of star age for each evolutionary sequence of first--overtone pulsators. Results of calculations agree with observational estimates of \dot\Pi recently obtained for the short--period Cepheids V532 Cyg, BG Cru and RT Aur.

Stellar evolution calculations to the stage of the cooling white dwarf were done for population~I stars with masses on the main sequence $1M_\odot\le M_0\le 1.5M_\odot$. The final helium flash LTP is shown to occur in post--AGB stars with initial masses $1.3M_\odot\le M_0\le 1.32M_\odot$ for the overshooting parameter $f=0.016$. In the case of more effective overshooting ($f=0.018$) the final helium flash occurs at initial masses $1.28M_\odot\le M_0\le 1.3M_\odot$. Fivefold variations of the parameter responsible for the mass loss rate during the post--AGB stage do not affect occurrence of the final helium flash but lead to perceptible changes of the evolutionary time. Selected models of two evolutionary sequences with initial mass $M_0 = 1.3M_\odot$ computed with overshooting parameters $f=0.016$ and $f=0.018$ were used as initial conditions in solution of the equations of hydrodynamics describing radial oscillations of stars on the stage of the final helium flash at effective temperatures less than $10^4$ K. The maximum pulsation period $\Pi=117$ day determined for the evolutionary sequence $M_0=1.3M_\odot$, $f=0.016$ is in a good agreement with observational estimates of the period of FG Sge. The mass, the radius and the effective temperature of the star are $M=0.565M_\odot$, $R=126R_\odot$ and $T_\mathrm{eff}=4445$ K, respectively. At the same time the average period change rate of FG Sge from 1960 to 1990 is nearly three time larger than its theoretical estimate.

Evolutionary calculations of population I stars with initial masses $M_0=1M_\odot$, $1.5M_\odot$ and $2M_\odot$ were carried out up to the stage of the proto--planetary nebula. Selected models of post--AGB evolutionary sequences with effective temperatures $3.6\times 10^3\,\mathrm{K}\lesssim T_\mathrm{eff}\lesssim 2\times 10^4\,\mathrm{K}$ were used as initial conditions in calculations of self--escited stellar oscillations. For the first time the sequences of hydrodynamic models of radially pulsating post--AGB stars were computed using the self--consistent solution of the equations of radiation hydrodynamics and time--dependent convection. Within this range of effective temperatures the post--AGB stars are the fundamental mode pulsators with period decreasing as the star evolves from $\Pi\approx 300$ day to several days. Period fluctuations are due to nonlinear effects and are most prominent at effective temperatures $T_\mathrm{eff} < 5000$K. The amplitude of bolometric light variations is $\Delta M_\mathrm{bol}\approx 1$ at $T_\mathrm{eff} \lesssim 6000$K and rapidly decreases with increasing $T_\mathrm{eff}$. The theoretical dependence of the pulsation period as a function of effective temperature obtained in the study can be used as a criterion for the evolutionary status of pulsating variables suspected to be post--AGB stars.

Evolutionary calculations of population II stars with chemical composition of the globular cluster M3 were carried out under various assumptions about the initial stellar mass ($0.809M_\odot\le M_\mathrm{ZAMS} \le 0.83M_\odot$) and the mass loss rate parameter in the Reimers formula ($0.45\le\eta_\mathrm{R}\le 0.55$). In general, 30 evolutionary tracks of the horizontal branch stars were computed. Selected models of evolutionary sequences were used as initial conditions for solution of the equations of hydrodynamics that describe radial stellar oscillations. Hydrodynamic models of RR Lyr type stars were computed for the core helium burning stage as well as for the preceding pre--ZAHB stage. Analytic relations for the effective temperature of the instability strip edges as a function of stellar luminosity are obtained. Theoretical histograms of the period distribution of RR Lyr type variables were produced for each evolutionary sequence using Monte--Carlo simulations based on the consistent stellar evolution and nonlinear stellar pulsation calculations. A satisfactory agreement with observations (i.e. the greater number of RRab variables) was found for the evolutionary sequence $M_\mathrm{ZAMS} = 0.811M_\odot$, $\eta_\mathrm{R}=0.55$ with the number fraction of fundamental mode pulsators $\approx 75\%$. At the same time the mean period of fundamental mode pulsators ($\langle\Pi\rangle_0=0.79$ day) is substantially greater compared to the observational estimate of $\langle\Pi\rangle_\mathrm{ab}$.

Hydrodynamic computations of nonlinear Cepheid pulsation models with periods from 20 to 100 day on the evolutionary stage of core helium burning were carried out. Equations of radiation hydrodynamics and time--dependent convection were solved with initial conditions obtained from selected models of evolutionary sequences of population I stars with initial masses from $8M_\odot$ to $12.5M_\odot$. For each crossing of the instability strip the pulsation period $\Pi$ and the rate of period change $\dot\Pi$ were derived as a function of evolutionary time. Comparing results of our computations with observational estimates of $\Pi$ and $\dot\Pi$ we determined fundamental parameters (the age, the mass, the luminosity and the radius) of seven long--period Cepheids. Theoretical estimates of the stellar radius are shown to agree with radius measurements by the Baade--Wesselink technique within 3\% for RS~Pup and GY~Sge whereas for SV~Vul the disagreement between theory and observations does not exceed 10\%.

The grid of evolutionary tracks of population II stars with initial masses from $0.81M_\odot$ to $0.85M_\odot$ and chemical composition of the globular cluster M 3 is computed. Selected models of horizontal branch stars were used as initial conditions for solution of the equations of radiation hydrodynamics and time--dependent convection describing radial stellar oscillations. The boundaries of the instability strip on the Herztsprung--Russel diagram were determined using nearly 100 hydrodynamic models of RR Lyr pulsating variables. For each evolutionary track crossing the instability strip the pulsation period was determined as a function of evolutionary time. The rate of period change of most variables is shown to range within $-0.02\le\dot\Pi\le 0.05~\mathrm{day}/10^6\mathrm{yr}$. Theoretical estimate of the mean period change rate obtained by the population synthesis method is $\langle\dot\Pi\rangle=6.0\times 10^{-3}~\mathrm{day}/10^6\mathrm{yr}$ and agrees well with observations of RR Lyr variables of the globular cluster M 3.

Stellar evolution calculations were carried out from the main sequence to the final AGB stage of stars with initial masses from 1 to $2M_\odot$ and metallicity $Z=0.01$. Selected models of evolutionary sequences were used as initial conditions for solution of the equations of radiation hydrodynamics and time-dependent convection describing radial stellar pulsations. The study was aimed to construct the hydrodynamic models of Mira-type stars that show the secular decrease in the pulsation period commenced in 1970--th at $\Pi=315$ day. We show that such a condition for the period change is satisfied with evolutionary sequences from 1 to $1.2M_\odot$ and the best agreement with observations is obtained for models with initial mass $1.2M_\odot$. The pulsation period reduction is due to both the stellar radius decrease during the thermal pulse of the helium burning shell and mode switch from the fundamental mode to the first overtone. Theoretical estimates of the fundamental parameters of the star at the onset of pulsation period reduction are as follows: the mass is $M=0.93M_\odot$, the luminosity is $L = 4080L_\odot$, and the radius is $R=220R_\odot$. The mode switch occurs 35 years after the onset of period reduction.

We performed hydrodynamic computations of nonlinear stellar pulsations of population I stars at the evolutionary stages of the ascending red giant branch and the following luminosity drop due to the core helium flash. Red giants populating this region of the Hertzsprung--Russel diagram were found to be the fundamental mode pulsators. The pulsation period is the largest at the tip of the red giant branch and for stars with initial masses from 1.1M_\odot to 1.9M_\odot ranges from 254 day to 33 day, respectively. The rate of period change during the core helium flash is comparable with rates of secular period change in Mira type variables during the thermal pulse in the helium shell source. The period change rate is largest (\dot\Pi/\Pi\approx -0.01 yr^{-1}) in stars with initial mass Mzams=1.1M_\odot and decreases to \dot\Pi/\Pi\sim -0.001\ yr^{-1} for stars of the evolutionary sequence Mzams=1.9M_\odot. Theoretical light curves of red giants pulsating with periods Pi > 200 day show the presence of the secondary maximum similar to that observed in many Miras.

Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2M_\odot <= Mzams <= 5M_\odot and metallicity Z=0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods from 30 to 400 day. The rate of period change gradually increases as the star evolves but is too small to be detected (d ln P/dt < 1.e-5 yr^-1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with intial masses Mzams=2M_\odot and 3M_\odot and 6 thermal pulses for Mzams=4M_\odot and 5M_\odot. TP-AGB stars with initial masses Mzams<=3M_\odot pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (Mzams>=4M_\odot) pulsate in the fundamental mode with periods as long as 1.e3 day. Most rapid pulsation period change with rate -0.02 < d ln P/dt < -0.01 yr^-1 occurs during decrease of the surface luminosity after the maximum helium luminosity. The rate of subsequent increase of the period is d ln P/dt <= 5.e-3 yr^-1.

Pulsation period changes in Mira type variables are investigated using the stellar evolution and nonlinear stellar pulsation calculations. We considered the evolutionary sequence of stellar models with initial mass 3 Msol and population I composition. Pulsations of stars in the early stage of the asymptotic giant branch are shown to be due to instability of the fundamental mode. In the later stage of evolution when the helium shell source becomes thermally unstable the stellar oscillations arise in either the fundamental mode (for the stellar luminosity L < 5.4e3 Lsol) or the first overtone (L > 7e3 Lsol). Excitation of pulsations is due to the kappa-mechanism in the hydrogen ionization zone. Stars with intermediate luminosities 5.4e3 < L/Lsol < 7e3 were found to be stable against radial oscillations. The pulsation period was determined as a function of evolutionary time and period change rates \dot\Pi were evaluated for the first ten helium flashes. The period change rate becomes the largest in absolute value (\dot\Pi/\Pi\approx -0.01 per year) between the helium flash and the maximum of the stellar luminosity. Period changes with rate |\dot\Pi/\Pi| > 1e-3 per year take place during approx 500 years, that is nearly one hundredth of the interval between helium flashes.

The study is devoted to radial pulsations in Population I Cepheids with masses from 5.4 to 6 Msun. Solution of the equations of radiation hydrodynamics and turbulent convection for nonlinear stellar pulsations was obtained with initial conditions as the core helium burning models of computed evolutionary sequences. For each value of the initial mass we considered stellar evolution from the zero age main sequence to central helium exhaustion with three values of the convective overshoot parameter: aov=0.1, 0.15 and 0.2. Models for the Cepheid SZ Tau with pulsation period 3.149 day can be constructed only for aov=0.1 and aov=0.15. The star is at the evolutionary stage of the second crossing of the instability strip and pulsates in the first overtone just near the boundary between the fundamental mode and the first overtone. The mass, radius and age of the star are in the ranges 5.46 <= M/Msun <= 5.75, 41.5 <= R/Rsun <= 42.3 and 6.9e7 yr <= t <= 8.0e7 yr, respectively. Predicted period change rates are of -0.4 s/yr.

Hydrodynamic models of short--period Cepheids were computed to determine the pulsation period as a function of evolutionary time during the first and third crossings of the instability strip. The equations of radiation hydrodynamics and turbulent convection for radial stellar pulsations were solved with the initial conditions obtained from the evolutionary models of population I stars (X=0.7, Z=0.02) with masses from 5.2 to 6.5 Msol and the convective core overshooting parameter 0.1 <= aov <= 0.3. In Cepheids with period of 4 d the rate of pulsation period change during the first crossing of the instability strip is over fifty times larger than that during the third crossing. Polaris is shown to cross the instability strip for the first time and to be the fundamental mode pulsator. The best agreement between the predicted and observed rates of period change was obtained for the model with mass of 5.4 Msol and the overshooting parameter aov=0.25. The bolometric luminosity and radius are L = 1.26e3 Lsol and R = 37.5 Rsol, respectively. In the HR diagram Polaris is located at the red edge of the instability strip.

Theoretical estimates of the rates of radial pulsation period change in Galactic Cepheids with initial masses 5.5M_\odot <= Mzams <= 13M_\odot, chemical composition X=0.7, Z=0.02 and periods 1.5 day <= P <= 100 day are obtained from consistent stellar evolution and nonlinear stellar pulsation computations. Pulsational instability was investigated for three crossings of the instability strip by the evolutionary track in the HR diagram. The first crossing occurs at the post-main sequence helium core gravitational contraction stage which proceeds in the Kelvin--Helmholtz timescale whereas the second and the third crossings take place at the evolutionary stage of thermonuclear core helium burning. During each crossing of the instability strip the period of radial pulsations is a quadratic function of the stellar evolution time. Theoretical rates of the pulsation period change agree with observations but the scatter of observational estimates of dP/dt noticeably exceeds the width of the band (\delta\log |dP/dt| <= 0.6) confining evolutionary tracks in the period - period change rate diagram. One of the causes of the large scatter with very high values of dP/dt in Cepheids with increasing periods might be the stars that cross the instability strip for the first time. Their fraction ranges from 2% for Mzams=5.5M_\odot to 9% for Mzams=13M_\odot. Variable stars \alpha UMi and IX Cas seem to belong to such objects.

Theoretical estimates of the pulsation period change rates in LMC Cepheids are obtained from consistent calculation of stellar evolution and nonlinear stellar pulsation for stars with initial chemical composition X=0.7, Z=0.008, initial masses 5M_\odot <= Mzams <= 9M_\odot and pulsation periods ranged from 2.2 to 29 day. The Cepheid hydrodynamical models correspond to the evolutionary stage of thermonuclear core helium burning. During evolution across the instability strip in the HR diagram the pulsation period of Cepheids is the quadratic function of the evolution time for the both fundamental mode and first overtone. Cepheids with initial masses Mzams >= 7M_\odot pulsate in the fundamental mode and the period change rate \dot\Pi varies nearly by a factor of two for both crossings of the instability strip. In the period -- period change rate diagram the values of the period and \dot\Pi concentrate within the strips, their slope and half--width depending on both the direction of the movement in the HR--diagram and the pulsation mode. For oscillations in the fundamental mode the half-widths of the strip are \delta\log\dot\Pi = 0.35 and \delta\log\dot\Pi = 0.2 for the first and the second crossings of the instability strip, respectively. Results of computations are compared with observations of nearly 700 LMC Cepheids. Within existing observational uncertainties of \dot\Pi the theoretical dependences of the period change rate on the pulsation period are in a good agreement with observations.

Excitation of radial pulsations in red supergiants of Magellanic Clouds is investigated using the stellar evolution calculations and the self-consistent solution of the equations of radiation hydrodynamics and turbulent convection. The stars with initial masses 6M_odot<=M_zams<=28 M_odot and the initial chemical composition X=0.7, 0.004<=Z<=0.008 are shown to be unstable against fundamental mode oscillations with periods from 17 to 1200 days as they become helium burning red supergiants. The period-luminosity relation slightly depends on the mass loss rate varying with a factor of three, whereas its dependence on the metal abundance is delta M_bol=0.89 delta log Z. In comparison with galactic red supergiants (Z=0.02) the low metal abundances in red supergiants of Magellanic Clouds are responsible for their higher effective temperatures and substantially narrower ranges of evolutionary stellar radius change during helium burning. Therefore on the period-mass diagram the red supergiants of Magellanic Clouds are located within the strip with width of delta log M = 0.09, so that the uncertainty of mass evaluation of the red supergiant with the known pulsational period is nearly 25%.

Excitation of radial oscillations in population I (X=0.7, Z=0.02) red supergiants is investigated using the solution of the equations of radiation hydrodynamics and turbulent convection. The core helium burning stars with masses 8M_odot <= M <= 20M_odot end effective temperatures T_eff < 4000K are shown to be unstable against radial pulsations in the fundamental mode. The oscillation periods range between 45 and 1180 days. The pulsational instability is due to the kappa-mechanism in the hydrogen and helium ionization zones. Radial pulsations of stars with mass M < 15M_odot are strictly periodic with the light amplitude Delta M_bol <= 0.5 mag. The pulsation amplitude increases with increasing stellar mass and for M > 15M_odot the maximum expansion velocity of outer layers is as high as one third of the escape velocity. The mean radii of outer mass zones increase due to nonlinear oscillations by <= 30% in comparison with the initial equilibrium. The approximate method (with uncertainty of a factor of 1.5) to evaluate the mass of the pulsating red supergiant with the known period of radial oscillations is proposed. The approximation of the pulsation constant Q as a function of the mass-to-radius ratio is given. Masses of seven galactic red supergiants are evaluated using the period-mean density relation.

Instability of population I (X=0.7, Y=0.02) massive stars against radial oscillations during the post-main sequence gravitational contraction of the helium core is investigated. Initial stellar masses are in the range from 65M_\odot to 90M_\odot. In hydrodynamic computations of self-exciting stellar oscillations we assumed that energy transfer in the envelope of the pulsating star is due to radiative heat conduction and convection. The convective heat transfer was treated in the framework of the theory of time-dependent turbulent convection. During evolutionary expansion of outer layers after hydrogen exhaustion in the stellar core the star is shown to be unstable against radial oscillations while its effective temperature is Teff > 6700K for Mzams=65M_\odot and Teff > 7200K for mzams=90M_\odot. Pulsational instability is due to the \kappa-mechanism in helium ionization zones and at lower effective temperature oscillations decay because of significantly increasing convection. The upper limit of the period of radial pulsations on this stage of evolution does not exceed 200 day. Radial oscillations of the hypergiant resume during evolutionary contraction of outer layers when the effective temperature is Teff > 7300K for Mzams=65M_\odot and Teff > 7600K for Mzams=90M_\odot. Initially radial oscillations are due to instability of the first overtone and transition to fundamental mode pulsations takes place at higher effective temperatures (Teff > 7700K for Mzams=65M_\odot and Teff > 8200K for Mzams=90M_\odot). The upper limit of the period of radial oscillations of evolving blueward yellow hypergiants does not exceed 130 day. Thus, yellow hypergiants are stable against radial stellar pulsations during the major part of their evolutionary stage.

In this study we consider the nonlinear radial oscillations exciting in LBV--stars with effective temperatures 1.5e4 K <= Teff <= 3e4 K, bolometric luminosities 1.2e6 L_odot <= L <= 1.9e6 L_odot and masses 35.7 M_odot <= M <= 49.1 M_odot. Hydrodynamic computations were carried out with initial conditions obtained from evolutionary sequences of population I stars (X=0.7, Z=0.02) with initial masses from 70M_odot to 90 M_odot. All hydrodynamical models show instability against radial oscillations with amplitude growth time comparable with dynamical time scale of the star. Radial oscillations exist in the form of nonlinear running waves propagating from the boundary of the compact core to the upper boundary of the hydrodynamical model. The velocity amplitude of outer layers is of several hundreds of km/s while the bolometric light amplitude does not exceed 0.2 mag. Stellar oscillations are not driven by the kappa-mechanism and are due to the instability of the gas with adiabatic exponent close to the critical value Gamma_1 = 4/3 due to the large contribution of radiation in the total pressure. The range of the light variation periods (6 day <= P <= 31 day) of hydrodynamical models agrees with periods of microvariability observed in LBV--stars.

Evolution of Population I stars with initial masses from 70M_\odot to 130M_\odot is considered under various assumptions on the mass loss rate \dot M. The mass-luminosity relation of W-R stars is shown to be most sensitive to the mass loss rate during the helium burning phase \dot M_{3\alpha}. Together with the mass-luminosity relation obtained for all evolutionary sequences several more exact relations are determined for the constant ratio f_{3\alpha}=\dot M/\dot M_{3\alpha} with 0.5 \le f_{3\alpha} \le 3. Evolutionary models of W-R stars were used as initial conditions in hydrodynamic computations of radial nonlinear stellar oscillations. The oscillation amplitude is larger in W-R stars with smaller initial mass or with lower mass loss rate due to higher surface abundances of carbon and oxygen. In the evolving W-R star the oscillation amplitude decreases with decreasing stellar mass M and for M < 10M_\odot the sufficiently small nonlinear effects allow us to calculate the integral of the mechanical work W done over the pulsation cycle in each mass zone of the hydrodynamical model. The only positive maximum on the radial dependence of W is in the layers with temperature of T\sim 2e5K where oscillations are excited by the iron Z--bump kappa-mechanism. Radial oscillations of W-R stars with mass of M > 10M_\odot are shown to be also excited by the kappa-mechanism but the instability driving zone is at the bottom of the envelope and pulsation motions exist in the form of nonlinear running waves propagating outward from the inner layers of the envelope.

We considered the structure of steady-state plane-parallel radiative shock waves propagating through the partially ionized hydrogen gas of temperature T_1 = 3000K and density 1e-12 gm/cm^3 <= \rho_1 <= 1e-9 gm/cm^3. The upstream Mach numbers range within 6 <= M_1 <= 14. In frequency intervals of hydrogen lines the radiation field was treated using the transfer equation in the frame of the observer for the moving medium, whereas the continuum radiation was calculated for the static medium. Doppler shifts in Balmer emission lines of the radiation flux emerging from the upstream boundary of the shock wave model were found to be roughly one-third of the shock wave velocity. The gas emitting the Balmer line radiation is located at the rear of the shock wave in the hydrogen recombination zone where the gas flow velocity in the frame of the observer is approximately one-half of the shock wave velocity. The ratio of the Doppler shift to the gas flow velocity of 0.7 results both from the small optical thickness of the shock wave in line frequencies and the anisotropy of the radiation field typical for the slab geometry. In the ambient gas with density of \rho_1 >= 1e-11 gm/cm^3 the flux in the H-alpha frequency interval reveals the double structure of the profile. A weaker H-beta profile doubling was found for \rho_1 >= 1e-10 gm/cm^3 and U_1 <= 50 km/s. The unshifted redward component of the double profile is due to photodeexcitation accompanying the rapid growth of collisional ionization in the narrow layer in front of the discontinuous jump.

We considered the structure of steady-state radiative shock waves propagating in the partially ionized hydrogen gas with density rho1 = 1e-10 gm/cm^3 and temperature 3000K <= T1 <= 8000K. The radiative shock wave models with electron thermal conduction in the vicinity of the viscous jump are compared with pure radiative models. The threshold shock wave velocity above of which effects of electron thermal conduction become perceptible is of U1=70 km/s and corresponds to the upstream Mach numbers from M1= 6 at T1=8000K to M1=11 at T1=3000K. In shocks with efficient electron heat conduction more than a half of hydrogen atoms are ionized in the radiative precursor, whereas behind the viscous jump the hydrogen gas undergoes the full ionization. The existence of the electron conductive precursor leads to the enhancement of the Lyman continuum flux trapped in the surroundings of the discontinuous jump. For upstream velocities ranged within 70 km/s <= U1 <= 85 km/s the partially ionized hydrogen gas of the radiative precursor undergoes the additional ionization (<= 5%), whereas the total radiave flux emerging from the shock wave increases by 10% <= delta(FRad) <= 25% .

The grid of the models of radiative shock waves propagating through partially ionized hydrogen gas with temperature 3000K <= T_1 <= 8000K and density 10^{-12} gm/cm^3 <= \rho_1 <= 10^{-9}gm/cm^3 is computed for shock velocities 20 km/s <= U_1 <= 90 km/s. The fraction of the total energy of the shock wave irreversibly lost due to radiation flux ranges from 0.3 to 0.8 for 20 km/s <= U_1 <= 70 km/s. The postshock gas is compressed mostly due to radiative cooling in the hydrogen recombination zone and final compression ratios are within 1 <\rho_N/\rho_1 \lesssim 10^2, depending mostly on the shock velocity U_1. The preshock gas temperature affects the shock wave structure due to the equilibrium ionization of the unperturbed hydrogen gas, since the rates of postshock relaxation processes are very sensitive to the number density of hydrogen ions ahead the discontinuous jump. Both the increase of the preshock gas temperature and the decrease of the preshock gas density lead to lower postshock compression ratios. The width of the shock wave decreases with increasing upstream velocity while the postshock gas is still partially ionized and increases as soon as the hydrogen is fully ionized. All shock wave models exhibit stronger upstream radiation flux emerging from the preshock outer boundary in comparison with downstream radiation flux emerging in the opposite direction from the postshock outer boundary. The difference between these fluxes depends on the shock velocity and ranges from 1% to 16% for 20 km/s <= U_1 <= 60 km/s. The monochromatic radiation flux transported in hydrogen lines significantly exceeds the flux of the background continuum and all shock wave models demonstrate the hydrogen lines in emission.

Models of steady-state plane-parallel shock waves propagating through the unperturbed hydrogen gas of temperature T=6000K and density rho = 1e-10 gm/cm^3 are computed for upstream velocities from 15 km/s to 70 km/s. The shock wave structure is considered in terms of the self-consistent solution of the radiation transfer, fluid dynamics and rate equations for 2 <= L <= 4 atomic bound levels with a continuum. The radiative flux Frad emergent from the shock wave was found to be independent of the lower limit nu_L of the frequency range provided that nu_L < nu_2, where nu_2 is the Balmer continuum head frequency, At the same time the decrease of nu_L is accompanied by decrease of the Lyman continuum flux and leads to smaller heating and weaker ionization of the hydrogen gas in the radiative precursor. For all models the size of the radiative precursor is of 10^4 cm and corresponds to several mean free paths of photons at the frequency of the Lyman continuum edge nu_1. The compression ratio at the discontinuous jump gradually increases with increasing upstream velocity U_1, reaches the maximum of \rho^+/\rho^- = 3.62 at U_1 = 55 km/s and slowly decreases for larger U_1 due to the strong rize of the preshock gas temperature. The radiative flux from the shock wave was determined as a function of the upstream velocity U_1 and its ratio to the total energy flux in the shock wave C_2 was found to range within 0.18 < Frad/C_2 < 0.92 for 15 km/s <= U_1 <= 65 km/s. Thus, at upstream velocities U_1 > 60 km/s the shock wave losses more than 90% of its total energy due to radiation. For all shock wave models the role of collisional processes in both bound-bound and bound-free atomic transitions was found to be negligible in comparison with corresponding radiative processes.

The structure of steady plane-parallel radiative shock waves propagating through the hydrogen gas undergoing partial ionization and excitation of bound atomic states is investigated in terms of the self-consistent solution of the equations of fluid dynamics, radiation transfer and atomic kinetics. The shock wave model is represented by a flat finite slab with no incoming radiation from external sources at both its boundaries. The self-consistent solution is obtained using the global iteration procedure each step of which involves (1) integration of the fluid dynamics and rate equations for the preshock and postshock regions, consecutively, both solutions being fitted by the Rankine-Hugoniot relations at the discontinuous jump; (2) solution of the radiation transfer equation for the whole slab. The global iteration procedure is shown to converge to the stable solution which allows for the strong coupling of the gas flow and the radiation field produced by this flow. Application of the method is demonstrated for the shock waves with upstream velocities from 15 km/s to 60 km/s (that is with upstream Mach numbers from 2.3 to 9.3) and the hydrogen gas of unperturbed temperature T=3000K and density rho = 1e-10 gm/cm^3.