• Nuclear Engineering and Technology
• /
• 제52권6호
• /
• pp.1137-1147
• /
• 2020

The discrete ordinates method (S_N) is one of the major shielding calculation method, which is suitable for solving deep-penetration transport problems. Our objective is to explore the available quadrature sets and to improve the accuracy in shielding problems involving strong anisotropy. The linear discontinuous finite-element (LDFE) quadrature sets based on the icosahedron (in short, ICLDFE quadrature sets) are developed by defining projected points on the surfaces of the icosahedron. Weights are then introduced in the integration of the discontinuous finite-element basis functions in the relevant angular regions. The multivariate secant method is used to optimize the discrete directions and their corresponding weights. The numerical integration of polynomials in the direction cosines and the Kobayashi benchmark are used to analyze and verify the properties of these new quadrature sets. Results show that the ICLDFE quadrature sets can exactly integrate the zero-order and first-order of the spherical harmonic functions over one-twentieth of the spherical surface. As for the Kobayashi benchmark problem, the maximum relative error between the fifth-order ICLDFE quadrature sets and references is only -0.55%. The ICLDFE quadrature sets provide better integration precision of the spherical harmonic functions in local discrete angle domains and higher accuracy for simple shielding problems.

• Structural Engineering and Mechanics
• /
• 제23권1호
• /
• pp.63-74
• /
• 2006

In this paper, the scattering of harmonic elastic anti-plane shear waves by two collinear cracks in functionally graded materials is investigated by means of nonlocal theory. The traditional concepts of the non-local theory are extended to solve the fracture problem of functionally graded materials. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. To make the analysis tractable, it is assumed that the shear modulus and the material density vary exponentially with coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularities are present at crack tips.

• 대한수학회보
• /
• 제50권1호
• /
• pp.263-273
• /
• 2013

A mean-value result, saying that the difference quotient of a differentiable function in a real interval is a mean value of its derivatives at the endpoints of the interval, leads to the functional equation $$\frac{f(x)-F(y)}{x-y}=M(g(x),\;G(y)),\;x{\neq}y$$, where M is a given mean and $f$, F, $g$, G are the unknown functions. Solving this equation for the arithmetic, geometric and harmonic means, we obtain, respectively, characterizations of square polynomials, homographic and square-root functions. A new criterion of the monotonicity of a real function is presented.

• 대한수학회논문집
• /
• 제17권1호
• /
• pp.165-173
• /
• 2002

We first review Ramaswami's find Apostol's identities involving the Zeta function in a rather detailed manner. We then present corrected, or generalized formulas, or a different method of proof for some of them. We also give closed-form evaluation of some series involving the Riemann Zeta function by an integral representation of ζ(s) and Apostol's identities given here.

• 충청수학회지
• /
• 제22권3호
• /
• pp.467-480
• /
• 2009

In 1859, Bernhard Riemann, in his epoch-making memoir, extended the Euler zeta function $\zeta$(s) (s > 1; $s{\in}\mathbb{R}$) to the Riemann zeta function $\zeta$(s) ($\Re$(s) > 1; $s{\in}\mathbb{C}$) to investigate the pattern of the primes. Sine the time of Euler and then Riemann, the Riemann zeta function $\zeta$(s) has involved and appeared in a variety of mathematical research subjects as well as the function itself has been being broadly and deeply researched. Among those things, we choose to make a further investigation of the following subjects: Evaluation of $\zeta$(2k) ($k {\in}\mathbb{N}$); Approximate functional equations for $\zeta$(s); Series involving the Riemann zeta function.

• 대한수학회지
• /
• 제44권5호
• /
• pp.1163-1184
• /
• 2007

In this presentation dedicated to the tricentennial birth anniversary of the great eighteenth-century Swiss mathematician, Leonhard Euler (1707-1783), we begin by remarking about the so-called Basler problem of evaluating the Zeta function ${\zeta}(s)$ [in the much later notation of Georg Friedrich Bernhard Riemann (1826-1866)] when s=2, which was then of vital importance to Euler and to many other contemporary mathematicians including especially the Bernoulli brothers [Jakob Bernoulli (1654-1705) and Johann Bernoulli (1667-1748)], and for which a fascinatingly large number of seemingly independent solutions have appeared in the mathematical literature ever since Euler first solved this problem in the year 1736. We then investigate various recent developments on the evaluations and representations of ${\zeta}(s)$ when $s{\in}{\mathbb{N}}{\backslash}\;[1],\;{\mathbb{N}}$ being the set of natural numbers. We emphasize upon several interesting classes of rapidly convergent series representations for ${\zeta}(2n+1)(n{\in}{\mathbb{N}})$ which have been developed in recent years. In two of many computationally useful special cases considered here, it is observed that ${\zeta}(3)$ can be represented by means of series which converge much more rapidly than that in Euler's celebrated formula as well as the series used recently by Roger $Ap\'{e}ry$ (1916-1994) in his proof of the irrationality of ${\zeta}(3)$. Symbolic and numerical computations using Mathematica (Version 4.0) for Linux show, among other things, that only 50 terms of one of these series are capable of producing an accuracy of seven decimal places.

• 한국측량학회지
• /
• 제29권5호
• /
• pp.543-550
• /
• 2011

중력이상값은 지구물리, 측지 및 국방 등 다양한 분야에서 활용되는 기초 지구물리 자료로서, 특정 위치에서의 중력이상값을 필요로 하는 경우 일반적으로 데이터베이스화 되어 있는 중력이상값으로부터 내삽하여 활용한다. 그러나 중력은 지형 및 지하광물 등에 의하여 다양하게 변할 수 있는 물리량으로, 내삽에서 가정한 선형성, 2차 곡선 등의 성질이 만족되지 않으면 그 결과로 계산된 중력이상값은 실제 중력값과 큰 차이를 나타내게 된다. 또한, 내삽을 통하여 계산되는 결과값은 이론적으로 조화함수를 만족하여야 한다는 중력의 물리적 성질을 반영하지 못한다. 본 연구에서는 이와 같은 문제점을 보완하기 위하여 필요에 따라 유연하게 중력이상값을 계산할 수 있도록 중력 모델링을 수행하였다. 모델링은 평면푸리에 시리즈와 point-mass 함수를 기저함수로 하는 두 방법을 기반으로 수행되었고, 구축된 모델은 내삽으로부터 산출된 결과와 비교하여 특성을 분석하였다. 또한 모델링의 결과와 내삽 방법을 중력참조항법에 적용하여 활용적인 측면을 검토하였다. 연구결과, 기복이 완만한 지역에서는 평면푸리에 시리즈와 point-mass 및 내삽으로부터 계산된 중력이상값이 유사하게 나타났으나, 중력의 기복이 큰 지역에서는 모델 및 내삽에 의한 결과가 큰 차이를 나타내었다. 특히 주변의 네 점을 이용하여 선형으로 계산하는 Bilinear 내삽함수를 이용한 경우가 가장 완만한 중력값을 보이는 반면 point-mass 함수로부터 산출된 결과가 고주파에서 가장 큰 값을 나타내었다. 또한, 모델링 및 내삽에 필요한 자료의 로딩 및 계산 시간을 비교한 결과, 중력참조항법의 경우 중력값의 계산은 모델링을 수행하는 경우가 데이터베이스에 기반을 둔 내삽보다 효율적임을 알 수 있었다. 본 연구에서는 중력모델링의 결과 및 특성을 분석하였으며, 향후 모델링은 중력참조항법과 같은 활용분야에 있어 가장 효율적인 신호의 특성과 해상도를 지닌 중력 자료를 제공할 수 있는 기술로 활용될 수 있을 것으로 사료된다.