• 천문학논총
• /
• 제25권4호
• /
• pp.177-186
• /
• 2010

We developed an efficient Monte-Carlo algorithm to solve dust-scattering radiative transfer problems for continuum radiation. The method calculates the scattered intensities for various anisotropic factors ($g_i$) all at once, while actual photon packets are tracked following a scattering phase function given by a single anisotropic factor ($g_0$). The algorithm was tested by applying the method to a dust cloud embedding a star at the cloud center and found to provide accurate results within the statistical fluctuation that is intrinsic in Monte-Carlo simulations. It was found that adopting $g_0$ = 0.4 - 0.5 in the algorithm is most efficient. The method would be efficient in estimating the anisotropic factor of the interstellar dust by comparing the observed data with radiative transfer models.

• 천문학회지
• /
• 제48권1호
• /
• pp.57-66
• /
• 2015

Monte-Carlo radiative models of the diffuse Galactic light (DGL) in our Galaxy are calculated using the dust radiative transfer code MoCafe, which is three-dimensional and takes full account of multiple scattering. The code is recently updated to use a fast voxel traversal algorithm, which has dramatically increased the computing speed. The radiative transfer models are calculated with the generally accepted dust scale-height of 0.1 kpc. The stellar scale-heights are assumed to be 0.1 or 0.35 kpc, appropriate for far-ultraviolet (FUV) and optical wavelengths, respectively. The face-on optical depth, measured perpendicular to the Galactic plane, is also varied from 0.2 to 0.6, suitable to the optical to FUV wavelengths, respectively. We find that the DGL at high Galactic latitudes is mostly due to backward or large-angle scattering of starlight originating from the local stars within a radial distance of r < 0.5 kpc from the Earth. On the other hand, the DGL measured in the Galactic plane is mostly due to stars at a distance range that corresponds to an optical depth of $${\sim_\sim}$$ 1 measured from the Earth. Therefore, the low-latitude DGL at the FUV wavelength band would be mostly caused by the stars located at a distance of $r{\leq}0.5$ kpc and the optical DGL near the Galactic plane mainly originates from stars within a distance range of $1{\leq}r{\leq}2kpc$. We also calculate the radiative transfer models in a clumpy two-phase medium. The clumpy two-phase models provide lower intensities at high Galactic latitudes compared to the uniform density models, because of the lower effective optical depth in clumpy media. However, no significant difference in the intensity at the Galactic plane is found.

• 천문학회지
• /
• 제33권1호
• /
• pp.29-36
• /
• 2000

We developed a Monte Carlo code that describes the resonant Ly$\alpha$ line transfer in an optically thick, dusty, and static medium. The code was tested against the analytic solution derived by Neufeld (1990). We explain the line transfer mechanism by tracing histories of photons in the medium. We find that photons experiences a series of wing scatterings at the moment of thier escape from the medium, during which polarization may develop. We examined the amount of dust extinction for a wide range of dust abundances, which are compared with the analytic solution. Brief discussions on the astrophysical application of our work are presented.

• 천문학회지
• /
• 제46권6호
• /
• pp.235-242
• /
• 2013

We compare mass-loss rates of OH/IR stars obtained from radio observations with those derived from the dust radiative transfer models and IR observations. We collect radio observational data of OH maser and CO line emission sources for a sample of 1533 OH/IR stars listed in Suh & Kwon (2011). For 1259 OH maser, 76 CO(J=1-0), and 55 CO(J=2-1) emission sources, we compile data of the expansion velocity and mass-loss rate. We use a dust radiative transfer model for the dust shell to calculate the mass-loss rate as well as the IR color indices. The observed mass-loss rates are in the range predicted by the theoretical dust shell models corresponding to $\dot{M}=10^{-8}M_{\odot}/yr-10^{-4}M_{\odot}/yr$. We find that the dust model using a simple mixture of amorphous silicate and amorphous $Al_2O_3$ (20% by mass) grains can explain the observations fairly well. The results indicate that the dust radiative transfer models for IR observations generally agree with the radio observations. For high mass-loss rate OH/IR stars, the mass-loss rates obtained from radio observations are underestimated compared to the mass-loss rates derived from the dust shell models. This could be because photon momentum transfer to the gas shell is not possible for the physical condition of high mass-loss rates. Alternative explanations could be the effects of different dust-to-gas ratios and/or a superwind.

• 한국연소학회지
• /
• 제2권1호
• /
• pp.1-8
• /
• 1997

Premixed combustion within porous ceramic media is numerically studied to understand burning characteristics and to find best configurations for burner implementations. Among many parameters, critical to burner performance, flame location and extinction coefficient are selected as major parameters for this study. The flame structure and burner performance with respect to these two parameters are observed. In the study, it is found that the location of flame is the most important in porous burner operation since it affects the rate of heat transfer and flame structure. Stability of the flame within the porous ceramic burner is discussed with respect to the flame location. It is found that to obtain high radiative output, the flame should be located downstream section of the burner. But the flame is to be unstable at most of downstream section except near the exit plane. To overcome this problem, new porous ceramic burner, using different ceramic properties in one burner instead of single property ceramic, is made and tested. With a combination of ceramics of high extinction coefficient at upstream and another material of low extinction coefficient at downstream of the burner, the flame can be stabilized at wider region of the burner with higher radiative output compared to the original burner configuration.

• 천문학논총
• /
• 제25권4호
• /
• pp.113-118
• /
• 2010

The Henyey-Greenstein (H-G) phase function, which is characterized by a single parameter, has been generally used to approximate the realistic dust-scattering phase function in investigating scattering properties of the interstellar dust. Draine (2003) proposed a new analytic phase function with two parameters and showed that the realistic phase function is better represented by his phase function. If the H-G and Draine's phase functions are significantly different, using the H-G phase function in radiative transfer models may lead to wrong conclusions about the dust-scattering properties. Here, we investigate whether the H-G and Draine's phase functions would indeed produce significant differences in radiative transfer calculations for two simple configurations. For the uniformly distributed dust with an illuminating star at the center, no significant difference is found. However, up to ~ 20% of difference is found when the central star is surrounded by a spherical-shell dust medium and the radiation of $\lambda$ < $2000\;{\AA}$ is considered. It would mean that the investigation of dust-scattering properties using the H-G phase function may produce errors of up to ~ 20% depending on the geometry of dust medium and the radiation wavelength. This amount of uncertainty would be, however, unavoidable since the configurations of dust density and radiation sources are only approximately available.

• 천문학논총
• /
• 제23권2호
• /
• pp.25-29
• /
• 2008

Scattering of incident light by the interstellar dust is usually approximated by Henyey-Greenstein scattering phase function. Recently, Draine (2003) proposed a new analytic phase function with two parameters. We describe an algorithm to generate random numbers distributed according to the Draine’s function, and compare two phase functions. It is also derived exact solutions of two parameters for given values ${\langle}cos{\theta}{\rangle}$ and ${\langle}cos^2{\theta}{\rangle}$. It is found that Henyey-Greenstein function with g = ${\langle}cos{\theta}{\rangle}$ provides a good approximation for ${\lambda}\;>\;2000{\AA}$. At shorter wavelengths, more realistic phase function may be needed for radiative transfer models.

• 천문학논총
• /
• 제23권2호
• /
• pp.31-35
• /
• 2008

We present analytical approximations for calculating the scattering and escape of non-ionizing photons from a plane-parallel medium with uniformly illuminated by external sources. We compare the results with the case of a spherical dust cloud. It is found that more scattering and absorption occur in the plane-parallel geometry than in the spherical geometry when the optical depth perpendicular to the plane and the radial optical depth of the sphere are the same. The results can provide an approximate way to estimate radiative transfer in a variety interstellar conditions and can be applied to the dust-scattered diffuse Galactic light.

• Journal of Astronomy and Space Sciences
• /
• 제25권4호
• /
• pp.329-334
• /
• 2008

We investigate the spectral energy distributions (SEDs) of Massive Young Stellar Objects (MYSOs) using the various infrared observational data including the Infrared Space Observatory (ISO) data. We model the dust envelopes around the stars using a radiative transfer model for spherically symmetric geometry. Comparing the model results with the observed SEDs of the two MYSOs (AFGL 4176 and AFGL 2591), we derive the relevant dust shell parameters including the dust opacity, the dust density distribution, and dust temperature distribution. We find that the spherical model can produce the SEDs roughly similar to the observations. We expect that the results would be helpful for making more realistic non-spherical dust envelope models for MYSOs.

• 천문학회지
• /
• 제42권4호
• /
• pp.81-91
• /
• 2009

We make a new catalog of AGB stars in our Galaxy from the sources listed in the Infrared Astronomical Satellite (IRAS) point source catalog (PSC) compiling the lists of previous works with verifying processes. We verify the class identification of AGB stars into oxygen-rich and carbon-rich stars using the information from recent investigations. For the large sample of AGB stars, we present infrared two-color diagrams from the observations at near infrared bands and IRAS PSC. On the two-color diagrams, we plot the tracks of theoretical radiative transfer model results with increasing dust shell optical depths. Comparing the observations with the theoretical tracks, we discuss the meaning of the infrared two-color diagrams.