• Water for future
• /
• v.25 no.1
• /
• pp.93-100
• /
• 1992

A study on the Derivation of the Unit Hydrograph using Multiple Regression Moe이. The purpose of this study is to deriver an optimal unit hydrograph suing the multiple regression model, particularly when only small amount of data is available. The presence of multicollinearity among the input data can cause serious oscillations in the derivation of the unit hydrograph. In this case, the oscillations in the unit hydrograph ordinate are eliminated by combining the data. The data used in this study are based upon the collection and arrangement of rainfall-runoff data(1977-1989) at the Soyang-river Dam site. When the matrix X is the rainfall series, the condition number and the reciprocal of the minimum eigenvalue of XTX are calculated by the Jacobi an method, and are compared with the oscillation in the unit hydrograph. The optimal unit hydrograph is derived by combining the numerous rainfall-runoff data. The conclusions are as follows; 1)The oscillations in the derived unit hydrograph are reduced by combining the data from each flood event. 2) The reciprocals of the minimum eigen\value of XTX, 1/k and the condition number CN are increased when the oscillations are active in the derived unit hydrograph. 3)The parameter estimates are validated by extending the model to the Soyang river Dam site with elimination of the autocorrelation in the disturbances. Finally, this paper illustrates the application of the multiple regression model to drive an optimal unit hydrograph dealing with the multicollinearity and the autocorrelation which cause some problems.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.8 no.5
• /
• pp.495-503
• /
• 1997

In this paper, the E-polarized electromagnetic scattering problems by a resistive strip grating with tapered resistivity on 3 dielectric layers are analyzed to find out the effects for the tapered resistivity of resistive strip and the relative permittivity and thickness of 3 die- lectric layers by applying the Fourier-Galerkin moment methods. The induced surface current density is expanded in a series of Jacobi-polynomial ${P^{(\chi,\beta)}}_p$(.) of the order $\alpha$= 0 and $\beta$=1 as a kind of orthogonal polyomians, and the tapered resistivity assumes to vary linearly from 0 at one edge to finite resistivity at the other edge. The normalized reflected and transmitted powers are obtained by varying the tapered resistivity and the relative permittivity and thickness of dielectric layers. The sharp variation points are observed when the higher order modes are transferred between propagating and evanescent modes, and in general the local minimum positions occur at less grating period for the more relative permittivity of dielectric layers. It should be noted that the patterns of the normalized reflected and transmitted powers for the tapered resistivity are very much different from those of the uniform resistivity and perfectly conducting cases. The proposed method of this paper cna solve the scattering problems for the tapered resistive, uniform resistive, and PEC strip cases.

• The Journal of Engineering Research
• /
• v.2 no.1
• /
• pp.113-123
• /
• 1997

The purpose of this paper is an optimum design of the shaped Cassegrainian antenna system for the base station. The process of the shaped Cassegrainian antenna design is as follows : 1) the aperture field distribution is determined so as to meet design specifications, 2) a proper design parameter is selected, 3) extracting of the dimension data for the main and sub-reflector antenna To do these, Hansen's distribution is chosen as the aperture field, and the far-field pattern from the aperture is predicted by the angular spectrum. Firstly, the aperture field distribution is designed to satisfy the specification for design frequency, it is confirmed if this distribution meet the specification for another frequency band. The main- and the sub-reflectors are synthesized so as for the given beamwaveguide feed pattern to be transformed into the prescribed aperture distribution. The designed system has circular aperture, left-right symmetry and no tilted structure. The continuous surface functions of reflectors are obtained by adopting the global interpolation technique to the discrete reflector profiles. Jacobi polynomial-sinusoidal is used as the basis function. A Ka-band Cassegrainian antenna operates over 17.7 – 20.2 GHz for down-link band and 27.5 – 30 GHz for up-link band is designed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• 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.

• International Journal of Control, Automation, and Systems
• /
• v.4 no.6
• /
• pp.689-696
• /
• 2006

This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of weakly coupled bilinear systems with time-varying parameter uncertainties and exogenous disturbance using the successive Galerkin approximation(SGA). By using weak coupling theory, the robust $H_{\infty}$ control can be obtained from two reduced-order robust $H_{\infty}$ control problems in parallel. The $H_{\infty}$ control theory guarantees robust closed-loop performance but the resulting problem is difficult to solve for uncertain bilinear systems. In order to overcome the difficulties inherent in the $H_{\infty}$ control problem, two $H_{\infty}$ control laws are constructed in terms of the approximated solution to two independent Hamilton-Jacobi-Isaac equations using the SGA method. One of the purposes of this paper is to design a closed-loop parallel robust $H_{\infty}$ control law for the weakly coupled bilinear systems with parameter uncertainties using the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

• Journal of the Korean Mathematical Society
• /
• v.38 no.3
• /
• pp.595-611
• /
• 2001

Let Q(n,1) be the set of even unimodular positive definite integral quadratic forms in n-variables. Then n is divisible by 8. For A[X] in Q(n,1), the theta series $\theta$(sub)A(z) = ∑(sub)X∈Z(sup)n e(sup)$\pi$izA[X] (Z∈h (※Equations, See Full-text) the complex upper half plane) is a modular form of weight n/2 for the congruence group Γ$_1$(8) = {$\delta$∈SL$_2$(Z)│$\delta$≡()mod 8} (※Equation, See Full-text). If n$\geq$24 and A[X], B{X} are tow quadratic forms in Q(n,1), the quotient $\theta$(sub)A(z)/$\theta$(sub)B(z) is a modular function for Γ$_1$(8). Since we identify the field of modular functions for Γ$_1$(8) with the function field K(X$_1$(8)) of the modular curve X$_1$(8) = Γ$_1$(8)＼h(sup)* (h(sup)* the extended plane of h) with genus 0, we can express it as a rational function of j(sub) 1,8 over C which is a field generator of K(X$_1$(8)) and defined by j(sub)1,8(z) = $\theta$$_3$(2z)/$\theta$$_3$(4z). Here, $\theta$$_3$ is the classical Jacobi theta series.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.1623-1628
• /
• 2004

In the 3-dimensional cyclic Lotka-Volterra equations, we show the solution on the invariant hyperplane. In addition, we show the existence of the invariant hyperplane by the center manifold theorem under the some conditions. With this result, we can lead the hyperplane of the n-dimensional cyclic Lotka-Volterra equaions. In other section, we study the 3- or 4-dimensional Hamiltonian Lotka-Volterra equations which satisfy the Jacobi identity. We analyze the solution of the Hamiltonian Lotka- Volterra equations with the functions called the split Liapunov functions by [4], [5] since they provide the Liapunov functions for each region separated by the invariant hyperplane. In the cyclic Lotka-Volterra equations, the role of the Liapunov functions is the same in the odd and even dimension. However, in the Hamiltonian Lotka-Volterra equations, we can show the difference of the role of the Liapunov function between the odd and the even dimension by the numerical calculation. In this paper, we regard the invariant hyperplane as the important item to analyze the motion of Lotka-Volterra equations and occur the chaotic orbit. Furtheremore, an example of the asymptoticaly stable and stable solution of the 3-dimensional cyclic Lotka-Volterra equations, 3- and 4-dimensional Hamiltonian equations are shown.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1510-1517
• /
• 1995

A theoretical method is described to simulate the propagation of turbulent premixed flames in a closed vessel. The objective is to develop and test an efficient technique to predict the propagation speed of flame as well as the geometric structure of the flame surfaces. Flame is advected by the statistically generated turbulent flow field and propagates as a wave by solving twodimensional Hamilton-Jacobi equation. In the simulation of the unburned gas flow field, following turbulence properties were satisfied: mean velocity field, turbulence intensities, spatial and temporal correlations of velocity fluctuations. It is assumed that these properties are not affected by the expansion of the burned gas region. Predictions were compared with existing experimental data for flames propagating in a closed vessel charged with hydrogen/air mixture with various turbulence intensities and Reynolds numbers. Comparisons were made in flame radius growth rate, rms flame radius fluctuations, and average perimeter and fractal dimensions of the flame boundaries. Two dimensional time dependent simulation resulted in correct trends of the measured flame data. The reasonable behavior and high efficiency proves the usefulness of this method in difficult problems of flame propagation such as in internal combustion engines.

• The Journal of Information Technology
• /
• v.8 no.2
• /
• pp.77-84
• /
• 2005

In this paper, electromagnetic scattering problems by a resistive strip grating with tapered resistivity on a grounded dielectric plane according to strip width and spacing, relative permittivity and thickness of dielectric layers, and incident angles of a electric wave are analyzed by applying the Fourier-Galerkin Moment Method known as a numerical procedure. The boundary conditions are applied to obtain the unknown field coefficients and the resistive boundary condition is used for the relationship between the tangential electric field and the electric current density on the strip. The resistivity of resistive strips in this paper varies from zeroes at one edge to infinite at the other edge, then the induced surface current density on the resistive strip is expanded in a series of Jacobi polynomials of the order ${\alpha}=0.2,\;{\beta}=-0.2$ as a orthogonal polynomials. The numerical results of the geometrically normalized reflected power in this paper are compared with those for the existing perfectly conducting strip. The numerical results of the normalized reflected power for conductive strips case with zero resistivity in this paper show in good agreement with those of existing papers.