• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.311-317
• /
• 2011

The present paper deals with the continuous works of extending the multi-dimensional limiting process (MLP) for compressible flows, which has been quite successful in finite volume methods, into discontinuous Galerkin (DG) methods. From the series of the previous, it was observed that the MLP shows several superior characteristics, such as an efficient controlling of multi-dimensional oscillations and accurate capturing of both discontinuous and continuous flow features. Mathematically, fundamental mechanism of oscillation-control in multiple dimensions has been established by satisfaction of the maximum principle. The MLP limiting strategy is extended into DG framework, which takes advantage of higher-order reconstruction within compact stencil, to capture detailed flow structures very accurately. At the present, it is observed that the proposed approach yields outstanding performances in resolving non-compressive as well as compressive flaw features. In the presentation, further numerical analyses and results are going to be presented to validate that the newly developed DG-MLP methods provide quite desirable performances in controlling numerical oscillations as well as capturing key flow features.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

• Composites Research
• /
• v.26 no.4
• /
• pp.235-244
• /
• 2013

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solid containing multiple isotropic elliptical fibers of arbitrary orientation subject to uniform stress at infinity. The fibers are assumed to be long parallel elliptical cylinders composed of isotropic elastic material perfectly bonded to the isotropic matrix. The solid is assumed to be under plane strain on the plane normal to the cylinders. A detailed analysis of the stress field at the matrix-fiber interface for square and hexagonal packing of the fibers is carried out for different values of the number, orientation angles and concentration of the elliptical fibers. The accuracy and efficiency of the method are examined through comparison with results obtained from analytical and finite element methods.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.2
• /
• pp.15-20
• /
• 2011

This study aimed at validating the adopted numerical methods to solve two-phase flow around a two-dimensional (2D) rectangular floating structure in regular waves. A structure with a draft equal to one half of its height was hinged at the center of gravity and free to roll with waves that had the same period as the natural roll period of a rectangular barge. In order to simulate the 2D incompressible viscous two-phase flow in a wave tank with the rectangular barge, the present study used the volume of fluid (VOF) method based on the finite volume method with a standard turbulence model. In addition, the sliding mesh technique was used to handle the motion of the rectangular barge induced by the fluid-structure interaction. Consequently, the present results for the flow field and roll motion of the structure had good agreement with those of the relevant previous experiment.

• Structural Engineering and Mechanics
• /
• v.50 no.1
• /
• pp.53-71
• /
• 2014

A semi-analytical method is developed to consider free vibrations of a functionally graded elastic plate resting on Winkler elastic foundation and in contact with a quiescent fluid. Material properties are assumed to be graded distribution along the thickness direction according to a power-law in terms of the volume fractions of the constituents. The fluid is considered to be incompressible and inviscid. In the analysis, the effect of an in-plane force in the plate due to the weight of the fluid is taken into account. By satisfying the compatibility conditions along the interface of fluid and plate, the fluid-structure interaction is taken into account and natural frequencies and mode shapes of the coupled system are acquired by employing energy methods. The results obtained from the present approach are verified by those from a finite element analysis. Besides, the effects of volume fractions of functionally graded materials, Winkler foundation stiffness and in-plane forces on the dynamic of plate are elucidated.

• Steel and Composite Structures
• /
• v.42 no.3
• /
• pp.375-388
• /
• 2022

In this work, limit elastic speed analysis of functionally graded porous rotating disks has been reported. The work proposes an effective approach for modeling the mechanical properties of a porous functionally graded rotating disk. Four different types of porosity models namely: uniform, symmetric, inner maximum, and outer maximum distribution are considered. The approach used is the variational principle, and the solution has been achieved using Galerkin's error minimization theory. The study aims to investigate the effect of grading indices, aspect ratio, porosity volume fraction, and porosity types on limit angular speed for uniform and variable disk geometries of constant mass. To validate the current study, finite element analysis has been used, and there is good agreement between the two methods. The study yielded a decrease in limit speed as grading indices and aspect ratio increase. The porosity volume fraction is found to be more significant than the aspect ratio effect. The research demonstrates a range of operable speeds for porous and non-porous disk profiles that can be used in industries as design data. The results show a significant increase in limit speed for an exponential disk when compared to other disk profiles, and thus, the study demonstrates a range of FG-based structures for applications in industries that will not only save material (lightweight structures) but also improve overall performance.

• Transactions of Materials Processing
• /
• v.16 no.1 s.91
• /
• pp.25-30
• /
• 2007

The present study covers an integrated simulation method to evaluate optical performance of an aspheric plastic lens by connecting a finite element (FE) analysis of injection molding with a ray tracing simulation. Traditional ray tracing methods have based on the assumption that the optical properties of a lens are homogeneous throughout the entire volume. This assumption is to a certain extent unrealistic for injection-molded plastic lenses because material properties vary at every point due to the injection molding effects. To take into account the effects of the inhomogeneous optical properties of the molded lens, a new.ay tracing scheme is proposed in conjunction with a FE analysis of the injection molding. A numerical scheme is developed to calculate ray paths on every element layer with more realistic information of the refractive indices which can be obtained through the FE analysis. This information is then used to calculate the ray paths based on the FE mesh of which nodal points have unique index values. The proposed tracing scheme is implemented on the tracing of an aspheric lens, and its validity is ascertained through experimental verification.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.29-34
• /
• 1999

To develop the efficient numerical optimization method for the design of an airfoil, an evaluation of various methods coupled with two-dimensional Naviev-Stokes analysis is presented. Simplex method and Hook-Jeeves method we used as direct search methods, and steepest descent method, conjugate gradient method and DFP method are used as indirect search methods and are tested to determine the search direction. To determine the moving distance, the golden section method and cubic interpolation method are tested. The finite volume method is used to discretize two-dimensional Navier-Stokes equations, and SIMPLEC algorithm is used for a velocity-pressure correction method. For the optimal design of two-dimensional airfoil, maximum thickness, maximum ordinate of camber line and chordwise position of maximum ordinate are chosen as design variables, and the ratio of drag coefficient to lift coefficient is selected as an objective function. From the results, it is found that conjugate gradient method and cubic interpolation method are the most efficient for the determination of search direction and the moving distance, respectively.