• Journal of Mechanical Science and Technology
• /
• 제18권4호
• /
• pp.671-680
• /
• 2004

We present numerical and experimental results of the periodic flows inside a rectangular container under a background rotation. In numerical computation, a parallel-computation technique with MPI is implemented. Flow visualization and PIV measurement are also performed to obtain velocity fields at the free surface. Through a series of numerical and experimental works, we aim to clarify the fundamental reasons of discrepancy between the two-dimensional computation and the experimental measurement, which was detected in the previous study for the same flow model. Specifically, we check if the various assumptions prerequisite for the validity of the classical Ekman pumping law are satisfied for periodic flows under a background rotation.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제12권2호
• /
• pp.119-131
• /
• 2008

The paper considers the computation method in the performance evaluation of cellular network in the phase-type distribution assumptions that the channel holding times induced from mobility are modeled by well-fitted distributions to reflect an actual situation. When ww consider a phase-type distribution model instead of exponential distribution, the complexity of the computation increase exponential even though the accuracy is improved. We consider an efficient numerical algorithm to compute the performance evaluations in cellular networks such as a handoff call dropping probability, new call blocking probability, and handoff arrival rate. Numerical experiment shows that numerical analysis results are well approximated to the results of simulation.

• 대한수학회보
• /
• 제50권6호
• /
• pp.2035-2051
• /
• 2013

In this paper, we study numerical schemes for solving multi-dimensional option pricing problem. We compare the direct solving method and the Operator Splitting Method(OSM) by using finite difference approximations. By varying parameters of the Black-Scholes equations for the maximum on the call option problem, we observed that there is no significant difference between the two methods on the convergence criterion except a huge difference in computation cost. Therefore, the two methods are compatible in practice and one can improve the time efficiency by combining the OSM with parallel computation technique. We show numerical examples including the Equity-Linked Security(ELS) pricing based on either two assets or three assets by using the OSM with the Monte-Carlo Simulation as the benchmark.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2004년도 춘계학술발표회
• /
• pp.565-573
• /
• 2004

This study is to present explicit analytical expressions for calculating bearing capacity factor $N_{\gamma}$, to provide results of the numerical computation instead of the graphical method. In this study, $N_{\gamma}$ is proposed in the critical failure surface on assumption that the center of log spiral in the radial shear zone can be located at the any points of around footing. The critical failure surface is one which yields minimum passive pressure $P_{\gamma}$ on the radial shear zone from the family of log spirals accoding to change of the center of log spiral. This study adoptes Terzaghi's bearing capacity principle(e.g., Prandtl's mechanism, limit equilibrium equation, superposition principle) but the soil wedge in an elastic zone makes angle $45^{\circ}+{\phi}/2$ with the horizontal and the location of the log spiral's center.

• Journal of applied mathematics & informatics
• /
• 제5권3호
• /
• pp.717-734
• /
• 1998

In the present paper is presented a new matrix pencil-based numerical approach achieving the computation of the elemen-tary divisors of a given matrix $A \in C^{n\timesn}$ This computation is at-tained without performing similarity transformations and the whole procedure is based on the construction of the Piecewise Arithmetic Progression Sequence(PAPS) of the associated pencil $\lambda I_n$ -A of matrix A for all the appropriate values of $\lambda$ belonging to the set of eigenvalues of A. This technique produces a stable and accurate numerical algorithm working satisfactorily for matrices with a well defined eigenstructure. The whole technique can be applied for the computation of the first second and Jordan canonical form of a given matrix $A \in C^{n\timesn}$. The results are accurate for matrices possessing a well defined canonical form. In case of defective matrices indications of the most appropriately computed canonical form. In case of defective matrices indication of the most appropriately computed canonical form are given.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.325-328
• /
• 2002

We present numerical and experimental result of the rotating flows inside a rectangular container under a background rotation. In numerical computation, a parallel-computation technique with MPI is implemented. Flow visualization and PIV measurement are also performed to obtain velocity fields at the free surface. Through a series of numerical and experimental works, we aim to clarify, if any, the fundamental reasons of discrepancy between the two-dimensional computation and the experimental measurement, which was detected in the previous study for the same flow model. Specifically, we check if the various assumptions prerequisite for the validity of the classical Ekman pumping law are satisfied for periodic flows under a background rotation.

• Nuclear Engineering and Technology
• /
• 제46권5호
• /
• pp.641-654
• /
• 2014

A new computer code, named CORONA (Core Reliable Optimization and thermo-fluid Network Analysis), was developed for the core thermo-fluid analysis of a prismatic gas cooled reactor. The CORONA code is targeted for whole-core thermo-fluid analysis of a prismatic gas cooled reactor, with fast computation and reasonable accuracy. In order to achieve this target, the development of CORONA focused on (1) an efficient numerical method, (2) efficient grid generation, and (3) parallel computation. The key idea for the efficient numerical method of CORONA is to solve a three-dimensional solid heat conduction equation combined with one-dimensional fluid flow network equations. The typical difficulties in generating computational grids for a whole core analysis were overcome by using a basic unit cell concept. A fast calculation was finally achieved by a block-wise parallel computation method. The objective of the present paper is to summarize the motivation and strategy, numerical approaches, verification and validation, parallel computation, and perspective of the CORONA code.

• 대한수학회보
• /
• 제60권1호
• /
• pp.1-22
• /
• 2023

Let L(s, χ) be the Dirichlet L-series associated with an f-periodic complex function χ. Let P(X) ∈ ℂ[X]. We give an expression for ∑^f_n=1 χ(n)P(n) as a linear combination of the L(-n, χ)'s for 0 ≤ n < deg P(X). We deduce some consequences pertaining to the Chowla hypothesis implying that L(s, χ) > 0 for s > 0 for real Dirichlet characters χ. To date no extended numerical computation on this hypothesis is available. In fact by a result of R. C. Baker and H. L. Montgomery we know that it does not hold for almost all fundamental discriminants. Our present numerical computation shows that surprisingly it holds true for at least 65% of the real, even and primitive Dirichlet characters of conductors less than 10⁶. We also show that a generalized Chowla hypothesis holds true for at least 72% of the real, even and primitive Dirichlet characters of conductors less than 10⁶. Since checking this generalized Chowla's hypothesis is easy to program and relies only on exact computation with rational integers, we do think that it should be part of any numerical computation verifying that L(s, χ) > 0 for s > 0 for real Dirichlet characters χ. To date, this verification for real, even and primitive Dirichlet characters has been done only for conductors less than 2·10⁵.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1572-1577
• /
• 2003

High-power pulsed laser ablation under atmospheric pressure is studied utilizing numerical and experimental methods with emphasis on recondensation ratio, and the dynamics of the laser induced vapor flow. In the numerical calculation, the temperature pressure, density and vaporization flux on a solid substrate are first obtained by a heat-transfer computation code based on the enthalpy method, and then the plume dynamics is calculated by using a commercial CFD package. To confirm the computation results, the probe beam deflection technique was utilized for measuring the propagation of a laser induced shock wave. Discontinuities of properties and velocity over the Knudsen layer were investigated. Related with the analysis of the jump condition, the effect of the recondesation ratio on the plume dynamics was examined by comparing the pressure, density, and mass fraction of ablated aluminum vapor. To consider the effect of mass transfer between the ablation plume and air, unlike the most previous investigations, the equation of species conservation is simultaneously solved with the Euler equations. Therefore the numerical model computes not only the propagation of the shock front but also the distribution of the aluminum vapor. To our knowledge, this is the first work that employed a commercial CFD code in the calculation of pulsed ablation phenomena.

• 대한조선학회논문집
• /
• 제53권1호
• /
• pp.1-9
• /
• 2016

It is important to research and understand the physical phenomenon around a semi-submersible offshore structure on waves and currents because the wave run-up and load occurs owing to the waves and currents. In this study, the numerical simulations are performed about flow around a fixed semi-submersible offshore structure. The Modified Marker-density method is adopted in the present computation procedure, this method is one of the various methods to define the free-surface. The present computation results are compared with existing experimental and numerical simulation(VOF method) results. And, the computation results are relatively coincident with the existing results of model test and numerical simulation by VOF method.