• Journal of Environmental Impact Assessment
• /
• v.8 no.4
• /
• pp.47-57
• /
• 1999

This study is aimed at the development of comprehensive environmental management system, which can be operated on the basis of world wide web, as a topic of G7 project. Even though there should be lots of works remaining to achieve this goal, preliminary products can be summarized as follows : 1) integrated environmental information management system, 2) web based control engine, 3) surface water environment management system, 4) subsurface water environment management system, 5) sewer and waterworks management system. The core methodology of the engine is the generalized multidimensional finite element matrices to depict the terms in the analysis of various partial differential equations. Spatial information management system (ArcView) and Visual Basic were extensively employed to construct GUI oriented web based engine. The developed systems were composed of very intense computer codes due to the necessity of combinatory management of environmental problems. The web based engine could be served as a decision tool for the integrated management of environmental projects in Korea.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.1 no.2
• /
• pp.14-19
• /
• 2010

Theories for advanced composite structures are too difficult for such design engineers for construction and some simple but accurate enough methods are necessary. The senior author has reported that some laminate orientations have decreasing values of $D_{16}$, $B_{16}$, $D_{26}$ and $B_{26}$ stiffnesses as the ply number increases. For such plates, the fiber orientations given above behave as specially orthotropic plates and simple formulas developed by the senior author. Most of the bridge and building slabs on girders have large aspect ratios. For such cases further simplification is possible by neglecting the effect of the longitudinal moment terms(Mx) on the relevant partial differential equations of equilibrium. In this paper. the influence of the aspect ratio on the natural frequency of the composite laminated plates is studied and it is concluded that the method used is sufficiently accurate for engineering purposes.

• Journal of Korean Society on Water Environment
• /
• v.23 no.5
• /
• pp.683-688
• /
• 2007

Water flow into unsaturated soils is most often modeled by Richards equation consisting of the mass conservation law and Darcy's law. Three standard forms of Richards equation are presented as the head (${\Psi}$)-based form, the moisture content (${\theta}$) based form, and the mixed form. Numerical solutions of these partial differential equations with highly nonlinear terms can cause poor results along with significant mass balance errors. The numerical solution based on the mixed form of Richards equation is known that the mass is perfectly conserved without any additional computational efforts. The aim of this study is to develop fully implicit numerical scheme of Richards equation for one-dimensional vertical unsaturated flow in homogeneous soils using the finite difference approximation, and then to perform sensitivity analysis of infiltration to the variations in the unsaturated soil properties and to different soil types.

• Fire Science and Engineering
• /
• v.26 no.3
• /
• pp.79-84
• /
• 2012

In this study thermal characteristics of heat-responsive element considering conduction, convection and rate of change of element using Response Time Index (RTI) applied to sensitivity test of sprinkler head at home and aborad are theoretically investigated. Analytic solution of temperature distributions with radial direction and time is obtained form energy transport equations, non-homogeneous 2th order partial differential equation, applying to constant wall temperature and symmetric condition in order to analyze thermal characteristics of heat-responsive element for circular cylindrical geometry. Base on the results, the analytic method of this study is fundamental data to practical use for sensitivity test of sprinkler head and design of heat-responsive element.

• Journal of applied mathematics & informatics
• /
• v.11 no.1_2
• /
• pp.305-316
• /
• 2003

Recently, an iterative algorithm for finding the interior eigenvalues of a definite matrix by CG-type method has been proposed. This method compares to the inverse power method. The given matrices A, and B are assumed to be large and sparse, and SPD( Symmetric Positive Definite) The CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for smallest eigenvalue. Also, it is very amenable to parallel computations, like the CG method for the linear systems. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. But for parallel computations we need to find an efficient parallel preconditioner. Our candidates we ILU(0) in the wave-front order, ILU(0) in the multi-coloring order, Point-SSOR(Symmetric Successive Overrelaxation), and Multi-Color Block SSOR preconditioner. Wavefront order is a simple way to increase parallelism in the natural order, and Multi-coloring realizes a parallelism of order(N), where N is the order of the matrix. Another choice is the Multi-Color Block SSOR(Symmetric Successive OverRelaxation) preconditioning. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test problem was drawn from the discretizations of partial differential equations by finite difference methods. The results show that for small number of processors Multi-Color ILU(0) has the best performance, while for large number of processors Multi-Color Block SSOR performs the best.

• Journal of Radiation Protection and Research
• /
• v.28 no.2
• /
• pp.97-108
• /
• 2003

A three-dimensional wind field model based on the variational technique has been developed for estimating the overall wind patterns over a complex terrain. The three-dimensional elliptic partial differential equations on Cartesian and terrain-following coordinates have been established to obtain the Lagrangian multiplier and the adjusted wind velocity. The simulations were performed to evaluate the variations of the velocity vectors on the hemisphere, half-cylinder, and saddle type obstacles. Also, the wind field model in the terrain-following coordinate has been applied for evaluating the characteristics of wind patterns according to the variations of Gauss precision moduli on the hemispheric topography. The results showed that the horizontal and vertical wind components were strongly governed by the selection of the values of Gauss precision moduli.

• Smart Structures and Systems
• /
• v.18 no.2
• /
• pp.355-374
• /
• 2016

This paper presents and compares a one-dimensional (1D) bending theory for piezoelectric thin beam-type structures with resistive-inductive electrodes to ANSYS$^{(R)}$ three-dimensional (3D) finite element (FE) analysis. In particular, the lateral deflections and vibrations of slender piezoelectric beams are considered. The peculiarity of the piezoelectric beam model is the modeling of electrodes in such a manner that is does not fulfill the equipotential area condition. The case of ideal, perfectly conductive electrodes is a special case of our 1D model. Two-coupled partial differential equations are obtained for the lateral deflection and for the voltage distribution along the electrodes: the first one is an extended Bernoulli-Euler beam equation (second-order in time, forth order in space) and the second one the so-called Telegrapher's equation (second-order in time and space). Analytical results of our theory are validated by 3D electromechanically coupled FE simulations with ANSYS$^{(R)}$. A clamped-hinged beam is considered with various types of electrodes for the piezoelectric layers, which can be either resistive and/or inductive. A natural frequency analysis as well as quasi-static and dynamic simulations are performed. A good agreement between the extended beam theory and the FE results is found. Finally, the practical relevance of this type of electrodes is shown. It is found that the damping capability of properly tuned resistive or resistive-inductive electrodes exceeds the damping performance of beams, where the electrodes are simply linked to an optimized impedance.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.6
• /
• pp.219-225
• /
• 2011

Advanced composite structures are too difficult for such design engineers for construction and some simple but accurate enough methods are necessary. The simply supported laminated plates are analyzed by the specially orthotropic laminates theory. This method, however, may be too difficult for some practising engineers. In this paper, the result of analysis for such plate by means of the beam theory with unit width is reported. The plate aspect ratio considered is from 1 : 1 to 1 : 5. Most of the bridge and building slabs on girders have large aspect ratios. For such cases further simplification is possible by neglecting the effect of the longitudinal moment terms($M_x$) on the relevant partial differential equations of equilibrium. In this paper. the influence of the aspect ratio on the natural frequency of the specially orthotropic laminated plates is studied and it is concluded that the method used is sufficiently accurate for engineering purposes. The result of this paper can be used for simply supported laminated plates analysis.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.17 no.2
• /
• pp.71-76
• /
• 1981

In order to make a numerical modeling for the one dimensional unsteady flow which expressed by Saint Venant partial differential equations, Preissmann's implicit scheme was used, and it's stability and accuracy was mentioned. By introducing recurrence relations make it possible to use double sweep algorithm. Effective parameters to the result were the values of the $\Delta$t/$\Delta$x, $\theta$ and the chezy coefficient. In oder to get numerical solutions with enough accuracy, $\Delta$t/$\Delta$x should not be far from the value of 1, and the criteria of the $\theta$ was 0.6<$\theta$<1.0 for the stability without condition. This model should be calibrated by real field data, and expected to be developed for the simulation of the river system and to the long wave analysis for one dimensional coastal zone problem.

• Journal of Aerospace System Engineering
• /
• v.14 no.3
• /
• pp.17-23
• /
• 2020

In this paper, the efficient periodic unsteady flow analysis is developed by using a Chebyshev collocation operator applied to the time differential term of the governing equations. The partial implicit time integration method was also applied in the governing equation for a fluid, which means flux terms were implicitly processed for a time integration and the time derivative terms were applied explicitly in the form of the source term by applying the Chebyshev collocation operator. To verify this method, we applied the 1D unsteady Burgers equation and the 2D oscillating airfoil. The results were compared with the existing unsteady flow frequency analysis technique, the Harmonic Balance Method, and the experimental data. The Chebyshev collocation operator can manage time derivatives for periodic and non-periodic problems, so it can be applied to non-periodic problems later.