• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1247-1256
• /
• 2009

Finite difference schemes are considered for a nonlinear $Ca^{2+}$ diffusion equations with stationary and mobile buffers. The scheme inherits mass conservation as for the classical solution. Stability and $L^{\infty}$ error estimates of approximate solutions for the corresponding schemes are obtained. using the extended Lax-Richtmyer equivalence theorem.

• Korean Journal of Mathematics
• /
• v.18 no.3
• /
• pp.299-309
• /
• 2010

In this paper, we study the first order nonlinear neutral difference equation: $${\Delta}(x(n)+px(n-{\tau}))+f(n,x(n-c),x(n-d))=r(n),\;n{\geq}n_0$$. Using the Banach fixed point theorem, we prove the existence of bounded positive solutions of the equation, suggest Mann iterative schemes of bounded positive solutions, and discuss the error estimates between bounded positive solutions and sequences generated by Mann iterative schemes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.299-306
• /
• 2008

Finite difference schemes are considered for a $Ca^{2+}$ diffusion equations, which discribe $Ca^{2+}$ buffering by using stationary and mobile buffers. Stability and $L^\infty$ error estimates of approximate solutions for the corresponding schemes are obtained using the extended Lax-Richtmyer equivalence theorem.

• East Asian mathematical journal
• /
• v.24 no.4
• /
• pp.407-414
• /
• 2008

Finite difference schemes are considered for a $Ca^{2+}$ diffusion equations with damping and convection terms, which describe $Ca^{2+}$ buffering by using stationary and mobile buffers. Stability and $L^{\infty}$ error estimates of approximate solutions for the corresponding schemes are obtained using the extended Lax-Richtmyer equivalence theorem.

• Journal of applied mathematics & informatics
• /
• v.23 no.1_2
• /
• pp.571-579
• /
• 2007

Finite difference schemes are considered for a generalization of the Cahn-Hilliard equation with Neumann boundary conditions and the Kuramoto-Sivashinsky equation with a periodic boundary condition, which is of the type $ut+\frac{{\partial}^2} {{\partial}x^2}\;g\;(u,\;u_x,\;u_{xx})=f(u,\;u_x,\;u_{xx})$. Stability and $L^{\infty}$ error estimates of approximate solutions for the corresponding schemes are obtained using the extended Lax-Richtmyer equivalence theorem.

• Journal of the Korean Mathematical Society
• /
• v.32 no.1
• /
• pp.115-139
• /
• 1995

Our purpose in this paper is to prove a new version of the nonlinear Chernoff theorem and to discuss the equivalence between resolvent consistency and converge nce for nonlinear algorithms acting on different Banach spaces. Such results are useful in the numerical treatment of partial differential equations via difference schemes.

• Nuclear Engineering and Technology
• /
• v.53 no.8
• /
• pp.2696-2707
• /
• 2021

The purpose of this study is to investigate the efficiency of various structural modeling schemes for evaluating seismic performances and fragility of the reactor containment building (RCB) structure in the advanced power reactor 1400 (APR1400) nuclear power plant (NPP). Four structural modeling schemes, i.e. lumped-mass stick model (LMSM), solid-based finite element model (Solid FEM), multi-layer shell model (MLSM), and beam-truss model (BTM), are developed to simulate the seismic behaviors of the containment structure. A full three-dimensional finite element model (full 3D FEM) is additionally constructed to verify the previous numerical models. A set of input ground motions with response spectra matching to the US NRC 1.60 design spectrum is generated to perform linear and nonlinear time-history analyses. Floor response spectra (FRS) and floor displacements are obtained at the different elevations of the structure since they are critical outputs for evaluating the seismic vulnerability of RCB and secondary components. The results show that the difference in seismic responses between linear and nonlinear analyses gets larger as an earthquake intensity increases. It is observed that the linear analysis underestimates floor displacements while it overestimates floor accelerations. Moreover, a systematic assessment of the capability and efficiency of each structural model is presented thoroughly. MLSM can be an alternative approach to a full 3D FEM, which is complicated in modeling and extremely time-consuming in dynamic analyses. Specifically, BTM is recommended as the optimal model for evaluating the nonlinear seismic performance of NPP structures. Thereafter, linear and nonlinear BTM are employed in a series of time-history analyses to develop fragility curves of RCB for different damage states. It is shown that the linear analysis underestimates the probability of damage of RCB at a given earthquake intensity when compared to the nonlinear analysis. The nonlinear analysis approach is highly suggested for assessing the vulnerability of NPP structures.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.7-12
• /
• 1998

Optimized high-order compact(OHOC) schemes were proposed, which have high spatial order of truncation and resolution to simulate the aeroacoustic problems due to unsteady compressible flows. Generally, numerical schemes are categorized explicit or implicit by time-marching method. In this research, OHOC differences which were developed with explicit time-marching method is used to have implicit formulation and the implicit OHOC differences result in block hepta-diagonal matrix. This paper presents the comparisons between the explicit and implicit OHOC schemes with a second order accuracy of time in the 1-d linear wave convection problem, and between the explicit OHOC scheme of 4th-order accuracy in time and the implicit OHOC scheme of 1st-order accuracy in tine for the 1-d nonlinear wave convection problem. With these comparisons, the characteristics of implicit OHOC scheme are shown in the point of CFL number.

• 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.