• Communications of the Korean Mathematical Society
• /
• v.26 no.4
• /
• pp.709-720
• /
• 2011

We obtain special type of differential equations which their solution are random variable with known continuous density function. Stochastic differential equations (SDE) of continuous distributions are determined by the Fokker-Planck theorem. We approximate solution of differential equation with numerical methods such as: the Euler-Maruyama and ten stages explicit Runge-Kutta method, and analysis error prediction statistically. Numerical results, show the performance of the Rung-Kutta method with respect to the Euler-Maruyama. The exponential two parameters, exponential, normal, uniform, beta, gamma and Parreto distributions are considered in this paper.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.4
• /
• pp.365-373
• /
• 2011

Damage of reinforced concrete panel subjected to blast load using parallel and domain decomposition is analyzed. The numerical results are sensitive to the mesh size because blast waves are generated during the extremely short term. In order to investigate the effect of mesh size on the blast wave, the analysis results from various wave mesh size using AUTODYN, the explicit finite element analysis program, were compared with existing experimental results. The smaller mesh size was, the higher accuracy was. However, in this case, the analysis was inefficient. Therefore, in order to increase numerical efficiency, the parallel analysis using decomposed method based on Euler and Lagrangian description was performed. Finally, the decomposed method using both the structure domain based on Lagrange description and the blast wave domain based on Euler description was more efficient than the decomposed method using only the Lagrange mesh on structure domain.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.1
• /
• pp.63-70
• /
• 1994

The three-dimensional incompressible clavier-Stokes equations are solved to simulate the flow field around a Wigley model with free-surface. The IAF(Implicit Approximate Factorization) method is used to show a good success in reducing the computing time. The CPU time is almost an half of that if the IAF method were used. The present method adopts the local linearization and Euler implicit scheme without the pressure-gradient terms for the artificial viscosity. Calculations are carried out at the Reynolds number of $10^6$ and the Froude numbers are 0.25, 0.289 and 0.316. For the approximations of turbulence, the Baldwin-Lomax model is used. The resulting free-surface wave configurations and the velocity vectors are compared with those by the explicit method and experiments.

• Kyungpook Mathematical Journal
• /
• v.52 no.3
• /
• pp.327-346
• /
• 2012

In this paper, we present a simple explicit type numerical method for discretizations in time for solving one dimensional Burgers' equations. The proposed method does not need an iteration process that may be required in most implicit methods and have good convergence and efficiency in computational sense compared to other known numerical methods. For evidences, several numerical demonstrations are also provided.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.49 no.4
• /
• pp.233-239
• /
• 2000

There are many difficult problems in analyzing the gas flow in puffer type circuit breaker such as complex geometry, moving boundary, shock wave and so on. To predict the interruption performance accurately, these should be considered in the simulation. In this paper, the analysis procedure of the cold gas flow in the circuit breaker is presented. Euler equation is solved by FVFLIC method which is an explicit time difference scheme for an unsteady flow computation. Moving boundaries are treated with a cell elimination-addition technique. The pressure and density in front of piston are calculated from the rate of the cell volume change. The presented method is applied to the real circuit breaker model and the pressure in front of the piston is good agreement with the experimental one.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.66-72
• /
• 1997

Optimum nozzle design exploiting the method of characteristic(M.O.C) has been in application as an efficient design methodology targeting a less weighted and short expansion nozzle. This paper treats the optimum nozzle design and the analysis of the inviscid compressible flow inside. Based on traditional Rao's method, the optimum nozzle design is coded with minor modifications for the identification of the control surface across which the mass flux should be conserved. Internal flow field is simulated numerically by M.O.C and implicit/explicit Taylor-Galerkin finite element method(F.E.M) with the aid of adaptive remeshing to capture the shock wave, hence improve the accuracy. Designed and calculated flow fields due to the separate analyses show that the mass flux predicted by optimum nozzle design with M.O.C is not conserved across the control surface and the sonic line should be located upstream of the nozzle throat. Rao's optimum nozzle design methodology exaggerates the momentum thrust and tends to overemphasize the engine performance loss.

• Convergence Security Journal
• /
• v.5 no.3
• /
• pp.93-97
• /
• 2005

MacCormack's method is an explicit, second order finite difference scheme that is widely used in the solution of hyperbolic partial differential equations. Apparently, however, it has shown entropy violations under small discontinuity. This non-physical shock grows fast and eventually all the meaningful information of the solution disappears. Some modifications of MacCormack's methods follow ideas of central schemes with an advantage of second order accuracy for space and conserve the high order accuracy for time step also. Numerical results are shown to perform well for the one-dimensional Burgers' equation and Euler equations gas dynamic.

• Journal of Energy Engineering
• /
• v.6 no.1
• /
• pp.52-57
• /
• 1997

Dynamic characteristics of compressible flow fields in super-circled constricted tube have been studied numerically. By applying MacCormack's explicit scheme, time marching method with predictor/corrector step, Euler equation is solved to find characteristics of fluid flow in a constricted tube where a two-dimensional inviscid compressible flow is assumed. The effects of tube diameter and aspect ratios on the pressure variations are discussed extensively. The results of the developed numerical schemes are compared with those of commercial FLUENT code, and show a good agreement.

• Bulletin of the Korean Mathematical Society
• /
• v.61 no.3
• /
• pp.671-697
• /
• 2024

In this paper, we formally introduce the notion of a general parametric digamma function Ψ(−s; A, a) and we find the Laurent expansion of Ψ(−s; A, a) at the integers and poles. Considering the contour integrations involving Ψ(−s; A, a), we present some new identities for infinite series involving Dirichlet type parametric harmonic numbers by using the method of residue computation. Then applying these formulas obtained, we establish some explicit relations of parametric linear Euler sums and some special functions (e.g. trigonometric functions, digamma functions, Hurwitz zeta functions etc.). Moreover, some illustrative special cases as well as immediate consequences of the main results are also considered.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.