the eddington limit corresponds to a situation where the radiative acceleration outwards equals the gravitational acceleration inwards .
eddington ( 1926 ) originally only considered electron scattering for the opacity ( the classical eddington limit ) , but in an ionization zone of an important element the true opacity can be much higher .
therefore , the opacity - modified eddington limit in stellar atmospheres can occur at significantly greater gravities than the classical eddington limit , as is evident for example from the calculations of lamers & fitzpatrick ( 1988 ) , gustafsson & plez ( 1992 ) and asplund & gustafsson ( 1996 ) .
if hydrostatic equilibrium is required , a stellar luminosity which locally exceeds the eddington luminosity automatically forces the development of a gas pressure ( @xmath2 ) inversion , as seen from the equation of hydrostatic equilibrium :    @xmath3    from rearranging the above expression one obtains    @xmath4,\ ] ]    with the radiative acceleration defined by :    @xmath5    ( e.g mihalas 1978 ) . here
@xmath6 denotes the physical flux and @xmath7 the total mass extinction coefficient ( with dimension @xmath8g@xmath1 ) .
obviously a positive @xmath2-gradient must occur when the eddington limit is encountered if hydrostatic equilibrium is valid .
it should be emphasized that a density inversion does not necessarily imply that the eddington limit is exceeded .
a density inversion predominantly occurs due to a changing molecular weight in the ionization zone of a dominant species , while a @xmath2-inversion reflects a super - eddington luminosity due to high opacity ( asplund et al .
since both are related to ionization they may often occur in the same atmospheric layers , which has probably caused the confusion found in the literature ( e.g. maeder 1989 ) .
it should be noted that @xmath2-inversions only seem possible to develop in optically thick conditions .
if the radiative acceleration exceeds gravity when lines dominate the opacity , the radiative force is highly unstable to perturbations induced by velocity gradients , since the spectral lines may then absorb unattenuated continuum flux .
therefore , rather than the development of a @xmath2-inversion , a stellar wind will be initiated , which is the case for radiatively driven winds of hot stars .
the aim of the present study is to investigate whether @xmath2-inversions in late - type stars , where the opacity instead is dominated by continuous opacity , are subject to a similar instability which would prevent their existence . before considering possible instabilities , in sect .
[ s : eddlimit ] the derived eddington limit for h - deficient stars is presented and a possible connection between the limit and the declines of the r coronae borealis ( rcrb ) stars is discussed .
the location of the eddington limit in a @xmath9-log@xmath10 diagram can be estimated by model atmosphere calculations ( lamers & fitzpatrick 1988 ; gustafsson & plez 1992 ) .
it should be emphasized that these models consistently include the radiative acceleration in the equation of hydrostatic equilibrium .
therefore , the location of the limit is not obtained by extrapolation towards lower gravities nor does it imply zero density , as has sometimes been claimed in the literature ( e.g. humphreys & davidson 1994 ; nieuwenhuijzen & de jager 1995 ) .
an interesting conclusion from the above studies is that the theoretical opacity - modified eddington limit coincides with the upper luminosity limit for stars ( lamers & fitzpatrick 1988 ) .
it is therefore possible that the eddington limit , possibly in connection with rotation ( langer 1997a , b ; owocki & gayley 1997 ) , prevents the evolution of stars into the super - eddington regime by the existence of radiative instabilities
. these may give rise to the eruptions of the luminous blue variables ( lbvs ) ( cf .
humphreys & davidson 1994 ) . here a similar investigation but for the h - deficient rcrb and related stars will be carried out ; a preliminary report can be found in asplund & gustafsson ( 1996 ) .
the presence of @xmath2-inversions in such model atmospheres was noted by asplund et al .
( 1997a ) , who speculated that it may be the trigger mechanism for the enigmatic visual declines of the rcrb stars .
the variability of these stars is likely due to the formation of an obscuring dust cloud close to the stellar surface in the line - of - sight , but how the dust condensation can proceed despite the relatively high equilibrium temperatures has long been an enigma ( cf .
clayton 1996 ) .
the presence of shocks seems to be a very promising mechanism for the gas to cool sufficiently and thereby trigger dust formation ( woitke et al .
1996 ) .
the models used here are recently constructed line - blanketed model photospheres with typical abundances of rcrb stars ( asplund et al .
1997a ) . in fig .
[ f : gamma ] the variation of @xmath11 with optical depth is displayed for a sequence of models corresponding to constant luminosity and mass ( e.g. @xmath12l@xmath13 and @xmath14m@xmath13 ) .
the super - eddington luminosities only occur in the deep atmospheric layers where hei ionizes and , as expected , where @xmath2-inversions occur .
since the rcrb stars presumably evolve towards higher @xmath9 ( kilkenny 1982 ) , the stars will move from the sub - eddington regime into having super - eddington luminosities , as is seen in fig .
[ f : eddlimit ] where the theoretical eddington limit in the @xmath9-log@xmath10 diagram is shown . compared with solar abundances the limit
is shifted towards higher @xmath9 due to the higher ionization potential of he compared with h , but also extends to slightly larger log@xmath10 .
the location of the eddington limit for higher @xmath9 has not yet been explored , since a self - consistent hydrodynamical treatment is needed for @xmath15k when @xmath16 occurs for @xmath17 and line opacity may determine @xmath18 .
the exact location of the computed eddington limit will naturally depend on the detailed properties of the model atmospheres .
for @xmath19k convection in the deeper layers of the hei ionization zone restricts the limit to lower log@xmath10 , since the radiative flux diminishes in the presence of a significant convective flux .
also , the lower temperatures decrease @xmath7 and @xmath18 further . including turbulent pressure
would therefore bring the limit towards greater gravities by typically 0.2dex , due to less efficient convection with lower densities .
a c / he ratio of 1% has been assumed on the basis of the measured ratio in extreme helium ( ehe ) stars ( jeffery 1996 ) , which presumably the rcrb stars are related to ( lambert et al .
a higher carbon abundance pushes the eddington limit towards higher gravities the same way as turbulent pressure does
. a lower overall metallicity by 0.5dex would decrease the limit by about 0.2dex by lowering the temperatures in the inner layers with a diminished line - blanketing .      in fig .
[ f : eddlimit ] the estimated parameters of the rcrb stars ( lambert et al . 1997 ; asplund et al .
1997c ) are also shown .
obviously , rather than being randomly located in the diagram as one a priori might have guessed , the rcrb stars fall close to the limit , which suggests that a connection between the declines of rcrb stars and the eddington limit exists .
thus , there is an alluring resemblance between the rcrb stars and the luminous blue variables in their proximity to the eddington limits and their eruptive behaviour ( asplund & gustafsson 1996 ) .
if the rcrb stars are evolving towards higher @xmath9 ( kilkenny 1982 ) and the ehe domain at constant luminosity , their immediate progenitors should be the h - deficient carbon ( hdc ) stars ( lambert et al .
the ehe and hdc stars experience no visual fadings despite being otherwise similar to the rcrb stars .
it is therefore plausible that the reason for the non - variability is that the hdc stars have not yet reached the eddington limit while the ehe stars are located on the stable side of the limit at higher @xmath9 .
if the declines of the rcrb stars are indeed triggered by the super - eddington luminosities , one could expect a correspondence between the maximum value of @xmath11 in the atmosphere and frequency of declines . in the pulsation - induced dust condensation model for the rcrb stars ( woitke et al .
1996 ) , in which the dust formation is triggered by propagating shock waves due to ordinary pulsations , one would sooner expect a correlation with pulsation strength .
intriguingly , v854cen , which is the rcrb star most frequent in decline , has very weak photospheric pulsations ( lawson & cottrell 1997 ) but by a wide margin the largest @xmath11 found of the 18 rcrb stars analysed sofar .
rysgr , on the other hand , has by far the strongest photospheric pulsations but only average amount of declines , which might reflect the fact that @xmath11 is not unusually high .
pulsation may , however , push some stars across the limit by periodically shifting them towards higher @xmath9 , e.g. @xmath20k for rcrb itself ( rao & lambert 1997 ) .
the proposed connection between the eddington limit and the declines of the rcrb stars is certainly not without problems , which deserve to be addressed further .
a few stars seem not to behave as expected if the proposal is correct , in particular this is true for some of the coolest ehe stars and the three hot rcrb stars .
as seen above , a @xmath2-inversion must develop at the eddington limit in hydrostatic equilibrium .
one should ask , however , whether strong convection may take place to prevent it , since the ionization zone causing the high opacity and thereby the super - eddington luminosities also favours convection .
in fact , a super - eddington luminosity in the diffusion regime automatically requires the presence of convection as shown by langer ( 1997a ) . generalizing langer s derivation for a non - negligible radiation pressure to also include the effects of ionization , the schwarzschild criterion @xmath21 for the onset of convection is equivalent to    @xmath22 % \right\ } } % { 1 } \frac{4\alpha \left\ { 2 + 8\alpha + x\left(1-x\right ) \left[5/2 + \epsilon
/ kt + 4\alpha\right ] \right\ } } { 5 + 40\alpha + 32\alpha^2 + x\left(1-x\right )
\left[5/2 + \epsilon / kt + 4\alpha\right]^2}\ ] ]    ( e.g. mihalas & weibel mihalas 1984 ) . here , while @xmath23 denotes the ionization fraction of the most abundant element with an ionization potential @xmath24 .
convection must be initiated as @xmath11(r ) approaches unity , since the right hand side is always less than 1 . thereby the radiative flux and hence @xmath11 diminishes when convection is efficient .
the eddington limit can therefore not be exceeded in adiabatic convection zones in the stellar interior and hence @xmath2-inversions will not develop .
thus , inversions can only be present in non - adiabatic convection zones close to the stellar surface .
it should be remembered that the location of the eddington limit is dependent on the convection treatment .
a @xmath2-inversion will for example disappear if the density scale height is applied instead of the normal pressure scale height for the calculation of the convective flux through the mixing length theory .
the mixing length theory is rarely a good description of stellar convection , and in particular for these extreme situations it may be very misleading .
ultimately , hydrodynamical simulations will be necessary to investigate the effects of convection on @xmath2-inversions .      using the momentum equation one can investigate under which conditions a @xmath2-inversion can exist in the presence of a stellar mass loss . or equivalently : which mass loss rates can one expect from a star if the super - eddington luminosities would instead drive a stellar wind ?
in the presence of a stellar outflow , the equation of hydrostatic equilibrium must be replaced with the momentum equation :    @xmath25    where we have restricted ourselves to the one - dimensional , time - independent case .
also other external forces , such as turbulent forces ( @xmath26 ) and centrifugal forces ( @xmath27 ) , can be contained in @xmath28 besides radiation .
it should therefore be kept in mind that the possibility of an instability may be underestimated when restricting the following discussion to only radiative forces , as other forces also tend to be destabilizing ( e.g. nieuwenhuijzen & de jager 1995 ; langer 1997a , b ; owocki & gayley 1997 ) .
together with a constant , spherically symmetric mass loss rate @xmath29 , and the equation of state , it is possible to rewrite the momentum equation into the so called  wind equation " :    @xmath30    where the isothermal sound speed defined by @xmath31 has been introduced . by an inspection of eq .
[ e : windeq ] , it is possible to determine under which circumstances @xmath32- or @xmath2-inversions are allowed . in normal cases ,
d@xmath33/dr@xmath34 in the surface layers and therefore the parenthesis on the right hand side of eq .
[ e : windeq ] is positive as long as a temperature inversion is not present , which is not possible in optically thick layers with a radiative flux .
therefore , if @xmath35 but @xmath36 , the velocity gradient @xmath37 must be negative . since a constant @xmath38 has been assumed , a density inversion will be present if @xmath39 and therefore a @xmath2-inversion is possible ( though not required ) .
if , on the other hand , @xmath40 while @xmath35 , then @xmath37 must be positive , and the density gradient can not be positive : a density inversion can not exist .
since a @xmath2-inversion requires a density inversion , neither inversion can occur in the stellar atmosphere despite the super - eddington luminosities . again
the possibility of a temperature inversion has been neglected .
a sufficient criterion for the disappearance of a @xmath2-inversion is therefore that the outflow velocity due to mass loss is high enough ( @xmath41 ) .
the critical mass loss rate can be estimated by comparing @xmath42 with @xmath43 .
the values for @xmath43 and @xmath32 will be taken from static models , which is justified for these rough estimates , since the atmosphere below the sonic point is normally little affected by mass loss . a model atmosphere with @xmath44k , log@xmath45 [ cgs ] and solar abundances requires a critical mass loss rate @xmath46 .
hence @xmath2-inversions can exist in stars despite high mass loss rates .
achmad & lamers ( 1997 ) have arrived at similar conclusions and estimates .
they also verify the finding numerically with dynamical steady state atmospheres , though not altogether self - consistently .
an even higher critical mass loss rate is found for the rcrb stars .
the lack of hydrogen makes the continuous opacity significantly lower than for solar abundances and hence the density at a given optical depth will be correspondingly greater , typically by a factor of 50 ( asplund et al .
since the density is higher , the outflow velocity will be lower in order to carry the same @xmath38 .
a similar model as above but with abundances typical for rcrb stars , will therefore require @xmath47 .
the mass loss rates of rcrb stars are poorly known , in particular since it may be episodic , but rough estimates suggest ( e.g. feast 1986 ) , orders of magnitude smaller than @xmath48 .    it is concluded that @xmath2-inversions are not an artifact of hydrostatic equilibrium but may also exist in dynamical steady state atmospheres , even with a high mass loss rate .
dynamical stability has been extensively studied , in particular in connection with the termination of the agb ( e.g. paczyski & zikowski 1968 ; wagenhuber & weiss 1994 ) .
radial adiabatic oscillations will grow exponentially in time for dynamical instability , which is formally encountered if the first generalized adiabatic exponent    @xmath49    is less than @xmath50 in a region effectively isolated from the rest of the star .
an ideal gas has @xmath51 while pure radiation has @xmath52 , but the combined effect of non - negligible radiation pressure and ionization may push @xmath53 below the limiting value . for dynamical instability to occur , the pressure weighted volumetric mean value of @xmath53 between the depth @xmath54 and the stellar surface @xmath55    @xmath56    must be reduced below @xmath50 in a significant fraction of the star .
stars with a pronounced core - envelope structure and high luminosity - to - mass ratios , like the supergiants investigated here , may show dynamical instabilities in ionization zones ( e.g. stothers & chin 1993 ; wagenhuber & weiss 1994 ) .    in fig .
[ f : gamma1 ] the run of @xmath53 and @xmath57 as functions of optical depth are shown for two model atmospheres corresponding to @xmath58 and log@xmath59 for a solar and a h - deficient composition . in particular , @xmath57 is reduced much below @xmath50 for the model with solar abundances in the hi and hei ionization zones . at the surface
@xmath60 dominates the total pressure and hence @xmath53 is very close to @xmath50 .
a more complete study for h - rich model atmospheres has been presented by lobel et al .
( 1992 ) , who arrive at the same conclusions . in the rcrb model
only the hei ionization zone exists , which explains the quite different depth variation of @xmath53 .
the higher temperatures at depth in the h - deficient model make hei ionization occur at smaller @xmath61 .
also seen is the minor effect of the ci ionization zone at @xmath62 .
these models only extend down to @xmath63 and hence it is not possibly to tell how @xmath57 varies further in , though the heii ionization zone will also reduce @xmath53 .
it seems that atmospheres of late - type supergiants are close to being dynamically unstable , or may even be so for sufficiently high luminosity - to - mass ratios .
violent instabilities due to this might lead to ejection of material , as found for the termination of the agb ( wagenhuber & weiss 1994 ) , as well as for the lbvs ( stothers & chin 1993 ) .
also rcrb stellar models seem to suffer from a similar dynamical instability as the stars on the tip of the agb ( a. weiss , private communication ) .
it should be remembered though that the above - mentioned evolutionary studies lack a proper hydrodynamical treatment , and therefore interpreting the instabilities found in the models as actually occurring in the stars must still be made with some caution .
it is clear from the discussion in sect . [
s : pgasinversion ] and [ s : dynamical ] that @xmath2-inversions may exist in stars with high luminosity - to - mass ratios and that such stars are close to being dynamically unstable .
thus , these stars may be vulnerable to radiative instabilities , since the additional destabilization from radiative forces may render them unstable . the existence of radiatively driven sound waves in hot stars due to high @xmath18 has been suggested by hearn ( 1972 , 1973 ) , and further elaborated by berthomieu et al .
( 1976 ) , spiegel ( 1976 ) and carlberg ( 1980 ) .
the instability is based on the increase in opacity upon compression , since the radiation force is in phase with the velocity perturbation and hence causes amplification . for these high temperatures it requires essentially isothermal perturbations , since    @xmath64    here @xmath65 , while @xmath66 is the flux weighted mean opacity .
thus only short wavenumbers will be amplified .
the growth rate is especially high as @xmath11 approaches unity ( berthomieu et al .
1976 ) but the super - eddington case has not been investigated .
carlberg ( 1980 ) found that the instability may provide rapid transfer of momentum from the radiation to the gas with typical growth rates of @xmath67s@xmath1 for @xmath68k .
it is possible that a similar instability is exhibited also in late - type supergiants , with the important differences being the super - eddington luminosities in the @xmath2-inversions and that isothermal perturbations are not necessary , since both ( @xmath69)@xmath70 and ( @xmath71)@xmath72 in the upper layers of the @xmath2-inversions ( fig .
[ f : opacity ] ) .    to study the stability of a @xmath2-inversion a local , linear , non - adiabatic stability analysis has been carried out .
small , sinusoidal perturbations around the equilibrium values are assumed for the temperature , density and vertical velocity :    @xmath73    where @xmath74 and @xmath75 are small .
the atmosphere is assumed to be static ( @xmath76 ) , since the outflow velocity must anyway be less than the isothermal sound speed for a @xmath2-inversion to exist , as demonstrated above .
the analysis will be restricted to vertical motion and therefore multi - dimensional instabilities of rayleigh - taylor - type will not be found .
carlberg ( 1980 ) investigated also such a gradient instability for hot stars but found slow growth rates compared to the radiatively driven sound waves
. such a gradient instability would resemble convective motion which , however , is already present in the layers of interest in the unperturbed model and is thus not specifically searched for here .
the equations to be linearized are the material equations for conservation of mass , momentum and energy :    @xmath77    ( e.g. mihalas & weibel mihalas 1984 ) , where d / dt , @xmath78 , @xmath79 and @xmath80 denote the lagrangean co - moving derivative , the specific internal energy of the gas , the wavelength integrated planck function ( @xmath81 ) and mean intensity ( @xmath82 ) , respectively .
the absorption and planck mean opacities have been replaced by the rosseland mean opacity for absorption @xmath83 , which is a good approximation for the relevant optical depths . the scattering contribution to the total extinction
is here denoted by @xmath84 .
furthermore lte ( @xmath85 ) has been assumed .
the material equations must be coupled to the equations for the radiation field to account for radiative transfer effects through exchanges of energy and momentum between the photons and the gas . in this exploratory investigation
the effects of radiation are included with a time - independent treatment , i.e. the radiation field is assumed to be quasi - static .
the radiation momentum and energy equations will then look like :    @xmath86    here @xmath87 and @xmath88 denote the wavelength integrated eddington flux ( @xmath89 ) and radiation pressure ( @xmath90 ) respectively . in the diffusion regime for the unperturbed atmosphere
, @xmath91 is valid .
the equilibrium radiation quantities are related by :    @xmath92    for the first relation radiative equilibrium is implicitly imposed , which is not strictly true in the cases which will be investigated below .
however , the analysis will be restricted only to atmospheric depths where the convective flux carries @xmath93 of the total flux ( i.e. in the upper part of the convection zone ) and radiative equilibrium can therefore be largely justified . to close the relations the eddington approximation ( @xmath94 ) will be made , which is , like lte , well justified in the atmospheric layers of interest . in the following
, the radiation energy and momentum equations will be combined to a single equation , which means that there is one more variable to perturb , for which @xmath80 is chosen :    @xmath95    together with the definitions of the isothermal sound speed @xmath43 , the density scale height @xmath96^{-1}$ ] and an ideal gas for @xmath2 , linearization to first order yields :    @xmath97 + v_1\left[ik - h_\rho^{-1}\right ] & = & 0 \\ % linearized conservation of momentum % \frac{\rho_1}{\rho_0}\left[g_{\rm eff } - g_\rho + \frac{\rho_1}{\rho_0}\left[g + c_{\rm t}^2\left(ik + \frac{1}{t_0}\frac{{\rm d}t}{{\rm dr } }   - \frac{1}{\mu_0}\frac{{\rm d}\mu}{{\rm dr } } \right)\right ]   & + & \nonumber \\ + \frac{\rho_1}{\rho_0}\left[- c_{\rm t}^2 \mu_\rho % \frac{\rho_0}{\mu_0 } % \frac{\partial \mu}{\partial \rho } \left(ik - h_\rho^{-1 } + \frac{1}{t_0}\frac{{\rm d}t}{{\rm dr } } - \frac{2}{\mu_0}\frac{{\rm d}\mu}{{\rm dr}}\right)\right ]   & + & \nonumber
\\ % + \frac{t_1}{t_0}\left[- g_{\rm t } + + \frac{t_1}{t_0}\left [ c_{\rm t}^2\left(ik - h_\rho^{-1 }   - \frac{1}{\mu_0}\frac{{\rm d}\mu}{{\rm dr } } \right)\right ]   & + & \nonumber \\ + \frac{t_1}{t_0}\left[- c_{\rm t}^2 \mu_t % \frac{t_0}{\mu_0 } % \frac{\partial \mu}{\partial t } \left(ik - h_\rho^{-1 } + \frac{1}{t_0}\frac{{\rm d}t}{{\rm dr } } - \frac{2}{\mu_0}\frac{{\rm d}\mu}{{\rm dr}}\right)\right ]   & + & \nonumber \\ + v_1\left [ - i\omega\right ]   % & + & \\ + \frac{j_1}{j_0}\left[- ik\frac{4\sigma_{\rm r } t_0 ^ 4}{3c\rho_0 } \right ] & = & 0 \\ % linearized conservation of energy \frac{\rho_1}{\rho_0}\left[-i\omega \left(\rho_0 \frac{\partial e}{\partial \rho }   % \frac{1}{t_0}\frac{\partial l}{\partial \rho } - c_{\rm t}^2 \right)\right ]   & + & \nonumber \\
+ \frac{t_1}{t_0}\left[-i\omega t_0 \frac{\partial e}{\partial t } + % \frac{1}{\rho_0}\frac{\partial l}{\partial t } 16\left(\kappa_0 - \sigma_0\right)\sigma_{\rm r } t_0 ^ 4\right ] & + & \nonumber \\ % + v_1\left[ik\gamma + \frac{1}{c_{\rm v}t_0}g_{\rm eff } -   % \gamma h_\rho^{-1 } +   % \gamma \frac{1}{t_0}\frac{{\rm d}t}{{\rm dr } } \right ] & = & 0 .
\\ % & + & \\
+ v_1\left[g_{\rm rad }   % + e_0\left ( ik - h_\rho^{-1}\right )   % \frac{h}{\rho_0}\frac{{\rm d}\rho}{{\rm dr } } + % - \rho_o\frac{\partial e}{\partial \rho } h_\rho^{-1 } + % \frac{\partial e}{\partial t}\frac{{\rm d}t}{{\rm dr } } + + \frac{{\rm d}e}{{\rm dr } } +   c_{\rm t}^2 h_\rho^{-1 }   \right ] & + & \nonumber \\
+ \frac{j_1}{j_0}\left[-4\left ( \kappa_0 - \sigma_0\right)\sigma_{\rm r } t_0 ^ 4\right ] & = & 0   \\ % linearized radiation energy + momentum \frac{\rho_1}{\rho_0}\left[-\frac{ik}{\rho_0(\kappa_0 - \sigma_0 ) } \frac{1}{4}\left(\frac{t_{\rm eff}}{t_0}\right)^4\left(1 +   \frac{\partial { \rm ln } \kappa}{\partial { \rm ln } \rho}\right)\right ]   & + & \nonumber \\ + \frac{t_1}{t_0}\left[-4 - \frac{ik}{\rho_0(\kappa_0 - \sigma_0 ) } % \frac{h_0}{b_0 } \frac{1}{4}\left(\frac{t_{\rm eff}}{t_0}\right)^4 \frac{\partial { \rm ln } \kappa}{\partial { \rm ln } t}\right ] & + & \nonumber \\ + \frac{j_1}{j_0}\left[1 +   \frac{k^2}{3\kappa_0(\kappa_0 - \sigma_0)\rho_0 ^ 2}\right ] & = & 0 .
\end{aligned}\ ] ]    the notation has been abbreviated according to and .
derivatives with respect to @xmath98 and @xmath80 have been ignored , which is justified in optically thick layers where line opacity is unimportant .
true enough , photoionization edges introduce some velocity dependence on @xmath18 but this minor effect is neglected .
lccc +   + variable & units & solar & rcrb +   + log@xmath61 & & 0.5 & 1.2 + @xmath18 & [ cms@xmath99 & @xmath100 & @xmath101 + @xmath102 & [ k ] & 9610 & 14950 + @xmath32 & [ g@xmath103 & @xmath104 & @xmath105 + @xmath2 & [ dyncm@xmath99 & @xmath107 & @xmath108 + @xmath109 & [ cms@xmath1 ] & @xmath111 & @xmath112 + @xmath113d@xmath102/dr & [ cm@xmath1 ] & @xmath115 & @xmath116 + @xmath117d@xmath32/dr & [ cm@xmath1 ] & @xmath118 & @xmath119 + @xmath120d@xmath33/dr & [ cm@xmath1 ] & @xmath121 & @xmath122 + @xmath123 & [ @xmath124g@xmath1 ] & @xmath125 & @xmath126 + @xmath84 & [ @xmath124g@xmath1 ] & @xmath127 & @xmath128 + @xmath129 & & @xmath130 & @xmath131 + @xmath132 & & @xmath133 & @xmath134 + d@xmath78/dr & [ @xmath135s@xmath99 & @xmath136 & @xmath137 + @xmath138 & [ @xmath8s@xmath99 & @xmath139 & @xmath140 + @xmath141 & [ @xmath8s@xmath99 & @xmath142 & @xmath143 +   +    in order for the above set of equations to have a non - trivial solution for the perturbations @xmath74 , @xmath75 and @xmath144 , the determinant of the equations must be equal to 0 .
the three complex eigenvalues @xmath145 for a given real wavenumber @xmath146 then correspond to the allowed perturbations .
two of the modes will represent radiation - modified gravity - acoustic waves while the third is a non - propagating , heavily damped , thermal wave . here
two different cases will be investigated numerically , one corresponding to a @xmath2-inversion in a supergiant with solar abundances and the second the same but for an rcrb star .
as before , both model atmospheres have @xmath147k and log@xmath148 .
values for some of the physical parameters at some atmospheric layer where the instability may occur are found in table [ t : values ] .
the calculated growth rates are shown in fig .
[ f : growthrate ] , which reveals amplification for @xmath149cm@xmath1 ( solar composition ) and @xmath151cm@xmath1 ( rcrb composition ) . for large @xmath146 ,
the perturbations are as expected damped due to the negligible optical depths of the disturbances ( @xmath152 ) , which make the sound waves little affected by the radiative forces . instead , the smoothing effects from radiation energy exchange due to large gradients in the perturbed variables dominate .
the maximum growth rates are slightly smaller for a solar composition than for a h - deficient composition : @xmath153 respectively @xmath154s@xmath1 , which corresponds to an e - folding time on the order of the sound crossing time of the atmosphere . for the largest wavenumbers with such growth rates ,
the real parts are significantly higher , which implies that the sound waves are only slightly amplified per wavelength . for smaller @xmath146 , however , the growth rate is comparable with the oscillation rate .
in fact , for the model with solar composition , the wave is essentially non - propagating for the smallest @xmath146 .
however , the use of a local analysis is not justified when the wavelengths of the perturbations are comparable with the scale heights of the variations for the variables . for @xmath155cm@xmath1 , the propagation speed is slightly less than the isothermal sound speed in both cases .
thus , radiation - modified sound waves can be amplified in late - type stars with significant radiative forces
. however , the perturbations are not very strongly excited with minimum growth times @xmath156s , contrary to the findings by carlberg ( 1980 ) for hot stars , mainly as a result of the much smaller radiative fluxes ( @xmath157 ) . for the instability to generate large amplification ,
the sound waves must be reflected at the surface and pass through the region with high radiative forces several times .
however , the inwards propagating modes are damped for the same reason as the outgoing are amplified , with the magnitude of the damping rate similar to the growth rate , which suggests that the instability will still not be very efficient . in order to estimate the behaviour of the instability in the global , non - linear regime ,
a more sophisticated study than presented here must be carried out .      in order to gain further theoretical insight into which atmospheric properties contribute to the instability presented above ,
the equations can be simplified by assuming the unperturbed atmosphere to be isothermal and neglecting ionization and scattering .
furthermore , only the limits of optically thin and optically thick disturbances will be studied and hence the combined radiative energy and momentum equation will not be needed .
the gas energy equation will thus instead be written as :    @xmath158    where @xmath159 denotes the net radiative cooling rate of the gas .
depending on the optical thickness of the disturbance , @xmath159 will be described by either newtonian cooling or equilibrium diffusion ( e.g. mihalas & weibel mihalas 1984 ) .
the rate of work done by @xmath18 has been ignored .
the linearized continuity equation remains the same but the linearized momentum equation reads as :    @xmath160 & + & \nonumber \\ + \frac{t_1}{t_0}\left [ c_{\rm t}^2\left(ik - h_\rho^{-1}\right ) - g_{\rm t}\right ] + v_1\left[-i\omega\right ] & = & 0,\end{aligned}\ ] ]    where @xmath161 and @xmath162 , while the energy equation looks like :    @xmath163 + \frac{t_1}{t_0}\left[\frac{i\omega}{\gamma -1 } -   \tilde{{\cal l}}\right ] + v_1\left[-h_\rho^{-1}\right ] = 0.\ ] ]    here @xmath164 represents the term arising from the net cooling rate @xmath159 and @xmath165 the ratio of specific heats . for newtonian cooling    @xmath166    and in the diffusion regime    @xmath167    @xmath164 is the inverse of the radiative cooling time - scale ( e.g. mihalas & weibel mihalas 1984 ) , where the factor @xmath168 takes into account the optical depth effects for optically thick conditions .
the simplified dispersion relation then reads :    @xmath169   & + & \nonumber \\ - \omega k\left[g_\rho + \left(\gamma-1\right)g_{\rm t } - g_{\rm eff } + \gamma c_{\rm t}^2 h_\rho^{-1}\right ]   & + & \nonumber \\ -\tilde{{\cal l}}\left(\gamma-1\right)\left [   k^2c_{\rm t}^2 + h_\rho^{-1}\left(g_{\rm eff } -   g_\rho\right)\right ]   & + & \nonumber \\ - ik\tilde{{\cal l}}\left(\gamma-1\right)\left [ g_\rho - g_{\rm eff } + c_{\rm t}^2 h_\rho^{-1}\right ]   & = & 0.\end{aligned}\ ] ]    for @xmath155cm@xmath1 , @xmath170 is much smaller than @xmath171 for the two investigated models and can be ignored .    in situations of rapid radiative cooling ( @xmath172 ) , the dispersion relation can be simplified to yield the solutions    @xmath173,\ ] ]    i.e. isothermally propagating , radiation - modified sound waves . in the limit of slow cooling ( @xmath174 ) , as expected adiabatic acoustic waves are recovered :    @xmath175,\ ] ]    which travel with the adiabatic sound speed . for @xmath176cm@xmath1
, the propagation speed will be modified by the neglected terms @xmath177 , though the growth / damping rates will remain the same . for both isothermal and adiabatic perturbations ,
the sound waves are amplified by the action of the radiative acceleration , as long as the opacity increases with temperature and density ( assuming constant radiative flux ) :    @xmath178    furthermore , in a @xmath2-inversion @xmath28 is negative and hence tends to excite an initial disturbance .
however , the normal momentum conservation counteracts this due to propagation in a density stratified medium ; the density inversion will tend to damp the amplitude for outwards propagating waves .
in fact , with this simplified atmosphere and the values in table [ t : values ] , in both models the last term dominates . hence , the amplitude of sound waves will decay as they propagate outwards with these assumptions , despite the super - eddington luminosities and the increase in opacity upon compression .
thus , since the numerical results in sect .
[ s : radinstab ] show amplification for the case without these simplifying assumptions , the effects of e.g. ionization and a temperature gradient tend to work together with the radiative force and its derivatives to overcome the damping contribution from the density gradient .
it has been shown that h - deficient model atmospheres may have radiative forces which exceed gravity , similar to the case of some luminous h - rich stars ( lamers & fitzpatrick 1988 ; gustafsson & plez 1992 ) .
such super - eddington luminosities manifest themselves as @xmath2-inversions in hydrostatic atmospheres .
the inversions are , however , not an artifact of this assumption but can also be present in dynamical atmospheres .
the inversions are not removed unless a very high mass loss rate is present , and is neither necessarily removed by convection if the large radiative forces occur close to the surface .
thus , super - eddington luminosities does _ not _ automatically initiate a stellar wind . the location of the h - deficient rcrb stars in the immediate proximity of the computed opacity - modified eddington limit ( fig .
[ f : eddlimit ] ) , is certainly suggestive of a connection between the enigmatic visual declines of the stars and the eddington limit .
it is therefore proposed that instabilities as the rcrb stars encounter the eddington limit during their evolution towards higher @xmath9 are the unknown trigger mechanism for their famous variability .
the eddington limit may thus be the `` smoking gun '' for the declines : gas is ejected due to the high radiative forces in the atmospheres , which then cools rapidly first radiatively and later from adiabatic expansion to reach sub - equilibrium temperatures and thus possibly enable dust condensation ( cf .
woitke et al .
1996 for the similar case of shocks instead induced by pulsations ) .
such gas ejections may be observable as absorption components in strong lines , as observed by rao & lambert ( 1997 ) in rcrb at maximum light .
if a connection indeed exists , it could explain why the ehe and hdc stars do not show rcrb - like variability as a result of being located on the stable side of the eddington limit .
it also points to an interesting similarity between the rcrb stars and the lbvs .
thus , the two types of stars known to be situated at the eddington limit both show eruptive behaviour , which suggests a similar underlying physical explanation for their variability .
a search for radiative and dynamical instabilities in the atmospheres has been carried out for both h - rich and h - deficient late - type supergiants , but only partly successfully .
sound waves are found to be amplified by the large radiative forces but the linear stability analysis reveal only rather slow growth rates compared to for early - type stars despite the super - eddington luminosities .
such radiation - modified sound waves may thus be partly responsible for the semi - regular pulsational variations of such supergiants , but it is doubtful whether the instability is efficient enough to eject gas clouds from the atmospheres as speculated above .
the atmospheres are also found to be close to dynamically unstable , which might give rise to increased mass loss , as previously suggested for lbvs ( stothers & chin 1993 ) and for the termination of the agb phase ( e.g. paczyski & zikowski 1968 ; wagenhuber & weiss 1994 ) .
this does not , however , explain in what way the rcrb stars are special compared to other similar stars . in fact , the h - rich models seem to be more vulnerable to such instabilities , which suggests that dynamical instabilities are not the answer to the declines of the rcrb stars .
another possibility is of course that the same dynamical instabilities occur in both types of stars but that the h - deficient and c - rich environment of the rcrb stars greatly favours dust condensation , and thus the observational manifestation , compared to for h - rich compositions .
one can also imagine that the eddington limit is still responsible for the behaviour of the rcrb stars though in a more indirect way .
( 1997a , b ) has shown that the effect of rotation coupled with strong radiative forces may lead to increased mass loss in the equatorial plane , and it is possible that such a mechanism is at work in these supergiants , despite the presumably relatively small rotational velocities .
the difference in observed variability between the rcrb stars and the hdc and ehe stars could then be the result of viewing angle .
there are some observational evidence for such bipolarity ( e.g. rao & lambert 1993 ; clayton et al .
1997 ) violent instabilities due to strange mode pulsations leading to dramatically increased mass loss is another possibility , which is partly related to the eddington limit , since both have their origin in ionization zones , as do the dynamical instabilities discussed above . whether any of these instabilities , or a combination thereof , could help explain the variability of the rcrb stars and the lbvs certainly deserves further investigation .
finally , it should be remembered that the stability analyses presented here assume the unperturbed model to be a realistic description of the stellar atmosphere .
at least for the rcrb stars , there are strong indications that this is in fact not the case ( gustafsson & asplund 1996 ; asplund et al .
1997a , b , c ; lambert et al .
1997 ) . if the supergiant atmospheres are indeed distinctly different from the predictions of standard models , the efficiency of radiative instabilities may have been underestimated .
one can suspect that models based on the normal assumptions , such as the mixing length theory for convection , near the eddington limit may be very unsuitable .
clearly , radiative hydrodynamical simulations of such atmospheres would be valuable , also since they would shed further light on possible instabilities .
the constructive suggestions by f.p .
pijpers and a. weiss have significantly improved the paper .
discussions with w. glatzel , b. gustafsson , l .- a . willson and p. woitke have been very helpful .
the author is grateful to l. achmad , h.j.g.l.m .
lamers , n. langer and s. owocki for communicating their work prior to publication .
achmad l. , lamers h.j.g.l.m .
, 1997 , a&a 320 , 196 asplund m. , gustafsson b. , 1996 , in : hydrogen - deficient stars , jeffery c.s .
, heber u. ( eds . ) .
series vol .
asplund m. , gustafsson b. , kiselman d. , eriksson k. , 1997a , a&a 318 , 521 asplund m. , gustafsson b. , lambert d.l .
, rao n.k . , 1997b , a&a 321 , l17 asplund m. , gustafsson b. , rao n.k .
, lambert d.l . , 1997c ,
a&a , in press berthomieu g. , provost j. , rocca a. , 1976 , a&a 47 , 413 carlberg r.g . , 1980 ,
apj 241 , 1131 clayton g.c . , 1996 ,
pasp 108 , 225 clayton g.c .
, bjorkman k.s . ,
nordsieck k.h . ,
zellner n.e.b . ,
schulte - ladbeck r.e . , 1997 ,
apj 476 , 870 eddington a.s . , 1926 , the internal constitution of the stars , cambridge university press , cambridge feast m.w . , 1986 , in : hydrogen - deficient stars and related objects , hunger k. , schnberner d. , rao n.k .
( eds . ) .
reidel , dordrecht , p. 151
gustafsson b. , asplund m. , 1996 , in : hydrogen - deficient stars , jeffery c.s .
, heber u. ( eds . ) .
series vol .
gustafsson b. , plez b. , 1992 , in : instabilities in evolved super- and hypergiants , de jager c. , nieuwenhuijzen h. ( eds . ) .
north holland , amsterdam , p. 86
hearn a.g . , 1972 ,
a&a 19 , 417 hearn a.g . , 1973 ,
a&a 23 , 97 humphreys r.m . , davidson k. , 1994 , pasp 106 , 1025 jeffery c.s . , 1996 , in : hydrogen - deficient stars , jeffery c.s . , heber u. ( eds . ) . asp conf
series vol .
96 , p. 152
kilkenny d. , 1982 , mnras 200 , 1019 lambert d.l . ,
, gustafsson b. , asplund m. , 1997 , submitted to a&a lamers h.j.g.l.m .
, fitzpatrick e.l . , 1988 ,
apj 324 , 279 langer n. , 1997a , in : luminous blue variables : massive stars in transition , nota a. , lamers h.j.g.l.m .
( eds . ) .
series , in press langer n. , 1997b , a&a , in press lawson w.a . ,
cottrell p.l . , 1997 ,
mnras 285 , 266 lobel a. , achmad l. , de jager c. , nieuwenhuijzen h. , 1992 , a&a 264 , 147 maeder a. , 1989 , in : physics of luminous blue variables , davidson k. , moffat a.f.j .
, lamers h.j.g.l.m .
( eds . ) .
kluwer , dordrecht , p. 15
mihalas d. , 1978 , stellar atmospheres , w.h .
freeman and company , san fransisco mihalas d. , weibel mihalas b. , 1984 , foundations of radiation hydrodynamics , oxford university press , oxford nieuwenhuijzen h. , de jager c. , 1995 , a&a 302 , 811 owocki s.p .
, gayley k.g . , 1997 , in : luminous blue variables : massive stars in transition , nota a. , lamers h.j.g.l.m .
( eds . ) .
series , in press paczyski b. , zikowski j. , 1968 , acta astr .
18 , 225 rao n.k .
, lambert d.l . , 1993 , aj 105 , 1915 rao n.k . , lambert d.l . , 1997 ,
mnras 284 , 489 seaton m.j . , yan y. , mihalas d. , pradhan a.k . , 1994 ,
mnras 266 , 805 spiegel e.a .
, 1976 , in : physique des mouvements dans les atmosphres stellaires , cayrel r. , steinberg m. ( eds ) , cnrs , paris , p. 19
stothers r.b .
, chin c .- w . , 1993 , apj 408 , l85 wagenhuber j. , weiss a. , 1994 , a&a 290 , 807 woitke p. , goeres a. , sedlmayr e. , 1996 , a&a 313 , 217


shock waves occur in many physical situations in nature . for example
, shock waves can be produced by solar winds ( bow shocks ) , supersonic or near sonic aircrafts ( transonic shocks around the body ) , explosions ( blast waves ) , and various natural processes .
when such a shock hits an obstacle ( steady or flying ) , shock reflection - diffraction phenomena occur .
one of the most important problems in mathematical fluid dynamics is the problem of shock reflection - diffraction by a wedge .
when the plane shock hits a wedge head on , it experiences a reflection - diffraction process , and then a fundamental question is what types of wave patterns of reflection - diffraction configurations it may form around the wedge .
the complexity of reflection - diffraction configurations was first reported by ernst mach @xcite in 1878 , who first observed two patterns of reflection - diffraction configurations : regular reflection ( two - shock configuration ) and mach reflection ( three - shock configuration ) ; also see @xcite .
the problem remained dormant until the 1940 s when von neumann , friedrichs , bethe , as well as many experimental scientists , among others , began extensive research into all aspects of shock reflection - diffraction phenomena , due to its importance in applications .
see von neumann @xcite and ben - dor @xcite ; also see @xcite and the references cited therein .
it has been found that the situation is much more complicated than what mach originally observed : the mach reflection can be further divided into more specific sub - patterns , and various other patterns of shock reflection - diffraction may occur such as the von neumann reflection and the guderley reflection ; see @xcite and the references cited therein .
then the fundamental scientific issues include :    \(i ) _ structure of the shock reflection - diffraction configurations _ ;    \(ii ) _ transition criteria between the different patterns of shock reflection - diffraction configurations _ ;    \(iii ) _ dependence of the patterns upon the physical parameters such as the wedge angle @xmath0 , the incident - shock - wave mach number @xmath1 , and the adiabatic exponent @xmath2_.    careful asymptotic analysis has been made for various reflection - diffraction configurations in lighthill @xcite , keller - blank @xcite , hunter - keller @xcite , morawetz @xcite , @xcite , and the references cited therein ; also see glimm - majda @xcite .
large or small scale numerical simulations have been also made ; see , e.g. @xcite , @xcite , @xcite , and the references cited therein .
however , most of the fundamental issues for shock reflection - diffraction phenomena have not been understood , especially the global structure and transition of different patterns of shock reflection - diffraction configurations .
this is partially because physical and numerical experiments are hampered by various difficulties and have not been able to select the correct transition criteria between different patterns . in particular , numerical dissipation or physical viscosity
smear the shocks and cause boundary layers that interact with the reflection - diffraction patterns and can cause _ spurious mach steams _ ; cf . woodward - colella @xcite .
furthermore , some difference between two different patterns are only fractions of a degree apart ( e.g. , see fig .
5 below ) , a resolution even by sophisticated modern experiments ( e.g. @xcite ) has been unable to reach , as pointed out by ben - dor in @xcite :  _ for this reason it is almost impossible to distinguish experimentally between the sonic and detachment criteria _ " ( cf .
section 5 below ) . in this regard
, it seems that the ideal approach to understand fully the shock reflection - diffraction phenomena , especially the transition criteria , is still via rigorous mathematical analysis . to achieve this ,
it is essential to establish first the global existence , regularity , and structural stability of solutions of the shock reflection - diffraction problem .    on the other hand ,
shock reflection - diffraction configurations are the core configurations in the structure of global solutions of the two - dimensional riemann problem for hyperbolic conservation laws ; while the riemann solutions are building blocks and local structure of general solutions and determine global attractors and asymptotic states of entropy solutions , as time tends to infinity , for multidimensional hyperbolic systems of conservation laws .
see @xcite and the references cited therein . in this sense , we have to understand the shock reflection - diffraction phenomena , in order to understand fully entropy solutions to multidimensional hyperbolic systems of conservation laws .    in this paper
, we first formulate the shock reflection - diffraction problem into an initial - boundary value problem in section 2 .
then we employ the essential feature of self - similarity of the initial - boundary value problem to reformulate the problem into a boundary value problem in the unbounded domain in section 3 . in section 4
, we present the unique solution of normal reflection for this problem when the wedge angle is @xmath3 . in section 5 ,
we exhibit the local theory of regular reflection - diffraction , introduce a stability criterion to determine state ( 2 ) at the reflection point on the wedge , and present the von neumann s detachment and sonic conjectures .
then we discuss the role of the potential flow equation in the shock reflection - diffraction problem even in the level of the full euler equations in section 6 .
based on the local theory , we reduce the boundary value problem into a free boundary problem in the context of potential flow in section 7 . in section 8 , we describe a global theory for regular reflection - diffraction for potential flow , established in chen - feldman @xcite and bae - chen - feldman @xcite . in section 9 ,
we discuss some open problems and new mathematics required for further developments , which are also essential for solving multidimensional problems in conservation laws and other areas in nonlinear partial differential equations .
the full euler equations for compressible fluids in @xmath4 , are of the following form : @xmath5 where @xmath6 is the density , @xmath7 the fluid velocity , @xmath8 the pressure , and @xmath9 the internal energy .
two other important thermodynamic variables are the temperature @xmath10 and the energy @xmath11 .
the notation @xmath12 denotes the tensor product of the vectors @xmath13 and @xmath14 .
choosing @xmath15 as the independent thermodynamical variables , then the constitutive relations can be written as @xmath16 governed by @xmath17    for a polytropic gas , @xmath18 or equivalently , @xmath19 where @xmath20 may be taken to be the universal gas constant divided by the effective molecular weight of the particular gas , @xmath21 is the specific heat at constant volume , @xmath22 is the adiabatic exponent , and @xmath23 is any constant under scaling .    when a flow is potential , that is , there is a velocity potential @xmath24 such that @xmath25 then the euler equations for the flow take the form : @xmath26 where @xmath27 especially , @xmath28
when @xmath29 , by scaling and @xmath30 is the bernoulli constant , which is usually determined by the boundary conditions if such conditions are prescribed . from the second equation in , we have @xmath31 then system can be rewritten as the following time - dependent potential flow equation of second order : @xmath32 with .    for a steady solution @xmath33 , i.e. , @xmath34 , we obtain the celebrated steady potential flow equation in aerodynamics : @xmath35    in applications in aerodynamics , or is used for discontinuous solutions , and the empirical evidence is that entropy solutions of or are fairly close to entropy solutions for , provided the shock strengths are small , the curvature of shock fronts is not too large , and the amount of vorticity is small in the region of interest .
furthermore , we will show in section 6 that , for the shock reflection - diffraction problem , the euler equations for potential flow is actually _ exact _ in an important region of the solution ( see theorem 6.1 below ) .
then the problem of shock reflection - diffraction by a wedge can be formulated as follows :    .
_ seek a solution of system satisfying the initial condition at @xmath36 : @xmath37 and the slip boundary condition along the wedge boundary : @xmath38 where @xmath39 is the exterior unit normal to the wedge boundary , and state @xmath40 and @xmath41 satisfy @xmath42 that is , given @xmath43 , and @xmath22 , the other variables @xmath44 and @xmath45 are determined by .
in particular , the mach number @xmath46 is determined by @xmath47 where @xmath48 is the sonic speed of fluid state @xmath41 . _
notice that the initial - boundary value problem ( * problem 2.1 * ) is invariant under the self - similar scaling : @xmath49 therefore , we seek self - similar solutions : @xmath50 then the self - similar solutions are governed by the following system : @xmath51 where @xmath52 , and @xmath53 is the pseudo - velocity .
the eigenvalues of system are @xmath54 where @xmath55 is the sonic speed .
when the flow is pseudo - subsonic , i.e. , @xmath56 , the eigenvalues @xmath57 become complex and thus the system consists of two transport equations and two nonlinear equations of hyperbolic - elliptic mixed type .
therefore , system is _ hyperbolic - elliptic composite - mixed _ in general .    since the problem is symmetric with respect to the axis @xmath58 , it suffices to consider the problem in the half - plane @xmath59 outside the half - wedge : @xmath60 then the initial - boundary value problem ( * problem 2.1 * ) in the @xmath61coordinates can be formulated as the following boundary value problem in the self - similar coordinates @xmath62 :    . _
seek a solution to system satisfying the slip boundary condition on the wedge boundary and the matching condition on the symmetry line @xmath58 : @xmath63 the asymptotic boundary condition as @xmath64 : @xmath65 _    , width=220,height=201 ]    it is expected that the solutions of * problem 3.1 * contain all possible patterns of shock reflection - diffraction configurations as observed in physical and numerical experiments ; cf .
@xcite and the references cited therein .
the simplest case of the shock reflection - diffraction problem is when the wedge angle @xmath0 is @xmath66 . in this case , the reflection - diffraction problem simply becomes the normal reflection problem , for which the incident shock normally reflects , and the reflected shock is also a plane . it can be shown that there exist a _ unique _
state @xmath67 and a _ unique _ location of the reflected shock @xmath68 such that state @xmath69 is subsonic inside the sonic circle with center at the origin and radius @xmath70 , and is supersonic outside the sonic circle ( see fig .
[ fig : nr ] ) .
that is , in this case , the normal reflection solution is unique .        in this case , @xmath71 and @xmath72 is the unique root of @xmath73 that is , @xmath74 in other words , given @xmath43 , and @xmath22 , state @xmath69 is uniquely determined through  and .
for a wedge angle @xmath75 , different reflection - diffraction patterns may occur .
various criteria and conjectures have been proposed for the existence of configurations for the patterns ( cf .
ben - dor @xcite ) .
one of the most important conjectures made by von neumann @xcite in 1943 is the _ detachment conjecture _ , which states that the regular reflection - diffraction configuration may exist globally whenever the two shock configuration ( one is the incident shock and the other the reflected shock ) exists locally around the point @xmath76 ( see fig .
the following theorem was rigorously shown in chang - chen @xcite ( also see sheng - yin @xcite , bleakney - taub @xcite , neumann @xcite ) .
[ local ] there exists @xmath77 such that , when @xmath78 , there are two states @xmath79 and @xmath80 such that @xmath81 where @xmath82 is the sonic speed .
then the conjecture can be stated as follows :    : _ there exists a global regular reflection - diffraction configuration whenever the wedge angle @xmath0 is in @xmath83_.    it is clear that the regular reflection - diffraction configuration is not possible without a local two - shock configuration at the reflection point on the wedge , so this is the weakest possible criterion . in this case
, the local theory indicates that there are two possible states for state ( 2 ) .
there had been a long debate to determine which one is more physical for the local theory ; see courant - friedrichs @xcite , ben - dor @xcite , and the references cited therein .    since
the reflection - diffraction problem is not a local problem , we take a different point of view that the selection of state ( 2 ) should be determined by the global features of the problem , more precisely , by the stability of the configuration with respect to the wedge angle @xmath0 , rather than the local features of the problem .    : _ since the solution is unique when the wedge angle @xmath84 , it is required that our global regular reflection - diffraction configuration should be stable and converge to the unique normal reflection solution when @xmath85 , provided that such a global configuration can be constructed .
_    we employ this stability criterion to conclude that our choice for state ( 2 ) must be @xmath86 . in general , @xmath86 may be supersonic or subsonic
. if it is supersonic , the propagation speeds are finite and state ( 2 ) is completely determined by the local information : state ( 1 ) , state ( 0 ) , and the location of the point @xmath76 .
this is , any information from the reflected region , especially the disturbance at the corner @xmath87 , can not travel towards the reflection point @xmath76 .
however , if it is subsonic , the information can reach @xmath76 and interact with it , potentially altering the reflection - diffraction type .
this argument motivated the second conjecture as follows :    : _ there exists a regular reflection - diffraction configuration when @xmath88 for @xmath89 such that @xmath90 at @xmath76 . _
this sonic conjecture is based on the following fact : if state ( 2 ) is sonic when @xmath91 , then @xmath92 for any @xmath93 .
this conjecture is stronger than the detachment one .
in fact , the regime between the angles @xmath94 and @xmath95 is very narrow and is only fractions of a degree apart ; see fig . 5 from sheng - yin @xcite .     when @xmath96.,width=211,height=192 ]
in this section , we discuss the role of the potential flow equation in the shock reflection - diffraction problem for the full euler equations .    under the hodge - helmoltz decomposition @xmath97 with @xmath98 ,
the euler equations become @xmath99 where @xmath100 is the vorticity of the fluid , @xmath101 is the entropy , and the gradient @xmath102 is with respect to the self - similar variables @xmath103 from now on .
when @xmath104 , and @xmath105 on a curve @xmath106 transverse to the fluid direction , we first conclude from that , in the domain @xmath107 determined by the fluid trajectories : @xmath108 past @xmath106 , @xmath109 this implies that @xmath110 since @xmath98 .
then we conclude that @xmath111 since @xmath112 , which yields that the right - hand side of equation vanishes .
furthermore , from , @xmath113 which implies that @xmath114 by scaling , we finally conclude that the solutions of system  in the domain @xmath107 is determined by the following system for self - similar solutions : @xmath115 or the potential flow equation for self - similar solutions : @xmath116 with @xmath117 then we have @xmath118    for our problem ( see fig .
4 ) , we note that , for state ( 2 ) , @xmath119 then , if our solution @xmath120 is @xmath121 and the gradient of the tangential component of the velocity is continuous across the sonic arc @xmath122 , we still have along @xmath122 on the side of @xmath123 .
thus , we have    [ p - dominate ] let @xmath120 be a solution of our * problem 3.1 * such that @xmath120 is @xmath121 in the open region @xmath124 and the gradient of the tangential component of @xmath125 is continuous across the sonic arc @xmath126 .
let @xmath107 be the subregion of @xmath123 formed by the fluid trajectories past the sonic arc @xmath122 , then , in @xmath107 , the potential flow equation with coincides with the full euler equations
, that is , equation with is exact in the domain @xmath107 for * problem 3.1*.    the regions such as @xmath107 also exist in various mach reflection - diffraction configurations . theorem [ p - dominate ] applies to such regions whenever the solution @xmath120 is @xmath121 and the gradient of the tangential component of @xmath125 is continuous .
in fact , theorem 8.3 indicates that , for the solutions @xmath127 of with , the @xmath128 regularity of @xmath127 and the continuity of the tangential component of the velocity field @xmath129 are optimal across the sonic arc @xmath122 .
the importance of the potential flow equation with in the time - dependent euler flows was also observed by hadamard @xcite through a different argument .
furthermore , when the wedge angle @xmath0 is close to @xmath66 , it is expected that the curvature of the reflected shock is small so that , in the other part @xmath130 of @xmath123 , the vorticity @xmath131 is small and the entropy is close to the constant .
then , in the reflection - diffraction domain @xmath132 , the potential flow equation with dominates , provided that the exact state along the reflected shock is given .    ,
width=230,height=163 ]    equation with is a nonlinear equation of mixed elliptic - hyperbolic type .
it is elliptic if and only if @xmath133 which is equivalent to @xmath134    the study of partial differential equations of mixed hyperbolic - elliptic type can date back 1940s ( cf .
linear models of partial differential equations of mixed hyperbolic - elliptic type include the lavrentyev - betsadze equation : @xmath135 the tricomi equation : @xmath136 and the keldysh equation : @xmath137 nonlinear models of mixed - type equations for with include the transonic small disturbance equation : @xmath138 or , for @xmath139 , @xmath140 which has been studied in @xcite and the references cited therein .
also see @xcite for the models for self - similar solutions from the pressure gradient system and nonlinear wave equations .
for the potential equation with , shocks are discontinuities in the pseudo - velocity @xmath141 .
that is , if @xmath142 and @xmath143 are two nonempty open subsets of @xmath144 and @xmath145 is a @xmath146curve where @xmath147 has a jump , then @xmath148 is a global weak solution of ( [ potential-1 ] ) with in @xmath149 if and only if @xmath127 is in @xmath150 and satisfies equation in @xmath151 and the rankine - hugoniot condition on @xmath11 : @xmath152_s=0,\ ] ] where the bracket @xmath153 $ ] denotes the difference of the values of the quantity along the two sides of @xmath11 .
then the plane incident shock solution in the @xmath154coordinates with states @xmath155 and @xmath156 corresponds to a continuous weak solution @xmath127 of ( [ potential - eq ] ) in the self - similar coordinates @xmath62 with the following form : @xmath157 respectively , where @xmath158 are the velocity of state ( 1 ) and the location of the incident shock , uniquely determined by @xmath159 through ( [ fbconditionselfsim-0 ] ) .
then @xmath160 in fig . 2 , and * problem 3.1 * in the context of the potential flow equation can be formulated as :    [ bvp ] ( boundary value problem ) ( see fig .
seek a solution @xmath127 of equation with in the self - similar domain @xmath161 with the boundary condition on @xmath162 : @xmath163 and the asymptotic boundary condition at infinity : @xmath164 where ( [ boundary - condition-2 ] ) holds in the sense that @xmath165 _    for our problem , since @xmath166 does not satisfy the slip boundary condition , the solution must differ from @xmath166 in @xmath167 , thus a shock diffraction - diffraction by the wedge vertex occurs .
in chen - feldman @xcite , we first followed the von neumann criterion and the stability criterion introduced in section 5 to establish a local existence theory of regular shock reflection near the reflection point @xmath76 in the level of potential flow , when the wedge angle is large and close to @xmath66 . in this case , the vertical line is the incident shock @xmath168 that hits the wedge at the point @xmath160 , and state ( 0 ) and state ( 1 ) ahead of and behind @xmath11 are given by @xmath169 and @xmath166 defined in and , respectively .
the solutions @xmath127 and @xmath166 differ only in the domain @xmath170 because of shock diffraction by the wedge vertex , where the curve @xmath171 is the reflected shock with the straight segment @xmath172 .
state ( 2 ) behind @xmath172 can be computed explicitly with the form : @xmath173 which satisfies @xmath174 on @xmath175 ; the constant velocity @xmath176 and the angle @xmath94 between @xmath172 and the @xmath177axis are determined by @xmath178 from the two algebraic equations expressing ( [ fbconditionselfsim-0 ] ) and continuous matching of state ( 1 ) and state ( 2 ) across @xmath172 , whose existence is exactly guaranteed by the condition on @xmath178 under which regular shock reflection - diffraction is expected to occur as in theorem 5.1 .
moreover , @xmath179 is the unique solution in the domain @xmath180 , as argued in @xcite .
denote @xmath181 the sonic arc of state @xmath182 with center @xmath183 and radius @xmath184 .
also we introduce the following notation for the other parts of @xmath185 : @xmath186    then * problem 7.1 * can be formulated as :    . _
seek a solution @xmath127 in @xmath123 to equation with subject to the boundary condition on @xmath187 , the rankine - hugoniot conditions on the shock @xmath188 : @xmath189_{{\gamma_{shock}}}=0 , \label{rh-1}\\ & & [ \rho(\nabla\varphi,\varphi,\rho_0)\nabla\varphi\cdot\nu]_{{\gamma_{shock}}}=0,\label{rh-2}\end{aligned}\ ] ] and the dirichlet boundary condition on the sonic arc @xmath122 : @xmath190 _    it should be noted that , in order that the solution @xmath127 in the domain @xmath123 is a part of the global solution to * problem 7.1 * , that is , @xmath127 satisfies the equation in the sense of distributions in @xmath161 , especially across the sonic arc @xmath122 , it is requaired that @xmath191 that is , we have to match our solution with state ( 2 ) , which is the necessary condition for our solution in the domain @xmath123 to be a part of the global solution . to achieve this , we have to show that our solution is at least @xmath146 with @xmath192 across @xmath122 .
then the problem can be reformulated as the following free boundary problem :    .    _
seek a solution @xmath127 and a free boundary @xmath193 such that _    * @xmath194 and @xmath195 * @xmath127 satisfies the free boundary condition along @xmath196 ; * @xmath197 solves in @xmath198 , is subsonic in @xmath198 , and satisfies @xmath199    the boundary condition on @xmath200 implies that @xmath201 and thus ensures the orthogonality of the free boundary with the @xmath177-axis .
formulation implies that the free boundary is determined by the level set @xmath202 , which is a convenient formulation to apply useful free boundary techniques .
the free boundary condition along @xmath196 is the conormal boundary condition on @xmath196 .
condition ensures that the solution of the free boundary problem in @xmath123 is a part of the global solution as pointed out earlier .
condition is the slip boundary condition .
* problem 7.3 * involves two types of transonic flow : one is a continuous transition through the sonic arc @xmath126 as a fixed boundary from the pseudo - supersonic region ( 2 ) to the pseudo - subsonic region @xmath123 ; the other is a jump transition through the transonic shock as a free boundary from the supersonic region ( 1 ) to the subsonic region @xmath123 .
in this section , we describe a global theory for regular reflection - diffraction established in chen - feldman @xcite and bae - chen - feldman @xcite .      in chen - feldman @xcite
, we have developed a rigorous mathematical approach to solve * problem 7.3 * and established a global theory for solutions of regular reflection - diffraction , which converge to the unique solution of the normal shock reflection when @xmath0 tends to @xmath66 .
introduce the polar coordinates @xmath203 with respect to the center @xmath204 of the sonic arc @xmath126 of state ( 2 ) , that is , @xmath205 then , for @xmath206 , we denote by @xmath207 the @xmath208-neighborhood of the sonic arc @xmath209 within @xmath123 ; see fig .
4 . in @xmath210
, we introduce the coordinates : @xmath211 then @xmath212 and @xmath213 .
there exist @xmath214 and @xmath215 such that , when @xmath216 , there exists a global self - similar solution @xmath217 with @xmath218 of * problem 7.1 * ( equivalently , * problem 7.2 * ) for shock reflection - diffraction by the wedge . the solution @xmath127 satisfies that , for @xmath219 , @xmath220 moreover ,    1 .
[ ellipticityinomega ] equation is elliptic in @xmath123 ; 2 .
[ phi - ge - phi2 ] @xmath221 ; 3 .
[ bdryreg ] the reflected shock @xmath171 is @xmath222 at @xmath223 and @xmath224 except @xmath223 ; 4
.   [ c11normestimate ] there exists @xmath225 such that @xmath226 ; in particular , in the coordinates ( [ coordnearsonic ] ) , @xmath227 5 .   [
psi - x - est ] there exists @xmath228 so that , in the coordinates ( [ coordnearsonic ] ) , @xmath229 6 .   [ fbnearsonic ] there exist @xmath230 and a function @xmath231 such that , in the coordinates ( [ coordnearsonic ] ) , @xmath232 and @xmath233)}<\infty,\qquad \frac{d\hat{f}}{dx}\!\ge\omega>0\!\;\;\ ; \text{for}\;\;0<x<{\varepsilon}_0.\ ] ]    furthermore , the solution @xmath127 is stable with respect to the wedge angle @xmath0 in @xmath234 and converges in @xmath234 to the unique solution of the normal reflection as @xmath235 .
the existence of a solution @xmath127 of * problem 7.1 * ( equivalently , * problem 7.2 * ) , satisfying ( [ phi - states-0 - 1 - 2 ] ) and property ( [ c11normestimate ] ) , follows from ( * ? ? ?
* main theorem ) .
property ( [ ellipticityinomega ] ) follows from lemma 5.2 and proposition 7.1 in @xcite .
property ( [ phi - ge - phi2 ] ) follows from proposition 7.1 and section 9 in @xcite which assert that @xmath236 , where the set @xmath237 is defined by ( 5.15 ) in @xcite .
property ( [ psi - x - est ] ) follows from propositions 8.18.2 and section 9 in @xcite .
property ( [ fbnearsonic ] ) follows from ( 5.7 ) and ( 5.25)(5.27 ) in @xcite and the fact that @xmath236 .
we remark that estimate above confirms that our solutions satisfy the assumptions of theorem 6.1 for the velocity field @xmath129 .
one of the main difficulties for the global existence is that the ellipticity condition for ( [ potential-1 ] ) with is hard to control , in comparison to our work on steady flow @xcite .
the second difficulty is that the ellipticity degenerates along the sonic arc @xmath126 .
the third difficulty is that , on @xmath126 , the solution in @xmath123 has to be matched with @xmath179 at least in @xmath146 , i.e. , the two conditions on the fixed boundary @xmath126 : the dirichlet and conormal conditions , which are generically overdetermined for an elliptic equation since the conditions on the other parts of boundary have been prescribed .
thus , one needs to prove that , if @xmath127 satisfies ( [ potential-1 ] ) in @xmath123 , the dirichlet continuity condition on the sonic arc , and the appropriate conditions on the other parts of @xmath238 derived from problem 7.3 , then the normal derivative @xmath239 automatically matches with @xmath240 along @xmath126 . indeed , equation ( [ potential-1 ] ) , written in terms of the function @xmath241 in the @xmath242coordinates defined near @xmath126 such that @xmath126 becomes a segment on @xmath243 , has the form : @xmath244 plus the _ small _ terms that are controlled by @xmath245 in appropriate norms .
equation is _ elliptic _ if @xmath246 .
hence , it is required to obtain the @xmath128 estimates near @xmath126 to ensure @xmath247 which in turn implies both the ellipticity of the equation in @xmath123 and the match of normal derivatives @xmath248 along @xmath126 .
taking into account the _ small _ terms to be added to equation , one needs to make the stronger estimate @xmath249 and assume that @xmath245 is suitably small to control these additional terms .
another issue is the non - variational structure and nonlinearity of this problem which makes it hard to apply directly the approaches of caffarelli @xcite and alt - caffarelli - friedman @xcite .
moreover , the elliptic degeneracy and geometry of the problem makes it difficult to apply the hodograph transform approach in chen - feldman @xcite and kinderlehrer - nirenberg @xcite to fix the free boundary .    for these reasons ,
one of the new ingredients in our approach is to develop further the iteration scheme in @xcite to a partially modified equation .
we modified equation ( [ potential-1 ] ) in @xmath123 by a proper shiffmanization ( i.e. a cutoff ) that depends on the distance to the sonic arc , so that the original and modified equations coincide when @xmath127 satisfies @xmath250 , and the modified equation @xmath251 is elliptic in @xmath123 with elliptic degeneracy on @xmath209 .
then we solved a free boundary problem for this modified equation : the free boundary is the curve @xmath196 , and the free boundary conditions on @xmath196 are @xmath202 and the rankine - hugoniot condition ( [ fbconditionselfsim-0 ] ) .    on each step , an _ iteration free boundary _
curve @xmath126 is given , and a solution of the modified equation @xmath252 is constructed in @xmath123 with the boundary condition ( [ fbconditionselfsim-0 ] ) on @xmath196 , the dirichlet condition @xmath253 on the degenerate arc @xmath126 , and @xmath254 on @xmath255 and @xmath256 .
then we proved that @xmath127 is in fact @xmath128 up to the boundary @xmath126 , especially @xmath257 , by using the nonlinear structure of elliptic degeneracy near @xmath126 which is modeled by equation ( [ degenerate - equation ] ) and a scaling technique similar to daskalopoulos - hamilton @xcite and lin - wang @xcite .
furthermore , we modified the _ iteration free boundary _ curve @xmath196 by using the dirichlet condition @xmath202 on @xmath196
. a fixed point @xmath127 of this _ iteration procedure _ is a solution of the free boundary problem for the modified equation .
moreover , we proved the precise gradient estimate : @xmath258 for @xmath259 , which implies that @xmath127 satisfies the original equation ( [ potential - eq ] ) .
this global theory for large - angle wedges has been extended in chen - feldman @xcite to the sonic angle @xmath260 , for which state ( 2 ) is sonic , such that , as long as @xmath261 $ ] , the global regular reflection - diffraction configuration exists .    the global existence result in theorem 8.1
can be extended up to the sonic wedge - angle @xmath94 for any @xmath2 and @xmath262 .
moreover , the solutions satisfy the properties ( i)(vi ) in theorem 8.1 .
the condition @xmath262 depends explicitly only on the parameters @xmath22 and @xmath263 . for the case @xmath264
, we have been making substantial progress as well , and the final detailed results can be found in chen - feldman @xcite .      by theorem 8.1([c11normestimate ] ) , the solution @xmath127 constructed there is at least @xmath128 near the sonic arc @xmath126 . the next question is to analyze the behavior of solutions @xmath265 to regular reflection - diffraction , especially the optimal regularity of the solutions .
we first define the class of regular reflection - diffraction solutions .
[ regreflsoldef ]    _ let @xmath266 , @xmath263 , and @xmath267 be constants , let @xmath44 and @xmath268 be defined by ( [ state-1-velocity ] ) and ( [ shocklocation ] ) . let the incident shock @xmath168 hits the wedge at the point @xmath160 , and let state @xmath40 and state @xmath41 ahead of and behind @xmath196 be given by and , respectively .
the function @xmath269 is a regular reflection - diffraction solution if @xmath270 is a solution to * problem 7.1 * such that _
\(a ) there exists state @xmath182 of form ( [ state2 ] ) with @xmath271 , satisfying the entropy condition @xmath272 and the rankine - hugoniot condition ( [ fbconditionselfsim-0 ] ) along the line @xmath273 which contains the points @xmath76 and @xmath223 , such that @xmath274 is on the sonic circle of state @xmath182 , and state @xmath182 is supersonic along @xmath172 ;    \(b ) there exists an open , connected domain @xmath275 such that ( [ phi - states-0 - 1 - 2 ] ) holds and equation is elliptic in @xmath123 ;    \(c ) @xmath276 on the part @xmath277 of the reflected shock .
the global solution constructed in @xcite is a regular reflection - diffraction solution , which is a part of the assertions in theorems 8.18.2 .
if state @xmath182 exists and is supersonic , then the line @xmath278 necessarily intersects the sonic circle of state @xmath182 ; see the argument in @xcite starting from ( 3.5 ) there .
thus , the only assumption regarding the point @xmath223 is that @xmath279 intersects the sonic circle within @xmath161 .
we note that , in the case @xmath280 , the regular reflection becomes the normal reflection , in which @xmath281 and the solution is smooth across the sonic line of state ( 2 ) ; see ( * ? ? ?
* section 3.1 ) .
condition @xmath282 in definition [ regreflsoldef ] rules out this case .
moreover , for @xmath283 , the property @xmath271 in part ( a ) of definition [ regreflsoldef ] is always true for state @xmath182 of form ( [ state2 ] ) , satisfying the entropy condition @xmath272 and the rankine - hugoniot condition ( 7.1 ) along the line @xmath273 which contains the point @xmath76 .
these are readily derived from the calculations in ( * ? ? ?
* section 3.2 ) .
there may exist a global regular reflection - diffraction configuration when state @xmath182 is subsonic which is a very narrow regime @xcite .
such a case does not involve the difficulty of elliptic degeneracy , which we are facing for the configurations in the class of solutions in the sense of definition 8.3 .
since @xmath202 on @xmath196 by , condition ( c ) in definition [ regreflsoldef ] is equivalent to @xmath284 that is , @xmath188 is below @xmath279 .    in bae - chen - feldman @xcite , we have developed a mathematical approach to establish the regularity of solutions of the regular reflection - diffraction problem in the sense of definition 8.3 .
first , we have shown that any regular reflection - diffraction solutions can not be @xmath222 across the sonic arc @xmath285 .
[ noc2acrosssoniclinethm ] there does not exist a global regular reflection - diffraction solution in the sense of definition 8.3 such that @xmath270 is @xmath222 across the sonic arc @xmath126 .
now we study the one - sided regularity up to @xmath126 from the elliptic side , i.e. , from @xmath123 . for simplicity of presentation
, we now use a localized version of @xmath286 : for a given neighborhood @xmath287 of @xmath122 and @xmath288 , define @xmath289 since @xmath287 will be fixed in the following theorem , we do not specify the dependence of @xmath286 on @xmath287 .
[ mainthm1 ] let @xmath270 be a regular reflection - diffraction solution in the sense of definition [ regreflsoldef ] and satisfy the following properties : there exists a neighborhood @xmath287 of @xmath122 such that    1 .
[ c11normestimate - a ] @xmath127 is @xmath128 across the sonic arc @xmath126 : @xmath290 ; 2 .
[ psi - x - est - a ] there exists @xmath228 so that , in the coordinates ( [ coordnearsonic ] ) , @xmath291 3 .   [ fbnearsonic - a ] there exists @xmath292 , @xmath230 , and a function @xmath231 such that , in the coordinates ( [ coordnearsonic ] ) , @xmath293 and @xmath294)}<\infty,\;\;\;\ ; \frac{d\hat{f}}{dx}\!\ge\omega>0\!\;\;\ ; \text{for}\;\;0<x<{\varepsilon}_0.\ ] ]    then we have    1 .   [ c2alpsonic ] @xmath270 is @xmath295 up to @xmath126 away from @xmath296 for any @xmath297 .
that is , for any @xmath298 and any given @xmath299 , there exists @xmath300 depending only on @xmath301 , and @xmath302 so that @xmath303 2 .   [ limitsonic ] for any @xmath304 , @xmath305 3 .
[ 2ndderjumpsonic ] @xmath306 has a jump across @xmath126 : for any @xmath307 , @xmath308 4 .
[ nolimatp1 ] the limit @xmath309 does not exist .
we remark that the solutions established in @xcite satisfy the assumptions of theorem 8.5 .
in particular , we proved that the @xmath128-regularity is _ optimal _ for the solution across the open part @xmath310 of the sonic arc ( the degenerate elliptic curve ) and at the point @xmath296 where the sonic circle meets the reflected shock ( as a free boundary ) .    to achieve the optimal regularity ,
one of the main difficulties is that the sonic arc @xmath126 is the transonic boundary separating the elliptic region from the hyperbolic region , near where the solution is governed by the nonlinear degenerate elliptic equation for @xmath241 .
we carefully analyzed the features of equation and established the @xmath311 regularity of solutions in the elliptic region up to the open sonic arc @xmath310 .
as a corollary , we showed that the @xmath128-regularity is actually optimal across the transonic boundary @xmath312 from the elliptic to hyperbolic region .
since the reflected shock @xmath312 is regarded as a free boundary connecting the hyperbolic region ( 1 ) with the elliptic region @xmath123 for the nonlinear second - order equation of mixed type , another difficulty for the optimal regularity of the solution is that the point @xmath296 is exactly the point where the degenerate elliptic arc @xmath310 meets a transonic free boundary for the nonlinear partial differential equation of second order . as far as we know , this is the first optimal regularity result for solutions to a free boundary problem of nonlinear degenerate elliptic equations at the point where an elliptic degenerate curve meets the free boundary . to achieve this ,
we carefully constructed two sequences of points on where the corresponding sequences of values of @xmath313 have different limits at @xmath296 ; this has been done by employing the one - sided @xmath311 regularity of the solution up to the open arc @xmath310 and by studying detailed features of the free boundary conditions on the free boundary @xmath312 , i.e. , the rankine - hugoniot conditions .
we remark that some efforts were also made mathematically for the reflection - diffraction problem via simplified models .
one of these models , the unsteady transonic small - disturbance ( utsd ) equation , was derived and used in keller - blank @xcite , hunter - keller @xcite , hunter @xcite , and morawetz @xcite for asymptotic analysis of shock reflection - diffraction . also see zheng @xcite for the pressure gradient equation and canic - keyfitz - kim @xcite for the utsd equation and the nonlinear wave system .
furthermore , in order to deal with the reflection - diffraction problem , some asymptotic methods have been also developed .
lighthill @xcite studied shock reflection - diffraction under the assumption that the wedge angle is either very small or close to @xmath66 .
keller - blank @xcite , hunter - keller @xcite , harabetian @xcite , and gamba - rosales - tabak @xcite considered the problem under the assumption that the shock is so weak that its motion can be approximated by an acoustic wave . for a weak incident shock and a wedge with small angle for potential flow , by taking the jump of the incident shock as a small parameter , the nature of the shock reflection - diffraction pattern was explored in morawetz @xcite by a number of different scalings , a study of mixed equations , and matching the asymptotics for the different scalings .
also see chen @xcite for a linear approximation of shock reflection - diffraction when the wedge angle is close to @xmath66 and serre @xcite for an apriori analysis of solutions of shock reflection - diffraction and related discussions in the context of the euler equations for isentropic and adiabatic fluids .
another related recent effort has been on various important physical problems in steady potential flow , as well as steady fully euler flow , for which great progress has been made .
the problems for global subsonic flow past an obstacle and for local supersonic flow past an obstacle with sharp head are classical , due to the works of shiffman @xcite , bers @xcite , finn - gilbarg @xcite , dong @xcite , and others for the pure elliptic case and to the works of gu @xcite , shaeffer @xcite , li @xcite , and others for the pure hyperbolic case .
recent progress has been made on transonic flow past nozzles ( e.g. @xcite ) , transonic flow past a wedge or conical body ( e.g. @xcite ) , and transonic flow past a smooth obstacle ( e.g. @xcite ) .
also see @xcite for the existence of global subsonic - sonic flow , @xcite for global supersonic flow past an obstacle with sharp head , and the references cited therein .
for some of other recent related developments , we refer the reader to chen @xcite for the local stability of mach configuration , elling - liu @xcite for physicality of weak prandtl - meyer reflection for supersonic potential flow around a ramp , serre @xcite for multidimensional shock interaction for a chaplygin gas , canic - keyfitz - kim @xcite for semi - global solutions for the shock reflection problem , glimm - ji - li - li - zhang - zhang - zheng @xcite for the formation of a transonic shock in a rarefaction riemann problem for polytropic gases , zheng @xcite for various solutions to some two - dimensional riemann problems , gues - mtivier - williams - zumbrun @xcite and benzoni - gavage and serre @xcite for revisits of the local stability of multidimensional shocks and phase boundaries , among many others .
as we have seen from the previous discussion , the shock reflection - diffraction problem involves several core challenging difficulties that we have to face in the study of nonlinear partial differential equations .
these nonlinear difficulties include _ free boundary problems , oblique derivative problems for nonsmooth domains , degenerate elliptic equations , degenerate hyperbolic equations , transport equations with rough coefficients , mixed and/or composite equations of hyperbolic - elliptic type , behavior of solutions when a free boundary meets an elliptic degenerate curve , and compressible vortex sheets_. more efficient numerical methods are also required for further understanding of shock reflection - diffraction phenomena .    furthermore , the wave patterns of shock reflection - diffraction configurations are the core patterns and configurations for the global solutions of the two - dimensional riemann problem ; these solutions are building blocks and local structure of general entropy solutions and determine global attractors and asymptotic states of entropy solutions , as time goes infinity , for two - dimensional systems of hyperbolic conservation laws .
therefore , a successful solution to the shock reflection - diffraction problem not only provides our understanding of shock reflection - diffraction phenomena and behavior of entropy solutions to multidimensional conservation laws , but also provides important new ideas , insights , techniques , and approaches for our developments of more efficient analytical techniques and methods to overcome the core challenging difficulties in multidimensional problems in conservation laws and other areas in nonlinear partial differential equations .
the shock reflection - diffraction problem is also an excellent test problem to examine our capacity and ability to solve rigorously various challenging problems for nonlinear partial differential equations and related applications .
bae , g .- q .
chen , and m. feldman , _ regularity of solutions to regular shock reflection for potential flow _ , inventiones mathematicae , 2009 ( to appear ) ; doi 10.1007/s00222 - 008 - 0156 - 4 ( online ) ; preprint arxiv:0804.2500 , april 2008 .
l.  a. caffarelli , _ a harnack inequality approach to the regularity of free boundaries , i. lipschitz free boundaries are @xmath314 _ , rev . mat .
iberoamericana , * 3 * ( 1987 ) , 139162 ; _ ii .
flat free boundaries are lipschitz _ , comm .
pure appl . math .
* 42 * ( 1989 ) , 5578 ; _ iii .
existence theory , compactness , and dependence on @xmath315 _ , ann . scuola norm .
sci . ( 4 ) , * 15 * ( 1989 ) , 583602 .
t. chang , g .- q .
chen , and s. yang , _
@xmath316-d riemann problem in gas dynamics and formation of spiral _ , in : nonlinear problems in engineering and science  numerical and analytical approach ( beijing , 1991 ) , 167179 , science press : beijing , 1992 .
t. chang , g .- q .
chen , and s. yang , _ on the riemann problem for two - dimensional euler equations i : interaction of shocks and rarefaction waves ; ii : interaction of contact discontinuities _ , discrete contin .
systems , * 1 * ( 1995 ) , 555584 ; * 6 * ( 2000 ) , 419430 .
chen , _ euler - poisson - darboux equations and hyperbolic conservation laws _ , in : nonlinear evolutionary partial differential equations , pp .
1125 , ams / ip studies in advanced mathematics , vol .
3 , international press : cambridge , massachusetts , * 1997*.    g .- q . chen , _ euler equations and related hyperbolic conservation laws _ , chapter 1 , handbook of differential equations , evolutionary equations , vol .
2 , eds . c. m. dafermos and e. feireisl , elsevier : amsterdam , the netherlands , 2005 .
chen , _ on nonlinear partial differential equations of mixed type _ , in : proceedings of the 4th international congress of chinese mathematicians , eds .
l. jin , k .- f .
liu , l. yang , and s .- t .
yau , isbn : 7 - 04 - 023267 - 7 , higher education press : beijing , jan
. 2008 .
a. yu . demyanov and a.  v. panasenko , _ numerical solution to the problem of the diffraction of a plane shock wave by a convex corner _ , fluid dynamics , * 16 * ( 1981 ) , 720725 ( translated from the original russian ) .                    i.  m. gamba and c.  s. morawetz , _ a viscous approximation for a @xmath316-d steady semiconductor or transonic gas dynamic flow : existence theorem for potential flow _ , comm .
pure appl . math .
* 49 * ( 1996 ) , 9991049 .
h.  m. glaz , p. colella , i.  i. glass , and r.  l. deschambault , _ a detailed numerical , graphical , and experimental study of oblique shcok wave reflection _ ,
lawrence berkeley laboratory report * lbl-20033 * ( 1985 ) .
j. glimm , x .-
ji , j .- q . li , x .-
li , p. zhang , t. zhang , and y .- x .
zheng , _ transonic shock formation in a rarefaction riemann problem for the 2-d compressible euler equations _ ,
siam j. appl . math .
2008 ( in press ) .
m.  s. ivanov , d. vandromme , v.  m. formin , a.  n. kudryavtsev , a. hadjadj , and d.  v. khotyanovsky , _ transition between regular and mach reflection of shock waves : new numerical and experimental results _
, shock waves , * 11 * ( 2001 ) , 199207 .                                          a. majda , _ one perspective on open problems in multi - dimensional conservation laws _ , in : multidimensional hyperbolic problems and computations ( minneapolis , mn , 1989 ) , 217238 , i m a vol . math . appl . *
29 * , springer : new york , 1991 .
p. woodward , _ simulation of the kelvin - helmholtz instability of a supersonic slip surface with the piecewise - parabolic method ( ppm ) , _ in : numerical methods for the euler equations of fluid dynamics ( rocquencourt , 1983 ) , 493508 , siam : philadelphia , pa , 1985 .