본 연구에서는 항아리 형태 젊은 초신성 잔해의 동력학을 설명하기 위해 분출물의 동력학적 효과를 고려하였다. 분출물과 성간 물질 사이에 존재하는 접촉불연속면에서 레일라-테일러(R-T)불안정에 기인한 자기장의 증폭호과가 고려되었다. 우주선 입자의 가속을 가정함으로서 합성전파 모형을 만들 수 있었으며 관측과의 비교를 위해 방위각을 따른 전파세기의 비(A)를 계산하였다. R-T불안정의 결과로 자란 가지의 경계면에서 압축, 전단, 인장의 결과로 충격파 후면의 자기장은 증폭되었다. 불안정의 시간에 따른 변화는 분출물의 감속에 민감하게 의존하며 초신성 잔해의 진화와 밀접히 관계됨을 볼 수 있었다. 자기장의 세기는 초기에 급격히 증가하며 역 충격파가 분출물의 등밀도지역으로 들어감에 따라 감소되었다. 그러나 이와 같은 초기 자기장 증폭의 효과는 초신성 잔해의 후기까지 남아 있음을 볼 수 있었다. 증폭된 자기장 영역에서 적도지역과 극지역의 자기장의 세기의 비는 최대 7.5까지 이를 수 있었다. 이와 같은 자기장의 증폭은 방위각에 따라 매우 큰 전파세기의 비를 만들 수 있었다(A=15). 증폭된 자기장의 세기는 수치계산의 분해능에 매우 민감함을 알 수 있었다. 본 연구에서는 우주선 입자의 가속효과가 자기장과 충격파 면이 이루는 각도에 의존한다는 가정 없이도 자기장의 증폭효과가 관측된 항아리 형태의 젊은 초신성 잔해를 만들 수 있음을 보였다. 그러나 이와 같은 기작이 효과적이기 위해서는 초신성 잔해 주위의 자기장이 잘 정열되어 있어야 한다. 항아리 형태의 젊은 초신성잔해의 수가 적게 관측되는 것은 이와 같은 조건이 성간에서 잘 이루어지지 않음을 의미한다.

We present the results of one-dimensional numerical simulations of SNR evolution in the in­homogeneous medium considering the effects of the evaporation of the cloud and the thermal conduction. We have included the effects of changing evaporation rate as a function of cloud size and the ambient temperature so that the clouds could be evaporated completely before they reach the center of the SNR. The heat conduction markedly changes the density distribution in the remnant interior. To explain the observed morphologies of the centrally peaked X-ray SNRs(for example W44), the maximal thermal conduction is required. However, this is unlikely due to the magnetic field and the turbulent motion. The effects of the evaporation of the cloud and the thermal conduction described here may explain the class of remnants observed to have centrally peaked X-ray emmision.

The structure and environments of the molecular clouds near the SNR HB3(G132.7±1.3) HB3(G132.7±1.3) are studied. The molecular complex which is located at the southern rim of HB3 was proposed by former investigators as the one interacting with HB3. This complex region of 2∘×2∘atl=133∘ 2∘×2∘atl=133∘ has been observed at 12CO,13CO,J=1−0ata1′ 12CO,13CO,J=1−0ata1′ , resolution with the 14-m radio telescope at Taeduk Radio Astronomy Observatory. We have reached to the following four conclusions. The possibility that these molecular complex and HB3 are interacting with each other cannot be supported with any of our data. The morphologies of the two show no similarities. Neither particular features for the interaction are found in the CO lines. The hypothetical 'Molecular wall' which was expected to exist on the northwestern rim of HB3 as a cause for the noncircular morphology of HB3 is turned out to be nonexistent in CO. The molecular complex which resembles a 'bar' at a low resolution is now resolved into a U-shaped shell. It seems that the U-shape is consist of two independent components. No peculiarities, such as unseen masses or bright stars capable of forming HlI regions, are found within the U-shape region. The total mass included in the complex is estimated to be Mtotal=2.9\~8.4×105M⨀ Mtotal=2.9\~8.4×105M⨀ , which is in good agreement with previous observations within errors. Considering about 12 clumps distinguishable within the complex, the total mass implies that masses of each of clumps are on the order of 104M⨀ 104M⨀ , which makes these good objects for further studies in relation to star-formation. Especially the clumps associated with W3 are worthy for more high resolution observations for better understanding of astrophysical phenomenon ongoing in them.

We have simulated the interaction of supernova remnants with constant ambient medium to explore the dynamics of Type Ia supernova remnant. We assumed the supernova ejecta density distribution of the central constant and the outer power-law density distribution(ρ∝γ−n) (ρ∝γ−n) . We have calculated four different cases with different n. By scaling the length and time scales from the initial parameters-ejecta mass, ejecta energy, the ambient density, we could compare effects of the different density distribution of the ejecta on the dynamics of the SNRs. The radius of the outer forward shock converges the Sedov-Talyor solution at t' = 2.3 when the swept-up mass is 8 times of the ejecta mass. On the other hand, the motion of the reverse shock are largely affected by n. The ejecta with smaller n takes comparably long time to thermalize the whole ejecta at t′≃5.3,Msw≃18Mej t′≃5.3,Msw≃18Mej . We have applied our calculated results to obtain the ejecta density distributions of Tycho and SN1006 with n≃6 n≃6 .

We have developed a spherical FCT code in order to simulate the interaction of supernova remnants with stellar wind bubbles. We assume that the density profile of the supernova ejecta follows the Chevalier mode1(1982) where the outer portion has a power-law density distribution(ρ∝γ−n ρ∝γ−n ) and the SN ejecta has a kinetic energy of 1051 1051 ergs. The structure of wind bubble has been calculated with the stellar mass loss rate ˙ M=5×10−6M⊙/yr M˙=5×10−6M⊙/yr and the wind velocity υ=2×103 υ=2×103 km/s We have simulated seven models with different initial conditions In the first two models we computed the evolution of SNRs with n=7 and n=14 in the uniform medium The numerical results agree with the Chevalier's similarity solution at early times. When all of the power-law portion of the ejecta is swept up by the reverse shock, the evolution slowly converges to the Sedov-Taylor stage. There is not much difference between the two cases with different n's The other five models simulate SNRs produced inside wind bubbles. In model III, we consider the SN ejecta of 1.4 M⊙ M⊙ and the radius of bubble ~2.76 pc so that ratio of the mass α(=MW.S/Mej α(=MW.S/Mej is 2. We follow the complex hydrodynamic flows produced by the interaction of SN shocks with stellar shocks and with the contact discontinuities, In the model III, the time scale for the SN shock to cross the wind shell τcross τcross is similar to the time scale for the reverse shock to sweep the power-law density profile τbend τbend . Hence the SN shock crosses the wind shell. At late times SN shock produces another shell in the ambient medium so that we have a SNR with double shell structure. From the numerical results of the remaining models, we have found that when τcross/τbend≤2 τcross/τbend≤2 , or equivalently when α≤50 α≤50 , the SNRs produced inside wind bubbles have double shell structure. Otherwise, either the SN shock does not cross the wind shell or even if it crosses at one time, the reverse shock reflected at the center accelerates the wind shell to merge into the SN shock Our results confirm the conclusion of Tenorio-Tagle et a1(1990).

We carried out high-resolution(FWHM=3' .3) HI 21 cm observations of the supernova remnant(SNR) PKS0607+17 and HII region S261 using Arecibo 305-m telescope. The observation was to investigate whether the high-velocity(HV) gas detected in the southern area of PKS0607+17 by Koo & Heiles(1991) is physically associated with the SNR or not. The velocity of the HV gas ranges from +64 km/s to +87 km/s, which is difficult to result from the Galactic rotation. The HV gas could be the gas accelerated by supernova blast wave. However, because the observation of Koo and Heiles(1991) was carried out using Hat Creek radio telescope(FWHM ≃ ≃ 36'), the association of the HV gas with the SNR could not be investigated. Using the Arecibo HI 21cm data, we have found that the HV gas appears m the southern part of the SNR and its velocity ranges from +61 km/s to +77 km/s. But the HV gas is scattered m the whole field, not only toward PKS0607+17 but also outside the SNR Accordingly the HV gas is probably not associated with the SNR, but is accidentally aligned along the same line of sight toward the SNR. Instead we have found that HI clouds at low velocities could be possibly associated with the SNR. In Arecibo HI 21cm channel maps the HI gas seems to surround the southern boundary of the SNR at VLSR VLSR = +19.6 ~ +40.2 km/s. But because the region of the Arecibo HI 21cm observation is not wide enough to examine the HI gas distribution, we investigated this area using the Berkely low-latitude HI survey data(Weaver & Williams 1974) too. There we found HI gas surrounding the radio continuum boundary of PKS0607+17 at VLSR VLSR = +21.6 ~ +258 km/s. It is possible that this HI gas is associated with the SNR, in which case, the velocity of the SNR Vo Vo ≃ ≃ +26 km/s, its distance d ≃ ≃ 12.5 kpc and its radius R ≃ ≃ 145 pc. If we assume that the expansion velocity is ~10 km/s, then the age of the SNR is ∼4.4×106 ∼4.4×106 years. PKS0607+17 could be one of the oldest SNRs in the Galaxy. We also studied HI propertities of the HII region S261, which is ∼1∘ ∼1∘ away from PKS0607+17. There has been no high-resolution m 21 cm observational study on S261. We discovered HI cloud located at the north-eastern part of S261 at VLSR VLSR = +5 km/s ~ +10 km/s, which is possibly associated with the HII region. The central velocity of the HI cloud VLSR VLSR = +7.2 km/s and the corresponding distance d = 1.5 kpc. This velocity is comparable to the radio recombination line velocities.

We have performed the high resolution computer simulation with 1D spherical hydrodynamic code in order to study the dynamical evolution of supernova ejecta interacting with a pre-existing fast wind structure. The fast wind structure has been calculated with Min=3×10−6M⊙yr−1 Min=3×10−6M⊙yr−1 and υin=1000km/sec υin=1000km/sec , which velocity is higher than the critical velocity relating to the initial radiative cooling. The fast wind becomes initially adiabatic. After a shell formation time of ∼4000yrs ∼4000yrs , the wind becomes radiative cooling at the shell zone, forming a thin dense radiative shell and an adiabatic wind bubble afterward. When supernova explodes in the wind center at 20,000yrs after the wind evolves, the supernova ejecta, which has a dense distribution of ρ∝r−n ρ∝r−n (here we have n = 9), interacts initially with, the understood wind zone, producing forward and reverse shocks. The reverse shock heats the supernova ejecta and its temperature increases. In this study, as the mass of the supernova ejecta is larger than that of the wind shell (Mej=5M⊙ Mej=5M⊙ , Msw=2M⊙ Msw=2M⊙ ), we can conform two shell structures: an outer shell by the supernova ejecta and a secondarily shocked wind shell by it. The secondarily shocked wind shell should accelerates in this case to be R-T unstable, consequently producing the knots

It has been recognized that the morphologies of the SNRs from the radio observation are "barrel shaped". To interpret the mechanism of the radiation and the physical state of the environments, we have analytically calculated the dynamical structure of the interacting region in the case where the ejectum has a steep power-law density profile(ρ∼r−n ρ∼r−n ) and the ambient medium has a shallow power-law density profile(ρ∼r−s ρ∼r−s ), assuming that the cosmic rays are isotropically accelerated in the shock wave and the magnetic fields are very weak. The calculated synchrotron radio maps show that the emission from the equator is intense and the emissions from the central and polar regions are less intense. Also the thicknesses of the shell are strongly dependent on s and weakly on n. The azimuthal intensity ratio α α increases as the efficiency of the cosmic ray acceleration increases and s decreases. We compared the results with the morphology of the SNR A. D. 1006(type I SNR). It does agree with the case of s = 0, w = 0.3 - 0.5. This value for w is consistent with the results by Eichler(1979). It provides us the evidence of the cosmic ray acceleration in the shock wave.