• 대한기계학회논문집B
• /
• 제29권8호
• /
• pp.903-910
• /
• 2005

Inverse radiation problems are solved for estimating the boundary conditions such as temperature distribution and wall emissivity in axisymmetric absorbing, emitting and scattering medium, given the measured incident radiative heat fluxes. Various regularization methods, such as hybrid genetic algorithm, conjugate-gradient method and finite-difference Newton method, were adopted to solve the inverse problem, while discussing their features in terms of estimation accuracy and computational efficiency. Additionally, we propose a new combined approach that adopts the hybrid genetic algorithm as an initial value selector and uses the finite-difference Newton method as an optimization procedure.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• Tribology and Lubricants
• /
• 제26권5호
• /
• pp.263-271
• /
• 2010

The rollers of cylindrical roller bearing are axially profiled to relieve high edge stress concentration caused by mainly their finite length and by misalignment. In this paper, a numerical analysis is carried to study the EHL of misaligned (tilted) rollers with axially profiled ends. Using a finite difference method with non-uniform grids and the Newton-Raphson method, the highly nonlinear EHL problems are systematically solved. Physically consistent solutions are obtained for moderate load, material parameters and very small misalignment. For different misalignment angles, contours and sectional plots of pressure and film shape near both edge regions are compared. The asymmetric pressure distributions and film shapes show that the EHL results of finite line contacts are highly dependent upon very small amounts of roller misalignment. Especially, the effect of misalignment on the EHL pressure distribution is much higher than the film shapes.

• 대한기계학회논문집
• /
• 제14권6호
• /
• pp.1382-1390
• /
• 1990

본 연구에서는 Yang의 이론을 바탕으로 여러가지 원형집중하중을 받는 원형평 판에 관한 최대하계해를 구하고자 한다.앞서 열거한 문헌에서는 여러가지 해석법이 제시되었으나, 문제의 성격에 비하여 그 해법이 매우 복잡하였다. 본 논문에서는 원 형평판에 관한 해법으로서 전산기를 응용하여 가능한 한 단순한 해법으로 최대화기법 에 의한 완전한 해를 얻기 위한 시도를 제시하고자 한다.

• 전산구조공학
• /
• 제10권2호
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• Tribology and Lubricants
• /
• 제28권4호
• /
• pp.153-159
• /
• 2012

Tapered roller bearings are widely used in high axial-load and radial-load applications. In this study, a numerical analysis is performed to study a finite line contacts EHL problem between a tapered roller and raceway in tapered roller bearings. Converged solutions are obtained for moderate load and material parameters using a finite difference method with non-uniform grids and the Newton-Raphson method. The contours and sectional plots of pressure distribution and film shape are compared. The pressure distribution and film shapes near both ends of the roller are very different from those in the central part and are transversely asymmetric. The maximum pressure and absolute minimum film thickness always occur at the small end of the roller.

• Tribology and Lubricants
• /
• 제29권6호
• /
• pp.353-359
• /
• 2013

The rollers and/or races in cylindrical and tapered roller bearings should be profiled to relieve high edge stress concentrations caused by their finite lengths and misalignment. In this study, a numerical analysis was performed to investigate the elastohydrodynamic lubrication (EHL) of a Lundberg profile-type cylindrical roller. A finite difference method with fully nonuniform grids and the Newton-Raphson method were used to present detailed EHL pressure distributions and film shapes, as well as the variations in the minimum and central film thicknesses with the profile modification coefficient. In the Lundberg profile, the maximum pressure and minimum film thickness always occurred near the edges. Proper modification of the Lundberg profile considerably increased the minimum film thickness.

• East Asian mathematical journal
• /
• 제33권1호
• /
• pp.53-65
• /
• 2017

Numerical solutions of the Rosenau-Burgers equation based on the cubic B-spline finite element method are introduced. The backward Euler method is used for discretization in time, and the obtained nonlinear algebraic system is changed to a linear system by the Newton's method. We show that those methods are unconditionally stable. Two test problems are studied to demonstrate the accuracy of the proposed method. The computational results indicate that numerical solutions are in good agreement with exact solutions.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2006년도 공동학술대회 논문집
• /
• pp.131-135
• /
• 2006

본 논문에서는 가우스 뉴튼법을 이용한 중합전 탄성파 자료의 파형역산에 관한 연구를 수행하였다. 탄성파 파형역산에 가우스 뉴튼법을 적용하는 방법은 80년대에 제시되었으나 최근 들어서야 활발히 연구가 진행되고 있는데 이는 연산 능력과 기억용량의 한계에 기인한 것이다. 이를 극복하기 위해 본 연구에서는, 파동 전파 수치모의와 역산과정에서 각각 다른 크기의 격자간격을 사용하고, 필요한 시간영역의 파동전파 모사와 가상 진원의 근사를 통해 편미분 파형을 계산하였으며, 효과적으로 슈퍼컴퓨터를 활용하기 위해 병렬처리 기법을 사용하였다. 수치모의를 통해, 가우스 뉴튼법을 이용한 파형 역산의 수렴속도가 빠르고 정확한 것을 알 수 있었으며, 이를 통해 본 연구에서 제시한 방법의 실제 탄성파 자료를 이용한 역산에의 적용가능성을 확인하였다.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1996년도 제23회 학술대회
• /
• pp.30-35
• /
• 1996

An elastohydrodynamic lubrication (EHL) analysis for an axially crown profiled cylindrical roller is carried out using a finite difference method and the Newton-Raphson method. Variations of the minimum film thickness with dimensionless parameters show considerably different from those of infinite solutions.