• 제목/요약/키워드: Thomas method

검색결과 180건 처리시간 0.032초

NUMERICAL METHOD FOR SINGULAR PERTURBATION PROBLEMS ARISING IN CHEMICAL REACTOR THEORY

  • Andargie, Awoke
    • Journal of applied mathematics & informatics
    • /
    • 제28권1_2호
    • /
    • pp.411-423
    • /
    • 2010
  • In this paper, a numerical method for singular perturbation problems arising in chemical reactor theory for general singularly perturbed two point boundary value problems with boundary layer at one end(left or right) of the underlying interval is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

NUMERICAL INTEGRATION METHOD FOR SINGULAR PERTURBATION PROBLEMS WITH MIXED BOUNDARY CONDITIONS

  • Andargie, Awoke;Reddy, Y.N.
    • Journal of applied mathematics & informatics
    • /
    • 제26권5_6호
    • /
    • pp.1273-1287
    • /
    • 2008
  • In this paper, the numerical integration method for general singularly perturbed two point boundary value problems with mixed boundary conditions of both left and right end boundary layer is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

  • PDF

A MULTISCALE MORTAR MIXED FINITE ELEMENT METHOD FOR SLIGHTLY COMPRESSIBLE FLOWS IN POROUS MEDIA

  • Kim, Mi-Young;Park, Eun-Jae;Thomas, Sunil G.;Wheeler, Mary F.
    • 대한수학회지
    • /
    • 제44권5호
    • /
    • pp.1103-1119
    • /
    • 2007
  • We consider multiscale mortar mixed finite element discretizations for slightly compressible Darcy flows in porous media. This paper is an extension of the formulation introduced by Arbogast et al. for the incompressible problem [2]. In this method, flux continuity is imposed via a mortar finite element space on a coarse grid scale, while the equations in the coarse elements (or subdomains) are discretized on a fine grid scale. Optimal fine scale convergence is obtained by an appropriate choice of mortar grid and polynomial degree of approximation. Parallel numerical simulations on some multiscale benchmark problems are given to show the efficiency and effectiveness of the method.

The Effect of Induced Multipoles on the Fifth-order Raman Response

  • Jansen, Thomas I.C.;Duppen, Koos;Snijders, Jaap G.
    • Bulletin of the Korean Chemical Society
    • /
    • 제24권8호
    • /
    • pp.1102-1106
    • /
    • 2003
  • In our previous work we developed the Finite Field method in order to calculate the fifth-order Raman response. The method was applied to calculate various polarization components of the two-dimensional response of liquid $CS_2$. So far, all calculations relied on the dipole-induced dipole. Accurate time-dependent density functional theory calculations have shown that this model has big discrepancies, when molecules are close together as in the liquid. We now report results of investigations on the importance of multipole and electron overlap effects on the polarizability and the fifth-order Raman response. It is shown that these collision effects, especially the induced multipoles, are crucial in the description of the fifth-order response. The impact is found to be especially pronounced for the χ_{mmzzzz}^{(5)}$response that is solely due to interaction induced effects. The calculated response will be compared with various experimental results.

A bivariate extension of the Hosking and Wallis goodness-of-fit measure for regional distributions

  • Kjeldsen, Thomas Rodding;Prosdocimi, Ilaria
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2015년도 학술발표회
    • /
    • pp.239-239
    • /
    • 2015
  • This study presents a bivariate extension of the goodness-of-fit measure for regional frequency distributions developed by Hosking and Wallis [1993] for use with the method of L-moments. Utilising the approximate joint normal distribution of the regional L-skewness and L-kurtosis, a graphical representation of the confidence region on the L-moment diagram can be constructed as an ellipsoid. Candidate distributions can then be accepted where the corresponding the oretical relationship between the L-skewness and L-kurtosis intersects the confidence region, and the chosen distribution would be the one that minimises the Mahalanobis distance measure. Based on a set of Monte Carlo simulations it is demonstrated that the new bivariate measure generally selects the true population distribution more frequently than the original method. An R-code implementation of the method is available for download free-of-charge from the GitHub code depository and will be demonstrated on a case study of annual maximum series of peak flow data from a homogeneous region in Italy.

  • PDF

A PRACTICAL THREE-DIMENSIONAL ESTIMATION TECHNIQUE FOR SPATIAL DISTRIBUTION OF GROUNDWATER CONTAMINANT CONCENTRATIONS

  • Richard Ewing;Kang, Sung-Kwon;Kim, Jeon-Gook;Thomas B.Stauffer
    • 대한수학회지
    • /
    • 제38권3호
    • /
    • pp.523-559
    • /
    • 2001
  • To predict the fate of groundwater contaminants, accurate spatially continuous information is needed. Because most field sampling of groundwater contaminants are not conducted spatially continuous manner, a special estimation technique is required to interpolate/extrapolate concentration distributions at unmeasured locations. A practical three-dimensional estimations method for in situ groundwater contaminant concentrations is introduced. It consistas of two general steps: estimation of macroscopic transport process and kriging. Using field data and nonlinear optimization techniques, the macroscopic behavior of the contaminant plume is estimated. A spatial distribution of residuals is obtained by subtracting the macroscopic transport portion from field data, then kriging is applied to estimate residuals at unsampled locations. To reduce outlier effects on obtaining correlations between residual data which are needed for determining variougram models, the R(sub)p-estimator is introduced. The proposed estimation method is applied to a field data set.

  • PDF

Iterative identification methods for ill-conditioned processes

  • Lee, Jietae;Cho, Wonhui;Edgar, Thomas F.
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 제어로봇시스템학회 1997년도 한국자동제어학술회의논문집; 한국전력공사 서울연수원; 17-18 Oct. 1997
    • /
    • pp.1762-1765
    • /
    • 1997
  • Some ill-conditioned processes are very sensitive to small element-wise uncertainties arising in classical element-by-element model identifications. For such processes, accurate identification of simgular values and right singular vectors are more important than theose of the elements themselves. Singular values and right singular vectors can be found by iteraive identification methods which implement the input and output transformations iteratively. Methods based on SVD decomposition, QR decomposition and LU decomposition are proposed and compared with the Kuong and Mac Gregor's method. Convergence proofs are given. These SVD and QR mehtods use normal matrices for the transformations which cannot be calculated analytically in general and so they are hoard to apply to dynamic processes, whereas the LU method used simple analyitc transformations and can be directly applied to dynamic processes.

  • PDF

풍력발전 시스템용 회전자 계통연계형 이중여자 권선형 유도발전기 (A Grid-connected to Rotor Type Doubly Fed Induction Generator for Wind Turbine Systems)

  • 유용민;;권병일
    • 대한전기학회:학술대회논문집
    • /
    • 대한전기학회 2011년도 제42회 하계학술대회
    • /
    • pp.800-801
    • /
    • 2011
  • This paper proposes a grid-connected to rotor type doubly fed induction generator (DFIG) in which the rotor winding is connected to the grid instead of the stator winding. The stator size and weight of the proposed grid-connected to rotor type DFIG can be reduced because the proposed type can use rotor core more efficiently compared to the stator type DFIG. In order to verify the size and weight reduction of the proposed type, the loading distribution method (LDM) is utilized. As a design result, the stator outer diameter and weight of the proposed type were decreased. The equivalent circuit analysis and finite element method also performed to verify the design results and to analyze characteristics of the novel DFIG.

  • PDF

Robust Controller Design Method for Systems with Parametric Uncertainties

  • Lee, Jietae;Koo, Doe-Gyoon;Thomas F. Edgar
    • Transactions on Control, Automation and Systems Engineering
    • /
    • 제2권2호
    • /
    • pp.140-148
    • /
    • 2000
  • This paper presents iterative schemes and continuation schemes for designing robust controllers which stabilize dynamic systems having bounded parametric uncertainties. Utilizing results of the cheap control problem, some existence conditions of the robust controller are obtained, which are different from the matching conditions. continuation schemes are used to overcome the divergence problem of iterative schemes. The roust controller design method is extended to nonlinear system ans easily implementable series solution is also obtained. Results are illustrated with simple examples.

  • PDF

A Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions

  • Most, Thomas;Bucher, Christian
    • Structural Engineering and Mechanics
    • /
    • 제21권3호
    • /
    • pp.315-332
    • /
    • 2005
  • The Element-free Galerkin Method has become a very popular tool for the simulation of mechanical problems with moving boundaries. The internally applied Moving Least Squares interpolation uses in general Gaussian or cubic weighting functions and has compact support. Due to the approximative character of this interpolation the obtained shape functions do not fulfill the interpolation conditions, which causes additional numerical effort for the application of the boundary conditions. In this paper a new weighting function is presented, which was designed for meshless shape functions to fulfill these essential conditions with very high accuracy without any additional effort. Furthermore this interpolation gives much more stable results for varying size of the influence radius and for strongly distorted nodal arrangements than existing weighting function types.