• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.4
• /
• pp.143-152
• /
• 2018

Numerical convergence and stability of TVD schemes have been applied in the FUNWAVE-TVD model were compared. The fourth order accurate MUSCL-TVD scheme using minmod limiter suggested by Yamamoto and Daiguji (1993), the fourth order accurate MUSCL-TVD scheme using van-Leer limiter suggested by Erduran et al. (2005) and the second order accurate MUSCL-TVD scheme using van-Leer limiter in Zhou et al. (2001) were compared. Comparisons of the numerical scheme were conducted with experimental data of Vincent and Briggs irregular wave experiments. In comparison with the fourth order accurate scheme using van-Leer limiter, the fourth order accurate scheme using minmod limiter is less dissipative but required lower CFL condition for stable numerical solution. On the other hand, the scheme using van-Leer limiter required smaller resolution spatial grid due to numerical dissipation, but relatively higher CFL condition can be used compared to the scheme using minmod limiter. In the breaking wave experiments which were conducted using high resolution spatial grid to reduce numerical dissipation, the characteristic of the schemes can be clearly observed. Numerical instabilities and blow-up of the numerical solutions were found in the irregular wave breaking simulation with the scheme using minmod limiter. However, the simulation can be completed with the scheme using van-Leer limiter, but required low CFL condition. Good agreements with the observed data were also observed in the results using van-Leer limiter.

• Interaction and multiscale mechanics
• /
• v.3 no.2
• /
• pp.111-121
• /
• 2010

Dynamic analysis of nano/micro bio-sensors based on a multiscale atomistic/continuum theory is introduced. We use a generalized atomistic finite element method (GAFEM) to analyze a bio-sensor which has $3{\times}N_a{\times}N_p$ degrees of freedom, where $N_p$ is the number of representative unit cells and $N_a$ is the number of atoms per unit cell. The stiffness matrix is derived from interatomic potential between pairs of atoms. This work contains two studies: (1) the resonance analysis of nano bio-sensors with different amount of target analyte and (2) the dependence of resonance frequency on finite element mesh. We also examine the Courant-Friedrichs-Lewy (CFL) condition based on the highest resonance frequency. The CFL condition is the criterion for the time step used in the dynamic analysis by GAFEM. Our studies can be utilized to predict the performance of micro/nano bio-sensors from atomistic perspective.

• Atmosphere
• /
• v.33 no.4
• /
• pp.331-341
• /
• 2023

A numerical forecasting models usually predict future states by performing time integration considering fixed static time-steps. A time-step that is too long can cause model instability and failure of forecast simulation, and a time-step that is too short can cause unnecessary time integration calculations. Thus, in numerical models, the time-step size can be determined by the CFL (Courant-Friedrichs-Lewy)-condition, and this condition acts as a necessary condition for finding a numerical solution. A static time-step is defined as using the same fixed time-step for time integration. On the other hand, applying a different time-step for each integration while guaranteeing the stability of the solution in time advancement is called an adaptive time-step. The adaptive time-step algorithm is a method of presenting the maximum usable time-step suitable for each integration based on the CFL-condition for the adaptive time-step. In this paper, the adaptive time-step algorithm is applied for the Korean Integrated Model (KIM) to determine suitable parameters used for the adaptive time-step algorithm through the monthly verifications of 10-day simulations (during January and July 2017) at about 12 km resolution. By comparing the numerical results obtained by applying the 25 second static time-step to KIM in Supercomputer 5 (Nurion), it shows similar results in terms of forecast quality, presents the maximum available time-step for each integration, and improves the calculation efficiency by reducing the number of total time integrations by 19%.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.25 no.9
• /
• pp.738-742
• /
• 2012

This paper conducted a study on how the heat radiation of light emitting diode(LED) luminaires affects the indoor temperature increase. The effect was compared with that of a 20 W compact fluorescent lamp(CFL) and a 50 W MR16 halogen lamp which are most widely used inside of cruises, a LED downlight and a 4W MR16 LED replacing each of them. We installed a luminarie inside a thermally shielded chamber, measuring the temperature changes under the same volume every 5 minutes and compared the result with theoretically calculated heat radiation. The temperature changes in the chamber was measured four times, on seven hours' period in order to keep sufficient time once the temperature reaches the thermal equilibrium state. The results showed that the temperature of the 20 W E26 CFL and the 10 W LED downlight increased by $21.1^{\circ}C$ and $10.4^{\circ}C$ respectively, while that of the 50 W halogen MR16 and the 4 W LED MR16 increased by $33.9^{\circ}C$ and $4.8^{\circ}C$ respectively. The experimental heat radiation were calculated from the results and the experimental heat radiation of the CFL and the LED downlight were 171.5 cal and 86.5 cal, and those of the halogen MR16 and the LED MR16 were 275.3 cal and 36.5 cal. Therefore, the heat radiation was reduced by 49.5% and 86.7%, respectively, by replacing conventional light source with LED. In conclusion, we can expect a reduction of power consumption in air condition system and the effect on indoor temperature increase by application of LED luminaires.

• Environmental Engineering Research
• /
• v.21 no.2
• /
• pp.171-179
• /
• 2016

This study analyzes guidelines to select optimum number of grids to represent behavior of a given water system appropriately. The EFDC model was chosen as a 3-D hydrodynamic and water quality model and salt was chosen as a surrogate variable of pollutant. The model is applied to an artificial canal that receives salt water from coastal area and fresh water from a river from respective gate according to previously developed gate operation rule. Grids are subdivided in vertical and horizontal (longitudinal) directions, respectively until no significant changes are found in salinity concentrations. The optimum grid size was determined by comparing errors in average salt concentrations between a test grid systems against the most complicated grid system. MSE (mean squared error) and MAE (mean absolute error) are used to compare errors. The CFL (Courant-Friedrichs-Lewy) number was used to determine the optimum number of grid systems for the study site though it can be used when explicit numerical method is applied only. This study suggests errors seem acceptable when both MSE and MAE are less than unity approximately.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.166-176
• /
• 2018

This paper presents a numerical sensitivity analysis for the simulation of the motion performance of an offshore structure in waves using computational fluid dynamics (CFD). Starting with 2D wave simulations with varying numerical parameters such as grid spacing and CFL value, proper numerical conditions were found for accurate wave propagation that avoids numerical diffusion problems. These results were mapped on 2D error distributions of wave amplitude and wave length against the numbers of grids per wave length and per wave height under a given CFL condition. Finally, the 2D numerical sensitivity result was validated through CFD simulation of the motion of a FPSO in waves showing good accuracy in motion RAOs compared with existing potential flow solutions.

• Journal of the Chungcheong Mathematical Society
• /
• v.28 no.4
• /
• pp.623-633
• /
• 2015

In this paper, a general Weighted Essentially Non-Oscillatory (WENO) interpolation is proposed and applied to a semi-Lagrangian method. The proposed method is based on the conservation law, and characteristic curves are used to complete the semi-Lagrangian method. Therefore, the proposed method satisfies conservation of mass and is free of the CFL condition which is a necessary condition for convergence. Using a several standard examples, the proposed method is compared with the third order Strong Stability Preserving (SSP) Runge-Kutta method to verify the high-order accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.116-120
• /
• 2009

A two-dimensional unsteady inviscid flow solver has been developed for the simulation of complex geometric configurations on adaptive Cartesian meshes. Embedded condition was used for boundary condition and a predictor-corrector explicit time marching scheme was used for time-accurate numerical simulation. The Cartesian mesh generator, which was previously developed for steady problem, was used grid generation for unsteady flow. The solver was based on ALE formulation for body motion. For diminishing the effects of cut-cells, the cell merging method was used. Using cell merging method, it was eliminated the CFL constraints. The conservation problem, which is caused cell-type variation around region swept by solid boundary, was also solved using cell merging method. The results are presented for 2D circular cylinder and missile launching problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).