• The Journal of Korean Society for Radiation Therapy
• /
• v.30 no.1_2
• /
• pp.97-105
• /
• 2018

Purpose : Proton Therapy using Bragg-peak, because it has distinct characteristics in providing maximum dosage for tumor and minimal dosage for normal tissue, a medical imaging system that can quantify changes in patient position or treatment area is of paramount importance to the treatment of protons. The purpose of this research is to evaluate the usefulness of the algorithm by comparing the image matching through the set-up and in-house code through the existing dips program by producing a Matlab-based in-house registration code to determine the error value between dips and DRR to evaluate the accuracy of the existing treatment. Materials and Methods : Thirteen patients with brain tumors and head and neck cancer who received proton therapy were included in this study and used the DIPS Program System (Version 2.4.3, IBA, Belgium) for image comparison and the Eclipse Proton Planning System (Version 13.7, Varian, USA) for patient treatment planning. For Validation of the Registration method, a test image was artificially rotated and moved to match the existing image, and the initial set up image of DIPS program of existing set up process was image-matched with plan DRR, and the error value was obtained, and the usefulness of the algorithm was evaluated. Results : When the test image was moved 0.5, 1, and 10 cm in the left and right directions, the average error was 0.018 cm. When the test image was rotated counterclockwise by 1 and $10^{\circ}$, the error was $0.0011^{\circ}$. When the initial images of four patients were imaged, the mean error was 0.056, 0.044, and 0.053 cm in the order of x, y, and z, and 0.190 and $0.206^{\circ}$ in the order of rotation and pitch. When the final images of 13 patients were imaged, the mean differences were 0.062, 0.085, and 0.074 cm in the order of x, y, and z, and 0.120 cm as the vector value. Rotation and pitch were 0.171 and $0.174^{\circ}$, respectively. Conclusion : The Matlab-based In-house Registration code produced through this study showed accurate Image matching based on Intensity as well as the simple image as well as anatomical structure. Also, the Set-up error through the DIPS program of the existing treatment method showed a very slight difference, confirming the accuracy of the proton therapy. Future development of additional programs and future Intensity-based Matlab In-house code research will be necessary for future clinical applications.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.499-504
• /
• 2001

A validation of a 3D CFD model for predicting local subcooling of moderator in the vicinity of calandria tubes in a CANDU reactor is performed. The small scale moderator experiments performed at Sheridan Park Experimental Laboratory(SPEL) in Ontario, Canada[1] is used for the validation. Also a comparison is made between previous CFD analyses based on 2DMOTH and PHOENICS, and the current model analysis for the same SPEL experiment. For the current model, a set of grid structures for the same geometry as the experimental test section is generated and the momentum, heat and continuity equations are solved by CFX-4.3, a CFD code developed by AEA technology. The matrix of calandria tubes is simplified by the porous media approach. The standard $k-\varepsilon$ turbulence model associated with logarithmic wall treatment and SIMPLEC algorithm on the body fitted grid are used and buoyancy effects are accounted for by the Boussinesq approximation. For the test conditions simulated in this study, the flow pattern identified is a buoyancy-dominated flow, which is generated by the interaction between the dominant buoyancy force by heating and inertial momentum forces by the inlet jets. As a result, the current CFD moderator analysis model predicts the moderator temperature reasonably, and the maximum error against the experimental data is kept at less than $2.0^{\circ}C$ over the whole domain. The simulated velocity field matches with the visualization of SPEL experiments quite well.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.4
• /
• pp.359-366
• /
• 2009

MATLAB Autocode generation is a feature that converts a block diagram model in Simulink to a c program. Utilizing this function makes MATLAB/Simulink an integrated developing environment, from controller design to implementation. It can reduce development cost and time significantly. However, this automated process requires high reliability on the software, especially the original Simulink block diagram model. And thus, the verification of the codes becomes important. In this study, a UAV flight program which is generated with Simulink is validated and modified according to DO-178B. As a result of applying the procedures, the final program not only satisfied the functional requirement but is also verified with structural point of view with Decision Coverage 93%, Condition Coverage 95% and MC/DC 90%.

• Computers and Concrete
• /
• v.19 no.1
• /
• pp.99-110
• /
• 2017

In this paper, the effect of normal load on the failure mechanism of echelon joint has been studied using PFC2D. In the first step, calibration of PFC was undertaken with respect to the data obtained from experimental laboratory tests. Then, six different models consisting various echelon joint were prepared and tested under two low and high normal loads. Furthermore, validation of the simulated models were cross checked with the results of direct shear tests performed on non-persistent jointed physical models. The simulations demonstrated that failure patterns were mostly influenced by normal loading, while the shear strength was linked to failure mechanism. When ligament angle is less than $90^{\circ}$, the stable crack growth length is increased by increasing the normal loading. In this condition, fish eyes failure pattern occur in rock bridge. With higher ligament angles, the rock bridge was broken under high normal loading. Applying higher normal loading increases the number of fracture sets while dilation angle and mean orientations of fracture sets with respect to ligament direction will be decreased.

• International Journal of Aeronautical and Space Sciences
• /
• v.11 no.2
• /
• pp.145-153
• /
• 2010

In this work, the comprehensive structural dynamics codes including DYMORE and CAMRAD II are used to validate the higher harmonic control aeroacoustic rotor test (HART) II data in descending flight condition. A total of 16 finite elements along with 17 aerodynamic panels are used for the CAMRAD II analysis; whereas, in the DYMORE analysis, 10 finite elements with 31 equally-spaced aerodynamic panels are utilized. To improve the prediction capability of the DYMORE analysis, the finite state dynamic inflow model is upgraded with a free vortex wake model comprised of near shed wake and trailed tip vortices. The predicted results on aerodynamic loads and blade motions are correlated with the HART II measurement data for the baseline, minimum noise and minimum vibration cases. It is found that an improvement of solution, especially for blade vortex interaction airloads, is achieved with the free wake method employed in the DYMORE analysis. Overall, fair to good correlation is achieved for the test cases considered in this study.

• Computers and Concrete
• /
• v.18 no.1
• /
• pp.39-51
• /
• 2016

This paper considers the tensile strength of concrete samples in direct, CTT, modified tension, splitting and ring tests using both of the experimental tests and numerical simulation (particle flow code 2D). It determined that which one of indirect tensile strength is close to direct tensile strength. Initially calibration of PFC was undertaken with respect to the data obtained from Brazilian laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, validation of the simulated models in four introduced tests was also cross checked with the results from experimental tests. By using numerical testing, the failure process was visually observed and failure patterns were watched to be reasonable in accordance with experimental results. Discrete element simulations demonstrated that the macro fractures in models are caused by microscopic tensile breakages on large numbers of bonded discs. Tensile strength of concrete in direct test was less than other tests results. Tensile strength resulted from modified tension test was close to direct test results. So modified tension test can be a proper test for determination of tensile strength of concrete in absence of direct test. Other advantages shown by modified tension tests are: (1) sample preparation is easy and (2) the use of a simple conventional compression press controlled by displacement compared with complicate device in other tests.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.4 s.148
• /
• pp.460-466
• /
• 2006

The three-dimensional RANS equations were applied to analyze the flow field of a submarine. To validate the code, the DARPA SUBOFF bare hull and an eliipsoid at angles of attack of $10^{\circ}\;and\;30^{\circ}$ were simulated and good agreement with experiments was obtained. After the code validation, the flows over the full configuration of DARPA SUBOFF model having a fairwater and four stern appendages were simulated at four angles of attack $(0^{\circ},\;10^{\circ},\;20^{\circ},\;30^{\circ})$ and three yaw angles $(10^{\circ},\;20^{\circ},\;30^{\circ})$ Specifically, the pressure contours and streamlines of fairwater and stern appendage were compared as the angle of attack and yaw angle changed. The variations of hydrodynamic forces were also calculated.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.2
• /
• pp.117-126
• /
• 2016

In this paper, a numerical procedure for the natural vibration analysis of plates with openings and carlings based on the assumed mode method is extended to assess their forced response. Firstly, natural response of plates with openings and carlings is calculated from the eigenvalue equation derived by using Lagrange's equation of motion. Secondly, the mode superposition method is applied to determine frequency response. Mindlin theory is adopted for plate modelling and the effect of openings is taken into account by subtracting their potential and kinetic energies from the corresponding plate energies. Natural and frequency response of plates with openings and carlings subjected to point excitation force and enforced acceleration at boundaries, respectively, is analysed by using developed in-house code. For the validation of the developed method and the code, extensive numerical results, related to plates with different opening shape, carlings and boundary conditions, are compared with numerical data from the relevant literature and with finite element solutions obtained by general finite element tool.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.4
• /
• pp.309-314
• /
• 2010

In this work, the helicopter price estimation formula from Harris & Scully is revised and updated for more accurate prediction of the current prices of helicopters which can be applied in the preliminary design stage. As many as 74 helicopter price and sizing data after the year 1997 are newly included to the out-dated coefficients in the cost estimation formula. In addition, the updated price estimation formula is added to Tishchenko's preliminary sizing method for price consideration in the preliminary design stage. Sikorsky's S-76C+ helicopter is used along with other various helicopters for comparison and validation of the sizing code as well as the new updated price estimation formula.

• The KSFM Journal of Fluid Machinery
• /
• v.17 no.6
• /
• pp.89-94
• /
• 2014

The 4th generation SFR is being designed with a milestone of construction by 2028. It is important to understand the subchannel flow characteristics in fuel assembly through the experimental investigations and to estimate the calculation uncertainties for insuring the confidence of the design code calculation results. The friction coefficient and the mixing coefficient are selected as primary parameters. The two parameters are related to the flow distribution and diffusion. To identify the flow distribution, an iso-kinetic method was developed based on the previous study. For the mixing parameters, a wire mesh system and a laser induced fluorescence methods were developed in parallel. The measuring systems were adopted on 37 rod bundle test geometry, which was developed based on the Euler number scaling. A scaling method for a design of experimental facility and the experimental identification techniques for the flow distribution and mixing parameters were developed based on the measurement requirement.