• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.635-648
• /
• 2015

We consider a weakly coupled system of singularly perturbed convection-diffusion equations with discontinuous source term. The diffusion term of each equation is associated with a small positive parameter of different magnitude. Presence of discontinuity and different parameters creates boundary and weak interior layers that overlap and interact. A numerical method is constructed for this problem which involves an appropriate piecewise uniform Shishkin mesh. The numerical approximations are proved to converge to the continuous solutions uniformly with respect to the singular perturbation parameters. Numerical results are presented which illustrates the theoretical results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.1
• /
• pp.11-18
• /
• 1983

For the steady-state solutions of vibratory systems where the dynamics involves nonlinearity and discontinuity, a method of numerical simulation has been normally used. This paper presents a new approach which may overcome some difficulties in the simulation method. This approach is based on the fundamental assumption that the steady-state forced vibration is periodic, so that the problem is formulated as a two-point boundary-value problem and can be solved by Waner's algorithm. This method is demonstrated through the solutions of a linear system, a system with Coulomb friction and an impact pair. It is found that the method gives true solutions well both for linear and nonlinear systems, which convinces us of the usefulness of the method.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.178-181
• /
• 2011

Recently, micro shock tube is being extensively used in various fields of engineering applications. The flow characteristics occurring in the micro shock tube may be significantly different from that of conventional macro shock tube due to very low Reynolds number and Knudsen number effects which are, in general, manifested in such flows of rarefied gas, solid-gas two-phase, etc. In these situations, Navier-Stokes equations cannot properly predict the micro shock tube flow. In the present study, a two-dimensional CFD method has been applied to simulate the micro shock tube, with slip velocity and temperature jump boundary conditions. The effects of wall thermal conditions on the unsteady flow in the micro shock tube were also investigated. The unsteady behaviors of shock wave and contact discontinuity were, in detail, analyzed. The results obtained show much more attenuation of shock wave, compared with macro-shock tubes.

• Journal of Welding and Joining
• /
• v.18 no.4
• /
• pp.104-110
• /
• 2000

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, an integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in plane finite bodies. The method was developed for in-plane (modes I and II) loadings only. In this paper, the method is formulated and applied to mode III problems involving smooth or kinked cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.implicity.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.48-58
• /
• 1999

The Characteristic Flux Difference Splitting (CFDS) scheme designed to adapt the characteristic boundary conditions at the wall and inflow/outflow boundary planes satisfies Roe's property U, although the CFDS Jacobian matrix is decomposed by a product of elaborate transformation matrices and explicit eigenvalue matrix. When the CFDS algorithm, thus a variant of Roe's scheme, is applied straightforwardly to hypersonic flows over a blunt body, the strong bow shock gradually breaks down near the stagnation point. This numerical instability is widely observed by many researchers employing flux-difference method, known in the literature as the carbuncle phenomenon. Many remedies have been proposed and resulted in partial cures. When the idea of Sanders et al. which identifies the minimum eigenvalues near the discontinuity present is applied to CFDS method, it is shown that the instability problem can be controlled successfully. A few flux splitting methods have also been tested and results are compared against the Nakamori's Mach 8 blunt body flow.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.05b
• /
• pp.643-650
• /
• 1995

A numerical method for treating transient diffusion Involving change of phase is presented. In other methods of dealing with this class of problems, the mass flux balance at the moving phase boundary requires explicit treatment of two distinct phases. The technique, originating from the apparent heat capacity method in transient heat conduction with the phase change, avoids the difficulty by transferring the concentration discontinuity at the boundary to smoothed physical property variations near the moving front. This technique accomodates the nonlinearities which preclude use of analytical solutions. It was tested against known analytical solutions for simple cases and turned out to be quite accurate.

• Journal of the Korean Society of Safety
• /
• v.23 no.5
• /
• pp.91-96
• /
• 2008

A reflected breakwater can be affected by wave pressure and power because it is to be concentrated by wave energy. The present study is to estimate hydraulic behavior affecting around a reflected breakwater, which is discontinuity cases and various angle of coner at the breakwater. The numerical model to investigate wave diffraction, which is important hydraulic factor in the ocean, is performed by using direct boundary element method. The present numerical results are compared with the solutions of approximate and absolute based on an eigenfunction, and the solution of analytical by Fresnel integral. The results of the present numerical simulation agreed well with those of the published numerical and analytical data. As a result of this study, wave height is high at the comer of breakwater, and it is to be high if angle of conner at the reflected breakwater is small.

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.598-610
• /
• 2022

Avoiding a fault zone would be a best practice for safety in underground construction, which is only sometimes possible because of many restrictions and other field conditions. For instance, there is an ongoing conception of Korea-Japan subsea tunnels that inevitably cross a massive fault system in the Korea Strait. Therefore it was deemed necessary to find an efficient way of predicting the likely behaviour of underground structures under fault slip. This paper presents the findings from simple numerical analysis for investigating the stress induced at a normal fault with a dip of 45 degrees. We used a boundary element software that assumed constant displacement discontinuity, which allowed the displacement to be estimated separately at both the fault's hangingwall and footwall sides. The results suggested that a principal stress rotation of 45 degrees occurred at the edges of the fault during the slip, which was in agreement with the phenomenon for fault plane suggested in the body of literature. A simple numerical procedure presented in this paper could be adopted to investigate other fault-related issues associated with underground structure construction.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.35 no.1
• /
• pp.55-62
• /
• 2017

The SRP (Solar Radiation Pressure) model has always been an issue in the dynamic GPS (Global Positioning System) orbit determination. The widely used CODE (Center for Orbit Determination in Europe) model and its variants have nine parameters to estimate the solar radiation pressure from the Sun and to absorb the remaining forces. However, these parameters show a very high correlation with each other and, therefore, only several of them are estimated at most of the IGS (International GNSS Service) analysis centers. In this study, we attempted to numerically verify the correlation between the parameters. For this purpose, a bi-directional, multi-step numerical integrator was developed. The correlation between the SRP parameters was analyzed in terms of post-fit residuals of the orbit. The integrated orbit was fitted to the IGS final orbit as external observations. On top of the parametric analysis of the SRP parameters, we also verified the capabilities of orbit prediction at later time epochs. As a secondary criterion for orbit quality, the positional discontinuity of the daily arcs was also analyzed. The resulting post-fit RMSE (Root-Mean-Squared Error) shows a level of 4.8 mm on average and there is no significant difference between block types. Since the once-per-revolution parameters in the Y-axis are highly correlated with those in the B-axis, the periodic terms in the D- and Y-axis are constrained to zero in order to resolve the correlations. The 6-hr predicted orbit based on the previous day yields about 3 cm or less compared to the IGS final orbit for a week, and reaches up to 6 cm for 24 hours (except for one day). The mean positional discontinuity at the boundary of two 1-day arcs is on the level of 1.4 cm for all non-eclipsing satellites. The developed orbit integrator shows a high performance in statistics of RMSE and positional discontinuity, as well as the separations of the dynamic parameters. In further research, additional verification of the reference frame for the estimated orbit using SLR is necessary to confirm the consistency of the orbit frames.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.23 no.5
• /
• pp.1299-1309
• /
• 1998

In this paper, we propose a new approach for reducing the blocking artifact that is one of drawbacks of the block-based Discrete Cosine Transform (DCT) without introducing additional information or significant blurring. We modify the inter-block discontinuity minimization technique to preserve edges within a block as well as to reduce visible block boundaries. The homogeneity of each block is decided by the threshold value reated to Q-factor, which is included in a JPEG as well as MPEG streams. The quantization error is estimated by minimizing the discontinuity, which is weighted in proportion to block discontinuity and added to each pixel in the block to compensate block artifacts. The proposed algorithm reconstructs images which have less noticeable block boundaries from a subjective viewpoit without anyconstraints.