• Publications of The Korean Astronomical Society
• /
• v.25 no.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.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.4
• /
• pp.1061-1066
• /
• 2014

Time-dependent formulations of the reactive scattering theory based on the wavepacket correlation functions with the M${\phi}$ller wavepackets for the electronically nonadiabatic reactions are presented. The calculations of state-to-state reactive probabilities for the quasi-Jahn-Teller scattering model system were performed. The conical intersection (CI) effects are investigated by comparing the results of the two-surface nonadiabatic calculations and the single surface adiabatic approximation. It was found that the results of the two-surface nonadiabatic calculations show interesting features in the reaction probability due to the conical intersection. Single surface adiabatic calculations with extended Born-Oppenheimer approximation using simple wavepacket phase factor was found to be able to reproduce the CI effect semi-quantitatively, while the single surface calculations with the usual adiabatic approximation cannot describe the scattering process for the Jahn-Teller model correctly.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.2
• /
• pp.291-301
• /
• 1999

A method for the synthesis of one-dimensional nonlinear distribution function is presented for the desired inverse scattering pattern. This method is based on the inverse transform of the solution of the Riccati equation derived from one-dimensional inverse scattering problem. Since the solution is analogous to the array factor or normalized space factor in collinear array antenna, the synthesis method for field pattern is applied for the construction of the involved line-source nonlinear distribution function. The suggested method is carried out under the optimization process, and is numerically verified by synthesizing the dispersive transmission line profile within the specified frequency band and control of scattered field on resistive strip.

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.324-329
• /
• 2023

Different types of deterministic solution methods were used to solve neutron transport equations corresponding to half-space and slab albedo problems. In these types of solution methods, in addition to the error of the numerical solutions, the obtained results contain truncation and discretization errors. In the present work, a non-analog Monte Carlo method is provided to simulate the half-space and slab albedo problems with Inönü and Anlı-Güngör strongly anisotropic scattering functions. For each scattering function, the sampling method of the direction of the scattered neutrons is presented. The effects of different beams with different angular dependencies and the effects of different scattering parameters on the reflection probability are investigated using the developed Monte Carlo method. The validity of the Monte Carlo method is also confirmed through the comparison with the published data.

• Atmosphere
• /
• v.23 no.1
• /
• pp.39-45
• /
• 2013

A parameterization for the scattering of longwave radiation by ice clouds has been developed based on spectral scattering property calculations with shapes and sizes of ice crystals. For this parameterization, the size distribution data by Fu (1996) and by Michell and Arnott (1994) are used. The shapes of ice crystal considered in this study are plate, solid column, hollow column, bullet-rosette, droxtal, aggregate, and spheroid. The properties of longwave scattering by ice crystals are presented as a function of the extinction coefficient, single-scattering albedo, and asymmetry factor. The heating rate and flux by the radiative parameterization model are calculated for wide range of ice crystal sizes, shapes, and optical thickness. The results are compared with the calculated results using a six-stream discrete ordinate scattering algorithm and Chou's method. The new method (with various habits and size distributions) provides a good simulation of the scattering properties and cooling rate in optically thin clouds (optical thickness < 5). Depending on the inclusion of scattering by ice clouds, the errors in the calculation of the cooling rates are significantly different.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 1999.11a
• /
• pp.552-554
• /
• 1999

Phonon scattering and impact ionization models have been presented to analyze hot carrier transport in high energy region, using full band model and Fermi's golden rule. We have investigated temperature dependent properties for impact ionization process of Si using realistic energy band structures at 77K and look. The realistic full band model, obtained from the empirical pseudopotential method with local from factors, is used to calculate scattering rate. The accurate calculation of impact ionization rate requires the use of a wavevector- and frequency-dependent dielectric function ξ ( q,$\omega$). The empirical phonon scattering rate P$\sub$ph/, is given by deriving from linear function for P$\sub$ph/ versus D(E) since the phonon scattering rate is linearly depended on density of states D(E). Impact ionization rate p,, is calculated from the first principle's theory. and fitted by modified Keldysh formula having power of above 2.

• Journal of The Korean Astronomical Society
• /
• v.51 no.4
• /
• pp.111-117
• /
• 2018

A numerical method is proposed to calculate the response of detectors measuring particle energies from incident isotropic fluxes of electrons and positive ions. The isotropic flux is generated by injecting particles moving radially inward on a hypothetical, spherical surface encompassing the detectors. A geometric projection of the field-of-view from the detectors onto the spherical surface allows for the identification of initial positions and momenta corresponding to the clear field-of-view of the detectors. The contamination of detector responses by particles penetrating through, or scattering off, the structure is also similarly identified by tracing the initial positions and momenta of the detected particles. The relative contribution from the contaminating particles is calculated using GEANT4 to obtain the geometric factor of the instrument as a function of the energy. This calculation clearly shows that the geometric factor is a strong function of incident particle energies. The current investigation provides a simple and decisive method to analyze the instrument geometric factor, which is a complicated function of contributions from the anticipated field-of-view particles, together with penetrating or scattered particles.

• Journal of the Korean Society of Radiology
• /
• v.3 no.3
• /
• pp.23-25
• /
• 2009

Gd-157 material has very large neutron capture cross section in the thermal region. So it is very useful to shield material for thermal neutrons. Futhermore, in the neutron capture experiment and calculation, the neutron absorption and scattering are very important. Especially these effects are conspicuous in the resonance energy region and below the thermal energy region. In the case of very narrow resonance, the effect of scattering is to be more considerable factor. In the present study, we obtained energy dependent neutron absorption ratios of natural indium in energy region from 0.003 to 100 keV by MCNP-4B Code. The coefficients for neutron absorption was calculated for circular type and 1 mm thickness. In the lower energy region, neutron absorption is larger than higher region, because of large capture cross section (1/v). Furthermore it seems very different neutron absorption in the large resonance energy region. These results are very useful to decide the thickness of sample and shielding materials.

• Economic and Environmental Geology
• /
• v.32 no.4
• /
• pp.379-384
• /
• 1999

In order to know the characteristics of attenuation of coda wave in the Kyungsang Sedimetary Basin, quality faclity factor for coda wave (coda Q) is estimated from the earthquake data recorded in the KIGAM local seismic network. Thesingle scattering model for coda wave generation is adopted is adopted in estimating coda Q. In the present study, coda Q(Qc)is estimated in the range of $\alpha$=1.5~3.0, where $\alpha$ denotes the normalized time to S-wave travel time and expressed in terms of frequency (f). The deduced function in the range of 1 to 25 Hz is Qc=36.8283$f^{1.15095}$ which represents the strong dependence of coda Q on frequency. It is found that the difference of Qc between U-D, N-S, and E-W components is negligible. This fact suports the back scattering therory that coda wave originates from scattered waves by randomly distributed heterogenities in the crust On the other hand, it is observed that the coda Q increases with increasing epicentral distence. This observation suggests that QC increases with depth.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1999.04a
• /
• pp.35-40
• /
• 1999

In order to know the characteristics of attenuation of coda wave in the Kyungsang Sedimentary Basin quality factor for coda wave or coda Q is estimated from the earthquake data recorded in the KIGAM microearthquake network. The single scattering model for coda wave generation is adopted in estimating coda Q. Coda Q appears to be largely dependent on the normalized time(a) which is the ratio of elapsed time to S-wave travel time. In the present study coda Q(Qc) is estimated in the range of a=1.5-3.Q and expressed in terms of frequency(f). The deduced function in the range of 1 to 25 Hz is Qc=36.8283 f1.15095 to represent the strong dependence of coda Q on frequency. It is found that the difference of Qc between U-D N-S and E-W components is negligible, This face supports the back-scattering theory that coda were originates from scattered waves by randomly distributed heterogeneities in the crust. On the other hand it is observed that the coda Q increases with depth.