• KSTLE International Journal
• /
• 제5권1호
• /
• pp.23-27
• /
• 2004

The object of the present work is the numerical analysis of the computer hard disk slider. The pressure between slider and disk surfaces is calculated using the Boundary Fitted Coordinate System and Divergence Formulation for the nonlinear Reynolds' equation solution. The optimization scheme is applied to search for the steady state position of the slider. The simplified method is given for the case of the fixed inclined pad. The film thickness ratios and pitching and rolling angles are considered as alternative choice of the slider's coordinates. The behavior of the objective function for the Negative Pressure slider is studied in details. Methods of conjugate directions and feasible directions are applied.

• Tribology and Lubricants
• /
• 제23권5호
• /
• pp.240-247
• /
• 2007

The design process of the dental handpiece is described. The parameters of the high speed air turbine are estimated. The effect of supply hole on the stiffuess and damping of the air bearing for handpiece is studied numerically. The Reynolds equation is solved by using the divergence formulation and the perturbation method. The test rig is built and the test procedure is developed for the turbine rotational speed measurement by using Fourier transform of noise generated by the turbine during steady operation.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.550-559
• /
• 2004

Theoretical studies have been carried out to examine the influence of the grain geometry-dependent driving forces, which control the internal flow pattern of solid rockets. Numerical studies have been executed with the help of a two-dimensional code. This code solves standard k-omega turbulence equations using the coupled second order implicit unsteady formulation. It has been concluded that the grain port divergence angles have significant leverage on the formation of recirculation bubbles leading for pressure oscillations, flow separation and reattachment. In solid rockets flow reattachment will favour secondary ignition and that will add to the complexity of the starting transient prediction.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Structural Engineering and Mechanics
• /
• 제48권6호
• /
• pp.879-914
• /
• 2013

In part I of the article, formulation and characteristics of the several well-known structural geometrical nonlinear solution techniques were studied. In the present paper, the efficiencies and capabilities of residual load minimization, normal plane, updated normal plane, cylindrical arc length, work control, residual displacement minimization, generalized displacement control and modified normal flow will be evaluated. To achieve this goal, a comprehensive comparison of these solution methods will be performed. Due to limit page of the article, only the findings of 17 numerical problems, including 2-D and 3-D trusses, 2-D and 3-D frames, and shells, will be presented. Performance of the solution strategies will be considered by doing more than 12500 nonlinear analyses, and conclusions will be drawn based on the outcomes. Most of the mentioned structures have complex nonlinear behavior, including load limit and snap-back points. In this investigation, criteria like number of diverged and complete analyses, the ability of passing load limit and snap-back points, the total number of steps and analysis iterations, the analysis running time and divergence points will be examined. Numerical properties of each problem, like, maximum allowed iteration, divergence tolerance, maximum and minimum size of the load factor, load increment changes and the target point will be selected in such a way that comparison result to be highly reliable. Following this, capabilities and deficiencies of each solution technique will be surveyed in comparison with the other ones, and superior solution schemes will be introduced.

• 한국전산응용수학회:학술대회논문집
• /
• 한국전산응용수학회 2003년도 KSCAM 학술발표회 프로그램 및 초록집
• /
• pp.9-9
• /
• 2003

Let $\Omega$ be a bounded, open, and polygonal domain in $R^2$ with re-entrant corners. We consider the following Partial Differential Equations: $$(I-\nabla\nabla\cdot+\nabla^{\bot}\nabla\times)u\;=\;f\;in\;\Omega$$, $$n\cdotu\;0\;0\;on\;{\Gamma}_{N}$$, $${\nabla}{\times}u\;=\;0\;on\;{\Gamma}_{N}$$, $$\tau{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$, $$\nabla{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$ where the symbol $\nabla\cdot$ and $\nabla$ stand for the divergence and gradient operators, respectively; $f{\in}L^2(\Omega)^2$ is a given vector function, $\partial\Omega=\Gamma_{D}\cup\Gamma_{N}$ is the partition of the boundary of $\Omega$; nis the outward unit vector normal to the boundary and $\tau$represents the unit vector tangent to the boundary oriented counterclockwise. For simplicity, assume that both $\Gamma_{D}$ and $\Gamma_{N}$ are nonempty. Denote the curl operator in $R^2$ by $$\nabla\times\;=\;(-{\partial}_2,{\partial}_1$$ and its formal adjoint by $${\nabla}^{\bot}\;=\;({-{\partial}_1}^{{\partial}_2}$$ Consider a weak formulation(WF): Find $u\;\in\;V$ such that $$a(u,v):=(u,v)+(\nabla{\cdot}u,\nabla{\cdot}v)+(\nabla{\times}u,\nabla{\times}V)=(f,v),\;A\;v{\in}V$$. (2) We assume there is only one singular corner. There are many methods to deal with the domain singularities. We introduce them shortly and we suggest a new Finite Element Methods by using Singular representation for the solution.

• 대한기계학회논문집
• /
• 제16권10호
• /
• pp.1971-1978
• /
• 1992

본 연구에서는 원형 파이프에 원추형 디퓨져가 연결된 덕트 내의 유동장에 대 하여 Launder-Sharma의 저 레이놀즈수 k-.epsilon. 난류모델을 이용하여 수치해석을 수행하였 으며, 수치해석 방법으로는 타원형 방법을 사용하였으며, 앞으로 일반적인 단면의 곡 관이나 스크롤 내부 유동 등의 연구 수행을 감안하여 지배방정식을 일반 비직교 좌표 계로 변환하여 계산을 수행하였다.

• 대한기계학회논문집A
• /
• 제36권10호
• /
• pp.1139-1146
• /
• 2012

P1 비순응 요소를 이용하여 정상 비압축성 Navier-Stokes 유동의 위상최적화 문제를 푸는 방법을 제시한다. 본 연구는 Stokes 유동의 위상최적화 문제에 P1 비순응 요소를 적용하여 그 수치적 효용성을 보인바 있는 이전 연구에 대한 후속 연구이다. 비압축성 물질 해석에서 잠김현상이 발생하지 않으며 선형형상함수를 가지는 P1 비순응 요소의 장점이 관성항을 가지는 유체 문제의 해석과 설계에도 유효한 지를 파악하고자 한다. 일반적으로 사용되는 혼합정식화법과 비교하여 P1 비순응 요소의 사용은 벌칙 함수를 이용하여 연속 방정식을 따로 사용하지 않고 운동방정식에 부과할 수 있기 때문에 자유도의 개수를 감소시킬 수 있다. 벌칙 파라미터가 해의 정확도에 주는 영향과 적정 범위는 수치적으로 검토하도록 한다. 또한 보통의 사각 비순응 요소들이 요소면의 중앙에 절점을 가지고 고차의 형상함수를 지니는데 비하여, 본 연구에서 제시하는 P1 비순응 요소는 요소의 꼭지점에 절점을 가지고 {1, x, y}의 P1 형상함수로 구성됨으로써 수치적인 구현의 용이함이 일반 선형 사절점 요소와 동일하다. 제안한 방법의 효용성을 다양한 레이놀즈수에 따른 유동최적화 문제들을 살펴봄으로써 검증하도록 한다.

• Journal of Mechanical Science and Technology
• /
• 제20권12호
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• 천문학회보
• /
• 제44권2호
• /
• pp.69.2-69.2
• /
• 2019

To construct a coronal force-free magnetic field, we must impose the boundary normal current density (or three components of magnetic field) as well as the boundary normal field at the photosphere as boundary conditions. The only method that is known to implement these boundary conditions exactly is the method devised by Grad and Rubin (1958). However, the Grad-Rubin method and all its variations (including the fluxon method) suffer from convergence problems. The magnetofrictional method and its variations are more robust than the Grad-Rubin method in that they at least produce a certain solution irrespective of whether the global solution is compatible with the imposed boundary conditions. More than often, the influence of the boundary conditions does not reach beyond one or two grid planes next to the boundary. We have found that the 2D solenoidal gauge condition for vector potentials allows us to implement the required boundary conditions easily and effectively. The 2D solenoidal condition is translated into one scalar function. Thus, we need two scalar functions to describe the magnetic field. This description is quite similar to the Chandrasekhar-Kendall representation, but there is a significant difference between them. In the latter, the toroidal field has both Laplacian and divergence terms while in ours, it has only a 2D Laplacian term. The toroidal current density is also expressed by a 2D Laplacian. Thus, the implementation of boundary normal field and current are straightforward and their effect can permeate through the whole computational domain. In this paper, we will give detailed math involved in this formulation and discuss possible lateral and top boundary conditions and their meanings.