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

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2006년도 제32회 KOSCO SYMPOSIUM 논문집
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제3권1호
• /
• pp.20-26
• /
• 2011

A weakly nonlinear seakeeping methodology for predicting motions and loads is presented in this paper. This methodology assumes linear radiation and diffraction forces, calculated in the frequency domain, and fully nonlinear Froude-Krylov and hydrostatic forces, evaluated in the time domain. The particular approach employed here allows to overcome numerical problems connected to the determination of the impulse response functions. The procedure is divided into three consecutive steps: evaluation of dynamic sinkage and trim in calm water that can significantly influence the final results, a linear seakeeping analysis in the frequency domain and a weakly nonlinear simulation. The first two steps are performed employing a three-dimensional Rankine panel method. Nonlinear Froude-Krylov and hydrostatic forces are computed in the time domain by pressure integration on the actual wetted surface at each time step. Although nonlinear forces are evaluated into the time domain, the equations of motion are solved in the frequency domain iteratively passing from the frequency to the time domain until convergence. The containership S175 is employed as a test case for evaluating the capability of this methodology to correctly predict the nonlinear behavior related to wave induced motions and loads in head seas; numerical results are compared with experimental data provided in literature.

• Coupled systems mechanics
• /
• 제6권4호
• /
• pp.501-522
• /
• 2017

This paper introduces a numerical procedure to incorporate elasto-plastic local deformation effects in the dynamic analysis of beams. The appealing feature is that simple beam type finite elements can be used for the global model which needs not to be altered by the localized elasto-plastic deformations. An overlapping local sophisticated 2D membrane model replaces the internal forces of the beam elements in the predefined region where the localized deformations take place. An iterative coupling technique is used to perform this replacement. Comparisons with full membrane analysis are provided in order to illustrate the accuracy and efficiency of the method developed herein. In this study, the membrane formulation is able to capture the elasto-plastic material behaviour based on the von Misses yield criterion and the associated flow rule for plane stress. The Newmark time integration method is adopted for the step-by-step dynamic analysis.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1997년도 추계 학술대회논문집
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• 대한기계학회논문집
• /
• 제14권6호
• /
• pp.1408-1416
• /
• 1990

본 연구에서는 과도 연성 열탄성문제의 해를 구할때 사용되는 직접시간적분방 법과 Laplace 변환방법은 상호 장 단점을 가지고 있다. 각 방법들의 장단점은 서로 배타적이므로 서로의 장점을 살리는 수치방법이 필요하다. 그런데, 대부분의 과도 열탄성문제는 급격한 온도변화로 인한 물체의 변형에 관심이 있기때문에 이 형태의 문 제를 효율적으로 다루는 데 주안점을 두고 본 연구를 수행하였다. 유도된 유한요소 방정식은 결국 열탄성 지배 방정식 중 열전달방정식인 에너지보존식은 Gurtin의 범함 수로부터 유도된 원래의 형태를 사용하나 수치적 안정성(numerical stability)을 보장 하기 위하여 운동방정식은 시간에 대한 2차미분 형태로 수정하였다. 에너지보존식은 시간에 대한 합성적분(convolution)형태로 표현되므로 온도의 시간미분항이 소거되므 로 경계에서의 급격한 온도변화로 인한 수치 해석적 문제점은 간단히 해결된다. 그 러므로, 제안된 수치해법은 직접시간적분방법의 일종이나 결과식인 유한요소방정식은 기존의 문헌들과 상당한 차이가 있다. 과도 연성 열탄성해석을 위한 새로운 근사수 치해법의 장점을 이론적으로 설명하기보다 수치계산면에서의 안전성, 정확성 및 효율 성이 있음을 증명하기 위하여 이미 발표된 문헌들에서 다룬 예제를 선정하여 해석결과 를 비교하였다.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2005년도 추계학술대회 논문집
• /
• pp.191-196
• /
• 2005

In this study, to describe vehicle-track-bridge dynamic interaction phenomena with 1/4 vehicle model, nonlinear Hertzian contact spring and nonlinear contact damper are introduced. In this approach external loads acting on 1/4 vehicle model are self weight of vehicle and geometry information of running surface. The constraint equation on contact surface is implemented by Penalty method. Also, to improve the numerical stability and to maintain accuracy of solution, the artificial damper and the reaction from constraint violation are introduced. A nonlinear time integration method, in this study, Newmark method is adopted for both equations of vehicles and structure. And to reduce the error caused by inadequate time step size, adaptive time-stepping technique is partially introduced. As the nonlinear Hertzian contact spring has no resistance to tensile force, the bouncing phenomena of wheelset can be described. Thus, it is expected that more versatile dynamic interaction phenomena can be described by this approach and it can be applied to various railway dynamic problems.

• 대한기계학회논문집
• /
• 제11권2호
• /
• pp.322-330
• /
• 1987

본 논문에서는 전자의 해석중에서 유한요소와 경계요소의 결합으로 생기는 문 제점 그리고 증가된 자요도로 인한 경제성의 손실 등을 고려하여 2차의 경계요소만을 사용한 경계요소법을 시도하였으며 다른 형태의 점탄성기본해 및 점탄성 경계적분식을 유도 사용하였다. 그리고 시도된 방법의 타당성을 보이기 위하여 사용된 예제들(평 면변형문제)의 이론해를 상응원리를 사용하여 유도하였으며 이를 수치결과와 비교하여 근사해의 정확함을 알 수가 있었다.

• 대한토목학회논문집
• /
• 제10권4호
• /
• pp.141-150
• /
• 1990

하구나 연안에서 해수의 순환형상을 모사(模寫)하게 위해 천수방정식(淺水方程式)을 여러 가지 경계조건 하에서 수치해석하였다. 공간영역은 Galerkin방법으로 이산화(離散化)하였으며 시간영역에 대해서는 유한차분법(Crank-Nicolson방법)을 사용하였다. 네 가지 검정실험이 해석적인 해가 있는 일차원 수로에서 행하여졌으며, 해석해를 구할 수 없는 이차원 모형에도 적용되었다. 해석해가 있는 경우 수치모사 결과가 이와 잘 일치하였으며, 이차원 모형에서의 결과도 매우 합당함을 알 수 있었다. 또 일차원 문제에서 4점 bilinear요소와 삼각형 요소를 사용한 결과를 각각 비교하였으며 시간적분도 2단계 Lax-Wendroff방법을 사용하여 결과를 비교하였다. 음해법을 사용할 경우 비교적 정확한 결과를 얻을 수 있으나 요소의 갯수가 많아지면 구성되는 대수방정식(代數方程式)이 커지기 때문에 각 시간마다의 계산량이 엄청나게 늘어나게 되며 양해법을 사용할 때는 원하는 만큼의 정확한 결과를 얻기 위하여 시간간격이나 공간격자 간격을 선정하는데 각별히 유의하여야 할 것이다.