High-order upwinding and the hydraulic jump

This work is concerned with the computation of incompressible axisymmetric and fall three-dimensional free-surface flows. In particular, the circular-hydraulic jump is simulated and compared with approximate analytic solutions. However, the principal thrust of this paper is to provide a real problem as a test bed for comparing the many existing convective approximations. Their performance is compared; SMART, HLPA and VONOS emerge as acceptable upwinding methods for this problem. Copyright (C) 2002 John Wiley Sons, Ltd.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier B.V. All rights reserved.

Mathematical modeling and control of population systems: Applications in biological pest control

The aim of this paper is to apply methods from optimal control theory, and from the theory of dynamic systems to the mathematical modeling of biological pest control. The linear feedback control problem for nonlinear systems has been formulated in order to obtain the optimal pest control strategy only through the introduction of natural enemies. Asymptotic stability of the closed-loop nonlinear Kolmogorov system is guaranteed by means of a Lyapunov function which can clearly be seen to be the solution of the Hamilton-Jacobi-Bellman equation, thus guaranteeing both stability and optimality. Numerical simulations for three possible scenarios of biological pest control based on the Lotka-Volterra models are provided to show the effectiveness of this method. (c) 2007 Elsevier B.V. All rights reserved.

Monitoring the mean vector and the covariance ­matrix of multivariate processes with sample means and sample ranges

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

Short term hydroelectric scheduling combining network flow and interior point approaches

In this paper, short term hydroelectric scheduling is formulated as a network flow optimization model and solved by interior point methods. The primal-dual and predictor-corrector versions of such interior point methods are developed and the resulting matrix structure is explored. This structure leads to very fast iterations since it avoids computation and factorization of impedance matrices. For each time interval, the linear algebra reduces to the solution of two linear systems, either to the number of buses or to the number of independent loops. Either matrix is invariant and can be factored off-line. As a consequence of such matrix manipulations, a linear system which changes at each iteration has to be solved, although its size is reduced to the number of generating units and is not a function of time intervals. These methods were applied to IEEE and Brazilian power systems, and numerical results were obtained using a MATLAB implementation. Both interior point methods proved to be robust and achieved fast convergence for all instances tested. (C) 2004 Elsevier Ltd. All rights reserved.

NATURAL PRODUCT SUBSTITUTION PATTERN AND SKELETAL RECOGNITION METHOD

This work describes a methodology for identification of skeletal types of diterpenes based on data base with 1500 compounds isolated from Asteraceae. One program named BOTOCSYS was built with the codification of the compounds and their botanical sources. An example of identification of a new substance is given.

TRICUT: a program to clip triangle meshes using the rapid and triangle libraries and the visualization toolkit

An efficient technique to cut polygonal meshes as a step in the geometric modeling of topographic and geological data has been developed. In boundary represented models of outcropping strata and faulted horizons polygonal meshes often intersect each other. TRICUT determines the line of intersection and re-triangulates the area of contact. Along this line the mesh is split in two or more parts which can be selected for removal. The user interaction takes place in the 3D-model space. The intersection, selection and removal are under graphic control. The visualization of outcropping geological structures in digital terrain models is improved by determining intersections against a slightly shifted terrain model. Thus, the outcrop line becomes a surface which overlaps the terrain in its initial position. The area of this overlapping surface changes with respect to the strike and dip of the structure, the morphology and the offset. Some applications of TRICUT on different real datasets are shown. TRICUT is implemented in C+ + using the Visualization Toolkit in conjunction with the RAPID and TRIANGLE libraries. The program runs under LINUX and UNIX using the MESA OpenGL library. This work gives an example of solving a complex 3D geometric problem by integrating available robust public domain software. (C) 2002 Elsevier B.V. Ltd. All rights reserved.

Strong and fast phase modulation for quantitative analysis of photorefractive gratings

A new approach for studying photorefractive gratings in two-wave mixing experiments by a phase modulation technique is presented. The introduction of a large-amplitude, high-frequency sinusoidal phase modulation in one of the input beams blurs the interference pattern and provides powerful harmonic signals for accurate measurements of the grating diffraction efficiency eta and the output phase shift rho between the transmitted and diffracted waves. The blurring of the light fringes can be used to suppress the higher spatial harmonics of the grating, allowing a space-charge field with sinusoidal profile to be recorded. Although the presence of such a strong phase modulation affects the beam coupling in a rather complicated way, it is shown that for the special case of equal intensity input beams, the effect of the phase modulation on eta and rho is reduced to a weakening of the coupling strength. The potentialities of the technique are illustrated in a study of refractive-index waves excited by running interference patterns in a Bi12TiO20 crystal. Expressions for the diffraction efficiency and the output phase shift are derived and used to match numerically calculated curves to the experimental data. The theoretical model is supported by the very good data fitting and allows the computation of important material parameters.

Analysing the new Saturnian rings, R/2004 S1 and R/2004 S2

The Cassini-Huygens arrival into the Saturnian system brought a large amount of data about the satellites and rings. Two diffuse rings were found in the region between the A ring and Prometheus. R/2004 S1 is coorbital to Atlas and R/2004 S2 is close to Prometheus. In this work we analysed the closest approach between Prometheus and both rings. As a result we found that the satellite removes particles from R/2004 S2 ring. Long-term numerical simulations showed that some particles can cross the F ring region . The well known region of the F ring, where small satellites are present and particles are being taking from the ring, gains a new insight with the presence of particles from R/2004 S2 ring. The computation of the Lyapunov Characteristic Exponent reveled that the R/2004 S2 ring lies in a chaotic region while R/2004 S1 ring and Atlas are in a stable region. Atlas is responsible for the formation of three regimes in the R/2004 S1 ring, as expected for a satellite embedded in a ring.

On the denominator values and barycentric weights of rational interpolants

We improve upon the method of Zhu and Zhu [A method for directly finding the denominator values of rational interpolants, J. Comput. Appl. Math. 148 (2002) 341-348] for finding the denominator values of rational interpolants, reducing considerably the number of arithmetical operations required for their computation. In a second stage, we determine the points (if existent) which can be discarded from the rational interpolation problem. Furthermore, when the interpolant has a linear denominator, we obtain a formula for the barycentric weights which is simpler than the one found by Berrut and Mittelmann [Matrices for the direct determination of the barycentric weights of rational interpolation, J. Comput. Appl. Math. 78 (1997) 355-370]. Subsequently, we give a necessary and sufficient condition for the rational interpolant to have a pole. (c) 2006 Elsevier B.V. All rights reserved.

X-ray powder data and bond valence of La0.65Sr0.35MnO3 after Rietveld refinement

Powder X-ray diffraction (XRD) data were collected for La0.65Sr0.35MnO3 prepared through an alternative method from a stoichiometric mixture of Mn2O3, La2O3, and SrO2, fired at 1300 degreesC for 16 h. XRD analysis using the Rietveld method was carried out and it was found that manganite has rhombohedral symmetry (space group R(3) over bar c). The lattice parameters are found to be a=5.5032 Angstrom and c=13.3674 Angstrom. The bond valence computation indicates that the initial inclusion of Sr occurs at higher temperature. (C) 2002 International Centre for Diffraction Data.

Using transient and steady state considerations to investigate the mechanism of loss of stability of a dynamical system

This work presents the complete set of features for solutions of a particular non-ideal mechanical system near the fundamental and near to a secondary resonance region. The system comprises a pendulum with a horizontally moving suspension point. Its motion is the result of a non-ideal rotating power source (limited power supply), acting oil the Suspension point through a crank mechanism. Main emphasis is given to the loss of stability, which occurs by a sequence of events, including intermittence and crisis, when the system reaches a chaotic attractor. The system also undergoes a boundary-crisis, which presents a different aspect in the bifurcation diagram due to the non-ideal supposition. (c) 2004 Published by Elsevier B.V.

Integration of polyharmonic functions

The results in this paper are motivated by two analogies. First, m-harmonic functions in R(n) are extensions of the univariate algebraic polynomials of odd degree 2m-1. Second, Gauss' and Pizzetti's mean value formulae are natural multivariate analogues of the rectangular and Taylor's quadrature formulae, respectively. This point of view suggests that some theorems concerning quadrature rules could be generalized to results about integration of polyharmonic functions. This is done for the Tchakaloff-Obrechkoff quadrature formula and for the Gaussian quadrature with two nodes.

Symbolic power, robotting, and surveilling

Symbolic power is discussed with reference to mathematics and formal languages. Two distinctions are crucial for establishing mechanical and formal perspectives: one between appearance and reality, and one between sense and reference. These distinctions include a nomination of what to consider primary and secondary. They establish the grammatical format of a mechanical and a formal world view. Such views become imposed on the domains addressed by means of mathematics and formal languages. Through such impositions symbolic power of mathematics becomes exercised. The idea that mathematics describes as it prioritises is discussed with reference to robotting and surveillance. In general, the symbolic power of mathematics and formal languages is summarised through the observations: that mathematics treats parts and properties as autonomous, that it dismembers what it addresses and destroys the organic unity around things, and that it simplifies things and reduces them to a single feature. But, whatever forms the symbolic power may take, it cannot be evaluated along a single good-bad axis.