• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.5
• /
• pp.47-60
• /
• 2007

In this paper, we investigate the natural frequencies and buckling loads of functionally graded material (FGM) plates and shells, using a quasi-conforming shell element that accounts for the transverse shear strains and rotary inertia. The eigenvalue of the FGM plates and shells are calculated by varying the volume fraction of the ceramic and metallic constituents using a sigmoid function, but their Poisson's ratios of the FGM plates and shells are assumed to be constant. The expressions of the membrane, bending and shear stiffness of FGM shell element are more complicated combination of material properties than a homogeneous element. In order to validate the finite element numerical solutions, the Navier's solutions of rectangular plates based on the first-order shear deformation theory are presented. The present numerical solutions of composite and sigmoid FGM (S-FGM) plates are proved by the Navier's solutionsand various examples of composite and FGM structures are presented. The present results are in good agreement with the Navier's theoretical solutions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.675-691
• /
• 2002

The focus of this work is on the development of large-scale numerical optimization methods for optimal control of steady incompressible Navier-Stokes flows. The control is affected by the suction or injection of fluid on portions of the boundary, and the objective function of fluid on portions of the boundary, and the objective function represents the rate at which energy is dissipated in the fluid. We develop reduced Hessian sequential quadratic programming. Both quasi-Newton and Newton variants are developed and compared to the approach of eliminating the flow equations and variables, which is effectively the generalized reduced gradient method. Optimal control problems we solved for two-dimensional flow around a cylinder. The examples demonstrate at least an order-of-magnitude reduction in time taken, allowing the optimal solution of flow control problems in as little as half an hour on a desktop workstation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.661-674
• /
• 2002

The objective of this study is to develop efficient numerical method to enable solution of optimal control problems of Navier-Stokes flows and to apply these technique to the problem of viscous drag minimization on a bluff body by controlling boundary velocities on the surface of the body. In addition to the industrial importance of the drag reduction problem, it serves as a model for other more complex flow optimization settings, and allows us to study, modify, and improve the behavior of the optimal control methods proposed here. The control is affected by the suction or injection of fluid on portions of the boundary, and the objective function represents the rate at which energy is dissipated in the fluid. This study shows how reduced Hessian successive quadratic programming method, which avoid converging the flow equations at each iteration, can be tailored to these problems.

• Journal of the Korean Chemical Society
• /
• v.40 no.6
• /
• pp.409-414
• /
• 1996

The limit of applicability of Navier-Stokes equation to stationary plane shock-waves is examined by using the principle of minimum entropy production of linear irreversible thermodynamics. In order to obtain analytic results, the equation is linearized near the equilibrium of downstream. Results show that the solution of Navier-Stokes equation which fits the boundary condition of far downstream flow is consistent with the thermodynamic requirement within the first order when the solution is expanded around the M=1, where M is the Mach number of upstream speed.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.63.1-63.1
• /
• 2010

수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two Step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험 결과와 비교 및 분석하였다.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.22 no.5
• /
• pp.13-22
• /
• 2018

In this study, a governing differential equation for the bending problem of orthotropic rectangular plate is drived. It's exact solution for various boundary conditions is presented. This solution follows traditional method like Navier's solution or Levy's solution that transforms the governing differential equation into an algebraic equation by using trigonometric series. To obtain a solution by Levy's method, it is required that two opposite edges of the plate be simply supported. And the boundary conditions, for which the Navier's method is applicable, are simply supported edge at all edges. In this study, it overcomes the limitations of the previous Navier's and Levy's methods.This solution is applicable for any combination of boundary conditions with simply supported edge and clamped edge in x, y direction. The plate could be subjected to uniform, partially uniform, and line loads. The advantage of the solution is that it is the exact solution as well as it overcomes the limitations of the previous Navier's and Levy's methods. Calculations are presented for orthotropic plates with nonsymmetric boundary conditions. Comparisons between the result of this paper and the result of Navier, Levy and Szilard solutions are made for the isotropic plates. The deflections were in excellent agreement.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.111-117
• /
• 2010

A three-dimensional numerical model is employed to investigate wave deformation due to a submerged structure. The three-dimensional numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES(large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS(sub-grid scale) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF(volume-of-fluid) method is used to track the distorted and broken free surface. A simple linear wave is generated on a constant depth and compared with analytical solutions. The model is then applied to study wave deformation due to a submerged structure and the predicted results are compared with available laboratory measurements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1581-1593
• /
• 2001

An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

• Tunnel and Underground Space
• /
• v.13 no.5
• /
• pp.389-396
• /
• 2003

It is shown that the cubic law can be modified regarding the steady-state Navier-Stokes equations by using perturbation approximation method for a sinusoidal aperture variation. In order to adopt the perturbation theory, the sinusoidal function needs to be non-dimensionalized for the amplitude and wavelength. Then, the steady-state Navier-Stokes equations can be solved by expanding the non-dimensionalized stream function with respect to the small value of the parameter (the ratio of the mean aperture to the wavelength), together with the continuity equation. From the approximate solution of the Navier-Stokes equations, the basic cubic law is successfully modified for the steady-state condition and a sinusoidal aperture variation. A finite difference method is adopted to calculate the pressure within a fracture model, and the results of numerical experiments show the accuracy and applicability of the modified cubic law. As a result, it is noted that the modified cubic law, suggested in this study, will be used for the analysis of fluid flow through aperture geometry of sinusoidal distributions.