• Communications for Statistical Applications and Methods
• /
• v.15 no.4
• /
• pp.563-585
• /
• 2008

This paper presents how to improve the efficiency and accuracy in the pricing and sensitivity evaluation for derivatives, since the need for the evaluation of complicated derivatives is increased. The Monte Carlo(MC) simulation using the quasi random number instead of pseudo random number can improve the elapsed time and accuracy for the valuation of European-type derivatives. However, the quasi MC simulation method has its limit for applying it in the multi-dimensional case such as American-type and path-dependent options due to the increased correlation between dimensions as the dimension of random numbers is increased. In order to complement this problem, we develop a modified method in which correlation values are controlled to be below a pre-specified value. Thus, this method is applicable for the pricing of either derivatives ill which underlying assets or risk factors are several or derivatives having path-dependent or early redemption property. Furthermore, we illustrate that it is important to take an appropriate grid interval for the use of finite difference method(FDM) by applying the FDM to one example of non-symmetrical butterfly spreads.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.36 no.4
• /
• pp.329-336
• /
• 2000

Welding is a important technological method in mechanical engineering. $CO_2$MAG(metal active gas) welding means that metal part in double capstan drum for the inshore and costal vessels are joined by melting(with or without a filler material) or that new material is added to a metal part by melting. The thermal stresses appear due to a non-uniform temperature field, inhomogeneous material properties, external restraint and volume changes during phase transformations. In this study analysis the elasto-plastic thermal stresses distribution of welding part in double capstan drum for the inshore and costal vessels using finite element method (FBM). Therefore it calculates the numerical value that can be applied to the optimum design of welding parts and the shapes. The significant results obtained in this study are summarized as fellows. At early stage of the cooling after welding process, the abrupt thermal stresses gradient has been shown in the vicinity of welding part. In the thermal stresses analysis due to temperature gradient and heat shocking maximum stress was occurred of welding part and stresses were distributed from 54MPa~48MPa.

• Journal of Sensor Science and Technology
• /
• v.8 no.2
• /
• pp.115-123
• /
• 1999

A micro-gas sensor with heater and sensing electrode on the same plane was fabricated on phosphosilicate glass(PSG, 800nm)/$Si_3N_4$ (150nm) dielectric membrane. PSG film was provided by atmospheric pressure chemical vapor deposition(APCVD), and $Si_3N_4$ film by low pressure chemical vapor deposition (LPCVD). Total area of the fabricated device was $3.78{\times}3.78mm^2$. The area of diaphragm was $1.5{\times}1.5mm^2$, and that of the sensing layer was $0.24{\times}0.24mm^2$. Finite-element simulation was employed to estimate temperature distribution for a square-shaped diaphragm. The power consumption of Pt heater was about 85mW at $350^{\circ}C$. Tin thin films were deposited on the silicon substrate by thermal evaporation at room temperature and $232^{\circ}C$, and tin oxide films($SnO_2$) were prepared by thermal oxidation of the metallic tin films at $650^{\circ}C$ for 3 hours in oxygen ambient. The film analyses were carried out by SEM and XRD techniques. Effects of humidity and ambient temperature on the resistance of the sensing layer were found to be negligible. The fabricated micro-gas sensor exhibited high sensitivity to butane gas.

• Korean Chemical Engineering Research
• /
• v.45 no.3
• /
• pp.242-248
• /
• 2007

For an optimal design of automotive electric system, it is important to have a reliable modeling tool to predict the charge-discharge behaviors of the automotive battery. In this work, a two-dimensional modeling was carried out to predict the charge-discharge behaviors of a 12-V automotive lead-acid battery. The model accounted for electrochemical kinetics and ionic mass transfer in a battery cell. In order to validate the modeling, modeling results were compared with the experimental data of the charge-discharge behaviors of a lead-acid battery. The discharge behaviors were measured with three different discharge rates of C/5, C/10, and C/20 at operating temperature of $25^{\circ}C$. The batteries were charged with constant current of 30A until the charging voltage reached to a predetermined value of 14.24 V and then the charging voltage was kept constant. The discharge and charge curves from the measurements and modeling were in good agreement. Based on the modeling, the distributions of the electrical potentials of the solid and solution phases, the porosity of the electrodes, and the current density within the electrodes as well as the acid concentration can be predicted as a function of charge and discharge time.

• Journal of the Korean Institute of Gas
• /
• v.2 no.2
• /
• pp.18-25
• /
• 1998

According to the requirements of ANSI B3l.8, the pipe thickness is determined with hoop stress resulted from internal pressure. And the other loads induced by soil, vehicle, thermal expansion, ground subsidence, etc shall be evaluated rationally. There are two ways of calculating stress of buried gas pipeline. The first is FEM. FEM can calculate the stress regardless of the complexity of pipeline shape and boundary conditions. But it needs high cost and long time. The second is the way to use equation. The reliable equations to calculate the stress of buried gas pipeline was developed and have been used in designing pipeline and evaluating pipeline safety, But these equation are very difficult to understand and use for non-specialist. For easy calculation of non-specialist, the new computer program to calculate stress of buried natural gas pipeline have been developed. The stress is calculated by the equations and extrapolation of the graph resulted from FEM. The full paper is consist of series I and II. In this paper, series I, the calculating equation of the program is explained in detail.

• Fire Science and Engineering
• /
• v.18 no.1
• /
• pp.24-30
• /
• 2004

In this work the dynamic heat transfer occurring in a cable penetration fire stop system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealants. Here was carried out an experiment to observe the heat transfer in the cable penetration fire stop system made of DOW CORNING products. The dynamic heat transfer occurring in the fire stop system is formulated in a parabolic partial differential equation subjected to a set of initial and boundary conditions. And it was modeled, simulated, and analyzed. The simulation results were illustrated in three-dimensional graphics and were compared with experimental data. Through the simulations, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable streams. It also was found that the dynamic heat transfer through the cable streams was one of the most dominant factors, and the feature of heat conduction could be understood as an unsteady-state process. It is certain that these numerical results are useful for making a performance-based design for the cable penetration fire stop system.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.1
• /
• pp.75-82
• /
• 2007

This paper addresses a method for shape optimization of a continuous elastic body considering stability, i.e., buckling behavior. The sensitivity formula for critical load is analytically derived and expressed in terms of shape variation, based on the continuum formulation of the stability problem. Unlike the conventional finite difference method (FDM), this method is efficient in that only a couple of analyses are required regardless of the number of design parameters. Commercial software such as ANSYS can be employed since the method requires only the result of the analysis in computation of the sensitivity. Though the buckling problem is more efficiently solved by structural elements such as a beam and shell, elastic solids have been chosen for the buckling analysis because solid elements can generally be used for any kind of structure whether it is thick or thin. Sensitivity is then computed by using the mathematical package MATLAB with the initial stress and buckling analysis of ANSYS. Several problems we chosen in order to illustrate the efficiency of the presented method. They are applied to the shape optimization problems to minimize weight under allowed critical loads and to maximize critical loads under same volume.

• Composites Research
• /
• v.13 no.5
• /
• pp.50-59
• /
• 2000

In this paper, two methods to suppress flutter of the composite panel are examined. First, in the active control method, a controller based on the linear optimal control theory is designed and control input voltage is applied on the actuators and a PZT is used as actuator. Second, a new technique, passive suppression scheme, is suggested for suppression of the nonlinear panel flutter. In the passive suppression scheme, a shunt circuit which consists of inductor-resistor is used to increase damping of the system and as a result the flutter can be attenuated. A passive damping technology, which is believed to be more robust suppression system in practical operation, requires very little or no electrical power and additional apparatuses such as sensor system and controller are not needed. To achieve the great actuating force/damping effect, the optimal shape and location of the actuators are determined by using genetic algorithms. The governing equations are derived by using extended Hamilton's principle. They are based on the nonlinear von Karman strain-displacement relationship for the panel structure and quasi-steady first-order piston theory for the supersonic airflow. The discretized finite element equations are obtained by using 4-node conforming plate element. A modal reduction is performed to the finite element equations in order to suppress the panel flutter effectively and nonlinear-coupled modal equations are obtained. Numerical suppression results, which are based on the reduced nonlinear modal equations, are presented in time domain by using Newmark nonlinear time integration method.

• Journal of Biomedical Engineering Research
• /
• v.23 no.6
• /
• pp.451-458
• /
• 2002

Osteoporosis is a systemic skeletal disease caused by low bone mass and the decrease of bone density in the microstructure of trabecular bone. Drug therapy(PTH Parathyroid hormone) may increase the trabecular thickness and thus bone strength. Vertebroplasty is a minimally invasive surgery foy the treatment of osteoporotic vertebral compression fracture. This Procedure includes Puncturing vertebrae and filling with Polymethylmethacrylate(PMMA). Although altering recommended monomer-to-Powder ratio affects material properties of bone cement, clinicians commonly alter the mixture ratio to decrease viscosity and increase the working time. The Purposes of this study were to analyze the effect of 4he monomer-to-powder ratio on the mechanical characteristics of trabecular. In this paper, the finite element model of human vertebral trabecualr bone was developed by modified Voronoi diagram, to analyze the relative effect of hormone therapy and vertebroplasty at the treatment of osteoporotic vertebrae. Trabeuclar bone models for vertebroplasty with varied monomer-to-Powder ratio(0.40∼1.07 ㎖/g) were analyzed. Effective modulus and strength of bone cement-treated models were approximately 60% of those of intact models and these are almost twice the values of hormone-treated models. The bone cement models with the ratio of 0.53㎖/g have the maximum modulus and strength. For the ratio of 1.07㎖/g, the modulus and strength were minimum(42% and 49% respectively) but these were greater than those for drug therapy. This study shows that bone cement treatment is more effective than drug therapy. It is found that in vertebroplasty, using a monomer-to-powder ratio different from that recommended by manufacturer nay significantly not only reduce the cement's material Properties but also deteriorate the mechanical characteristics of osteoporotic vertebrae.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.141-154
• /
• 1993

In this paper, a method to calculate the ultimate compressive strength of welded one-sided stiffened plates simply supported along all edges is proposed. At first initial imperfections such as distortions and residual stresses due to welding are predicted by using simplified methods. Then, the collapse modes of the stiffened plate are assumed and collapse loads for each mode are calculated. Among these loads, the lowest value is selected as the ultimate strength of the plate. Collapse modes are assumed as follows ; (1) Overall buckling of the stiffened plate$\rightarrow$Overall collapse due to stiffener bending (2) Local buckling of the plate part$\rightarrow$Local collapse of the plate part$\rightarrow$Overall collapse due to stiffener yielding (3) Local buckling of the plate part$\rightarrow$Overall collapse due to stiffener berthing (4) Local buckling of the plate part$\rightarrow$Local collapse of the plate part$\rightarrow$Overall collapse due to stiffener tripping. The elastic large deflection analysis based on the Rayleigh-Ritz method is carried out, and plastic analysis assuming hinge lines is also carried out. Collapse load is defined as the cross point of the two analysis curves. This method enables the utimate strength to be calculated with small computing time and a good accuracy. Using the present method, characteristics of the stiffener including torsional rigidity, bending and tripping can also be clarified.