• 천문학논총
• /
• 제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.

• Bulletin of the Korean Chemical Society
• /
• 제35권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.

• 한국전자파학회논문지
• /
• 제10권2호
• /
• pp.291-301
• /
• 1999

원하는 역산란 패턴을 만족시키는 일차원 비선형 분포함수를 합성하는 방법을 제안하였다. 이 방법은 일차 원 역산란 문제로부터 유도된 Riccati 방정식의 해의 역변환에 기본을 두고 있다. 그 해는 collinear 배열 안테 나에서의 배열인자 혹은 정규화된 공간인자와 유사한 형태를 가지므로 연관된 비선형 분포함수를 구하기 위 해 line-source 법에 입각한 패턴 합성법을 사용하였다. 이 과정들은 최적화 기법에 의해 수행되며, 분산특성 을 갖는 매질내에서 임의의 전송선로를 특정 주파수 대역내에서 합성하는 예와 저항성 스트립의 산란패턴 제어 등의 예로부터 제안한 방법의 타당성을 수치적으로 입증하였다.

• Nuclear Engineering and Technology
• /
• 제55권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.

• 대기
• /
• 제23권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.

• 한국정보통신학회:학술대회논문집
• /
• 한국해양정보통신학회 1999년도 추계종합학술대회
• /
• pp.552-554
• /
• 1999

포논산란과 임팩트이온화 모델은 풀밴드 모델과 페르미 황금법칙을 이용한 고에너지 영역에서 hot carrier 전송을 해석하기 위해서 제시되어 왔다. 본 연구에서는 77K와 300k에서 실제에너지 밴드 구조를 이용한 Si의 임팩트이온화 과정에 대해서 온도 의존성을 조사했다. 풀밴드 모델은 local form factors을 이용한 의사포텐셜방법에서 얻어지고, 산란율 계산에 이용된다. 정확한 임팩트이온화율 계산은 파동벡터와 주파수에 의존하는 유전함수를 필요로 한다. 포논 산란율은 상태밀도에 선형적으로 의존하기 때문에 포논산란율과 상태밀도의 비에 대한 선형함수 유도에 의해서 구해진다 임팩트이온화율 $P_{ii}$ 는 first principle's theory로부터 계산되어지고, 수정된 Keldysh 공식에 의해서 최적화되었다.

• 천문학회지
• /
• 제51권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.

• 한국방사선학회논문지
• /
• 제3권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.

• 자원환경지질
• /
• 제32권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.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 1999년도 춘계 학술발표회 논문집 Proceedings of EESK Conference-Spring
• /
• 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.