• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Journal of electromagnetic engineering and science
• /
• v.17 no.2
• /
• pp.91-97
• /
• 2017

An accurate and efficient evaluation for hypersingular integrals (HIs), strongly singular integrals (SSIs), and weakly singular integrals (WSIs) plays an essential role in the numerical solutions of 3D electromagnetic scattering problems. We derive analytic formulas for WSIs based on Stokes' theorem, which can be expressed in elementary functions. Several numerical examples are presented to validate these analytic formulas. Then, to show the feasibility of the proposed formulations for numerical methods, these formulations are used with the existing analytical expressions of HIs and SSIs to correct the near-field interaction in an iterative physical optics (IPO) scheme. Using IPO, the scattering caused by a dihedral reflector is analyzed and compared with the results of the method of moments and measurement data.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3498-3506
• /
• 1996

A loading device to be used in fracture experiment is presented. It's loading angle can be adjusted from $-45^{\circ}$ to $105^{\circ}$ at intervals of $15^{\circ}$ for a CTS ( compact tension-shear) specimen, so that it is to be useful to measure mixed mode toughness. The equations to give the $K_ I$, $K_ I1$ and J-integral for the experiment are evluated though finite elemetn analysis in which the loading procedure is simulated and the behaviors of the specimen such as load-displacement curve are estimated. In the course of the evaluation the values $K_ I$, $K_ I1$ and J-integral calculated through recentrly released numerical methods are employed as the reference ones.

• Journal of applied mathematics & informatics
• /
• v.17 no.1_2_3
• /
• pp.73-91
• /
• 2005

The application of Green's function in calculation of flow characteristics around submerged and floating bodies due to a regular wave is presented. It is assumed that the fluid is homogeneous, inviscid and incompressible, the flow is irrotational and all body motions are small. Two methods based on the boundary integral equation method (BIEM) are applied to solve associated problems. The first is a low order panel method with triangular flat patches and uniform distribution of velocity potential on each panel. The second method is a high order panel method in which the kernels of the integral equations are modified to make it nonsingular and amenable to solution by the Gaussian quadrature formula. The calculations are performed on a submerged sphere and some floating spheroids of different aspect ratios. The excellent level of agreement with the analytical solutions shows that the second method is more accurate and reliable.

• Structural Engineering and Mechanics
• /
• v.36 no.6
• /
• pp.669-678
• /
• 2010

In general means, the stress concentration problem of elastic plate with a rectangular hole can be investigated by numerical methods, and only approximative results are derived. This paper deduces an analytical study of the stress concentration due to a rectangular hole in an elastic plate under bending loads. Base on classical elasticity theory and FEM applying the U-transformation technique, the uncoupled governing equations with 3-DOF are established, and the analytical displacement solutions of the finite element equations are derived in series form or double integral form. Therefore, the stress concentration factor can then be discussed easily and conveniently. For the plate subjected to unidirectional bending loads, the non-conforming plate bending element with four nodes and 12-DOF is taken as examples to demonstrate the application of the proposed method. The inner force distribution is obtained. The solutions are adequate for the condition when the hole is far away from the edges and the thin plate subjected to any transverse loadings.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.3 no.3
• /
• pp.185-190
• /
• 1974

Floor panel heating system is popular in Korea as dwelling house heating system. There are two methods for keeping floor surface warm. One method is delivering warm air under the floor such as Roman Hypocaust and Korean traditional Ondol. The other method is imbedding hot water pipes into the concrete floor slab. This paper gives basic equations for steady and transient solutions of heat conduction from hurried pipes. For steady-state solution, fin Efficiency Method and Sink and Source Method were introduced. Sink and Source Method is applied to transient state and basic solution is given in the form of Exponential Integral Function. Numerical solutions can be solved easily by digital computer from these equations.

• Journal of the Korean Geophysical Society
• /
• v.5 no.1
• /
• pp.63-71
• /
• 2002

The integral equation method is a powerful tool for electromagnetic numerical modeling. But the difficulty of this technique is the size of their linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. To overcome this limitation, a lot of approximation and series methods, such as conventional Born, modifed Born and extended Born, were developed. But all the methods need volume integration of Green tensor, which is very time consuming. In electromagnetic theory, Green tensor rapidly decreases as the distance between source and field cell increases. Therefore, the source cell which are far away from the field cell does not make an effect on the electric field of the field cell. Consequently, by ignoring the effect of Green tensor due to far away source cells, computing time for electromagnetic numerical modeling can be reduced dramatically. Comparisons of this new method against a full integral equation, extended Born approximation and series code show that the method is accurate enough much less time consuming.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.1010-1015
• /
• 2002

Over the last few decades, noise has played a major role in the development of aircraft engines. The dominant noise is generated by the wake interactions of fan and downstream stator. Engine inlet and exhaust ducts are being fitted with liner materials that aid in damping fan related noise. In this paper, the radiation of duct internal noise from duct open ends with liners is studies via numerical methods. The linearized Euler's equations in generalized curvilinear coordinates are solved by the DRP scheme. The far field sound pressure levels are computed by the Kirchhoff integration method. Through comparison of sound directivity from bell-mouth duct with and without liners, it is shown that radiation from engine inlet is affected by liner effects or a soft wall boundary condition.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1408-1416
• /
• 1990

A powerful and efficient method for finding approximate solutions to initial-boundary-value problems in the transient coupled thermoelasticity is formulated in time domain using the finite element technique with time-marching strategy. The final system equations can be derived by the Guritin's variational principle using the definition of convolution integral. But, the finite element formulation for the equations of motion is modified by differentiating in time. Numerical results to some test problems are compared with analytical and other sophisticated approximate solutions. Stable responces are observed in all the given examples irrespective of incremental time steps and mesh shapes. In addition, it is shown that good numerical results are obtained even in coarser mesh or larger time step comparing to other numerical methods.

• Interaction and multiscale mechanics
• /
• v.4 no.2
• /
• pp.85-105
• /
• 2011

The influence of surface elasticity and surface residual stress on the elastic field of an isotropic nanoscale elastic layer of finite thickness bonded to a rigid material base is considered by employing the Gurtin-Murdoch continuum theory of elastic material surfaces. The fundamental solutions corresponding to buried vertical and horizontal line loads are obtained by using Fourier integral transform techniques. Selected numerical results are presented for the cases of a finite elastic layer and a semi-infinite elastic medium to portray the influence of surface elasticity and residual surface stress on the bulk stress field. It is found that the bulk stress field depends significantly on both surface elastic constants and residual surface stress. The consideration of out-of-plane terms of the surface stress yields significantly different solutions compared to previous studies. The solutions presented in this study can be used to examine a variety of practical problems involving nanoscale/soft material systems and to develop boundary integral equations methods for such systems.