• Title/Summary/Keyword: Galerkin approximation

Search Result 90, Processing Time 0.034 seconds

Parallel Robust $H_{\infty}$ Control for Weakly Coupled Bilinear Systems with Parameter Uncertainties Using Successive Galerkin Approximation

  • Kim, Young-Joong;Lim, Myo-Taeg
    • International Journal of Control, Automation, and Systems
    • /
    • v.4 no.6
    • /
    • pp.689-696
    • /
    • 2006
  • This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of weakly coupled bilinear systems with time-varying parameter uncertainties and exogenous disturbance using the successive Galerkin approximation(SGA). By using weak coupling theory, the robust $H_{\infty}$ control can be obtained from two reduced-order robust $H_{\infty}$ control problems in parallel. The $H_{\infty}$ control theory guarantees robust closed-loop performance but the resulting problem is difficult to solve for uncertain bilinear systems. In order to overcome the difficulties inherent in the $H_{\infty}$ control problem, two $H_{\infty}$ control laws are constructed in terms of the approximated solution to two independent Hamilton-Jacobi-Isaac equations using the SGA method. One of the purposes of this paper is to design a closed-loop parallel robust $H_{\infty}$ control law for the weakly coupled bilinear systems with parameter uncertainties using the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

$H_{\infty}$ Composite Control for Singularly Perturbed Nonlinear Systems via Successive Galerkin Approximation

  • Kim, Young-Joong;Kim, Beom-Soo;Shin, Eun-Chul;Yoo, Ji-Yoon;Lim, Myo-Taeg
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 2004.08a
    • /
    • pp.407-412
    • /
    • 2004
  • This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of singularly perturbed nonlinear systems with a exogenous disturbance, using the successive Galerkin approximation(SGA). The singularly perturbed nonlinear system is decomposed into two subsystems of a slow-time scale and a fast-time scale via singular perturbation theory, and two $H_{\infty}$ control laws are obtained to each subsystem by using the SGA method. The composite control law that consists of two $H_{\infty}$ control laws of each subsystem is designed. One of the purposes of this paper is to design the closed-loop $H_{\infty}$ composite control law for the singularly perturbed nonlinear systems via the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

  • PDF

ON THE GALERKIN-WAVELET METHOD FOR HIGHER ORDER DIFFERENTIAL EQUATIONS

  • Fukuda, Naohiro;Kinoshita, Tamotu;Kubo, Takayuki
    • Bulletin of the Korean Mathematical Society
    • /
    • v.50 no.3
    • /
    • pp.963-982
    • /
    • 2013
  • The Galerkin method has been developed mainly for 2nd order differential equations. To get numerical solutions, there are some choices of Riesz bases for the approximation subspace $V_j{\subset}L^2$. In this paper we shall propose a uniform approach to find suitable Riesz bases for higher order differential equations. Especially for the beam equation (4-th order equation), we also report numerical results.

Stability Improved Split-step Parabolic Equation Model

  • Kim, Tae-Hyun;Seong, Woojae
    • The Journal of the Acoustical Society of Korea
    • /
    • v.21 no.3E
    • /
    • pp.105-111
    • /
    • 2002
  • The parabolic equation technique provides an excellent model to describe the wave phenomena when there exists a predominant direction of propagation. The model handles the square root wave number operator in paraxial direction. Realization of the pseudo-differential square root operator is the essential part of the parabolic equation method for its numerical accuracy. The wide-angled approximation of the operator is made based on the Pade series expansion, where the branch line rotation scheme can be combined with the original Pade approximation to stabilize its computational performance for complex modes. The Galerkin integration has been employed to discretize the depth-dependent operator. The benchmark tests involving the half-infinite space, the range independent and dependent environment will validate the implemented numerical model.

FULLY DISCRETE MIXED FINITE ELEMENT METHOD FOR A QUASILINEAR STEFAN PROBLEM WITH A FORCING TERM IN NON-DIVERGENCE FORM

  • Lee, H.Y.;Ohm, M.R.;Shin, J.Y.
    • Journal of applied mathematics & informatics
    • /
    • v.24 no.1_2
    • /
    • pp.191-207
    • /
    • 2007
  • Based on a mixed Galerkin approximation, we construct the fully discrete approximations of $U_y$ as well as U to a single-phase quasilinear Stefan problem with a forcing term in non-divergence form. We prove the optimal convergence of approximation to the solution {U, S} and the superconvergence of approximation to $U_y$.

Energy extraction from the motion of an oscillating water column

  • Wang, Hao;Falzarano, Jeffrey M.
    • Ocean Systems Engineering
    • /
    • v.3 no.4
    • /
    • pp.327-348
    • /
    • 2013
  • An Oscillating Water Column (OWC) is a relatively practical and convenient device that converts wave energy to a usable form, which is electricity. The OWC is kept inside a fixed truncated vertical cylinder, which is a hollow structure with one open end submerged in the water and with an air turbine at the top. This research adopts potential theory and Galerkin methods to solve the fluid motion inside the OWC. Using an air-water interaction model, OWC design for energy extraction from regular wave is also explored. The hydrodynamic coefficients of the scattering and radiation potentials are solved for using the Galerkin approximation. The numerical results for the free surface elevation have been verified by a series of experiments conducted in the University of New Orleans towing tank. The effect of varying geometric parameters on the response amplitude operator (RAO) of the OWC is studied and modification of the equation for evaluating the natural frequency of the OWC is made. Using the model of air-water interaction under certain wave parameters and OWC geometric parameters, a computer program is developed to calculate the energy output from the system.

Development of Galerkin Finite Element Method Three-dimensional Computational Code for the Multigroup Neutron Diffusion Equation with Unstructured Tetrahedron Elements

  • Hosseini, Seyed Abolfazl
    • Nuclear Engineering and Technology
    • /
    • v.48 no.1
    • /
    • pp.43-54
    • /
    • 2016
  • In the present paper, development of the three-dimensional (3D) computational code based on Galerkin finite element method (GFEM) for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the GFEM-3D computational code through a developed interface. The forward/adjoint multiplication factor, forward/adjoint flux distribution, and power distribution in the reactor core are calculated using the power iteration method. Criticality calculations are benchmarked against the valid solution of the neutron diffusion equation for International Atomic Energy Agency (IAEA)-3D and Water-Water Energetic Reactor (VVER)-1000 reactor cores. In addition, validation of the calculations against the $P_1$ approximation of the transport theory is investigated in relation to the liquid metal fast breeder reactor benchmark problem. The neutron fixed source calculations are benchmarked through a comparison with the results obtained from similar computational codes. Finally, an analysis of the sensitivity of calculations to the number of elements is performed.

Sensitivity Analysis of the Galerkin Finite Element Method Neutron Diffusion Solver to the Shape of the Elements

  • Hosseini, Seyed Abolfazl
    • Nuclear Engineering and Technology
    • /
    • v.49 no.1
    • /
    • pp.29-42
    • /
    • 2017
  • The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

PARAMETRIC INVESTIGATIONS ON THE DOUBLE DIFFUSIVE CONVECTION IN TRIANGULAR CAVITY

  • Kwon, SunJoo;Oh, SeYoung;Yun, Jae Heon;Chung, Sei-Young
    • Journal of the Chungcheong Mathematical Society
    • /
    • v.20 no.4
    • /
    • pp.419-432
    • /
    • 2007
  • Double-diffusive convection inside a triangular porous cavity is studied numerically. Galerkin finite element method is adopted to derive the discrete form of the governing differential equations. The first-order backward Euler scheme is used for temporal discretization with the second-order Adams-Bashforth scheme for the convection terms in the energy and species conservation equations. The Boussinesq-Oberbeck approximation is used to calculate the density dependence on the temperature and concentration fields. A parametric study is performed with the Lewis number, the Rayleigh number, the buoyancy ratio, and the shape of the triangle. The effect of gravity orientation is considered also. Results obtained include the flow, temperature, and concentration fields. The differences induced by varying physical parameters are analyzed and discussed. It is found that the heat transfer rate is sensitive to the shape of the triangles. For the given geometries, buoyancy ratio and Rayleigh numbers are the dominating parameters controlling the heat transfer.

  • PDF