• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.34-39
• /
• 1996

A consistent general order nodal method for solving the 3-D neutron diffusion equation in (x-y-z) geometry has ben derived by using a weighted integral technique and expanding the spatial variables by the Legendre orthogonal series function. The equation set derived can be converted into any order nodal schemes. It forms a compact system for general order of nodal moments. The method utilizes the analytic solutions of the transverse-integrated quasi -one dimensional equations and a consistent expansion for the spatial variables so that it renders the use of an approximation for the transverse leakages no necessary. Thus, we can expect extremely accurate solutions and the solution would converge exactly when the mesh width is decreased or the approximation order is increased since the equation set is consistent mathematically.

• Kyungpook Mathematical Journal
• /
• v.50 no.2
• /
• pp.177-193
• /
• 2010

Working with the various special functions of mathematical physics and applied mathematics we often encounter inverse relations of the type $F_n(x)=\sum\limits_{k=0}^{n}A^n_kG_k(x)$ and $ G_n(x)=\sum\limits_{k=0}^{n}B_k^nF_k(x)$, where 0, 1, 2,$\cdots$. Here $F_n(x)$, $G_n(x)$ denote special polynomial functions, and $A_k^n$, $B_k^n$ denote coefficients found by use of the orthogonal properties of $F_n(x)$ and $G_n(x)$, or by skillful series manipulations. Typically $G_n(x)=x^n$ and $F_n(x)=P_n(x)$, the n-th Legendre polynomial. We give a collection of inverse series pairs of the type $f(n)=\sum\limits_{k=0}^{n}A_k^ng(k)$ if and only if $g(n)=\sum\limits_{k=0}^{n}B_k^nf(k)$, each pair being based on some reasonably well-known special function. We also state and prove an interesting generalization of a theorem of Rainville in this form.

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.430-438
• /
• 2021

The spatial distribution of neutron flux or reaction rate was calculated by cell or mesh tally in traditional Monte Carlo simulation. However, either cell or mesh tally leads to the increase of memory consumption and simulation time. In this paper, the function expansion tally (FET) method was developed in Reactor Monte Carlo code RMC to solve this problem. The FET method was applied to the tallies of neutron flux distributions of uranium block and PWR fuel rod models. Legendre polynomials were used in the axial direction, while Zernike polynomials were used in the radial direction. The results of flux, calculation time and memory consumption of different expansion orders were investigated, and compared with the mesh tally. Results showed that the continuous distribution of flux can be obtained by FET method. The flux distributions were consistent with that of mesh tally, while the memory consumption and simulation time can be effectively reduced. Finally, the convergence analysis of coefficients of polynomials were performed, and the selection strategy of FET order was proposed based on the statistics uncertainty of the coefficients. The proposed method can help to determine the order of FET, which was meaningful for the efficiency and accuracy of FET method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.28 no.2
• /
• pp.147-154
• /
• 2017

The impedance approximation has been widely used to model an earth surface such as ocean surface. In calculation of the dyadic Green's function for the impedance half plane, Sommerfeld integral and its partial derivatives are required. It is known that two far-field approximation of the Sommerfeld integral can be represented in terms of Legendre or Laguerre polynomials. Hence, a criterion is required to choose one of two far-field approximations for a given application, which can be expressed in a complex plane of the surface impedance. Also, we approximate the required partial derivatives of Sommerfeld integral and numerically verify the accuracy of the approximation.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.78.2-78.2
• /
• 2019

We develop an Alcock-Paczynski (AP) test method that uses the evolution of redshift-space two-point correlation function (2pCF) of galaxies. The method improves the AP test proposed by Li et al. (2015) in that it uses the full two-dimensional shape of the correlation function. Similarly to the original method, the new one uses the 2pCF in redshift space with its amplitude normalized. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF. This is because the 2pCF should not change apart from the expected small non-linear evolution if galaxy clustering is not distorted by incorrect choice of cosmology used to convert redshift to comoving distance. Our new method decomposes the redshift difference of the 2-dimensional correlation function into the Legendre polynomials whose amplitudes are modelled by radial fitting functions. The shape of the normalized 2pCF suffers from small intrinsic time evolution due to non-linear gravitational evolution and change of type of galaxies between different redshifts. It can be accurately measured by using state of the art cosmological simulations. We use a set of our Multiverse simulations to find that the systematic effects on the shape of the normalized 2pCF are quite insensitive to change of cosmology over \Omega_m=0.21 - 0.31 and w=-0.5 - -1.5. Thanks to this finding, we can now apply our method for the AP test using the non-linear systematics measured from a single simulation of the fiducial cosmological model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.481-486
• /
• 2007

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.1
• /
• pp.61-68
• /
• 2006

The p-version nonlinear finite element model has been developed to analyze the nonlinear behavior of simply supported RC slabs strengthened with carbon fiber reinforced plastic sheets. The shape function is adopted with integral of Legendre polynomials. The compression model of concrete is based on the Kupfer's yield criterion, hardening rule, and crushing condition. The cracking behavior is modeled by a smeared crack model. In this study, the fixed crack approach is adopted as being geometrically fixed in direction once generated. Each steel layer has a uniaxial behavior resisting only the axial force in the bar direction. Identical behavior is assumed fur tension and compression of steel according to the elastic modulus. The carbon fiber reinforced plastic sheets are considered as reinforced layers of equivalent thickness with uniaxial strength and rigidity properties in the present model. It is shown that the proposed model is able to adequately predicte the displacement and ultimate load of nonlinear simply supported RC slabs by a patch with respect to reinforcement ratio, thickness and angles of CFRP sheets.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.319-326
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted for in the sense of von Karman hypothesis. The material model Is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The Integrals of Legendre Polynomials we used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone

• The Bulletin of The Korean Astronomical Society
• /
• v.37 no.2
• /
• pp.210.2-210.2
• /
• 2012

Using Fresnel reflection and Hapke BRDF model with Apollo 10084 soil sample's scattering properties, we constructed a real scale optical lunar model and used it to simulate lunar apparent albedo and moonshine. For Fresnel reflection, the refractive index of $1.68{\pm}0.5$ was used. For Hapke BRDF parameters from BUGs BRDF measurement, the single scattering with w=0.33, hot spot width h=0.017, average phase angle ${\zeta}$=-0.086 and Legendre polynomial coefficients b=0.308, c=0.425 in wavelength 700nm with two types of Henyey-Greenstein phase function was applied. The computation model includes the Sun as a Lambertian scattering sphere, emitting 1.5078 W/m2 at 700nm in wavelength. The Sun and Moon models were then imported into the IRT based radiative transfer computation. The trial simulation of the irradiance levels of moonshine lights shows that they agree well with the ROLO measurement data. We then estimate the lunar apparent albedo to 0.11. The results are to be compared with the measurement data.

• Kyungpook Mathematical Journal
• /
• v.63 no.2
• /
• pp.155-166
• /
• 2023

In 2018, Andrews introduced the partition functions ${\mathfrak{EO}}$(n) and ${\bar{\mathfrak{EO}}}$(n). The first of these denotes the number of partitions of n in which every even part is less than each odd part, and the second counts the number of partitions enumerated by the first in which only the largest even part appears an odd number of times. In 2021, Pore and Fathima introduced a new partition function ${\mathfrak{EO}}_e$(n) which counts the number of partitions of n which are enumerated by ${\bar{\mathfrak{EO}}}$(n) together with the partitions enumerated by ${\bar{\mathfrak{EO}}}$(n) where all parts are odd and the number of parts is even. They also proved some particular congruences for ${\bar{\mathfrak{EO}}}$(n) and ${\mathfrak{EO}}_e$(n). In this paper, we establish infinitely many families of congruences modulo 2, 4, 5 and 8 for ${\bar{\mathfrak{EO}}}$(n) and modulo 4 for ${\mathfrak{EO}}_e$(n). For example, if p ≥ 5 is a prime with Legendre symbol $({\frac{-3}{p}})=-1$, then for all integers n ≥ 0 and α ≥ 0, we have ${\bar{\mathfrak{EO}}}(8{\cdot}p^{2{\alpha}+1}(pn+j)+{\frac{19{\cdot}p^{2{\alpha}+2}-1}{3}}){\equiv}0$ (mod 8); 1 ≤ j ≤ (p - 1).