• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.113-123
• /
• 1990

The paper is concerned with the elasto-plastic and geometrically nonlinear analysis of shell structures using an improved degenerated shell element. In the formulation of the element stiffness, the combined use of three different techniques was made. They are; 1) an enhanced interpolation of transverse shear strains in the natural coordinate system to overcome the shear locking problem ; 2) the reduced integration technique in in-plane strains to avoid the membrane locking behavior ; and 3) selective addition of the nonconforming displacement modes to improve the element performances. This element is free of serious shear/membrane locking problems and undesirable compatible/commutable spurious kinematic deformation modes. In the formulation for plastic deformation, the concept of a layered element model is used and the material is assumed von Mises yield criterion. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements and rotations. The resulting non-linear equilibrium equations are solved by the Netwon-Raphson method combined with load or displacement increment. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Journal of Astronomy and Space Sciences
• /
• v.4 no.1
• /
• pp.11-23
• /
• 1987

In dward nova we assume the primary star as a white dwarf and the secondary as the late type star which filled Roche lobe. Mass flow from the secondary star leads to the formation of thin accretion disk around the white dwarf. We use the $\alpha$ parameter as viscosity to maintain the disk form and propose that the outburst in dwarf nova cause the step increase of source term. With these assumptions we solve the basic equations of stellar structure using New-Raphson method. We show the physical parameters like temperature, density, pressure, opacity, surface density, height and flux to the radius of disk. Changing the value of $\alpha$, we compare several parameters when mass flow rate is constant with those of when luminosity of disk is brightest. At the same time, we obtain time-dependent variations of luminosity and mass of disk. We propose the suitable range of $\alpha$is 0.15-0.18 to the difference of luminosity. We compare several parameters of disk with those of the normal late type stars which have the same molecular weight of disk. These show the temperature and pressure of disk are similar to those of normal stars but the density of disk is lower. Maybe the outburst in dwarf nova is due to the variation of the $\alpha$ value instead of increment of mass flow from the secondary star.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.1-15
• /
• 1986

The algorithm Proposed utilizes the two-levels technique. In the first level which consists of two phases, the cross-sectional area of the truss member is optimized by transforming the nonlinear problem into SUMT, and solving SUMT utilizing the modified Newton-Rahson method. In the second level, the geometric shape is optimized utilizing the unindirectional search technique of the Powell method which make it possible to minimize only the objective function. The algorithm Proposed in this study is numerically tested for several truss structures with various shapes, loading conditions and design criteria, and compared with the results of the other algorithms to examine its applicability and stability. The numerical comparisons show that the two-Levels algorithm Proposed in this study is safely applicable to any design criteria, and the convergency rate is relathely fast and stable compared with other iteration methods for the geometric optimization of truss structures.

• Journal of Mechanical Science and Technology
• /
• v.17 no.6
• /
• pp.844-853
• /
• 2003

In this paper, a dynamic analysis of the reciprocating compression mechanism of a small refrigeration compressor is performed. In the problem formulation of the mechanism dynamics, the viscous frictional force between the piston and the cylinder wall is considered in order to determine the coupled dynamic behaviors of the piston and the crankshaft. Simultaneous solutions are obtained for the equations of motion of the reciprocating mechanism and the time-dependent Reynolds equations for the lubricating film between the piston and the cylinder wall and for the oil films on the journal bearings. The hydrodynamic forces of the journal bearings are calculated by using a finite bearing model along with the Gumbel boundary condition. A Newton-Raphson procedure is employed in solving the nonlinear equations for the piston and crankshaft. The developed computer program can be used to calculate the complete trajectories of the piston and the crankshaft as functions of the crank angle under compressor-running conditions. The results explored the effects of the radial clearance of the piston, oil viscosity, and mass and mass moment of inertia of the piston and connecting rod on the stability of the compression mechanism.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.3 s.36
• /
• pp.487-495
• /
• 1998

Recently, the space truss has been attracted by many designers because of their ability to support significant loads with a minimum material. And it is relatively flexible to design the configuration of structures. This paper presents a volume optimization for the space truss on the basis of result evaluated from nonlinear analysis. The optimization of the truss is done by nonlinear optimum GINO(General Interactive Nonlinear Optimizer) program. The objective function considered is the volume of the steel bars. The constraints for optimum design are the design limits, such as the axial force strength, maximum slenderness, minimum thickness, allowable deflection and ratio of the external diameter to thickness of the circular tube bars.

• Journal of the Korean Solar Energy Society
• /
• v.25 no.1
• /
• pp.87-96
• /
• 2005

본 논문은 압력차로 인한 유체의 유동장에서 서스펜션의 입자 밀도를 분포 규명하기 위해 적용할 수 있는 Electric Impedance Tomography (EIT)의 새로운 기법에 대한 효율성을 다루고 있다. Regularized Newton-Raphson iterative method를 근간으로 inverse problem의 해를 구하는데, 이는 곧 목적 함수(object function)를 몇 가지의 제한조건(constraints) 하에서 최소화시키는 과정이라 할 수 있다. 한편, 관련 forward problem은 유한요소법(FEM)을 이용하여 해결하며, 기존의 연구와는 달리 선형 형상 함수(linear shape function)를 이용하여 전도도가 연속적인 물성치로 유동장에 분포되어 있는 것으로 가정하였다. 여러 경우의 test run에 대한 결과는 본 논문에서 적용한 방법론의 타당성을 보여 주고 있다. 태양에너지의 추출을 위해 직접촉식 열교환기가 종종 이용되고 있는데, 본 연구는 열교환기 내부의 분산 유체에 대한 해석에 일조를 할 수 있을 것으로 기대된다.

• Smart Structures and Systems
• /
• v.9 no.3
• /
• pp.273-286
• /
• 2012

As civil structures are exposed to various external loads, it is essential to assess the structural condition, especially the structural displacement, in every moment. Therefore, a visually servoed paired structured light system was proposed in the previous study. The proposed system is composed of two screens facing with each other, each with a camera, a screen, and one or two lasers controlled by a 2-DOF manipulator. The 6-DOF displacement can be calculated from the positions of three projected laser beams and the rotation angles of the manipulators. In the estimation process, one of well-known iterative methods such as Newton-Raphson or extended Kalman filter (EKF) was used for each measurement. Although the proposed system with the aforementioned algorithms estimates the displacement with high accuracy, it takes relatively long computation time. Therefore, an incremental displacement estimation (IDE) algorithm which updates the previously estimated displacement based on the difference between the previous and the current observed data is newly proposed. To validate the performance of the proposed algorithm, simulations and experiments are performed. The results show that the proposed algorithm significantly reduces the computation time with the same level of accuracy compared to the EKF with multiple iterations.

• Structural Engineering and Mechanics
• /
• v.47 no.4
• /
• pp.559-573
• /
• 2013

Geometrically non-linear axisymmetric bending of a shallow spherical shell with a clamped or a simply supported edge under axisymmetric load was investigated numerically. The partial load was introduced by the Heaviside step function, and the solution was obtained by the finite difference and the Newton-Raphson methods. The thickness of the shell was considered to be uniform and the material was assumed to be homogeneous and isotropic. Sensitivity analysis was made for three geometrical parameters. The accuracy of the algorithm was checked by comparing the central deflection, the radial membrane stress at the edge, or the transverse shear force with the solutions of plates and shells in the literature and good agreement was obtained. The main findings of the study can be outlined as follows: (i) If the shell is fully loaded the central deflection of a clamped shell is larger than that of a simply supported shell provided that the shell is not very shallow, (ii) if the shell is partially loaded the central deflection of the shell is sensitive to the parameters of thickness, depth, and partial loading but the influence of the boundary conditions is negligible.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.