• Title/Summary/Keyword: discrete time approximation method

Search Result 51, Processing Time 0.027 seconds

ERROR ESTIMATIES FOR A FREQUENCY-DOMAIN FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS WITH A NEUMANN BOUNDARY CONDITION

  • Lee, Jong-Woo
    • Bulletin of the Korean Mathematical Society
    • /
    • v.35 no.2
    • /
    • pp.345-362
    • /
    • 1998
  • We introduce and anlyze a naturally parallelizable frequency-domain method for parabolic problems with a Neumann boundary condition. After taking the Fourier transformation of given equations in the space-time domain into the space-frequency domain, we solve an indefinite, complex elliptic problem for each frequency. Fourier inversion will then recover the solution of the original problem in the space-time domain. Existence and uniqueness of a solution of the transformed problem corresponding to each frequency is established. Fourier invertibility of the solution in the frequency-domain is also examined. Error estimates for a finite element approximation to solutions fo transformed problems and full error estimates for solving the given problem using a discrete Fourier inverse transform are given.

  • PDF

Robust Intelligent Digital Redesign of Nonlinear System with Parametric Uncertainties (불확실성을 갖는 비선형 시스템의 강인한 지능형 디지털 재설계)

  • Sung, Hwa-Chang;Joo, Young-Hoon;Park, Jin-Bae
    • Journal of the Korean Institute of Intelligent Systems
    • /
    • v.16 no.2
    • /
    • pp.138-143
    • /
    • 2006
  • This paper presents intelligent digital redesign method for hybrid state space fuzzy-model-based controllers. For effectiveness and stabilization of continuous-time uncertain nonlinear systems under discrete-time controller, Takagi-Sugeno(TS) fuzzy model is used to represent the complex system. And global approach design problems viewed as a convex optimization problem that we minimize the error of the norm bounds between nonlinearly interpolated linear operators to be matched. Also, by using the bilinear and inverse bilinear approximation method, we analyzed nonlinear system's uncertain parts more precisely. When a sampling period is sufficiently small, the conversion of a continuous-time structured uncertain nonlinear system to an equivalent discrete-time system have proper reason. Sufficiently conditions for the global state-matching of the digitally controlled system are formulated in terms of linear matrix inequalities (LMIs). Finally, a TS fuzzy model for the chaotic Lorentz system is used as an . example to guarantee the stability and effectiveness of the proposed method.

Ray Effect Analysis Using the Discrete Elements Method in X-Y Geometry (2차원 직각좌표계에서 DEM을 이용한 ray effect의 해석)

  • Choi, Ho-Sin;Kim, Jong-Kyung
    • Journal of Radiation Protection and Research
    • /
    • v.17 no.1
    • /
    • pp.43-56
    • /
    • 1992
  • As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.

  • PDF

Prediction of initiation time of corrosion in RC using meshless methods

  • Yao, Ling;Zhang, Lingling;Zhang, Ling;Li, Xiaolu
    • Computers and Concrete
    • /
    • v.16 no.5
    • /
    • pp.669-682
    • /
    • 2015
  • Degradation of reinforced concrete (RC) structures due to chloride penetration followed by reinforcement corrosion has been a serious problem in civil engineering for many years. The numerical simulation methods at present are mainly finite element method (FEM) and finite difference method (FDM), which are based on mesh. Mesh generation in engineering takes a long time. In the present article, the numerical solution of chloride transport in concrete is analyzed using radial point interpolation method (RPIM) and element-free Galerkin (EFG). They are all meshless methods. RPIM utilizes radial polynomial basis, whereas EFG uses the moving least-square approximation. A Galerkin weak form on global is used to attain the discrete equation, and four different numerical examples are presented. MQ function and appropriate parameters have been proposed in RPIM. Numerical simulation results are compared with those obtained from the finite element method (FEM) and analytical solutions. Two case of chloride transport in full saturated and unsaturated concrete are analyzed to test the practical applicability and performance of the RPIM and EFG. A good agreement is obtained among RPIM, EFG, and the experimental data. It indicates that RPIM and EFG are reliable meshless methods for prediction of chloride concentration in concrete structures.

Analysis of Transient Scattering from Arbitrarily Shaped Three-Dimensional Conducting Objects Using Combined Field Integral Equation (결합 적분방정식을 이용한 삼차원 임의형태 도체 구조물의 전자파 지연산란 해석)

  • Jung, Baek-Ho
    • The Transactions of the Korean Institute of Electrical Engineers C
    • /
    • v.51 no.11
    • /
    • pp.551-558
    • /
    • 2002
  • A time-domain combined field integral equation (CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time-domain electric field integral equation (EFIE) with the magnetic field integral equation (MFIE). The time derivative of the magnetic vector potential in EFIE is approximated using a central finite difference approximation and the scalar potential is averaged over time. The time-domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. The incident spectrum of the field may contain frequency components, which correspond to the internal resonance of the structure. For the numerical solution, we consider both the explicit and implicit scheme and use two different kinds of Gaussian pulses, which may contain frequencies corresponding to the internal resonance. Numerical results for the EFIE, MFIE, and CFIE are presented and compared with those obtained from the inverse discrete Fourier transform (IDFT) of the frequency-domain CFIE solution.

Bayesian Multiple Change-Point Estimation and Segmentation

  • Kim, Jaehee;Cheon, Sooyoung
    • Communications for Statistical Applications and Methods
    • /
    • v.20 no.6
    • /
    • pp.439-454
    • /
    • 2013
  • This study presents a Bayesian multiple change-point detection approach to segment and classify the observations that no longer come from an initial population after a certain time. Inferences are based on the multiple change-points in a sequence of random variables where the probability distribution changes. Bayesian multiple change-point estimation is classifies each observation into a segment. We use a truncated Poisson distribution for the number of change-points and conjugate prior for the exponential family distributions. The Bayesian method can lead the unsupervised classification of discrete, continuous variables and multivariate vectors based on latent class models; therefore, the solution for change-points corresponds to the stochastic partitions of observed data. We demonstrate segmentation with real data.

Instability of (Heterogeneous) Euler beam: Deterministic vs. stochastic reduced model approach

  • Ibrahimbegovic, Adnan;Mejia-Nava, Rosa Adela;Hajdo, Emina;Limnios, Nikolaos
    • Coupled systems mechanics
    • /
    • v.11 no.2
    • /
    • pp.167-198
    • /
    • 2022
  • In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.

Analysis of Stress Concentration Problems Using Moving Least Squares Finite Difference Method(I) : Formulation for Solid Mechanics Problem (이동최소제곱 유한차분법을 이용한 응력집중문제 해석(I) : 고체문제의 정식화)

  • Yoon, Young-Cheol;Kim, Hyo-Jin;Kim, Dong-Jo;Liu, Wing Kam;Belytschko, Ted;Lee, Sang-Ho
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.20 no.4
    • /
    • pp.493-499
    • /
    • 2007
  • The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

Robust Digital Redesign for Observer-based System (관측기 기반 시스템에 대한 강인 디지털 재설계)

  • Sung, Hwa-Chang;Joo, Young-Hoon;Park, Jin-Bae
    • Journal of the Korean Institute of Intelligent Systems
    • /
    • v.17 no.3
    • /
    • pp.285-290
    • /
    • 2007
  • In this paper, we presents robust digital redesign (DR) method for observer-based linear time-invariant (LTI) system. The term of DR involves converting an analog controller into an equivalent digital one by considering two condition: state-matching and stability. The design problems viewed as a convex optimization problem that we minimize the error of the norm bounds between interpolated linear operators to be matched. Also, by using the bilinear and inverse bilinear approximation method, we analyzed the uncertain parts of given observer-based system more precisely, When a sampling period is sufficiently small, the conversion of a analog structured uncertain system to an equivalent discrete-time system have proper reason. Sufficiently conditions for the state-matching of the digitally controlled system are formulated in terms of linear matrix inequalities (LMIs).

Development of High Erosion Resistant Fe-based Alloy for Continuous Hot Dipping Line (연속용융 도금라인 용 고내침식 Fe계 합금 개발)

  • Baek, Min-Sook;Kim, Yong-Cheol;Baek, Kyeong-Cheol;Kwak, Joon-Seop;Yoon, Dong-Joo
    • Journal of the Korean institute of surface engineering
    • /
    • v.53 no.3
    • /
    • pp.95-103
    • /
    • 2020
  • In this study, the material used in the hot dip galvanizing equipment was poorly corrosion-resistant, so it was performed to solve the cost and time problems caused by equipment replacement. The theoretical calculation was performed using the DV-Xα method(Discrete Variational Local-density approximation method). The alloy (STS4XX series) of the equipment currently used has a martensite phase. Therefore, the theoretical calculation was performed by applying P4 / mmm, which is a tetragonal structure. The new alloy was chosen by designing theoretical values close to existing materials. Considering elements that contribute to corrosion, most have high prices. Therefore, the design was completed by adjusting the content using only the components of the reference material in the theoretical design. The final design alloys were chosen as D6 and D9. Designed D6 and D9 were dissolved and prepared using an induction furnace. After the heat treatment process was completed, the corrosion rate of the alloys was confirmed by using the potentiodynamic polarization test. The surface of the prepared alloys were processed horizontally and then polished to # 1200 using sand paper to perform potentiodynamic polarization test. Domestic products: 4.735 mpy (mils / year), D6: 0.9166 mpy, D9: 0.3372 mpy, alloys designed than domestic products had a lower corrosion rate. Therefore, the designed alloy was expected to have better erosion resistance.