• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.43-43
• /
• 2015

천수방정식을 사용하는 초기 수치모형은 프로드수($F_4$)가 변화하는 흐름 즉, 상류방향과 하류방향으로 전파하는 홍수파를 동시에 해석하기 위해 중앙 차분기법이 필요한 상류(sub-critical flow)와 흐름방향에 따른 상류이송(upwinding)기법이 필요한 사류(super-critical flow)가 나타나는 흐름해석에서 어려움이 있었다. 하지만, 근사 Riemann 해법의 등장으로 흐름방향에 관계없이 특성선을 따라 정확한 상향가중기법의 적용이 가능하게 되어, 천수방정식을 지배방정식으로 하는 수치모형이 더욱 실용적으로 적용될 수 있도록 하였다. 따라서, 현재 근사 Riemann 해법은 Godunov 형 유한체적 기법, 불연속 Galerkin 혹은 Petrov-Galerkin 유한요소기법 그리고 Boussinesq 기법에도 적용되고 있으며, 특히 Godunov 형 유한체적기법과 결합한 근사 Riemann 해법은 댐 붕괴, 하천 범람 그리고 도시 및 해안지역 침수에 이르기까지 여러 가지 문제에 폭넓게 적용되고 있다. 지금까지 홍수 모델링에 적용된 Godunov형 유한체적모형은 정형 사각격자나 비정형 삼각격자 중에서 한가지의 격자 종류만을 적용한 연구가 주로 수행되었으며, 유한요소모형과 같이 이 두 가지 격자를 동시에 적용한 연구는 거의 이루어지지 않고 있다. 일반적으로, 삼각격자는 사각격자와 는 달리 연구유역의 경계나 지형이 복잡한 경우에도 큰 노력없이 격자의 생성이 가능하나, 격자와 노드의 수가 사각격자보다 많아 계산시간이 많이 소요되는 단점이 있다. 반면, 사각격자는 하천과 같이 선형으로 변하는 지형에 대해서는 표현하기가 용이하며 계산시간의 효율성도 뛰어나다. 본 연구에서는 하천, 도시 그리고 해안지역에서의 효율적이고 정확한 홍수 모델링을 위해 삼각 및 사각격자 그리고 이 두 격자를 동시에 고려한 하이브리드 격자의 적용이 가능한 Godunov형 2차원 유한체적 모형을 개발하였다. 그리고 개발모형을 정확해가 있는 댐 붕괴 문제, 실측치가 존재하는 실험하도 및 실제하도에 삼각, 사각 그리고 혼합격자를 생성하여 모의를 수행하고, 각 적용 격자에 따른 정확성과 효율성 및 장점과 단점을 연구하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.38-52
• /
• 1991

In this paper, a semi-Lagrangian method is used to solve the nonlinear hydrodynamics of a three-dimensional body beneath the free surface in the time domain. The boundary value problem is solved by using the boundary integral method. The geometries of the body and the free surface are represented by the curved panels. The surfaces are discretized into the small surface elements using a bi-cubic B-spline algorithm. The boundary values of $\phi$ and $\frac{\partial{\phi}}{\partial{n}}$ are assumed to be bilinear on the subdivided surface. The singular part proportional to $\frac{1}{R}$ are subtracted off and are integrated analytically in the calculation of the induced potential by singularities. The far field flow away from the body is represented by a dipole at the origin of the coordinate system. The Runge-Kutta 4-th order algorithm is employed in the time stepping scheme. The three-dimensional form of the integral equation and the boundary conditions for the time derivative of the potential Is derived. By using these formulas, the free surface shape and the equations of motion are calculated simultaneously. The free surface shape and fille forces acting on a body oscillating sinusoidally with large amplitude are calculated and compared with published results. Nonlinear effects on a body near the free surface are investigated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1999.04a
• /
• pp.22-22
• /
• 1999

교란에 대한 고체 추진제의 연소율의 반등에 대한 이해는 고성능 추진제를 설계하는데 매우 중요한 요소이다. 그 동안의 연구는 고체 추진제의 표면에서 발생하는 교란이 매우 작은 크기로 발생한다는 선형적인 가정을 사용하여 이론적인 응답 함수를 구하였다. 특히 실험실에서 행해지는 교란에 대한 추진제의 응답 함수를 구하기 위하여 이용한 비집촉식 교란 방법을 사용하였다. 이 경우 추진제 표면으로 전달되는 복사열 전달의 크기는 레이저에 의한 복사 일전달과 기체 영의 화염에 의한 열전달을 동시에 고려하여야 한다. 그러나 언급하였던 것처럼 대부분의 이론적 연구는 추진제 표면의 온도 구매가 단열인 것으로 가정하여 진행하였다. 이러한 가정을 기체 영역으로부터 추진제로 전달되는 열전달 량이 작은 점소화초기 등에서 타당한 가정이나, 기체 영역에서 연소가 활발하게 진행되는 경우에는 비합리적인 가정이다. 본 연구에서는. 추진제의 응축 영역에서 분포 화학 반응이 발생하여, 기체 영역에서 화학반응에 의한 연소가 진행되는 경우, 복사 열전달의 교란에 대한 추진제의 응답함수를 수치적으로 계산하였다. 이때 기체 영역에서 발생하는 연소 반응은 De Luca 등에 의하여 제안된 실험적 모델인 $\alpha$ $\beta$ ${\gamma}$ 화염 모델을 사용하였으며, 추진제 표면에서의 열전달 균형에 의한 경계 조건을 사용하였다. 그러나 외부로부터 입사되는 복사광 레이저와 기체 영역의 상호 간섭은 고려하지 않았다. 수치 계산에 의한 응답 함수의 특징은 단열 조건이 사용된 이론적 응답 함수에 비하여 낮은 값을 나타내었으며, 최대치를 보이는 주파수 영역도 이론 함수에 비하여 다른 값을 보여주고 있다.연구 분석 결과 기술적 문제점으로는 배기 가스온도가 낮은데 따른 출구 부분의 Bearing, Sealing이 문제가 될 수 있다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14

• Proceedings of the KSR Conference
• /
• 2007.11a
• /
• pp.1141-1151
• /
• 2007

In this paper through Honam high-speed railroad which is planned with the north and south axis, we will verify the feasibility of the coordinate conversion using railroad control points after regarding current planned-railroad as the linear central axises. From analysis, distortion of Y axis varies 21cm to 40cm diminishing to a gentle straight line, distortion of X axis varies 14cm to 29cm. Through a revision, the deviation value between the coordinates were 6mm to 9mm and it satisfied the allowable error of national geographic information institute which is following ITRF (International Terrestrial Reference Frame) and cadastral boundary survey(10cm). consequently the coordinate conversion is possible using railroad control points as common control points.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.139-145
• /
• 2009

Reduction scheme is remarkably useful in the case requiring the repeated calculation procedure. Recently, the efficiency of the reduction scheme has been improved by combining scheme of sub-domain method. But, when the global domain is partitioned into a few sub-domains, sub-domains without constraints can be produced. it is needed to extract the ritz vector from each sub-domain to construct the reduced system of each sub-domain. it is easy to extract the ritz vector from sub-domain with constraint. on the other hand, pseudo inverse method should be employed to extract the ritz vector from sub-domain without constraint. generally, the pseudo inverse takes a large number of computing time to obtain a reduced system of a sub-domain without boundary condition. This trouble can be overcome by the reduced pseudo inverse scheme which proposed in this study. This scheme is based on the static condensation that is not related with selection of the primary degrees of freedom. Numerical examples demonstrate that present method saves computational cost effectively and predicts the accurate eigenvalues.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.9
• /
• pp.29-38
• /
• 2015

In this study, track dynamic interaction characteristics caused by the vehicle running through transitional section such as bridge abutments were studied using the finite element analysis program. The geometric condition of track was generated by trigonometric function and allowable maximum track irregularity is determined by KORAIL track maintenance criteria. The sub-infrastructure under rail fastener system was modelled by 3D solid elements. To reduce computational cost only half track line is numerically considered and the roller boundary condition was applied to each side of model. In this study, the vehicle-track dynamic interaction analysis was carried out for standard Korean transition section of concrete track and the dynamic behaviors were investigated. The dynamic characteristics considered are wheel load variation, vertical acceleration at body, and maximum Mises stress at each part of transitional section.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.91-100
• /
• 2004

A Green's function approach is adopted for analyzing the thermoelastic deformations and stresses of a plate made of functionally graded materials(FGMs). The solution to the 3-dimensional unsteady temperature is obtained by using the laminate theory. The fundamental equations for thermoelastic problems are derived in terms of out-plane deformation and in-plane force, separately. The thermoelastic deformation and the stress distributions due to the bending and in-plane forces are analyzed by using a Green's function based on the Galerkin method. The eigenfunctions of the Galerkin Green's function for the thermoelastic deformation and the stress distributions are approximated in terms of a series of admissible functions that satisfy the homogeneous boundary conditions of the rectangular plate. Numerical analysis for a simply supported plate is carried out and effects of material properties on unsteady thermoclastic behaviors are discussed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.173-179
• /
• 2009

Reduction scheme is remarkably useful in the case requiring the repeated calculation procedure. Recently, the efficiency of the reduction scheme has been improved by combining scheme of sub-domain method. But, when the global domain is partitioned into a few sub-domains, sub-domains without constraints can be produced. it is needed to extract the ritz vector from each sub-domain to construct the reduced system of each sub-domain. it is easy to extract the ritz vector from sub-domain with constraint. on the other hand, pseudo inverse method should be employed to extract the ritz vector from sub-domain without constraint. generally, the pseudo inverse takes a large number of computing time to obtain a reduced system of a sub-domain without boundary condition. This trouble can be overcome by the reduced pseudo inverse scheme which proposed in this study. This scheme is based on the static condensation that is not related with selection of the primary degrees of freedom. Numerical examples demonstrate that present method saves computational cost effectively. In addition, it is shown that the reduced system based on the proposed scheme predicts the accurate eigenvalues of global system.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.12
• /
• pp.1811-1820
• /
• 2010

Most components in the real world show nonlinear response. The nonlinearity may arise because of contact between the parts, nonlinear material, or large deformation of the components. Structural optimization considering nonlinearities is fairly expensive because sensitivity information is difficult to calculate. To overcome this difficulty, the equivalent load method was proposed for nonlinear response optimization. This method was originally developed for size and shape optimization. In this study, the equivalent load method is modified to perform topology optimization considering all kinds of nonlinearities. Equivalent load is defined as the load for linear analysis that generates the same response field as that for nonlinear analysis. A simple example demonstrates that results of the topology optimization using equivalent load are very similar to the numerical results. Nonlinear response topology optimization is performed with a practical example and the results are compared with those of conventional linear response topology optimization.