• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• 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
• /
• v.54 no.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
• /
• v.26 no.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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1382-1390
• /
• 1990

Accurate load-carrying capacities and moment distributions of thin circular plates are obtained for clamped or simply-supported boundary condition under various concentrated circular loadings. The material is regarded as perfectly-plastic based on an arbitrary yield function such as the Tresca yield function, the Johansen yield function, and the farmily of .betha.-norms which possesses the von Mises yield function and the Frobenius norm. To obtain the lower bound solutions, a maximization formulation is derived and implemented for efficient numerical calculation with a finite difference method and the modified Newton's method. The numerical results demonstrate plastic collapse behavior of circular plates and provide their design criteria.

• Computational Structural Engineering
• /
• v.10 no.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
• /
• v.28 no.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
• /
• v.29 no.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
• /
• v.33 no.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.06a
• /
• pp.131-135
• /
• 2006

A seismic waveform inversion for prestack seismic data based on the Gauss-Newton method is presented. The Gauss-Newton method for seismic waveform inversion was proposed in the 80s but has rarely been studied. Extensive computational and memory requirements have been principal difficulties. To overcome this, we used different sizes of grids in the inversion stage from those of grids in the wave propagation simulation, temporal windowing of the simulation and approximation of virtual sources for calculating partial derivatives, and implemented this algorithm on parallel supercomputers. We show that the Gauss-Newton method has high resolving power and convergence rate, and demonstrate potential applications to real seismic data.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1996.05a
• /
• 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.