• 제목/요약/키워드: Numerical Computation

검색결과 1,368건 처리시간 0.025초

접면포착법에 의한 수중익 주위의 이층류 유동계산 (Computation of Two-Fluid Flows with Submerged hydrofoil by Interface Capturing Method)

  • 곽승현
    • 한국항만학회지
    • /
    • 제13권1호
    • /
    • pp.167-174
    • /
    • 1999
  • Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

  • PDF

수직평판 주위를 흐르는 느린 점성류의 수치해석 (Computation of the Slow Viscous Flow about a Normal Plate)

  • 인기문;최도형;김문언
    • 대한기계학회논문집
    • /
    • 제17권9호
    • /
    • pp.2329-2338
    • /
    • 1993
  • An accurate analysis procedure to solve the flow about a flat plate at various incidences has been developed. The Navier-Stokes equations of stream function and vorticity form are solved in a sufficiently large computational domain, in which the grid lines are mutually orthogonal. The details of the flow near the singularity at the tip of the plate is well captured by the analytic solution which is asymptotically matched to the numerically generated outer solution. The solution for each region is obtained iteratively : the solution of one (inner or outer) region uses that of the other as the boundary condition after each cycle. The resulting procedure is accurate everywhere and also computationally efficient as the singularity has been removed. It is applied to the flat plate for a wide range of Re : the results agree very well with the existing computation and experiment.

기호계산 기법을 이용한 현가장치의 기구학적 민감도 해석 (Kinematic Design Sensitivity Analysis of Suspension System Using a Symbolic Computation Method)

  • 송성재;탁태오
    • 한국자동차공학회논문집
    • /
    • 제4권6호
    • /
    • pp.247-259
    • /
    • 1996
  • Kinematic design sensitivity analysis for vehicle in suspension systems design is performed. Suspension systems are modeled using composite joins to reduce the number of the constraint equations. This allows a semi-analytical approach that is computerized symbolic manipulation before numerical computations and that may compensate for their drawbacks. All the constraint equations including design variables are derived in symbolic equations for sensitivity analysis. By directly differentiating the equations with respect to design variables, sensitivity equations are obtained. Since the proposed method only requires the hard point data, sensitivity analysis is possible in suspension design stage.

  • PDF

COMPUTATION OF THE DYNAMIC FORCE COMPONENT ON A VERTICAL CYLINDER DUE TO SECOND ORDER WAVE DIFFRACTION

  • Bhatta, Dambaru
    • Journal of applied mathematics & informatics
    • /
    • 제26권1_2호
    • /
    • pp.45-60
    • /
    • 2008
  • Here we consider the evaluation of the the dynamic component of the second order force due to wave diffraction by a circular cylinder analytically and numerically. The cylinder is fixed, vertical, surface piercing in water of finite uniform depth. The formulation of the wave-structure interaction is based on the assumption of a homogeneous, ideal, incompressible, and inviscid fluid. The nonlinearity in the wave-structure interaction problem arises from the free surface boundary conditions, namely, dynamic and kinematic free surface boundary conditions. We expand the velocity potential and free surface elevation functions in terms of a small parameter and then consider the second order diffraction problem. After deriving the pressure using Bernoulli's equation, we obtain the analytical expression for the dynamic component of the second order force on the cylinder by integrating the pressure over the wetted surface. The computation of the dynamic force component requires only the first order velocity potential. Numerical results for the dynamic force component are presented.

  • PDF

다관절 보의 진동해석 알고리즘의 개발에 관한 연구 (On Development of Vibration Analysis Algorithm of Beam with Multi-Joints)

  • 문덕홍;여동준;최명수
    • Journal of Advanced Marine Engineering and Technology
    • /
    • 제18권5호
    • /
    • pp.68-77
    • /
    • 1994
  • The authors apply the transfer influence coefficient method to the 3-dimensional vibration analysis of beam with multi-joints and formulate a general algorithm to analysis the longitudinal, flexural and torsional coupled free vibration. In this paper, the structure, which is mainly founded in the robot arms, cranes and so on, has some crooked parts, subsystems and joints but has no closed loop in this system. It is modeled as the beam of a distributed mass system with massless translational, rotational and torsional springs in each node, and joint elements of release or roll at which node the displacement vector is discontinuous. The superiorty of the present method to the transfer matrix method in the computation accuracy was confirmed by the numerical computation results. Moreover, we confirmed that boundary and intermediate conditions could ve controlled by varying the values o the spring constants.

  • PDF

Reflections of shocks in nonequilibrium flow of air

  • Park, Tae-Hoon
    • 대한수학회논문집
    • /
    • 제10권3호
    • /
    • pp.767-781
    • /
    • 1995
  • In this paper we present computation of a reflected shock in the hypersonic flow of air with chemical reactions. We consider two dimensional steady inviscid hypersonic flow of air around bodies including chemical reaction effects. At a high Mach number, a strong shock is formed in front of the body when a wedge is placed against the flow. In front of the shock, temperature and pressure increase greatly and the flow is in nonequilibrium state. If the shock hits a wall, then a reflected shock is formed in the nonequilibrium flow region. Behind this reflected shock, the temperature and pressure are very high. We carry out the computation of the reflected shock and the flow behind it. The jump conditions at the reflected shock are presented. A technique combining smooth transforms of domain and implicit difference methods is used to overcome numerical difficulties associated with the lack of resolution behind the shock and near the body.

  • PDF

정적-외연적 강소성 유한요소법의 개발 및 Osakada방법에서 변형율 증분에 따른 영향분석 (Development of Static-explicit rigid-plastic finite Element Method and Investigate the offset of strain increment in Osakada method)

  • 정동원;이승훈
    • 한국정밀공학회지
    • /
    • 제21권2호
    • /
    • pp.116-121
    • /
    • 2004
  • In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study, static-explicit rigid-plastic finite element method will be introduced. This method is the way that restrict the convergence interval. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming process.

유압펌프용 실린더 블록의 윤활 및 동특성 해석 (Analysis of Lubrication and Dynamic Characteristics of a Cylinder Block for Hydraulic Pump)

  • 안성용;임윤철;홍예선
    • Tribology and Lubricants
    • /
    • 제20권4호
    • /
    • pp.209-217
    • /
    • 2004
  • Lubrication characteristics between a cylinder block and a valve plate for high speed bent-axis type hydraulic pump play an important role in volumetric efficiency and durability of pump. In this paper, a finite element method is presented for the computation of the pressure distribution between a cylinder block and a valve plate for high speed bent-axis type hydraulic pump. Also, a Runge-Kutta method is applied to simulate the cylinder block dynamics of three-degrees of freedom motion. From the results of computation, we can draw two major conclusions. One is related to the fluid film characteristics between a cylinder block and a valve plate and the other is related to the average leakage that is determined by the pressure gradient and the clearance near the discharge port. The numerical results of cylinder block dynamics were compared with the experimental results using eddy-current type gap sensors those are fixed at a pump housing.

VARIABLE TIME-STEPPING HYBRID FINITE DIFFERENCE METHODS FOR PRICING BINARY OPTIONS

  • Kim, Hong-Joong;Moon, Kyoung-Sook
    • 대한수학회보
    • /
    • 제48권2호
    • /
    • pp.413-426
    • /
    • 2011
  • Two types of new methods with variable time steps are proposed in order to valuate binary options efficiently. Type I changes adaptively the size of the time step at each time based on the magnitude of the local error, while Type II combines two uniform meshes. The new methods are hybrid finite difference methods, namely starting the computation with a fully implicit finite difference method for a few time steps for accuracy then performing a ${\theta}$-method during the rest of computation for efficiency. Numerical experiments for standard European vanilla, binary and American options show that both Type I and II variable time step methods are much more efficient than the fully implicit method or hybrid methods with uniform time steps.

SPECTRAL-COLLOCATION METHOD FOR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS

  • Yang, Yin;Chen, Yanping;Huang, Yunqing
    • 대한수학회지
    • /
    • 제51권1호
    • /
    • pp.203-224
    • /
    • 2014
  • We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Fredholm-Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in $L^{\infty}$ norm and weighted $L^2$-norm. The numerical examples are given to illustrate the theoretical results.