• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• Journal of Bio-Environment Control
• /
• v.27 no.1
• /
• pp.20-26
• /
• 2018

This study was conducted to develop a device for measuring the air flow by space variation through monitoring program, which acquires data by each point from each environmental sensor located in the greenhouse. The distribution of environmental factors(air temperature, humidity, wind speed, etc.) in the greenhouse is arranged at 12 points according to the spatial variation and a large number of measurement points (36 points in total) on the X, Y and Z axes were selected. Considering data loss and various greenhouse conditions, a bit rate was at 125kbit/s at low speed, so that the number of sensors can be expanded to 90 within greenhouse with dimensions of 100m by 100m. Those system programmed using MATLAB and LabVIEW was conducted to measure distributions of the air flow along the greenhouse in real time. It was also visualized interpolated the spatial distribution in the greenhouse. In order to verify the accuracy of CFD modeling and to improve the accuracy, it will compare the environmental variation such as air temperature, humidity, wind speed and $CO_2$ concentration in the greenhouse.

• Nuclear Engineering and Technology
• /
• v.52 no.9
• /
• pp.2092-2099
• /
• 2020

Many damages due to flow-accelerated corrosion and cracking have been observed during recent in-service inspections of nuclear power plants. To determine the operability or repair for damaged pipes, an integrity evaluation related to the damaged piping system should be performed by using already proven code and standards. One of them, the ASME Code Case is most popularly used to integrity assessment in nuclear power plants. However, the recent version of CC N-513 still recommends the simplified method which means a damaged elbow is assumed as an equivalent straight pipe. In addition, to enhance the accuracy integrity assessment in elbow, several previous studies recommend that the SIF and elastic COD values for an elbow with relatively large crack could be predicted by an interpolation technique. However, those estimates for elbow with relatively large crack might be derived to inaccurate results for crack growth analysis, such as for the allowable crack size and life estimation. Therefore, in this paper, the effect of crack length (0.3≤θ₁/π≤0.5) on SIF and elastic COD for elbow is systematically investigated. Then, for large crack in elbow, accurate estimates for SIF and elastic COD, which are widely used to assess the integrity of elbows, are proposed. Those proposed solutions are expected to be the technical basis for revisions of CC N-513-4 through the validation.

• Journal of Environmental Science International
• /
• v.6 no.5
• /
• pp.541-550
• /
• 1997

A least-square regression analysis is applied for the estimation of velocity streamfunction field based on discretely sampled current meter data. The coefficients of a streamfuunction that is expanded in terms of trigonometric basis function are obtained by enforcing the horizontal non-divergence of two-dimensional flow field. This method avoids Interpolation and gives a root-mean-square (rms) residual of fit which Includes the divergent part and noisiness of oceanic data. The Implementation of the method Is done by employing a boundary-fitted, curvilinear orthogonal coordinate which facilitates the specification of boundary conditions. An application is successfully made to the Texas-Louisiana shelf using the 32 months current meter data (31 moorings) observed as a part of the Texas-Louisiana Shelf and Transport Processes Study (LATEX). The rms residual of the fitting is relatively small for the shelf, which indicates the field Is Ivell represented by the streamnfunction.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.197-201
• /
• 2007

Vibration characteristics of an annular cylindrical fuel rod, which was proposed as a candidate design of fuel's cross section for the ultra-high burnup nuclear fuel, according to the cross-sectional dimensions and the number of supports or the span length were analytically studied. Finite element(FE) modeling for the annular cross sectional fuel was based on the methodology, that have been proven by the test verification, for the conventional PWR nuclear fuel rod. A commercial FEA code, ABAQUS, was used for the FE modeling and analysis. A planar beam element (B21) that uses a linear interpolation was used for the fuel rod and a linear spring element for the spring and dimple of the SG. Natural frequencies and mode shape were calculated according to the preliminary design candidates for the fuel's cross sectional dimension and the number of span. From the analysis results, the design scheme of the annular fuel compatible to the present PWR nuclear reactor core was discussed in terms of the number of supports and fuel's cross section.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.18 no.1
• /
• pp.47-53
• /
• 1982

An analysis of convection, in a fully developed laminar steady flow through the vertical square duct under the condition of variational symmetric heat flux, is considered. Finite element solution algorithm by Galerkin's method with triangular elements and linear interpolation polynominals for the temperature and velocity profiles are derived for the vertical square duct. The comparison of temperature distribution due to variational symmetric heat flux in the duct were made with available the other data when the condition of peripheral heat flux were uniform and zero. Numerical values for the dimensionless temperatures and Nusselt numbers at selected Rayleigh numbers and pressure gradient parameters were obtained at a few nodal points for the vertical square ducts and effects of corner in the duct were investigated.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.1
• /
• pp.15-20
• /
• 2017

The main purpose of the ball valve ball is to be moved by the rotation of the stem when fully open or completely closed. In this study the heat of the initial model, which used a structure interaction analysis technique, tried to examine the structural safety of the high temperature for the ball valve. In the initial model the stress of the exiting sheet was more than the yield strength. We selected two design shapes with variables of length and thickness for the optimization of the sheet. The Kriging interpolation method was applied to a meta-model-based optimization technique. As a result, it was possible to find a thickness and length for the sheet within the yield strength. This was done by measuring the value of the capacity coefficient of the valve and evaluating the performance of the ball valve.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.2
• /
• pp.17-25
• /
• 2006

In the present study, a 2-D chimera technique for overlapped bodies in relative motion is developed using unstructured triangular meshes. The solid boundary nodes located next to the intersecting point between bodies are merged to the intersecting point to assure accurate representation of the intersecting region. In order to assign proper value of flow variables at the nodes located out of the computational field, interpolation is conducted for non-active nodes. For validation, the motions of a NACA64A006 airfoil and a NACA0012 airfoil with a plane flap are computed and the results are compared with other simulations. The motion of a launching missile ejected from a NACA0012 airfoil is also simulated.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.4 no.3
• /
• pp.133-140
• /
• 2002

Site-specific minimum temperature forecasts are critical in a short-term decision making procedure for preventive measures as well as a long-term strategy such as site selection in fruits industry. Nocturnal cold air pools frequently termed in mountainous areas under anticyclonic systems are very dangerous to the flowering buds in spring over Korea, but the spatial resolution to detect them exceeds the current weather forecast scale. To supplement the insufficient spatial resolution of official forecasts, we developed a GIS - assisted frost risk assesment scheme for using in mountainous areas. Daily minimum temperature data were obtained from 6 sites located in a 2.1 by 2.1 km area with complex topography near the southern edge of Sobaek mountains during radiative cooling nights in spring 2001. A digital elevation model with a 10 m spatial resolution was prepared for the entire study area and the cold air inflow was simulated for each grid cell by counting the number of surrounding cells coming into the processing cell. Primitive temperature surfaces were prepared for the corresponding dates by interpolating the Korea Meteorological Administration's automated observational data with the lapse rate correction. The cell temperature values corresponding to the 6 observation sites were extracted from the primitive temperature surface, and subtracted from the observed values to obtain the estimation error. The errors were regressed to the flow accumulation at the corresponding cells, delineating a statistically significant relationship. When we applied this relationship to the primitive temperature surfaces of frost nights during April 2002, there was a good agreement with the observations, showing a feasibility of site-specific frost warning system development in mountainous areas.

• Journal of the Korea Society for Simulation
• /
• v.29 no.3
• /
• pp.41-48
• /
• 2020

In this paper, a dynamic prediction algorithm using the cumulative distribution function for traffic volume is presented as a new method for predicting highway traffic rate more accurately, where an approximation function of the cumulative distribution function is obtained through numerical methods such as natural cubic spline interpolation and Levenberg-Marquardt method. This algorithm is a new structure of random number generation algorithm using the cumulative distribution function used in financial mathematics to be suitable for predicting traffic flow. It can be confirmed that if the highway traffic rate is simulated with this algorithm, the result is very similar to the actual traffic volume. Therefore, this algorithm is a new one that can be used in a variety of areas that require traffic forecasting as well as highways.