• Korean Journal of Soil Science and Fertilizer
• /
• v.36 no.4
• /
• pp.181-192
• /
• 2003

We developed a mathematical simulation model to portray the vertical distribution of soil water from the measured weather data and the known soil hydraulic properties, and then compared simulation results with the periodically measured soil water profiles obtained on Jungdong sandy loam to verify the model, In this model, we solved potential-based Richards' equation by the implicit finite difference method superimposed on the predictor-corrector scheme. We presumed that: soil hydraulic properties are homogeneous; soil water flows isothermally; hysteresis is not considered; no vapor flows; no heat transfers into the soil profiles; and water added to soil surface is distributed along the soil profile following partial displacement principle. The input data were broadly classified into two groups: (1) daily weather data such as rainfall, maximum and minimum air temperatures, relative humidity and solar radiation and (2) soil hydraulic data to approximate unsaturated hydraulic conductivity and water retention. Each hydraulic polynomial function approximated using the Chebyshev polynomial and least square difference technique in tandem showed a fairly good fit of the given set of data. Vertical distribution of soil water as approximations to the Richards' equation subject to changing surface and phreatic boundaries was solved numerically during 53 days with a comparatively large time increment, and this pattern agreed well with field neutron scattering data, except for the surface 0.1 m slab.

• Proceeding of KASS Symposium
• /
• 2008.05a
• /
• pp.107-112
• /
• 2008

This paper investigates the variation of post-buckling behaviours of 2-dimensional shallow arch type truss after sizing optimization. The mathematical programming technique is used to produce the optimum member size of 2D arch truss against a central point load. Total weight of structure is considered as the objective function to be minimized and the displacement occurred at loading point and member stresses of truss are used as the constraint functions. The finite difference method is used to calculate the design sensitivity of objective function with respect to design variables. The postbuckling analysis carried out by using the geometrically nonlinear finite element analysis code ISADO-GN. It is found to be that there is a huge change of post-buckling behaviour between the initial structure and optimum structure. Numerical results can be used as useful information for future research of large spatial structures.

• Journal of the Korean Institute of Gas
• /
• v.17 no.3
• /
• pp.20-26
• /
• 2013

In this research, an empirical study on the designed burst pressure of Type3 composite cylinder was performed by hydrostatic burst test equipment. The designed burst pressure of Type3 composite cylinders, which are 6.8 liter and 31 MPa of service pressure, was estimated with the analysis using the finite element method. In order to confirm its accuracy, the burst test of small Type3 composite cylinders was perfomed through three times. The burst test equipment can pressurize to 400 MPa. As a result of comparison between the designed burst pressure and actual burst pressure, the difference was less than 4 percentage. With a test result, the analysis accuracy was verified. This technique will be applied to both qualification and inspection for the composite cylinder.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.16 no.4 s.95
• /
• pp.373-379
• /
• 2005

Papers on improvement of electromagnetic field uniformity in a reverberation chamber with 1D Quadratic Residue Diffuser of Schroeder method has been published several times. In this paper, to obtain improved electromagnetic field characteristics and field uniformity in a reverberation chamber, cubical residue diffuser sets of Schroeder type are designed for a chamber in $2.3\;\cal{GHz}\~3\;\cal{GHz}$. The FDTD(Finite-Difference Time-Domain) technique is used to analyze the field characteristics in a chamber. Cubical residue algorithm and 2D arrangement show more randomness than the previous study results. The characteristics of tolerance, polarity, deviations, as well as power efficency, are improved with cubical residue diffuser sets in a chamber.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.8
• /
• pp.67-73
• /
• 2019

An ultra-high pressure treatment device is a device used for increasing the shelf life of food by sterilizing it by applying hydrostatic pressure to solid or liquid food. The ultrahigh pressure treatment system developed in this study is a pressure vessel with a processing capacity of 100 L and a maximum pressure of 700 MPa. Pressure vessels for ultrahigh-pressure processing equipment are manufactured using wire-winding techniques. The design formula for making ultra-high pressure vessels with wire windings is given in ASME Section VIII - Division 3. In this study, the ratio of the cylinder to the winding area that can be applied in a wire-winding application was analyzed using a finite element analysis. Furthermore, the relationship between the variation of the residual stress in the vessel and the ratio of the winding area due to the variation of the winding tension was analyzed, and a design guide applicable to the actual product design was developed. Finally, the design equation was modified by presenting the coefficients to correct the difference between the finite element analysis and the design equation.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.27 no.5
• /
• pp.120-129
• /
• 2023

Hanger cable tension is a major response that can determine the integrity and safety of suspension bridges. In general, the vibration method is used to estimate hanger cable tension on operational suspension bridges. It measures natural frequencies from hanger cables and indirectly estimates tension using the geometry conditions of the hanger cables. This study estimated the hanger cable tension of the Palyeong Bridge using a vision-based system. The vision-based system used digital camcorders and tripods considering the convenience and economic efficiency of measurement. Measuring the natural frequencies for high-order modes required for the vibration method is difficult because the hanger cable response measured using the vision-based system is displacement-based. Therefore, this study proposed a back analysis technique for estimating tension using the natural frequencies of low-order modes. Optimization for the back analysis technique was performed by defining the difference between the natural frequencies of hanger cables measured in the field and those calculated using finite element analysis as the objective function. The direct search method that does not require the partial derivatives of the objective function was applied as the optimization method. The reliability and accuracy of the back analysis technique were verified by comparing the tension calculated using the method with that estimated using the vibration method. Tension was accurately estimated using the natural frequencies of low-order modes by applying the back analysis technique.

• Journal of applied mathematics & informatics
• /
• v.25 no.1_2
• /
• pp.435-443
• /
• 2007

Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.309-313
• /
• 2004

Prediction of green water loads acting an the bow deck is an essential part for the design of bow structures against the green water impact. Proper technique of the green water simulation is highly required for the prediction of green water loads. In this paper, a new numerical method for green water simulation, which is based an predictor-corrector-upwind finite difference scheme of the 2nd kind, is introduced. Through the comparisons between computed' results and experimental measurements, it is verified that the present numerical tool is adequate as a practical calculation tool for the green water problem.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.7
• /
• pp.1211-1217
• /
• 2000

The three-dimensional finite-difference time-domain (FDTD) method and the two-dimensional quasi-static formulation have been used to calculate the characteristic impedance and the microwave effective index of coplanar waveguide structures on Lithium Niobate ($LiNBO_3$) single crystal substrates with a yttria-stabilized zirconia (YSZ) or $SiO_2$ buffer layer. The results shown can be a good source to predict the modulator characteristics. The effects of the thin buffer layer and anisotropy of the $LiNBO_3$ crystal (x-cut and z-cut) are discussed. The comparison between the FDTD and quasi-static results shows good agreement. In this paper, the efficient modeling technique of the FDTD method for the coplanar waveguide (CPW) structures based on an anisotropic substrate with a thin buffer layer is developed.