• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.114-125
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.359-366
• /
• 2008

Current state and perspective of DNS of turbulence and turbulent combustion are discussed with feature trend of the fastest supercomputer in the world. Based on the perspective of DNS of turbulent combustion, possibility of perfect simulations of IC engine is shown. In 2020, the perfect simulation will be realized with 30 billion grid points by 1EXAFlops supercomputer, which requires 4 months CPU time. The CPU time will be reduced to about 4 days if several developments were achieved in the current fundamental researches. To shorten CPU time required for DNS of turbulent combustion, two numerical methods are introduced to full-explicit full-compressible DNS code. One is compact finite difference filter to reduce spatial resolution requirements and numerical oscillations in small scales, and another is well-known point-implicit scheme to avoid quite small time integration of the order of nanosecond for fully explicit DNS. Availability and accuracy of these numerical methods have been confirmed carefully for auto-ignition, planar laminar flame and turbulent premixed flames. To realize DNS of IC engine with realistic kinetic mechanism, several DNS of elemental combustion process in IC engines has been conducted.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.4
• /
• pp.651-661
• /
• 2017

A truss-braced wing (TBW) aircraft has recently received increasing attention due to higher aerodynamic efficiency compared to conventional cantilever wing aircraft. For conceptual TBW aircraft design, we developed a propulsion-and-airframe integrated design environment by replacing a semi-empirical turbofan engine model with a thermodynamic cycle-based one built upon the numerical propulsion system simulation (NPSS). The constructed NPSS model benefitted TBW aircraft design study, as it could handle engine installation effects influencing engine fuel efficiency. The NPSS model also contributed to broadening TBW aircraft design space, for it provided turbofan engine design variables involving a technology factor reflecting progress in propulsion technology. To effectively consolidate the NPSS propulsion model with the TBW airframe model, we devised a rapid, approximate substitute of the NPSS model by reduced-order modeling (ROM) to resolve difficulties in model integration. In addition, we formed an artificial neural network (ANN) that associates engine component attributes evaluated by object-oriented weight analysis of turbine engine (WATE++) with engine design variables to determine engine weight and size, both of which bring together the propulsion and airframe system models. Through propulsion-andairframe design space exploration, we optimized TBW aircraft design for fuel saving and revealed that a simple engine model neglecting engine installation effects may overestimate TBW aircraft performance.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.4
• /
• pp.582-591
• /
• 2001

It is known that the reduced integration, modified shape function, anisoparametric and non-conforming element can reduce the error induced by stiffness locking phenomenon in the finite element analysis. In this study, we propose new anisoparametric curved beam element. The new element based on reduced minimization theory is composed of different shape functions in each displacement field. By the substitution of this modified shape function, the unmatched coefficient that cause stiffness locking in the constraint energy is eliminated. To confirm the availability of this new model, we performed numerical tests for a simple model. As a result of numerical test, the undulate stress patterns are disappeared in static analysis, and displacements and stresses are close to exact solution. Not only in the static analysis but also in the eigen analysis of free vibrated curved beam model, this element shows successful convergent results.

• Journal of Institute of Control, Robotics and Systems
• /
• v.14 no.12
• /
• pp.1270-1277
• /
• 2008

A Haar wavelet based numerical method for solving singularly perturbed linear time invariant system is presented in this paper. The reduced pure slow and pure fast subsystems are obtained by decoupling the singularly perturbed system and differential matrix equations are converted into algebraic Sylvester matrix equations via Haar wavelet technique. The operational matrix of integration and its inverse matrix are utilized to reduce the computational time to the solution of algebraic matrix equations. Finally a numerical example is given to demonstrate the validity and applicability of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary element.

• Structural Engineering and Mechanics
• /
• v.57 no.4
• /
• pp.733-762
• /
• 2016

Based on a reduced displacement field, a layer-wise (LW) formulation is developed for analysis of thick shell panels which is subjected to axial tension. Employing the principle of minimum total potential energy, the local governing equations of thick panel which is subjected to axial extension are obtained. An analytical method is developed for solution of the governing equations for various edge conditions. The governing equations are solved for free and simply supported edge conditions. The interlaminar stresses in the panel are investigated by means of Hooke's law and also by means of integration of the equilibrium equations of elasticity. Dependency of the result upon the number of numerical layers in the layerwise theory (LWT) is studied. The accuracy of the numerical results is validated by comparison with the results of the finite element method and with other available results in the open literature and good agreement is seen between the results. Numerical results are then presented for the distribution of interlaminar normal and shear stresses within the symmetric and un-symmetric cross-ply thick panels with free and simply supported boundaries. The effects of the geometrical parameters such as radius to thickness and width to thickness ratio are investigated on the distribution of the interlaminar stresses in thick panels.

• Structural Engineering and Mechanics
• /
• v.38 no.4
• /
• pp.429-441
• /
• 2011

In large-scale problem, a huge size of computational resources is needed for a reliable solution which represents the detailed description of dynamic behavior. Recently, eigenvalue reduction schemes have been considered as important technique to resolve computational resource problems. In addition, the efforts to advance an efficiency of reduction scheme leads to the development of the multi-level system condensation (MLSC) which is initially based on the two-level condensation scheme (TLCS). This scheme was proposed for approximating the lower eigenmodes which represent the global behavior of the structures through the element-level energy estimation. The MLSC combines the multi-level sub-structuring scheme with the previous TLCS for enhancement of efficiency which is related to computer memory and computing time. The present study focuses on the implementation of the MLSC on the direct time response analysis and the frequency response analysis of structural dynamic problems. For the transient time response analysis, the MLSC is combined with the Newmark's time integration scheme. Numerical examples demonstrate the efficiency of the proposed method.