• 한국정보통신학회논문지
• /
• 제12권12호
• /
• pp.2253-2258
• /
• 2008

산화물 반도체 감지막이 동작온도에 따라 감응특성을 가지는 마이크로 흐름센서를 설계하기 위해서 통계적 수법에 기초한 Markov 체인 MCM을 이용하여 기초방정식을 정식화하고 마이크로 소자의 열 전달특성을 해석하였다. 계산 결과를 통하여 기존 유한차분법이 가지는 계산 정밀도와 차이가 없음을 확인하였다. 본 논문에서 제안한 Markov 체인 MCM을 활용하면 다양한 마이크로 소자의 열전달 특성과 같은 물리적 특성을 해석하고 설계하는데 유용할 것으로 판단된다.

• 대한전자공학회논문지TC
• /
• 제39권4호
• /
• pp.201-208
• /
• 2002

본 논문에서는 3 차원 임의 형태 도체의 지연 산란 응답을 얻기 위한, 새로운 시간영역 전장 적분방정식(Time-Domain Electric Field Integral Equation: TD-EFIE)을 제안한다. 자기 벡터 전위의 시간 미분항은 중앙 차분으로, 전기 스칼라 전위는 시간에 대한 평균을 취한 두 개의 항으로 근사하였다. 이로부터 도체에 의한 산란 지연 응답 해의 산출시, 기존의 방법보다 정확하고 더욱 안정된 해를 얻을 수 있었다. 제안된 방법의 자세한 정식화 과정을 보였으며, 주파수 영역에서의 이산 푸리에 역변환 (Inverse Discrete Fourier Transform: IDFT) 결과치와 제안된 방법에 의한 수치해를 각각 비교하였다.

• 대한전자공학회논문지
• /
• 제26권1호
• /
• pp.48-61
• /
• 1989

본 논문에서는 게이트의 길이가 0.7${\mu}m$인 n형 GaAs MESFET를 2차원적으로 수치 해석하였으며, 이동도를 국부 전계의 함수로 취하는 드리프트 -확산 모델을 사용하였다. 이산화 방법으로는 종래에 사용되던 FDM(finite difference method), FEM(finite element method)을 사용치 아낳고 Control-Volume Formulation을 사용하였으며, numerical scheme으로는 기존의 hybrid scheme이나 upwind scheme 대신에 exponential scheme과 거의 근사한 power-law scheme을 사용하였다. 이때 드리프트 속도와 확산 속도의 비율을 나타내는 Peclet number의 개념을 사용하였으며, 이 개념을 사용하여 control volume의 경계에서 numerical scheme을 고려한 전류식을 제안하였다. 앞에서 고려한 모델들과 수치해석 방법을 사용하여 시뮬레이션한 I-V 특성은 기존 노문의 결과와 일치하였다. 따라서 본 논문의 결과가 GaAs MESFET를 위한 유용한 2차원 시뮬레이터가 될 수 있음을 확인하였다. 또한 I-V 특성외에 채널 밑바닥에서이 속도 및 전계 분포를 통해 드리프트-확산 모델을 고려한 경우에 발생하는 속도 포화의 메카니즘을 제시했고, Dipole의 발생위치 및 발생 원인과 드레인 전류와의 관계 등에 대해서도 제시했다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• 한국정밀공학회지
• /
• 제20권8호
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제51권11호
• /
• pp.551-558
• /
• 2002

A time-domain combined field integral equation (CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time-domain electric field integral equation (EFIE) with the magnetic field integral equation (MFIE). The time derivative of the magnetic vector potential in EFIE is approximated using a central finite difference approximation and the scalar potential is averaged over time. The time-domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. The incident spectrum of the field may contain frequency components, which correspond to the internal resonance of the structure. For the numerical solution, we consider both the explicit and implicit scheme and use two different kinds of Gaussian pulses, which may contain frequencies corresponding to the internal resonance. Numerical results for the EFIE, MFIE, and CFIE are presented and compared with those obtained from the inverse discrete Fourier transform (IDFT) of the frequency-domain CFIE solution.

• 한국주조공학회지
• /
• 제13권4호
• /
• pp.323-332
• /
• 1993

An implicit finite difference formulation with three methods of latent heat treatment, such as equivalent specific heat method, temperature recovery method and enthalpy method, was applied to solidification analysis. The Neumann problem was solved to compare the numerical results with the exact solution. The implicit solutions with the equivalent specific heat method and the temperature recovery method were comparatively consistent with the Neumann exact solution for smaller time steps, but its error increased with increasing time step, especially in predicting the solidification beginning time. Although the computing time to solve energy equation using temperature recovery method was shorter than using enthalpy method, the method of releasing latent heat is not realistic and causes error. The implicit formulation of phase change problem requires enthalpy method to treat the release of latent heat reasonably. We have modified the enthalpy formulation in such a way that the enthalpy gradient term is not needed, and as a result of this modification, the computation stability and the computing time were improved.

• Geomechanics and Engineering
• /
• 제23권4호
• /
• pp.365-379
• /
• 2020

One of the leading causes for the damage of various underground structures during an earthquake is soil liquefaction, and among this liquefaction-induced uplift of these structures is a major concern. In this study, finite-difference modelling is carried out to study the liquefaction-induced uplift of an underground structure of 5 m diameter (D) with and without the replacement of the in-situ fine sand around the structure with the coarse sand. Soil replacements are carried out by three methods: replacement of soil above the structure, around the structure, and below the structure. The soil behaviour is represented using the elastic-perfectly plastic Mohr-Coulomb model, where the pore pressures were computed using Finn-Byrne formulation. The predicted pore pressure and uplift of the structure due to sinusoidal input motion were validated with the centrifuge tests reported in the literature. Based on numerical studies, an empirical equation is developed for the determination of liquefaction-induced maximum uplift of the underground structure without replacement of the in-situ sand. It is found that the replacement of soil around the structure with 2D width and spacing of D can reduce the maximum uplift by 50%.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• 한국정밀공학회지
• /
• 제22권1호
• /
• pp.89-96
• /
• 2005

This paper is concerned with a balancing motion formulation and control of the ZMP (Zero Moment Point) for a biped-walking robot that has a prismatic balancing weight or a revolute balancing weight. The dynamic stability equation of a walking robot which have a prismatic balancing weight is conditionally linear but a walking robot's stability equation with a revolute balancing weight is nonlinear. For a stable gait, stabilization equations of a biped-walking robot are modeled as non-homogeneous second order differential equations for each balancing weight type, and a trajectory of balancing weight can be directly calculated with the FDM (Finite Difference Method) solution of the linearized differential equation. In this paper, the 3dimensional graphic simulator is developed to get and calculate the desired ZMP and the actual ZMP. The operating program is developed for a real biped-walking robot IWRⅢ. Walking of 4 steps will be simulated and experimented with a real biped-walking robot. This balancing system will be applied to a biped humanoid robot, which consist legs and upper body, as a future work.