• Journal of Electrical Engineering and Technology
• /
• v.1 no.3
• /
• pp.381-386
• /
• 2006

In dielectric waveguide analysis and synthesis, we often encounter an awkward task of solving the eigenvalue equation to find the value of propagation constant. Since the dispersion equation is an irrational equation, we cannot solve it directly. Taking advantage of approximated calculation, we attempt here to solve this irrational dispersion equation. A new type of eigenvalue equation, in which guide index is expressed as a function of frequency, has been developed. In practical optical waveguide designing and in calculating the propagation mode, this equation will be used more conveniently than the previous one. To expedite the design of the waveguide, we then solve the eigenvalue equation of a slab waveguide, which is sufficiently accurate for practical purpose.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.34-39
• /
• 1996

A consistent general order nodal method for solving the 3-D neutron diffusion equation in (x-y-z) geometry has ben derived by using a weighted integral technique and expanding the spatial variables by the Legendre orthogonal series function. The equation set derived can be converted into any order nodal schemes. It forms a compact system for general order of nodal moments. The method utilizes the analytic solutions of the transverse-integrated quasi -one dimensional equations and a consistent expansion for the spatial variables so that it renders the use of an approximation for the transverse leakages no necessary. Thus, we can expect extremely accurate solutions and the solution would converge exactly when the mesh width is decreased or the approximation order is increased since the equation set is consistent mathematically.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.161-165
• /
• 2004

In the vortex particle method based on the vorticity-velocity formulation for solving the Wavier-Stokes equations, the unsteady, incompressible, viscous laminar flow over a NACA 0012 foil is simulated. By applying an operator-splitting method, the 'convection' and 'diffusion' equations are solved sequentially at each time step. The convection equation is solved using the vortex particle method, and the diffusion equation using the particle strength exchange(PSE) scheme which is modified to avoid a spurious vorticity flux. The scheme is improved for variety body shape using one image layer scheme. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsively started NACA 0012 foil for Reynolds number 550.

• Journal of the Korean Society of Safety
• /
• v.12 no.4
• /
• pp.180-190
• /
• 1997

A new method in the fault tree analysis (FTA) for the reliability calculation is suggested. Two steps are necessary in traditional method in evaluation of the occurrence probability of top event in fault tree (FT). The first step is to find the minimal outsets, and the second one is to substitute the result into the poincare equation. In order to reduce the enormous computing time of this method, lots of rapid algorithms have been developed. Almost of all achievements were, however, based on the partial structural properties of FT. In this paper, the FT is transformed to a non-linear graph G which has the same minimal outsets of original n, and then the reliability is calculated using the domination theory. In this new method, the required number of equation terms are at most $2^n$ (n is node number of graph G), while $2^m$-1 (m is the number of minimal cutsets) calculation terms are required in the poincare equation in traditional method. Since m>>n in general. our new method reduces the calculation time significantly.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.364-370
• /
• 2005

A stable walking pattern generation method for a biped robot is presented in this paper. In general, the ZMP (zero moment point) equations, which are expressed as differential equations, are solved to obtain a stable walking pattern. However, the number of differential equations is less than that of unknown coordinates in the ZMP equations. It is impossible to integrate the ZMP equations directly since one or more constraint equations are involved in the ZMP equations. To overcome this difficulty, DAE (differential and algebraic equation) solution method is employed. The proposed method has enough flexibility for various kinematic structures. Walking simulation for a virtual biped robot is performed to demonstrate the effectiveness and validity of the proposed method. The method can be applied to the biped robot for stable walking pattern generation.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.13 no.6
• /
• pp.480-489
• /
• 1988

Although CDP stacking of common depth gathering is used to get the zero-offset-section, the exploding reflector concept is examined for the modeling of zero source to receiver offset sections in this paper. The acoustic wave equation is compared with a one way wave equation which represents the upgoing wave field only. The one way wave equation used is not derived through an expansion and, therefore, can represent dips up do 90b degrees and may not lost the signals by the dipping angles. There is apparently no simple counterpart of this equation is the space domain and it can be conveniently implemented only by a Fourier method. This paper compares their modeling technique with ray tracing and wave method for over thrust structure which is one of the geological structures are dificult to process and interpret. As a result of modeling much clean and accurate signals, especially, diffractions form the corner and dipping angles can be gathered.

• Structural Engineering and Mechanics
• /
• v.57 no.5
• /
• pp.823-835
• /
• 2016

Using the Complementary Functions Method (CFM), a general solution for the one-dimensional steady-state thermal and mechanical stresses in a hollow thick sphere made of functionally graded material (FGM) is presented. The mechanical properties are assumed to obey the exponential variations in the radial direction, and the Poisson's ratio is assumed to be constant, with general thermal and mechanical boundary conditions on the inside and outside surfaces of the sphere. In the present paper, a semi-analytical iterative technique, one of the most efficient unified method, is employed to solve the heat conduction equation and the Navier equation. For different values of inhomogeneity constant, distributions of radial displacement, radial stress, circumferential stress, and effective stress, as a function of radial direction, are obtained. Various material models from the literature are used and corresponding temperature distributions and stress distributions are computed. Verification of the proposed method is done using benchmark solutions available in the literature for some special cases and virtually exact results are obtained.

• The Journal of the Acoustical Society of Korea
• /
• v.15 no.2
• /
• pp.72-78
• /
• 1996

This thesis deals with a method to solve a two-and-one-half-dimensional ($2\frac12$ D) problem, which means that the ocean environment is two-dimensional whereas the source is fully three-dimensionally propagating, including three-dimensional refraction phenomena and three-dimensional back-scattering, using two-dimensional two-way parabolic equation method combined with Fourier synthesis. Two dimensional two-way parabolic equation method uses Galerkin's method for depth and Crank-Nicolson method and alternating direction for range and provides a solution available to range-dependent problem with wave-field back-scattered from discontinuous interface. Since wavenumber, k, is the function of depth and vertical or horizontal range, we can reduce a dimension of three-dimensional Helmholtz equation by Fourier transforming in the range direction. Thus transformed two-dimensional Helmholtz equation is solved through two-way parabolic equation method. Finally, we can have the $2\frac12$ D solution by inverse Fourier transformation of the spectral solution gained from in the last step. Numerical simulation has been carried out for a canonical ocean environment with stair-step bottom in order to test its accuracy using the present analysis. With this spectral parabolic equation method, we have examined three-dimensional acoustic propagation properties in a specified site in the Korean Straits.

• Journal of Energy Engineering
• /
• v.16 no.1 s.49
• /
• pp.46-52
• /
• 2007

The piecewise constant angular approximation is developed to replace the conventional angular quadrature sets in the solution of the second-order, multi-dimensional $S_{N}$ neutron transport equations. The newly generated quadrature sets by this method substantially mitigate ray effects and can be used in the same manner as the conventional quadrature sets are used. The discrete-ordinates and the piecewise-constant approximations are applied to both the first-order Boltzmann and the second-order form of neutron transport equations in treating angular variables. The result is that the mitigation of ray effects is only achieved by the piecewise-constant method, in which new angular quadratures are generated by integrating angle variables over the specified region. In other sense, the newly generated angular quadratures turn out to decrease the contribution of mixed-derivative terms in the even-parity equation that is one of the second-order neutron transport equation. This result can be interpreted as the entire elimination or substantial mitigation of ray effect are possible in the simplified even-parity equation which has no mixed-derivative terms.