• Journal of Radiation Protection and Research
• /
• v.47 no.1
• /
• pp.50-57
• /
• 2022

Background: There are several proton therapy facilities in operation or planned in Taiwan, and these facilities are anticipated to not only treat cancer but also provide beam services to the industry or academia. The simplified approach based on the Monte Carlo-based data sets (source terms and attenuation lengths) with the point-source line-of-sight approximation is friendly in the design stage of the proton therapy facilities because it is intuitive and easy to use. The purpose of this study is to expand the Monte Carlo-based data sets to allow the simplified approach to cover the application of proton beams more widely. Materials and Methods: In this work, the MCNP6 Monte Carlo code was used in three simulations to achieve the purpose, including the neutron yield calculation, Monte Carlo-based data sets generation, and dose assessment in simple cases to demonstrate the effectiveness of the generated data sets. Results and Discussion: The consistent comparison of the simplified approach and Monte Carlo simulation results show the effectiveness and advantage of applying the data set to a quick shielding design and conservative dose assessment for proton therapy facilities. Conclusion: This study has expanded the existing Monte Carlo-based data set to allow the simplified approach method to be used for dose assessment or shielding design for beam services in proton therapy facilities. It should be noted that the default model of the MCNP6 is no longer the Bertini model but the CEM (cascade-exciton model), therefore, the results of the simplified approach will be more conservative when it was used to do the double confirmation of the final shielding design.

• 제어로봇시스템학회:학술대회논문집
• /
• 1996.10a
• /
• pp.64-68
• /
• 1996

In this note, we propose the discrete model reduction method over disc-type analytic domains. We define Hankel singular value over the disc that is mapped by standard bilinear mapping. And GSPA(generalized singular perturbation approximation) and DT(direct truncation) are generalized to GSPA and DT over a disc. Furthermore we show that the reduced order model over a smaller domain has a smaller L$_{\infty}$ norm error bound..

• Journal of Institute of Control, Robotics and Systems
• /
• v.5 no.4
• /
• pp.481-488
• /
• 1999

A thermal model based on the chamber geometry of the industry-standard AST SHS200MA rapid thermal processing system has been developed for the study of thermal uniformity and process repeatability thermal model combines radiation energy transfer directly from the tungsten-halogen lamps and the steady-state thermal conducting equations. Because of the difficulties of solving partial differential equation, calculation of wafer temperature was performed by using finite-difference approximation. The proposed thermal model was verified via titanium silicidation experiments. As a result, we can conclude that the thermal model show good estimation of wafer surface temperature distribution.

• Journal of Korean Institute of Industrial Engineers
• /
• v.32 no.1
• /
• pp.9-17
• /
• 2006

This paper deals with a discrete simulation optimization method for designing a complex probabilistic discrete event simulation. The developed algorithm uses the configuration algorithm that can change decision variables and the stopping algorithm that can end simulation in order to satisfy the given objective value during single run. It tries to estimate an auto-regressive model for evaluating correctly the objective function obtained by a small amount of output data. We apply the proposed algorithm to M/M/s model, (s, S) inventory model, and known-function problem. The proposed algorithm can't always guarantee the optimal solution but the method gives an approximate feasible solution in a relatively short time period. We, therefore, show the proposed algorithm can be used as an initial feasible solution of existing optimization methods that need multiple simulation run to search an optimal solution.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.34-40
• /
• 1999

Drag reduction on vehicles are the main concern for the body shape designers in order to lower the fuel consumption rate and to aid the driving stability. The drag of bluff bodies like transportation vehicles is mostly pressure drag due to the flow separation, which can be minimized by controlling the location and size of the separation bubble. In the present study, the TURBO-3D code is incorporated with optimal algorithm based on analytical approximation method to obtain an optimal afterbody shape of the MIRA Model corresponding to the lowest drag coefficient. For this purpose three mutually independent afterbody angles are chosen as design variables, while the drag coefficient is chosen as an objective function. It is demonstrated in the present study that an optimal body shape having the lowest drag coefficient which is about 6% lower than that of the original shape has been successfully obtained within number of iterations of tile optimal design loop.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10b
• /
• pp.1858-1863
• /
• 1991

The theory of stock option pricing has, recently, attracted attention of many researchers interested not only in finance but also in statistics and control theory. In this field, the problem of estimating stock return volatility is, above all, of great importance in calculating actual stock option value. In this paper, we assume that the stock market is represented by the stochastic volatility model which is the same as that of Hull and White. Then, we propose an approximation function of option value. It is a type of Black-Sholes option formula in which the first and the second order moments of logarithmic stock value are modified in a special form from the original model. Finally, an algorithm of estimating the parameters of the stochastic volatility model is given, and parameters are estimated by using Nikkei 225 index option data.

• Publications of The Korean Astronomical Society
• /
• v.15 no.spc2
• /
• pp.65-71
• /
• 2000

We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

• Journal of the Korean Data and Information Science Society
• /
• v.27 no.6
• /
• pp.1653-1660
• /
• 2016

Smoothly clipped absolute deviation (SCAD) penalty is known to satisfy the desirable properties for penalty functions like as unbiasedness, sparsity and continuity. In this paper, we deal with the regression function estimation and variable selection based on SCAD penalized censored regression model. We use the local linear approximation and the iteratively reweighted least squares algorithm to solve SCAD penalized log likelihood function. The proposed method provides an efficient method for variable selection and regression function estimation. The generalized cross validation function is presented for the model selection. Applications of the proposed method are illustrated through the simulated and a real example.

• Journal of applied mathematics & informatics
• /
• v.31 no.5_6
• /
• pp.711-732
• /
• 2013

In this paper, we are interested in finding explicit numerical formulas to evaluate defaultable bonds prices of firms. For this purpose, we use a default intensity whose values depend on the credit rating of these firms. Each credit rating corresponds to a state of the default intensity. Then, this regime switches as soon as one of the credit rating of a firm also changes. Moreover, this regime switching default intensity model allows us to capture well some market features or economics behaviors. Thus, we obtain two explicit different formulas to evaluate the conditional Laplace transform of a regime switching Cox Ingersoll Ross model. One using the property of semi-affine of the model and the other one using analytic approximation. We conclude by giving some numerical illustrations of these formulas and real data estimation results.

• 한국정보통신설비학회:학술대회논문집
• /
• 2007.08a
• /
• pp.27-30
• /
• 2007

This paper deals with an equivalent circuit model for RF passive components. Rational functions are obtained from the frequency responses of EM simulation by using Foster canonical partial fraction expressions. The Vector Fitting(VF) and the Adaptive Frequency Sampling(AFS) scheme are also implemented to obtain the rational functions. A passivity enforcement algorithm is applied to ensure the stability of the equivalent circuit model. In order to verify the schemes, S parameters of the equivalent circuit model is compared to those of EM simulation in case of the microstrip line structure with 3 slots in ground.