• 제어로봇시스템학회:학술대회논문집
• /
• 1997.10a
• /
• pp.1750-1753
• /
• 1997

This paper describes how the directivity pattern of the back-scattered sound pressure is distributed when a plane acoustic wave is incident on a righid spherical shell underwater. A coupled Finite Element-Boundary Element mehtod is developed as numerical technique. The result of the coupled FE-BE method is agreed with theoretical solution for algorithmic confirmation.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.5
• /
• pp.57-68
• /
• 1995

A general procedure for the numerical solution of coupled, nonlinear, differential two-point boundary-value problems, solutions of which are crucial to the controller design, has been developed and demonstrated. A fixed-end-points, free-terminal-time, optimal-control problem, which is derived from Pontryagin's Maximum Principle, is solved by an extension of Davidenko's method, a differential form of Newton's method, for algebraic root finding. By a discretization process like finite differences, the differential equations are converted to a nonlinear algebraic system. Davidenko's method reconverts this into a pseudo-time-dependent set of implicitly coupled ODEs suitable for solution by modern, high-performance solvers. Another important advantage of Davidenko's method related to the time-optimal problem is that the terminal time can be computed by treating this unkown as an additional variable and sup- plying the Hamiltonian at the terminal time as an additional equation. Davidenko's method uas used to produce optimal trajectories of a single-degree-of-freedom problem. This numerical method provides switching times for open-loop control, minimized terminal time and optimal input torque sequences. This numerical technique could easily be adapted to the multi-point boundary-value problems.

• Journal of the Korean Mathematical Society
• /
• v.37 no.1
• /
• pp.31-44
• /
• 2000

We study a free boundary problem satisfying Bernoulli type boundary condition along which the gradient of a piecewise harmonic solution jumps zero to a given constant value. In such problem, the free boundary splits the domain into two regions, the zero set and the harmonic region. Our main interest is to identify the global shape and the location of the zero set. In this paper, we find the lower and the upper bound of the zero set. In a convex domain, easier estimation of the upper bound and faster disk test technique are given to find a rough shape of the zero set. Also a simple proof on the convexity of zero set is given for a connected zero set in a convex domain.

• Korean Journal of Materials Research
• /
• v.19 no.11
• /
• pp.588-595
• /
• 2009

The knowledge of grain growth of carbide particles is very important for manufacturing micrograined cemented carbides. In the present study, continuous and discontinuous grain growth in WC-Co and WC-VC-Co cemented carbides is investigated using the Monte Carlo computer simulation technique. The Ostwald ripening process (solution/re-precipitation) and the grain boundary migration process are assumed in the simulation as the grain growth mechanism. The effects of liquid phase fraction, grain boundary energy and implanted coarse grain are examined. At higher liquid phase content, mass transfer via solid/liquid interfaces plays a major role in grain growth. Growth rate of the implanted grain was higher than that of the matrix grains through solution/re-precipitation and coalescence with neighboring grains. The results of these simulations qualitatively agree with experimental ones and suggest that distribution of liquid phase and carbide particle/carbide grain boundary energy as well as contamination by coarse grain are important factors controlling discontinuous grain growth in WC-Co and WC-VC-Co cemented carbides. The contamination by coarse grains must by avoided in the manufacturing process of fine grain cemented carbides, especially with low Co.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.9 no.1
• /
• pp.43-54
• /
• 1997

This paper presents an approximate analytical solution to a two-region one-dimensional model for the charging process of stratified thermal storage tanks with variable inlet temperature in the presence of momentum-induced mixing. Based on the superposition principle, an arbitrary-varying inlet temperature is decomposed into inherent discontinuous steps and continuous intervals approximated as a finite number of piecewise linear functions. This approximation allows the temperature of the upper perfectly-mixed layer to be expressed in terms of constant, linear and exponential functions with respect to time. Applying the Laplace transform technique to the model equation for the lower thermocline layer subject to each of three representative interfacial conditions yields compact-form solutions, a linear combination of which constitutes the final temperature profile. A systematic method for deriving solutions to the plug-flow problem having polynomial-type boundary conditions is also established. The effect of adiabatic exit boundary on solution behaviors proves to be negligible under the actual working conditions, which justifies the assumption of semi-infinite domain introduced in the solution procedure. Finally, the approximate solution is validated by comparing it with an exact solution obtained for a specific variation of inlet temperature. Excellent agreements between them suffice to show the necessity and utility of this work.

• East Asian mathematical journal
• /
• v.26 no.3
• /
• pp.403-413
• /
• 2010

In this paper, we consider a doubly nonlinear parabolic equation related to the Leray-Lions operator with Dirichlet boundary condition and initial data given. By exploiting a suitable implicit time-discretization technique, we obtain the existence of global strong solution.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.2 s.245
• /
• pp.149-156
• /
• 2006

A boundary-based design sensitivity analysis(DSA) technique is proposed for addressing shape optimization issues in the elastostatics problems. Sensitivity formula is derived based on the continuum formulation in a boundary integral form, which consists of the boundary solutions and shape variation vectors. Though the boundary element method(BEM) has been mainly used to obtain the boundary solution, the FEM is used in this paper because this is much more popular, and has greatly improved meshing and computing power recently. The advantage of the boundary DSA is that the shape variation vectors, which are also known as design velocity fields, are needed only on the boundary. Then, the step for determining the design velocity field over the whole domain, which was necessary in the domain-based DSA, is eliminated, making the process easy to implement and efficient. Problem of fillet design is chosen to illustrate the efficiency of the proposed method. Accuracy of the sensitivity is good with this method even by employing the free mesh for the FE analysis.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.05a
• /
• pp.707-712
• /
• 1995

Grain boundary diffusion plays significant role in the fission gas release, which is one of the crucial processes dominating nuclear fuel performance. Gaseous fission products such as Xe and Kr generated inside fuel pellet have to diffuse in the lattice and in the grain boundary before they reach open space in the fuel rod. In the mean time, the grains in the fuel pellet grow and shrink according to grain growth kinetics, especially at elevated temperature at which nuclear reactors are operating. Thus the boundary movement ascribed to the grain growth greatly influences the fission gas release rate by lengthening or shortening the lattice diffusion distance, which is the rate limiting step. Sweeping fission gases by the moving boundary contributes to the increment of the fission gas release as well. Lattice and grain boundary diffusion processes in the fission gas release can be studied by 'tracer diffusion' technique, by which grain boundary diffusion can be estimated and used directly for low burn-up fission gas release analysis. However, even for tracer diffusion analysis, taking both the intragranular grain growth and the diffusion processes simultaneously into consideration is not easy. Only a few models accounting for the both processes are available and mostly handle them numerically. Numerical solutions are limited in the practical use. Here in this paper, an approximate analytical solution of the lattice and stationary grain boundary diffusion in a polycrystalline solid is developed for the tracer diffusion techniques. This short closed-form solution is compared to available exact and numerical solutions and turns out to be acceptably accurate. It can be applied to the theoretical modeling and the experimental analysis, especially PIE (post irradiation examination), of low burn up fission. gas release.

• Tunnel and Underground Space
• /
• v.12 no.3
• /
• pp.171-178
• /
• 2002

A dynamic analysis of a horseshoe_shaped tunnel near to cavity was performed to study the effect of the cavity on the dynamic behavior of the tunnel. In order to obtain the dynamic response of the tunnel embedded in a semi-infinite domain, a hybrid numerical technique was primarily developed. A dynamic fundamental solution in frequency domain for multi-layered half planes was derived and subsequently incorporated in the boundary element method. Coupling of the boundary element method for the far field with the finite element method for the near field is made by imposing compatibility condition of a displacement at the interface. The boundary element method is then coupled with the finite element method, which is utilized to model the near field including the tunnel and the cavity. In order to demonstrate the validity of the proposed technique, dynamic responses of single and multiply-layered semi-infinite structural systems are obtained by using the Kicker waveform and investigated in the limestone layer to find how the being and the location of the cavity affect the dynamic characteristics of the system.

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.