• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3181-3204
• /
• 2022

The variational nodal method for solving the neutron transport equation has evolved over 40 years. Based on a functional form of the Boltzmann neutron transport equation, the method now comprises a complete set of variants that can be employed for different problems. This paper presents an extensive review of the development of the variational nodal method. The emphasis is on summarizing the whole theoretical system rather than validating the methodologies. The paper covers the variational nodal formulation of the Boltzmann neutron transport equation, the Ritz procedure for various application purposes, the derivation of boundary conditions, the extension for adjoint and perturbation calculations, and treatments for anisotropic scattering sources. Acceleration approaches for constructing response matrices and solving the resulting system of algebraic equations are also presented.

• Bulletin of the Korean Mathematical Society
• /
• v.59 no.5
• /
• pp.1215-1235
• /
• 2022

In this paper we study Szegö kernel projections for Hardy spaces in quaternionic Clifford analysis. At first we introduce the matrix Szegö projection operator for the Hardy space of quaternionic Hermitean monogenic functions by the characterization of the matrix Hilbert transform in the quaternionic Clifford analysis. Then we establish the Kerzman-Stein formula which closely connects the matrix Szegö projection operator with the Hardy projection operator onto the Hardy space, and we get the matrix Szegö projection operator in terms of the Hardy projection operator and its adjoint. At last, we construct the explicit matrix Szegö kernel function for the Hardy space on the sphere as an example, and get the solution to a Diriclet boundary value problem for matrix functions.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.9
• /
• pp.137-147
• /
• 1995

The unbalance response analysis is one of the essential area in the forced vibration analysis of rotor-bearing systems. Local bearing parameters in rotor-bearing systems are the major sources which give rise to a difficulty in unbalance response computation due to the complicated dynamic properties such as rotational speed dependency and anisotropy. In the present paper, an efficient method for unbalance responses is proposed so as to easily take into account bearing parameters in computation. An exact matrix condensation procedure is proposed which enables the present method to compute unbalance responses by dealing with condensed, small matrices. The proposed method causes no errors even though the computation procedure is based on the small matrices condensed from the full matrices. The present method is illustrated through a numerical example and compared with the conventional method.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.23 no.3
• /
• pp.397-400
• /
• 1986

In this paper, the orthogonal properties of the eigenmodes of optical resonators which have phase conjugate mirrors at both ends are derived. The modes which propagate in resonators are descdribed by the Huygens integral. Then, the eigeneuqation which is needed in order to prove the orthogonality of the eigenmodes of the resonator is obtained from this representation. When the kernel being a part of the eigenequation is Hermitian as in conventional resonators and in optical resonator with only one phase conjugate mirror, one can show that the eigenmodes have essentially biorthogonal relations, which are satisfied the set of modes propagating in one direcdtion around the resonator and the adjoint set of modes propagating in the reversed direction.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.337-343
• /
• 2014

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• Journal of Energy Engineering
• /
• v.16 no.4
• /
• pp.155-163
• /
• 2007

To insure nuclear reactor safety, the reactivity of control rods should be calculated by measuring the criticality of reactor core and it is regularly performed during the annual physics test period. Also, the core criticality should be monitored during the start-up operation to avoid reactivity induced accidents. Many research works on control rod reactivity measurement and subcriticality measurement have been accomplished throughout the world for decades and recently a new method named "Modified Neutron Source Multiplication Method (MNSM)" was proposed in Japan which is known to be improved overcoming limitations of traditional Neutron Source Multiplication Method (NSM). In this study, MNSM was tested in calculation of subcriticalities and in evaluation of application validity using the educational reactor in Kyung Hee University, AGN-201. For this study, a revised nuclear data library and a neutron transport code system TRANSX - PARTISN were established. Correction factors for various control rod positions were produced using the k-effective values and the corresponding flux distributions and adjoint flux distributions. Experimental values of the core criticality were obtained using the neutron count rates of the BF3 proportional counters. The results showed that the expected reactivity worth of control rods by MNSM agreed well with the theoretical values and the correction factors contributed much for this purpose.

• Computational Structural Engineering
• /
• v.1 no.1
• /
• pp.87-94
• /
• 1988

The linerization method as one of the recursive quadratic programming method is applied for the optimal design of a large-scale structure supported by Pshenichny's proof of global convergence of the algorithm and convergence rate estimates. The linearization method transforms all constants of the design problem into an equivalent linearized constraint and employs the active-set strategy. This results in substantial computational savings by reducing the number of sate and adjoint to be solved at every design iteration. The illustrative example of plates with beams supported by columns is the typical one of a large-scale structure to give successful optimum solutions with satisfactory convergence criteria. Hopefully, the method may be applicable to all classes of optimization problems.

• Journal of the Korean Mathematical Society
• /
• v.34 no.4
• /
• pp.1029-1036
• /
• 1997

Let ${A_t)}_{t>0}$ be a dilation group given by $A_t = exp(-P log t)$, where P is a real $n \times n$ matrix whose eigenvalues has strictly positive real part. Let $\nu$ be the trace of P and $P^*$ denote the adjoint of pp. Suppose that $K$ is a function defined on $R^n$ such that $$\mid$K(x)$\mid$ \leq k($\mid$x$\mid$_Q)$ for a bounded and decreasing function $k(t) on R_+$ satisfying $k \diamond $\mid$\cdot$\mid$_Q \in \cup_{\varepsilon >0}L^1((1 + $\mid$x$\mid$)^\varepsilon dx)$ where $Q = \int_{0}^{\infty} exp(-tP^*) exp(-tP)$ dt and the norm $$\mid$\cdot$\mid$_Q$ stands for $$\mid$x$\mid$_Q = \sqrt{}, x \in R^n$. For $f \in L^1(R^n)$, define $mf(x) = sup_{t>0}$\mid$K_t * f(x)$\mid$$ where $K_t(X) = t^{-\nu}K(A_{1/t}^* x)$. Then we show that $m$ is a bounded operator of $L^1(R^n) into L^{1, \infty}(R^n)$.