• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.164-175
• /
• 2011

Using natural electromagnetic (EM) fields at low frequencies, magnetotelluric (MT) surveys can investigate conductivity structures of the deep subsurface and thus are used to explore geothermal energy resources and investigate proper sites for not only geological $CO_2$ sequestration but also enhanced geothermal system (EGS). Moreover, marine MT data can be used for better interpretation of marine controlled-source EM data. In the interpretation of MT data, MT modeling schemes are important. This study improves a three dimensional (3D) MT modeling algorithm which uses edge finite elements. The algorithm computes magnetic fields by solving an integral form of Faraday's law of induction based on a finite difference (FD) strategy. However, the FD strategy limits the algorithm in computing vertical magnetic fields for a topographic model. The improved algorithm solves the differential form of Faraday's law of induction by making derivatives of electric fields, which are represented as a sum of basis functions multiplied by corresponding weightings. In numerical tests, vertical magnetic fields for topographic models using the improved algorithm overcome the limitation of the old algorithm. This study recomputes induction vectors and tippers for a 3D hill and valley model which were used for computation of the responses using the old algorithm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.347-352
• /
• 2007

Eigenvalue reduction schemes approximate the lower eigenmodes that represent the global behavior of the structures. In the previous study, we proposed a two-level condensation scheme (TLCS) for the construction of a reduced system. In the first step, the selection of candidate elements by energy estimation, Rayleigh quotient, through Ritz vector calculation. In the second step, the primary degrees of freedom are selected by the sequential elimination method from the degrees of freedom connected to the candidate elements in the first step. In the present study, we propose TLCS combined with iterative improved reduced system (IIRS) to increase accuracy of the higher modes in the intermediate range. Also, it is possible to control the accuracy of the eigenvalues and eigenmodes of the reduced system. Finally, numerical examples demonstrate the performance of the proposed method.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.3 no.1
• /
• pp.33-44
• /
• 1995

As the occupant safety receives more attention from automobile industries. protection systems have been developed quite well. Developed protection systems must be evaluated through real tests in crash environment Since the real tests are extremely expensive. computer simulations are replaced for some prediction of the real test In the computer simulation. it is very crucial to express the real environment precisely in the modeling precess. The energy absorbing(EA) steering system has a very important rote in vehicle crashes because the occupant can hit the system directly. In this study. the EA steering system is modeled precisely. analyzed for the safely and designed by an optimization technology. First. the EA steering system is disassembled by parts and modeled by segments and joints. The segments are modeled by rigid bodies in motion and they have resistances in contact. Spring-damper elements and force-deflection curves are utilized to represent the joints. The body block test is cal lied out to validate. the modeling. When the test results are not enough for the detailed modeling. the differences between tests and simulations are minimized to calculate unknown parameters using optimization. The established model is applied to a crash simulation of a full-car model and tuned again. After the modeling is finished. components of the steering system are designed by an optimization algorithm. In the optimization process. the compound injury of a driver is defined and minimized to determine the chracteristics of the components. The second. order approximation algorithm has been adopted for the optimization.

• Journal of Korean Society of Steel Construction
• /
• v.14 no.6
• /
• pp.783-791
• /
• 2002

An improved formulation for multi-objective optimization was proposed. This formulation was applied to steel seismic loads. The multi-objective optimization problem was formulated with minimum structural weight, maximum strstability. The global criterion method was employed to find a rational solution closest to the ideal solution for the optimization problem using standard steel profile, To efficiently solve the optimization problem, the decomposition meth both system-level and element-level was used. In addition, various techniques including efficient reanalysis technique intermediate variables and sensitivity analysis using an automatic differentiation(AD) were incorporated. Moreover the reamong section properties fitted to the section profile used in order to link the system level and the element level. From numerical investigation, it could be stated that the proposed method will lead to the more rational design compared with one.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.1
• /
• pp.85-92
• /
• 2013

This paper deals with the resonant type of a WEC (wave energy converter) and the determination method of its geometric parameters which were obtained to construct the robust and optimal structure, respectively. In detail, the optimization problem is formulated with the constraints composed of the response surfaces which stand for the resonance period(heave, pitch) and the meta center height of the buoy. Use of a signal-to-noise ratio calculated from normalized multi-objective results with the weight factor can help to select the robust design level. In order to get the sample data set, the motion responses of the power buoy were analyzed using the BEM (boundary element method)-based commercial code. Also, the optimization result is compared with a robust design for a feasibility study. Finally, the power efficiency of the WEC with the optimum design variables is estimated as the captured wave ratio resulting from absorbed power which mainly related to PTO (power take off) damping. It could be said that the resultant of the WEC design is the economical optimal design which satisfy the given constraints.

• Fire Science and Engineering
• /
• v.15 no.1
• /
• pp.7-15
• /
• 2001

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP (Smoke Movement Estimating Program) codo to the simulation of fire induced flows in the atrium space (SIVANS atrium at Japan) containing smoke radiation effect. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown a better prediction than the result calculated by only convection effect in comparison with the experimental data. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire should be necessary in order to get more realistic result. Also the numerical results indicated that the smoke layer is developing at a rate of about 0.1 m/s. It would take about 450 seconds after starting the ultra fast fire of 560 kW that the smoke layer move down to 1.5m above the escape level.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2011.08a
• /
• pp.273-273
• /
• 2011

현재 사용되고 있는 플로팅 게이트를 이용한 플래시 메모리 소자는 비례축소에 의해 발생하는 단 채널 효과, 펀치스루 효과 및 소자간 커플링 현상과 같은 문제로 소자의 크기를 줄이는데 한계가 있다. 이러한 문제를 해결하기 위하여 silicon nitride와 같은 절연체를 전자의 트랩층으로 사용하는 charge trap flash (CTF) 메모리 소자에 대한 연구가 활발히 진행되고 있다. CTF 메모리 소자의 전기적 특성에 대한 연구는 활발히 진행 되었지만, 수치 해석 모델을 사용하여 메모리 소자의 전하수송 메커니즘을 분석한 연구는 매우 적다. 본 연구에서는 수치 해석 모델을 적용하여 개발한 시뮬레이터를 사용하여 CTF 메모리 소자의 프로그램 동작 시 전하수송 메커니즘에 대한 연구를 하였다. 시뮬레이터에 사용된 모델은 연속방정식, 포아송 방정식과 Shockley-Read-Hall 재결합 모델을 수치해석적 방법으로 계산하였다. 또한 CTF 소자 프로그램 동작 시 트랩 층으로 주입되는 전자의 양은 Wentzel-Kramers-Brillouin 근사 법을 이용하여 계산하였다. 트랩 층에 트랩 되었던 전자의 방출 모델은 이온화 과정을 사용하였다. 게이트와 트랩 층 사이의 터널링은 Fowler-Nordheim (FN) tunneling 모델, Direct tunneling 모델, Modified FN tunneling 모델을 적용하였다. FN tunneling 만을 적용했을때 보다 세가지 모델을 적용했을 때가 더 실험치와의 오차가 적었다. 그 이유는 시뮬레이션 결과를 통해 인가된 전계에 의해 Bottom Oxide 층의 에너지 밴드 구조가 변화하여 세가지 tunneling 모델의 구역이 발생하는 것을 확인 할 수 있었다. 계산된 결과의 전류-전압 곡선을 통해 CTF 메모리 소자의 프로그램 동작 특성을 관찰하였다. 트랩 층의 전도대역과 트랩 층 내부에 분포하는 전자의 양을 시간에 따라 계산하여 트랩 밀도가 시간이 지남에 따라 일정 값에 수렴하고 많은 전하가 트랩 될 수록 전하 주입이 줄어듬을 관찰 하였다. 이와 같은 시뮬레이션 결과를 통해 CTF 메모리의 트랩층에서 전하의 이동에 대해 더 많이 이해하여 CTF 소자가 가진 문제점 해결에 도움을 줄 것이다.

• Journal of the Korean Chemical Society
• /
• v.40 no.3
• /
• pp.187-195
• /
• 1996

Ab initio SCF calculations have been carried out for the five conceivable trimers formed between one hydrogen cyanide and two hydrogen fluorides using a basis set of TZ＋P quality. Several ground state properties of these trimeric complexes have been evaluated, and compared with those of isolated monomers and appropriate dimers. Computed equilibrium geometries, stabilization energies, and dipole moments are given in order to suppliment the available experimental data. At this level of approximation, intramolecular bond distances are consistently shorter than experimental ones. However, intermolecular distances upon complex formation, and dipole moments are overestimated compared with experimental ones. HCN$(HF)_2$ trimer appears to be the most favourable among the five kinds of trimer complex, and also more stable than $(HCN)_2HF$. The typical features of the non-additivity of intermolecular interaction are relatively strong in the HCN$(HF)_2$ trimer.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.4 s.107
• /
• pp.351-357
• /
• 2006

In this paper, a direct $TM_{01}-TE_{11}$ mode converter of circular waveguide which is used the approximated design method of mode conversion between $TM_{01}$ and $TE_{11}$ modes is designed for radiation with the maximum directivity of high power energy into air (rom high power microwave antenna system. The proposed direct $TM_{01}-TE_{11}$ mode converter is calculated and designed with an approximated method for non-linear function about an extended non-constant serpentine type of circular waveguide. Then, the designed mode converter is optimized by length of 200 mm and efficiency above 95% by FIM simulation and achieved short length and high efficiency of the antenna by results of fabricated and measured characteristics.