• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.568-571
• /
• 2008

The computation of moving interface by the level set method typically requires reinitializations of level set function. An inaccurate estimation of level set function ${\phi}$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, accurate and robust reinitialization process is essential to the free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates directly level set function ${\phi}$ using a normal vector in the interface without solving the re-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1splitting FEM are adopted to discretize advection equation of the level set function and the Navier-Stokes equation, respectively. Advection equation of free surface and re-initialization process are validated with benchmark problems, i.e., a broken dam flow and time-reversed single vortex flow. The simulation results are in good agreement with the existing results.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.2
• /
• pp.160-168
• /
• 2008

This paper presents the study on the stability of nonlinear oscillations of a thin cantilever beam subject to harmonic base excitation in vertical direction. Two partial differential governing equations under combined parametric and external excitations were derived and converted into two-degree-of-freedom ordinary differential Mathieu equations by using the Galerkin method. We used the method of multiple scales in order to analyze one-to-one combination resonance. From these, we could obtain the eigenvalue problem and analyze the stability of the system. From the thin cantilever experiment using foamax, we could observe the nonlinear modes of bending, twisting, sway, and snap-through buckling. In addition to qualitative information, the experiment using aluminum gave also the quantitative information for the stability of combination resonance of a thin cantilever beam under parametric excitation.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.2 s.35
• /
• pp.317-326
• /
• 1998

This paper describes a reduced finite element formulation for nonlinear transient heat transfer analysis based on Lanczos Algorithm. In the proposed reduced formulation all material nonlinearities of irradiation boundary element are included using the pseudo force method and numerical time integration of the reduced formulation is conducted by Galerkin method. The results of numerical examples demonstrate the applicability and the accuracy of the proposed method for the nonlinear transient heat transfer analysis.

• Proceedings of the KIEE Conference
• /
• 1994.07a
• /
• pp.147-149
• /
• 1994

In designing high efficiency electrical machines, calculation of iron loss is very important. And it is reported that in the induction motor and in the T-joint of 3 phase transformer, there occurred rotational magnetic field and much iron loss is generated owing to this field. In this paper, rotational power loss in the electrical machine under rotational magnetic field is discussed. Until now, loss analysis is based on the magnetic properties under alternating field. And with this one dimensional magnetic propertis, it is difficult to express iron loss under rotational field. In this paper, we used two dimensional magnetic property data for the numerical calculation of rotational power loss. We used finite element method for calculation and the analysis model is two dimensional magnetic property measurement system. We used permeability tensor instead of scalar permeability to present two dimensional magnetic properties. And in this case, we cannot uniquely define energy functional because of the asymmetry of the permeability tensor, so Galerkin method is used for finite element analysis.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.19 no.5
• /
• pp.857-871
• /
• 1994

In this paper, six kinds of spiral antenna, a combination of two types of spiral arm-width and three types of spiral curvature are analyzed by using moment method. Dividing spiral arms into N sections, the current distribution is calculated by Galerkin`s method. The radiation pattern and the antenna gain are derived from antenna currents. All os the six spiral antenna have amni-dirctional and wide-band characteristics, although the antenna gain changes within +_ 5dB bound for operating range(600MHz-2GHz). The variation of antenna`s gain is caused by the return loss in connection the Balun to the antenna. Simulation and experimental results on the radiation pattern also show spiral antennas have omni-directional and wide-band characteristics.

• Wind and Structures
• /
• v.17 no.1
• /
• pp.1-19
• /
• 2013

A multiscale finite element method is applied to the Spalart-Allmaras turbulence model based detached-eddy simulation (DES). The multiscale arises from a decomposition of the scalar field into coarse (resolved) and fine (unresolved) scales. It corrects the lack of stability of the standard Galerkin formulation by modeling the scales that cannot be resolved by a given spatial discretization. The stabilization terms appear naturally and the resulting formulation provides effective stabilization in turbulent computations, where reaction-dominated effects strongly influence near-wall predictions. The multiscale DES is applied in the context of high-Reynolds flow over the Commonwealth Advisory Aeronautical Council (CAARC) standard tall building model, for both uniform and turbulent inflows. Time-averaged pressure coefficients on the exterior walls are compared with experiments and it is demonstrated that DES is able to resolve the turbulent features of the flow and accurately predict the surface pressure distributions under atmospheric boundary layer flows.

• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.295-303
• /
• 2018

The semi-analytical method to study the nonlinear dynamic behavior of simply supported spiral stiffened functionally graded (FG) cylindrical shells subjected to an axial compression is presented. The FG shell is surrounded by damping and linear/nonlinear elastic foundation. The proposed linear model is based on the two-parameter elastic foundation (Winkler and Pasternak). A three-parameter elastic foundation with hardening/softening cubic nonlinearity is used for nonlinear model. The material properties of the shell and stiffeners are assumed to be FG. Based on the classical plate theory of shells and von $K{\acute{a}}rm{\acute{a}}n$ nonlinear equations, smeared stiffeners technique and Galerkin method, this paper solves the nonlinear vibration problem. The fourth order Runge-Kutta method is used to find the nonlinear dynamic responses. Results are given to consider effects of spiral stiffeners with various angles, elastic foundation and damping coefficients on the nonlinear dynamic response of spiral stiffened simply supported FG cylindrical shells.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.6 no.1
• /
• pp.19-31
• /
• 1996

To invetigate the impurity distribution in GaAs crystal grown by horizontal Bridgman method, we constructd the mathematical model describing heat transfer, mass transfer and fluid flow n transient growth of GaAs. Galerkin finite element method and implicit time integration were used to solve the equations and simulate the transient growth. The concentration distribution is similar to the case of diffusion controlled growth when Gr - 0. With the increase of Gr the concentration profile is distroted and the minimum solute concentration appears near the interface. As solidification prosceeds, interface deflection increases steadily and transverse segregation increases until mixing by flow becomes steady. The axial segregation increases with solidification. But, with high intensity of flow axial segregation becomes steady after short transient. At small and large Gr the result showed a good agreememt with the prediction Smith and Scheil.

• Advances in aircraft and spacecraft science
• /
• v.7 no.3
• /
• pp.215-228
• /
• 2020

As is shown in the paper, the Koltunov-Rzhanitsyn singular kernel of heredity (when constructing mathematical models of the dynamics problem of the hereditary theory of viscoelasticity) adequately describes real mechanical processes, best approximates experimental data for a long period of time. A mathematical model of the problem of the flutter of viscoelastic plates moving in a gas with a high supersonic velocity is given. Using the Bubnov-Galerkin method, discrete models of the problem of the flatter of viscoelastic plates flowed over by supersonic gas flow are obtained. A numerical method is developed to solve nonlinear integro-differential equations (IDE) for the problem of the hereditary theory of viscoelasticity with weakly singular kernels. A general computational algorithm and a system of application programs have been developed, which allow one to investigate the nonlinear dynamic problems of the hereditary theory of viscoelasticity with weakly singular kernels. On the basis of the proposed numerical method and algorithm, nonlinear problems of the flutter of viscoelastic plates flowed over in a gas flow at an arbitrary angle are investigated. In a wide range of changes in various parameters of the plate, the critical velocity of the flutter is determined. It is shown that the singularity parameter α affects not only the oscillations of viscoelastic systems, but the critical velocity of the flutter as well.

• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.505-512
• /
• 2020

This paper presents an error estimation technique for 2-D crack analysis by an enriched natural element (more exactly, enriched Petrov-Galerkin NEM). A bare solution was approximated by PG-NEM using Laplace interpolation functions. Meanwhile, an accurate quasi-exact solution was obtained by a combined use of enriched PG-NEM and the global patch recovery. The Laplace interpolation functions are enriched with the near-tip singular fields, and the approximate solution obtained by enriched PG-NEM was enhanced by the global patch recovery. The quantitative error amount is measured in terms of the energy norm, and the accuracy (i.e., the effective index) of the proposed method was evaluated using the errors which obtained by FEM using a very fine mesh. The error distribution was investigated by calculating the local element-wise errors, from which it has been found that the relative high errors occurs in the vicinity of crack tip. The differences between the enriched and non-enriched PG-NEMs have been investigated from the effective index, the error distribution, and the convergence rate. From the comparison, it has been justified that the enriched PG-NEM provides much more accurate error information than the non-enriched PG-NEM.