• Journal of applied mathematics & informatics
• /
• 제9권2호
• /
• pp.561-577
• /
• 2002

In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집B
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회(1)
• /
• pp.406-406
• /
• 2005

• Composites Research
• /
• 제34권2호
• /
• pp.129-135
• /
• 2021

본 논문에서는 필라멘트 와인딩 공정에서 발생하는 두께 방향 섬유체적비 불균일이 원통형 복합재 격자 구조의 좌굴하중에 미치는 영향을 확인하기 위해서 Vasiliev가 제안한 원통형 복합재 격자 구조 좌굴하중 이론식을 변형하여 섬유체적비에 따른 좌굴하중 저하를 확인하였다. 섬유체적비에 따라 격자 구조 리브의 각 층의 두께를 달리하였으며, 혼합법칙을 사용하여 각 층별로 물성치를 다르게 적용하였다. 구조물 크기, 두께, 섬유체적비 평균값을 달리한 유한요소모델에 대한 선형좌굴해석을 수행하였다. 최종적으로 이론식을 사용한 등가모델의 좌굴 하중 계산 결과와 유한요소해석 결과를 비교하여 두께 방향 섬유체적비 불균일이 원통형 복합재 격자 구조의 좌굴하중 저하의 원인이 될 수 있음을 확인하였다.

• Nuclear Engineering and Technology
• /
• 제41권8호
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• Nuclear Engineering and Technology
• /
• 제45권7호
• /
• pp.839-846
• /
• 2013

This paper describes a finite element model of the microstructure of dispersion type nuclear fuels, which can be used to determine the effective thermal conductivity of the fuels during irradiation. The model simulates a representative region of the fuel as a prism shaped unit cell made of brick elements. The elements within the unit cell are assigned material properties of either the fuel or the matrix depending on position, in such a way as to represent randomly distributed fuel particles with a size distribution similar to that of the as manufactured fuel. By applying an appropriate heat flux across the unit cell it is possible to determine the effective thermal conductivity of the unit cell as a function of the volume fraction of the fuel particles. The presence of a fuel/matrix interaction layer is simulated by the addition of a third set of material properties that are assigned to the finite elements that surround each fuel particle. In this way the effective thermal conductivity of the material may also be determined as a function of the volume fraction of the interaction layer. Work is on going to add fission gas bubbles in the fuel as a fourth phase to the model.

• 한국수자원학회논문집
• /
• 제44권3호
• /
• pp.249-262
• /
• 2011

본 연구에서는 댐 붕괴파와 같이 연속 및 불연속 흐름해석에 적용되고 있는 HLLC 기법을 불규칙한 하상지형에서의 흐름해석에 적용할 때 생성항과 흐름률항의 사이의 수치적 불균형으로 인한 수치진동을 감소시키기 위해 well-balanced HLLC 기법과 천수방정식에 근거한 비구조적 유한체적모형을 개발하였으며, 이를 댐 붕괴파 문제에 적용하였다. 적용된 well-balanced HLLC 기법은 단순히 흐름률항을 계산할 때 하상지형경사를 직접 포함시키는 것으로 정상상태의 천이류에 적용하였을 때 생성항과 흐름률항 사이에 수치적 균형이 이루어짐을 확인하였다. 수치모형의 검증을 위해 댐 붕괴파 문제와 관련된 세 가지 서로 다른 수리모형실험과 프랑스 Mapasset 댐 붕괴에 대한 현장사례에 적용하였으며, 적용결과 본 연구에서 개발된 모형은 수리모형실험 그리고 현장에서 관측된 결과와 비교적 잘 일치하는 것으로 나타났으며, 기존의 모형으로부터 계산된 모의결과에 비해 비교적 정확한 결과를 나타내었다.

• 한국수자원학회논문집
• /
• 제36권1호
• /
• pp.1-11
• /
• 2003

본 연구에서는 2차원 부정류 해석을 위한 유한체적 수치모형을 개발하였다. 개발된 모형에서는 천수방정식의 수치해를 구하기 위하여 Lax-Friedrichs 해법에 근간을 둔 ENO (essentially non-oscillatory) 기법을 적용하였으며, 공간적으로 2차 정도 모의를 위하여 수정 MUSCL 기법을 적용하였다. 개발된 모형의 적용성을 평가하기 위하여 해석해가 있는 1차원 댐 붕괴파 모의에 적용한 결과, 해석해와 유사한 결과를 나타냈으며, 수정 MUSCL 기법의 모의결과를 minmod 경사제한자를 사용한 모의결과와 비교한 결과, 수정 MUSCL 기법을 사용하는 경우 해석해와 더 근접한 결과를 나타냈다. 댐붕괴파의 2차원 해석에 대한 모형의 적용성을 평가하기 위하여 대칭형, 비대칭형 댐 붕괴에 적용한 결과, 저수지에서 수로 구간으로 충격파 전파를 잘 모의했으며, 기존의 연구결과와 유사한 모의결과를 나타냈다. 개수로의 부정류 해석에 대한 모형의 적용성을 평가하기 위하여 말단이 폐쇄된 수로의 입구에서 정상파가 주어졌을 때 수위와 유속을 모의하여 해석해와 비교하였으며, 입구와 말단이 폐쇄된 수로에서 파의 진동을 모의하여 모형의 거동을 평가하였다.

• 한국해양공학회지
• /
• 제11권4호
• /
• pp.141-151
• /
• 1997

This paper deals with a numerical study of flow-guider applied to controlling current in a bay. Two dimensional numerical model for tidal currents based on the depth averaged equation is developed and standard k-.epsilon. model is adopted to determine the turbulence diffusion. Equations are described in a generalized coordinate system to be implemented by non-staggered grid system and discretized by using finite volume method. Unsteady flow is simulated by fully implicit scheme. Hybrid scheme and central differencing are used to compute the convective terms and source terms, respectively. The tidal current in a rectangular bay is simulated and it gives satisfactory results. The realistic and distinct models of a large structure placed in bay are also exemplified with or without flow-guiders. The simulation results show that the flow-guider gives the residual tidal current in the bay by the different flux with respect to the direction of tidal current.

• 설비공학논문집
• /
• 제11권6호
• /
• pp.733-739
• /
• 1999

A modified low-Re $k-\varepsilon$ model is used for the calculation of drag-reducing turbulent flow by polymer injection in a pipe. With the viscoelastic model, molecular viscosity in the definition of turbulent viscosity is related to elongations viscosity of the solution to account for the effects of drag reduction. Finite volume method is used for the discretization, and power-law scheme is used as a numerical scheme. Computed dimensionless velocity profiles are in good agreements with the experimental data in case of low drag reductions. However, in case of high drag reductions, they deviate largely from the measurements in the central zone of the flow field.