Acessibilidade / Reportar erro
Computational & Applied Mathematics, Volume: 30, Número: 2, Publicado: 2011
  • A study of axi-symmetric waves through an isotropic thermoelastic diffusive medium

    Kumar, Rajneesh; Panchal, Meenakshi

    Resumo em Inglês:

    The paper deals with the propagation of axial symmetric cylindrical surface waves in a cylindrical bore through a homogeneous isotropic thermoelastic diffusive medium of infinite extent. The three theories of thermoelasticity namely, Coupled Thermoelasticity (CT), Lord and Shulman (L-S) and Green and Lindsay (G-L) are used to study the problem. The frequency equations, connecting the phase velocity with wave number, radius of bore and other material parameters, for empty and liquid filled bore are derived. The numerical results obtained have been illustrated graphically to understand the behaviour of phase velocity and attenuation coefficient versus wave number of a wave. A particular case of interest has also been deduced from the present investigation.
  • A semi-analytical computation of the Kelvin kernel for potential flows with a free surface

    D'elía, Jorge; Battaglia, Laura; Storti, Mario

    Resumo em Inglês:

    A semi-analytical computation of the three dimensional Green function for seakeeping flow problems is proposed. A potential flow model is assumed with an harmonic dependence on time and a linearized free surface boundary condition. The multiplicative Green function is expressed as the product of a time part and a spatial one. The spatial part is known as the Kelvin kernel, which is the sum of two Rankine sources and a wave-like kernel, being the last one written using the Haskind-Havelock representation. Numerical efficiency is improved by an analytical integration of the two Rankine kernels and the use of a singularity subtractive technique for the Haskind-Havelock integral, where a globally adaptive quadrature is performed for the regular part and an analytic integration is used for the singular one. The proposed computation is employed in a low order panel method with flat triangular elements. As a numerical example, an oscillating floating unit hemisphere in heave and surge modes is considered, where analytical and semi-analytical solutions are taken as a reference.
  • Comparing stochastic optimization methods to solve the medium-term operation planning problem

    Gonçalves, Raphael E. C.; Finardi, Erlon C.; Silva, Edson L. da; Santos, Marcelo L. L. dos

    Resumo em Inglês:

    The Medium-Term Operation Planning (MTOP) of hydrothermal systems aims to define the generation for each power plant, minimizing the expected operating cost over the planning horizon. Mathematically, this task can be characterized as a linear, stochastic, large-scale problem which requires the application of suitable optimization tools. To solve this problem, this paper proposes to use the Nested Decomposition, frequently used to solve similar problems (as in Brazilian case), and Progressive Hedging, an alternative method, which has interesting features that make it promising to address this problem. To make a comparative analysis between these two methods with respect to the quality of the solution and the computational burden, a benchmark is established, which is obtained by solving a single Linear Programming problem (the Deterministic Equivalent Problem). An application considering a hydrothermal system is carried out.
  • A family of uniformly accurate order Lobatto-Runge-Kutta collocation methods

    Yakubu, D. G.; Manjak, N. H.; Buba, S. S.; Maksha, A. I.

    Resumo em Inglês:

    We consider the construction of an interpolant for use with Lobatto-Runge-Kutta collocation methods. The main aim is to derive single symmetric continuous solution(interpolant) for uniform accuracy at the step points as well as at the off-step points whose uniform order six everywhere in the interval of consideration. We evaluate the continuous scheme at different off-step points to obtain multi-hybrid schemes which if desired can be solved simultaneously for dense approximations. The multi-hybrid schemes obtained were converted to Lobatto-Runge-Kutta collocation methods for accurate solution of initial value problems. The unique feature of the paper is the idea of using all the set of off-step collocation points as additional interpolation points while symmetry is retained naturally by integration identities as equal areas under the various segments of the solution graph over the interval of consideration. We show two possible ways of implementing the interpolant to achieve the aim and compare them on some numerical examples.
  • Initial values for Riccati ODEs from variational PDEs

    Costanza, Vicente; Rivadeneira, Pablo S.

    Resumo em Inglês:

    The recently discovered variational PDEs (partial differential equations) for finding missing boundary conditions in Hamilton equations of optimal control are applied to the extended-space transformation of time-variant linear-quadratic regulator (LQR) problems. These problems become autonomous but with nonlinear dynamics and costs. The numerical solutions to the PDEs are checked against the analytical solutions to the original LQR problem. This is the first validation of the PDEs in the literature for a nonlinear context. It is also found that the initial value of the Riccati matrix can be obtained from the spatial derivative of the Hamiltonian flow, which satisfies the variational equation. This last result has practical implications when implementing two-degrees-of freedom control strategies for nonlinear systems with generalized costs.
  • Encoding through generalized polynomial codes

    Shah, T.; Khan, A.; Andrade, A. A.

    Resumo em Inglês:

    This paper introduces novel constructions of cyclic codes using semigroup rings instead of polynomial rings. These constructions are applied to define and investigate the BCH, alternant, Goppa, and Srivastava codes. This makes it possible to improve several recent results due to Andrade and Palazzo [1].
  • Solving of time varying quadratic optimal control problems by using Bézier control points

    Gachpazan, Mortaza

    Resumo em Inglês:

    In this paper, linear quadratic optial control probles are solved by applying least square method based on Bézier control points. We divide the interval which includes t, into k subintervals and approximate the trajectory and control functions by Bézier curves. We have chosen the Bézier curves as piacewise polynomials of degree three, and determined Bézier curves on any subinterval by four control points. By using least square ethod, e introduce an optimization problem and compute the control points by solving this optimization problem. Numerical experiments are presented to illustrate the proposed method.
  • Dynamical behavior of a pest management model with impulsive effect and nonlinear incidence rate

    Wang, Xia; Guo, Zhen; Song, Xinyu

    Resumo em Inglês:

    In this paper, we consider the pest management model with spraying microbial pesticide and releasing the infected pests, and the infected pests have the function similar to the microbial pesticide and can infect the healthy pests, further weaken or disable their prey function till death. By using the Floquet theory for impulsive differential equations, we show that there exists a globally asymptotically stable pest eradication periodic solution when the impulsive period τ < τmax, we further prove that the system is uniformly permanent if the impulsive period τ > τmax. Finally, by means of numerical simulation, we showthatwith the increaseof impulsive period, the system displays complicated behaviors.
  • Relaxation approaches to the optimal control of the Euler equations

    Ngnotchouye, Jean Medard T.; Herty, Michael; Steffensen, Sonja; Banda, Mapundi K.

    Resumo em Inglês:

    The treatment of control problems governed by systems of conservation laws poses serious challenges for analysis and numerical simulations. This is due mainly to shock waves that occur in the solution of nonlinear systems of conservation laws. In this article, the problem of the control of Euler flows in gas dynamics is considered. Numerically, two semi-linear approximations of the Euler equations are compared for the purpose of a gradient-based algorithm for optimization. One is the Lattice-Boltzmann method in one spatial dimension and five velocities (D1Q5 model) and the other is the relaxation method. An adjoint method is used. Good results are obtained even in the case where the solution contains discontinuities such as shock waves or contact discontinuities.
  • A global linearization approach to solve nonlinear nonsmooth constrained programming problems

    Vaziri, A.M.; Kamyad, A.V.; Jajarmi, A.; Effati, S.

    Resumo em Inglês:

    In this paper we introduce a new approach to solve constrained nonlinear non-smooth prograing probles ith any desirable accuracy even hen the objective function is a non-smooth one. In this approach for any given desirable accuracy, all the nonlinear functions of original problem (in objective function and in constraints) are approximated by a piecewise linear functions. We then represent an efficient algorithm to find the global solution of the later problem. The obtained solution has desirable accuracy and the error is completely controllable. One of the main advantages of our approach is that the approach can be extended to problems with non-smooth structure by introducing a novel definition of Global Weak Differentiation in the sense of L1 norm. Finally some numerical examples are given to show the efficiency of the proposed approach to solve approximately constraints nonlinear non-smooth programming problems.
  • Accuracy of analytical-numerical solutions of the Michaelis-Menten equation

    Gonzalez-Parra, Gilberto; Acedo, Luis; Arenas, Abraham

    Resumo em Inglês:

    It is the aim of this paper to investigate a suitable approach to compute solutions of the powerful Michaelis-Menten enzyme reaction equation with less computational effort. We obtain analytical-numerical solutions using piecewise finite series by means of the differential transformation method, DTM. In addition, we compute a global analytical solution by a modal series expansion. The Michaelis-Menten equation considered here describes the rate of depletion of the substrate of interest and in general is a powerful approach to describe enzyme processes. A comparison of the numerical solutions using DTM, Adomian decomposition and Runge-Kutta methods is presented. The numerical results show that the DTM is accurate, easy to apply and the obtained solutions retain the positivity property of the continuous model. It is concluded that the analytic form of the DTM and global modal series solutions are accurate, and require less computational effort than other approaches thus making them more convenient.
  • A universal constant for semistable limit cycles

    Artés, Joan C.; Llibre, Jaume; Teixeira, Marco Antonio

    Resumo em Inglês:

    We consider one-parameter families of 2-dimensional vector fields Xµ having in a convenient region R a semistable limit cycle of multiplicity 2m when µ = 0, no limit cycles if µ < 0, and two limit cycles one stable and the other unstable if µ > 0. We show, analytically for some particular families and numerically for others, that associated to the semistable limit cycle and for positive integers n sufficiently large there is a power law in the parameter µ of the form µn ≈ Cnα< 0 with C, α ∈ R, such that the orbit of Xµn through a point of p ∈ R reaches the position of the semistable limit cycle of X0 after given n turns. The exponent α of this power law depends only on the multiplicity of the semistable limit cycle, and is independent of the initial point p ∈ R and of the family Xµ. In fact α = -2m/(2m - 1). Moreover the constant C is independent of the initial point p ∈ R, but it depends on the family Xµ and on the multiplicity 2m of the limit cycle Γ.
Sociedade Brasileira de Matemática Aplicada e Computacional Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC, Rua Maestro João Seppe, nº. 900 , 16º. andar - Sala 163, 13561-120 São Carlos - SP Brasil, Tel./Fax: 55 16 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br