• East Asian mathematical journal
• /
• v.40 no.1
• /
• pp.109-118
• /
• 2024

This paper proposes a numerical method based on domain decomposition to find approximate solutions for one-dimensional convection-diffusion equations with Neumann boundary conditions. First, the equations are transformed into convection-diffusion equations with Dirichlet conditions. Second, the author introduces the Prediction/Correction Domain Decomposition (PCDD) method and estimates errors for the interface prediction scheme, interior scheme, and correction scheme using known error estimations. Finally, the author compares the PCDD algorithm with the fully explicit scheme (FES) and the fully implicit scheme (FIS) using three examples. In comparison to FES and FIS, the proposed PCDD algorithm demonstrates good results.

• The Pure and Applied Mathematics
• /
• v.17 no.4
• /
• pp.289-298
• /
• 2010

In this paper, we propose a second-order prediction/correction (SPC) domain decomposition method for solving one dimensional linear hyperbolic partial differential equation $u_{tt}+a(x,t)u_t+b(x,t)u=c(x,t)u_{xx}+{\int}(x,t)$. The method can be applied to variable coefficients problems and singular problems. Unconditional stability and error analysis of the method have been carried out. Numerical results support stability and efficiency of the method.

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

A continuum-based design sensitivity analysis (DSA) method fur non-shape problems is developed for geometrically nonlinear elastic structures. The non-shape problem is characterized by the design variables that are not associated with the domain of system like sizing, material property, loading, and so on. Total Lagrangian formulation with the Green-Lagrange strain and the second Piola-Kirchhoff stress is employed to describe the geometrically nonlinear structures. The spatial domain is discretized using the 4-node isoparametric plane stress/strain elements. The resulting nonlinear system is solved using the Newton-Raphson iterative method. To take advantage of the derived analytical sensitivity In topology optimization, a fast and efficient design sensitivity analysis method, adjoint variable method, is employed and the material property of each element is selected as non-shape design variable. Combining the design sensitivity analysis method and a gradient-based design optimization algorithm, an automated design optimization method is developed. The comparison of the analytical sensitivity with the finite difference results shows excellent agreement. Also application to the topology design optimization problem suggests a very good insight for the layout design.

• Communications of the Korean Mathematical Society
• /
• v.9 no.2
• /
• pp.439-448
• /
• 1994

It is known that the steepest descent(SD) method and the conjugate gradient(CG) method [1, 2, 5, 6] converge when these methods are applied to solve linear systems of the form Ax = b, where A is symmetric and positive definite. For some finite difference discretizations of elliptic problems, one gets positive definite matrices that are almost symmetric. Practically, the SD method and the CG method work for these matrices. However, the convergence of these methods is not guaranteed theoretically. The SD method is also called Orthores(1) in iterative method papers. Elman [4] states that the convergence proof for Orthores($\kappa$), with $\kappa$ a positive integer, is not heard. In this paper, we prove that the SD method and the CG method converge when the $\iota$$^2$ matrix norm of the non-symmetric part of a positive definite matrix is less than some value related to the smallest and the largest eigenvalues of the symmetric part of the given matrix.(omitted)

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.107-119
• /
• 2022

The aim of the work presented in this paper is development of numerical model for prediction of temperature distribution in pavement according to the measured meteorological parameters, with introduction of non-linear heat transfer coefficient which is a function of temerature difference between the air and the pavement. Developed model calculates heat radiated from the pavement back in the air, which is an important part of the heat trasfer process in the open air surfaces. Temperature of the pavement surface, heat radiation together with many meteorological parameters were measured in series during two years in order to validate the model and calibrate model parameters. Special finite element method for temperature heat transfer towards the soil together with the time integration scheme are used to solve the governing equation. It is proved that non-linear heat transfer coefficient, which is a function of time and temperature difference between the air and the pavement, is required to decribe this phenomena. Proposed model includes heat tranfer coefficient callibration for specific climate region, through the iterative inverse procedure.

• The Pure and Applied Mathematics
• /
• v.13 no.4 s.34
• /
• pp.281-294
• /
• 2006

Many partial differential equations defined on a rectangular domain can be solved numerically by using a domain decomposition method. The most commonly used decompositions are the domain being decomposed in stripwise and rectangular way. Theories for non-overlapping domain decomposition(in which two adjacent subdomains share an interface) were often focused on the stripwise decomposition and claimed that extensions could be made to the rectangular decomposition without further discussions. In this paper we focus on the comparisons of the two ways of decompositions. We consider the unconditionally stable scheme, the MIP algorithm, for solving parabolic partial differential equations. The SOR iterative method is used in the MIP algorithm. Even though the theories are the same but the performances are different. We found out that the stripwise decomposition has better performance.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.2
• /
• pp.93-99
• /
• 1998

Parabolic approximation and sponge layer are applied as open boundary condition for elliptic finite difference wave model. Generalized conjugate gradient method is used as a solution procedure. Using parabolic approximation a large part of spurious reflection is removed at the spherical shoal experiment and sponge layer boundary condition needs more than 2 wave lengths of sponge layer to give similar results. Simulating the propagation of waves on a rectangular harbor, it is identified that iterative scheme can be applied easily for the non-rectangular computational region.

• Geomechanics and Engineering
• /
• v.4 no.4
• /
• pp.263-279
• /
• 2012

The present paper deals with the analysis of combined footings resting on geosynthetic reinforced granular fill overlying stone column improved poor soil. An attempt has been made to study the influence of inclusion of geosynthetic layer on the deflection of the footing. The footing has been idealized as a beam having finite flexural rigidity. Granular fill layer has been represented by Pasternak shear layer and stone columns and poor soil have been represented by nonlinear Winkler springs. Nonlinear behavior of granular fill layer, stone columns and the poor soil has been considered by means of hyperbolic stress strain relationships. Governing differential equations for the soil-foundation system have been derived and solution has been obtained employing finite difference scheme by means of iterative Gauss Elimination method. Results of a detailed parametric study have been presented, for a footing supporting typically five columns, in non-dimensional form in respect of deflection with and without geosynthetic inclusion. Geosynthetic layer has been found to significantly reduce the deflection of the footing which has been quantified by means of parametric study.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.