• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.414-418
• /
• 2005

A finite element model is developed for the numerical simulation of bed elevation change in a curved channel. The SU/PG (Streamline-Upwind/Petrov-Galerkin) method is used to solve 2D shallow water equations and the BG (Bubnov-Galerkin) method is used for the Exner equation. For the time derivative terms, the Crank-Nicolson scheme is used. The developed model is a decoupled model in a sense that the bed elevation does not change simultaneously with the flow during the computational time step. The total load formula with is used for the sediment transport model. The slip conditions are described along the lateral boundaries. The effects of gravity force due to geometry change and the secondary flows in a curved channel are considered in the model. For the verification, the model is applied to two laboratory experiments. The first is $140^{\circ}$ bended channel data at Delft Hydraulics Laboratory and the second is $140^{\circ}$ bended channel data at Laboratory of Fluid Mechanics of the Delft University of Technology. The finite element grid is constructed with linear quadrilateral elements. It is found that the computed results are in good agreement with measured data, showing a point bar at the inner bank and a pool at the outer bank.

• 대한기계학회논문집B
• /
• 제21권11호
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• 한국수산과학회지
• /
• 제15권2호
• /
• pp.154-160
• /
• 1982

1) 수-단백질-지질의 3성분계 식품으로서 '두부'를 선택하여 일차원 비정상 열전도에 관한 실험을 행하였다. 두부는 불균질계 물질이지란 '유효'한 개념을 도입하여 유효열확산율을 측정함으로써 수학적으로 균질계 물질로 취급이 가능했다. 2) 저자가 고안 제작한 장치의 성능(performance)및 정도를 검토하기 위하여 표준물질로서 toluene은 물론 얼음(빙)의 열확산율을 측정한 결과, $-5\sim-25^{\circ}C$의 온도역에서 문헌치와 비교하여 재현성 $\pm1.5\%$, 오차 $-2\%$ 정도로 그 신뢰도는 매우 높았다. 3) 경계조건이 경시적으로 변화하는 경우라 찰지라도 Crank-Nicolson's implicit method를 이용하여 열확산율을 구할 수 있는 computer program을 개발했다. 4) 성분조성이 각각 상이한 3종류의 두부의 열확산율을 측정한 결과, 그 양들은 수분, 단백질 및 지질의 각 성분 함량의 분율에 정확하게 의존하고 있음을 보여 주었다.

• Journal of Biosystems Engineering
• /
• 제31권4호
• /
• pp.355-362
• /
• 2006

A one-dimensional heat conduction problem in cylindrical coordinate system was solved using Crank-Nicolson finite difference method to predicting the temperature distribution in rice storage bin with wall insulator. The model can simulate the grain temperatures in insulated round bins using the input data of initial grain temperature. ambient air temperature, wind velocity, solar radiation on a horizontal surface, and thermal properties of grain, bin wall, wall insulator, insulator cover, and air. Temperatures were collected at the bin center, 0.65m in radial direction from the center, and near the bin wall in 2.7m diameter bin filled with rough rice to depth of 3.0m were used to validate the simulation model. Grain temperatures predicted by the model were in very good agreement with the measured temperatures. The residual mean square error between measured and predicted grain temperatures at the bin center was $1.38^{\circ}C$.

• 설비공학논문집
• /
• 제24권3호
• /
• pp.218-223
• /
• 2012

This study is focused on the evaluation of thermal output of TABS (Thermally Activated Building System). The aim of this study is to evaluate TABS in terms of the temperature difference between heating medium supply temperature ($T_s$) and return temperature ($T_r$), thermal output and the surface temperature distribution according to the design flow rate and the design flow temperature. Through the transient heat transfer simulation using temperature calculation using Crank-Nicolson FDM using Physibel Voltra 6.0 W, the temperature difference between $T_s$ and ��$T_r$, thermal output and the surface temperature distribution of specific TABS was calculated and evaluated. The results show that specific thermal output and temperature difference at $60^{\circ}C$ of supply water temperature were about 162 $W/m^2$, $13.6^{\circ}C$ respectively.

• Journal of the Optical Society of Korea
• /
• 제12권2호
• /
• pp.112-117
• /
• 2008

Single-mode propagation conditions of X-ray waveguides are investigated by numerical calculations in order to understand the importance of waveguide design parameters, such as core thickness and the optical constants of waveguide materials, on the transmission and coherence properties of the waveguide. The simulation code for mode analyzing is developed based on a numerical solution of the parabolic wave equation. The initial boundary value problem is solved numerically using a finite-difference scheme based on the Crank-Nicolson scheme. The E-field intensities in a core layer are calculated at an X-ray energy of 8.0 keV for air and beryllium(Be) core waveguides with different cladding layers such as Pt, Au, W, Ni and Si to determine the dependence on waveguide materials. The highest E-field intensity radiated at the exit of the waveguide is obtained from the Pt cladded beryllium core with a thickness of 20 nm. However, the intensity from the air core waveguide with Pt cladding reaches 64% of the Be-Pt waveguide. The dependence on the core thickness, which is the major parameter used to generate a single mode in the waveguide, is investigated for the air-Pt, and Be-Pt waveguides at an X-ray energy of 8.0 keV. The mode profiles at the exit are shown for the single mode at a thickness of up to 20 nm for the air-Pt and the Be-Pt waveguides.

• International Journal of Fuzzy Logic and Intelligent Systems
• /
• 제3권1호
• /
• pp.13-17
• /
• 2003

A continuous solution of the Dirichlet boundary value problem for the heat equation $u_t$＝$a2u_{xx}$ using a fuzzy system is described. We first apply the Crank-Nicolson method to obtain a discrete solution at the grid points for the heat equation. Then we find a continuous function to represent approximately the discrete values at the grid points in the form of a bicubic spline function (equation omitted) that can in turn be represented exactly by a fuzzy system. We show that the computed values at non-grid points using the bicubic spline function is much smaller than the ones obtained by linear interpolations of the values at the grid points. We also show that the fuzzy rule table in the fuzzy system representation of the bicubic spline function can be viewed as a gray scale image. Hence, the fuzzy rules provide a visual representation of the functions of two variables where the contours of different levels for the function are shown in different gray scale levels

• 한국재료학회지
• /
• 제29권9호
• /
• pp.547-552
• /
• 2019

In the present study, the thermal conductivity of a silicon nitride($Si_3N_4$) thin-film is evaluated using the dual-wavelength pump-probe technique. A 100-nm thick $Si_3N_4$ film is deposited on a silicon (100) wafer using the radio frequency plasma enhanced chemical vapor deposition technique and film structural characteristics are observed using the X-ray reflectivity technique. The film's thermal conductivity is measured using a pump-probe setup powered by a femtosecond laser system of which pump-beam wavelength is frequency-doubled using a beta barium borate crystal. A multilayer transient heat conduction equation is numerically solved to quantify the film property. A finite difference method based on the Crank-Nicolson scheme is employed for the computation so that the experimental data can be curve-fitted. Results show that the thermal conductivity value of the film is lower than that of its bulk status by an order of magnitude. This investigation offers an effective way to evaluate thermophysical properties of nanoscale ceramic and dielectric materials with high temporal and spatial resolutions.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2011년도 정기 학술발표대회
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.