• Bulletin of the Korean Mathematical Society
• /
• v.40 no.3
• /
• pp.385-397
• /
• 2003

In this paper, we introduce and study a system of nonlinear implicit variational inclusions (SNIVI) in real Banach spaces: determine elements $x^{*},\;y^{*},\;z^{*}\;\in\;E$ such that ${\theta}\;{\in}\;{\alpha}T(y^{*})\;+\;g(x^{*})\;-\;g(y^{*})\;+\;A(g(x^{*}))\;\;\;for\;{\alpha}\;>\;0,\;{\theta}\;{\in}\;{\beta}T(z^{*})\;+\;g(y^{*})\;-\;g(z^{*})\;+\;A(g(y^{*}))\;\;\;for\;{\beta}\;>\;0,\;{\theta}\;{\in}\;{\gamma}T(x^{*})\;+\;g(z^{*})\;-\;g(x^{*})\;+\;A(g(z^{*}))\;\;\;for\;{\gamma}\;>\;0,$ where T, g : $E\;{\rightarrow}\;E,\;{\theta}$ is zero element in Banach space E, and A : $E\;{\rightarrow}\;{2^E}$ be m-accretive mapping. By using resolvent operator technique for n-secretive mapping in real Banach spaces, we construct some new iterative algorithms for solving this system of nonlinear implicit variational inclusions. The convergence of iterative algorithms be proved in q-uniformly smooth Banach spaces and in real Banach spaces, respectively.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.39-48
• /
• 2004

In the present study, an efficient implementation technique of the van Leer's implicit operator is suggested in accordance with the Roe's explicit operator. By using an efficient treatment of the off-diagonal terms, which occupy most of the memory requirement for the linear system of equations, it is shown that the improved scheme only requires less than 30% of memory and is approximately 10-20% faster than the baseline scheme.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.909-920
• /
• 2011

A regularized Newton method for nonlinear ill-posed problems is considered. In each Newton step an implicit iterative method with an appropriate stopping rule is proposed and analyzed. Under certain assumptions on the nonlinear operator, the convergence of the algorithm is proved and the algorithm is stable if the discrepancy principle is used to terminate the outer iteration. Numerical experiment shows the effectiveness of the method.

• East Asian mathematical journal
• /
• v.26 no.3
• /
• pp.403-413
• /
• 2010

In this paper, we consider a doubly nonlinear parabolic equation related to the Leray-Lions operator with Dirichlet boundary condition and initial data given. By exploiting a suitable implicit time-discretization technique, we obtain the existence of global strong solution.

• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• Journal of Aerospace System Engineering
• /
• v.14 no.3
• /
• pp.17-23
• /
• 2020

In this paper, the efficient periodic unsteady flow analysis is developed by using a Chebyshev collocation operator applied to the time differential term of the governing equations. The partial implicit time integration method was also applied in the governing equation for a fluid, which means flux terms were implicitly processed for a time integration and the time derivative terms were applied explicitly in the form of the source term by applying the Chebyshev collocation operator. To verify this method, we applied the 1D unsteady Burgers equation and the 2D oscillating airfoil. The results were compared with the existing unsteady flow frequency analysis technique, the Harmonic Balance Method, and the experimental data. The Chebyshev collocation operator can manage time derivatives for periodic and non-periodic problems, so it can be applied to non-periodic problems later.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.17-22
• /
• 1999

Convergence characteristics and efficiency of three implicit approximate factorization schemes(ADI, DDADI and MAF) are examined using 2-Dimensional compressible upwind Navier-Stokes code. Second-order CSCM(Conservative Supra Characteristic Method) upwind flux difference splitting method with Fromm scheme is used for the right-hand side residual evaluation, while generally first-order upwind differencing is used for the implicit operator on the left-hand side. Convergence studies are performed using an example of the flow past a NACA0012 airfoil at steady transonic flow condition, i. e. Mach number 0.8 at $1.25^{\circ}$ angle of attack. The results were compared with other computational results in order to validate the current numerical analysis. The results from the implicit AF algorithms were compared well in low surface with the other computational results; however, not well in upper surface. It might be due to lack of the grid around the shock position. Because the algorithm minimizes the errors of the approximate decomposition, the improved convergence rate with MAF were observed.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.49-60
• /
• 2022

Motivated by the high-resolution properties of high-order Weighted Essentially Non-Oscillatory (WENO) and flux limiter (FL) for steep-gradient problems and the robust convergence of Jacobian-free Newton-Krylov (JFNK) methods for nonlinear systems, the preconditioned JFNK fully implicit high-order WENO and FL schemes are proposed to solve the transient two-phase two-fluid models. Specially, the second-order fully-implicit BDF2 is used for the temporal operator and then the third-order WENO schemes and various flux limiters can be adopted to discrete the spatial operator. For the sake of the generalization of the finite-difference-based preconditioning acceleration methods and the excellent convergence to solve the complicated and various operational conditions, the random vector instead of the initial condition is skillfully chosen as the solving variables to obtain better sparsity pattern or more positions of non-zero elements in this paper. Finally, the WENO_JFNK and FL_JFNK codes are developed and then the two-phase steep-gradient problem, phase appearance/disappearance problem, U-tube problem and linear advection problem are tested to analyze the convergence, computational cost and efficiency in detailed. Numerical results show that WENO_JFNK and FL_JFNK can significantly reduce numerical diffusion and obtain better solutions than traditional methods. WENO_JFNK gives more stable and accurate solutions than FL_JFNK for the test problems and the proposed finite-difference-based preconditioning acceleration methods based on the random vector can significantly improve the convergence speed and efficiency.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.177-180
• /
• 2007

Numerical simulations on the suspended load in the Do jang fish port are carried out. Suspended load is analysed by using the two-dimensional advection-diffusion equation. To describe behaviors of a pollutant in costal zone, a split-operator method is applied to the numerical model. The advection part is first solved by SOWMAC and then the diffusion part is solved by a three-level locally implicit scheme.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.