• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.36-41
• /
• 1999

Shock wave propagating in the particle suspension has important applications. Examples are shock waves occurring in the solid rocket plume and detonation of dusty particles by shock waves. Experimental and numerical investigations on this subject have drawn much attention. More recently, Sivier et al. numerically simulated the experiment of Sommerfeld using the unstructured adaptive grid. They used the Eulerian-Eulerian approach based on the continuum assumption for both gas and particles. In the present paper, a new numerical method using the Lagrangian particle tracing technique and unstructured particle-adaptive grid for the polydisperse system is presented. It is explained why the existing numerical calculation has showed discrepancy with the experimental results by Sommerfeld.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.137-142
• /
• 1996

A new wall impaction model for spray and its assessment are described in this paper. The gas phase is modelled in terms of the Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled following the discrete droplet model approach. The droplet parcel contains many thousands of drops assumed to have the same size, temperature and velocity components. The droplet parcel equations of trajectory, momentum, mass and energy are written in Lagrangian form. The new drop-wall interaction model is proposed, which is based on experimental investigations on individual drops, and it is applied for the general non-orthogonal gird. The model is then assessed through comparison with experiments over a wide range of test conditions of sprays. The results are in good agreement with experimental data.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.7 s.238
• /
• pp.990-997
• /
• 2005

The objective of this paper is to present the element connectivity parameterization (ECP) fur three dimensional problems. In the ECP method, a continuum structure is viewed as discretized finite elements connected by zero-length elastic links whose stiffness values control the degree of inter-element connectivity. The ECP method can effectively avoid the formation of the low-density unstable elements. These elements appear when the standard element density method is used for geometrical nonlinear problems. In this paper, this ECP method developed fur two-dimensional problems is expanded to the design of three-dimensional geometrical nonlinear structures. Among others, the automatic procedure converting standard finite element models to the models suitable for the ECP approach is developed and applied for optimization problems defined on general three-dimensional design domains.

• Journal of Advanced Marine Engineering and Technology
• /
• v.23 no.1
• /
• pp.47-53
• /
• 1999

Small compression-ignition direct injection engines have been developed as a measure to improve a fuel efficiency and reduce harmful exhaust gases. Those small engines generally employ high injection pressure increase on the spray impacting on a wall is discussed in this paper. The gas phase is modelled by the Eulerian continuum conservation equations of mass momentum energy and fuel vapour fraction. The liquid phases is modelled following the discrete droplet model approach in Lagrangian form and the droplet wall interaction is modelled as a func-tion of the velocity normal to impaction lands. The droplet distributions vapor fractions and gas flows are analyzed in various injection pres-sure cases. The penetrations of wall spray and vapor increase and the Sauter mean diameter decreases with increasing injection pressure.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.17 no.3
• /
• pp.101-107
• /
• 2008

A composite mechanics for discontinuous fiber reinforced polymer matrix composites(PMC) is analysed in order to predict fiber axial stresses. In continuum approach. frictional slip which usually takes place between fibers and polymers is accounted to derive PMC equations. The interfacial friction stress is treated by the product of the coefficient of friction and the compressive stress norma1 to the fiber/matrix interface. The residual stress and the Poisson's contraction implemented by the rule of mixture(ROM) are considered for the compressive stress normal to the fiber/matrix interface. In addition. the effects of fiber aspect ratio and fiber volume fraction on fiber axial stresses are evaluated using the derived equations. Results are illustrated numerically using the present equations with reasonable materials data. It is found that the fiber axial stress in the center region shows no great discrepancy for different fiber aspect ratios and fiber volume fractions while some discrepancies are shown in the fiber end region.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.9 no.4
• /
• pp.18-26
• /
• 2001

OSKA engine was developed to remove the dense core of injection sprays. The engine uses impinging spray on a small pip, which spray after impinging is broken into smaller drops and disperses into fee space in chamber. In this paper the pip size is analyzed to give more dispersion of spray and fuel vapor. The gas phase is modelled by the Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled following the discrete droplet model approach in Lagrangian form, and the droplet wall interaction is modelled as a function of the velocity normal to impaction lands. The droplet distributions, vapor fractions and gas flows are analyzed for various injection pressure cases. Numerical results indicate that the land diameter of 5.6mm has the best performance of spray dynamics and vaporization in the test sizes.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.8
• /
• pp.224-232
• /
• 1999

Topology optimization has evolved into a very efficient concept design tool and has been incorporated into design engineering processes in many industrial sectors. In recent years, topology optimization has become the focus of structural design community and has been researched and applied widely both in academia and industry. There are mainly tow approaches for topology optimization of continuum structures ; homogenization and density methods. The homogenization method is to compute is to compute an optimal distribution of microstructures in a given design domain. The sizes of the micro-calvities are treated as design variables for the topology optimization problem. the density method is to compute an optimal distribution of an isotropic material, where the material densities are treated as design variables. In this paper, the density method is used to formulate the topology optimization problem. This optimization problem is solved by using an optimality criteria method. Several example problems are solved to show the usefulness of the present approach.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.3
• /
• pp.215-221
• /
• 2013

The basic quantitative tool for the study of hypersonic rarefied flows is the direct simulation Monte Carlo method (DSMC). The DSMC method requires a large amount of computer memory and performance and is unreasonably expensive at the first stage of spacecraft design and trajectory analysis. A possible solution to this problem is approximate engineering methods. However, the Monte Carlo method remains the most reliable approach to compare to the engineering methods that provide good results for the global aerodynamic coefficients of various geometry designs. This paper presents the calculation results of aerodynamic characteristics for spacecraft vehicles in the free molecular, the transitional and the continuum regimes using the local engineering method. Results and methods would be useful to calculate aerodynamics for new-generation hypersonic vehicle designs.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.10a
• /
• pp.393-400
• /
• 2001

A simple numerical modelling technique is proposed for estimating the shear stress distribution in beams of framed tube structures with multiple internal tubes. The structures are analysed using a continuum approach in which each tube is individually modelled by a tube beam that accounts for the flexural and shear deformations, as well as the shear lag effects. The method idealises the discrete tubes-in-tube structures as an assemblage of equivalent multiple beams, each composed of orthotropic plate panels. The numerical analysis of shear stress is based on the elastic theory in conjunction with the minimum potential energy principle. By simplifying assumptions regarding the form of strain distributions in external and internal tubes, the shear stress distributions are expressed in terms of a series of linear functions of the second moments of area of the structures and the corresponding geometric and material properties, as well as the applied loads. The simplicity and accuracy of the proposed method are demonstrated through the solutions of three numerical examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.10a
• /
• pp.45-52
• /
• 2001

A series of unreinforced masonry (URM) walls were analytically investigated for a limited version of seismic in-plane performance. For this URM walls were assumed to be an elastic continuum and modeled as isotropic plane stress elements within which the nature of cracking was propagated. Accordingly, cracking mode of behavior in URM was modeled by smeared-crack approach. Total of 70 cases were considered for various parameters such as axial load ratio, aspect ratio and effective section area ratio due to the existence of opening, etc. The analysis results indicated a general tendency in base shear coefficient and deformability of URM walls for these variables.