• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• 한국수학교육학회지시리즈B:순수및응용수학
• /
• 제19권2호
• /
• pp.147-169
• /
• 2012

We introduce a Heaviside-function formulation of the interaction between incompressible two-phase fluid and a non-deformable solid. Fluid and solid interact in two ways : fluid satises the Dirichlet boundary condition imposed by the velocity field of solid, and solid is accelerated by the surface traction exerted by fluid. The two-way couplings are formulated by the Heaviside function to the interface between solid and fluid. The cumbersome treatment of interface is taken care of by the Heaviside function, and the interaction is discretized in a simple manner. The discretization results in a stable and accurate projection method.

• Structural Engineering and Mechanics
• /
• 제36권3호
• /
• pp.279-299
• /
• 2010

A methodology on application of the discrete singular convolution (DSC) technique to the free vibration analysis of thin plates with curvilinear quadrilateral platforms is developed. In the proposed approach, irregular physical domain is transformed into a rectangular domain by using geometric coordinate transformation. The DSC procedures are then applied to discretization of the transformed set of governing equations and boundary conditions. For demonstration of the accuracy and convergence of the method, some numerical examples are provided on plates with different geometry such as elliptic, trapezoidal having straight and parabolic sides, sectorial, annular sectorial, and plates with four curved edges. The results obtained by the DSC method are compared with those obtained by other numerical and analytical methods. The method is suitable for the problem considered due to its generality, simplicity, and potential for further development.

• Structural Engineering and Mechanics
• /
• 제27권6호
• /
• pp.683-696
• /
• 2007

This article is concerned with the presentation of a time-domain BEM approach applied to the solution of the scalar wave equation for 2D problems. The basic idea is quite simple: the basic variables of the problem at time $t_n$ (potential and flux) are computed with the results related to the potential and to its time derivative at time $t_{n-1}$ playing the role of "initial conditions". This time-marching scheme needs the computation of the potential and its time derivative at all boundary nodes and internal points, as well as the entire discretization of the domain. The convolution integrals of the standard time-domain BEM formulation, however, are not computed; the matrices assembled, only at the initial time interval, are those related to the potential, flux and to the potential time derivative. Two examples are presented and discussed at the end of the article, in order to verify the accuracy and potentialities of the proposed formulation.

• Steel and Composite Structures
• /
• 제24권5호
• /
• pp.523-536
• /
• 2017

A layerwise (LW) formulation based on the Galerkin method is presented to investigate the three-dimensional stress state in long sandwich plate which is subjected to tension force and pure bending moment. Based on the Galerkin method and the LW discretization approach, the equilibrium equations of elasticity for the long plate are written in the weak form and discretized through the thickness of the plate. The discretized equations are written in terms of displacement components of the numerical layers. The governing equations of the plate are solved analytically for the free edge boundary conditions. The distribution of stress state especially the 3D stress state in the vicinity of the edges of the sandwich plate which is subjected to tension and pure bending is studied. In order to increase the accuracy, the out of plane stresses are obtained by integrating the equilibrium equations of elasticity. The convergence and accuracy of the predictions are studied and various numerical results are presented for distribution of the in-plane and out of plane stresses in symmetric and un-symmetric sandwich plates.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2012년도 춘계학술대회 논문집
• /
• pp.748-754
• /
• 2012

In this paper free vibration analysis of symmetric and cross-ply elastic laminated shells based on FSDT was performed through discretization of equations of motion and boundary condition. Model of laminated composite shells subjected to a combination of magnetic and thermal fields is developed. These coupled equations of motion are based on the electromagnetic equations (Faraday, Ampere, Ohm, and Lorenz equations) and thermal equations which are involved in constitutive equations. Dynamic characteristic of composite shells for change of magnetic fields is investigated.

• Journal of applied mathematics & informatics
• /
• 제7권1호
• /
• pp.139-159
• /
• 2000

In this manuscript we study perturbed Newton-like methods for the solution of nonlinear operator equations in a Banach space and their discretized versions in connection with the mesh independence principle. This principle asserts that the behavior of the discretized process is asymptotically the same as that for the original iteration and consequently, the number of steps required by the two processes to converge to within a given tolerance is essentially the same. So far this result has been proved by others using Newton's method for certain classes of boundary value problems and even more generally by considering a Lipschitz uniform discretization. In some of our earlierpapers we extend these results to include Newton-like methods under more general conditions. However, all previous results assume that the iterates can be computed exactly. This is mot true in general. That in why we use perturbed Newton-like methods and even more general conditions. Our results, on the one hand, extend, and on the other hand, make more practical and applicable all previous results.

• Journal of applied mathematics & informatics
• /
• 제29권3_4호
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.

• 대한전기학회논문지
• /
• 제40권9호
• /
• pp.859-865
• /
• 1991

This paper presents a new shaoe optimal design method using Desing Sensitivity Analysis. Design Sensitivity, defined as the variation of the object function with respect to the design variables, is based on the total differentiation of the matrix equation obtained from discretization of governing equations with respect to design variables. Analysis of the various characteristics adn calculation of Design Sensitivity of optimization model are achieved by using finite or boundary element methods. The proposee algorithm is applied to the optimal shape design of high voltage electrode under specified conditions. It is shown, from the numerical results, that the algorithm is very useful for the optimal shape design of electric apparatus.

• Nuclear Engineering and Technology
• /
• 제32권1호
• /
• pp.46-56
• /
• 2000

Numerical results are presented for the 2-dimensional turbulent transient heat transfer of the shell/tube heat exchanger with a step change of the inlet temperature in the primary side. Heat transfer boundary conditions outside the pipe are given partially by the convection heat transfer conditions and partially by insulated conditions. Calculation results were obtained by solving the unsteady two-dimensional elliptic forms for the Reynolds-averaged governing equations for the mass, momentum and energy. Finite-difference method was used to obtain discretization equations, and the SIMPLER solution algorithm was employed for the calculation procedure. Turbulent model used is the algebraic model proposed by Cebeci-Smith. Results presented include the time variant Nusselt number distribution, average temperature distribution and outlet temperatures for the various inlet temperatures and flow rates.