Integration Methods Used in Numerical Simulations of Electromagnetic Transients

This paper made an analysis of some numerical integration methods that can be used in electromagnetic transient simulations. Among the existing methods, we analyzed the trapezoidal integration method (or Heun formula), Simpson's Rule and Runge-Kutta. These methods were used in simulations of electromagnetic transients in power systems, resulting from switching operations and maneuvers that occur in transmission lines. Analyzed the characteristics such as accuracy, computation time and robustness of the methods of integration.

Comparison of methodologies for degree-day estimation using numerical methods

O desenvolvimento de projetos relacionados ao desempenho de diversas culturas tem recebido aperfeiçoamento cada vez maior, incorporado a modelos matemáticos sendo indispensável à utilização de equações cada vez mais consistentes que possibilitem previsão e maior aproximação do comportamento real, diminuindo o erro na obtenção das estimativas. Entre as operações unitárias que demandam maior estudo estão aquelas relacionadas com o crescimento da cultura, caracterizadas pela temperatura ideal para o acréscimo de matéria seca. Pelo amplo uso dos métodos matemáticos na representação, análise e obtenção de estimativas de graus-dia, juntamente com a grande importância que a cultura da cana-de-açúcar tem para a economia brasileira, foi realizada uma avaliação dos modelos matemáticos comumente usados e dos métodos numéricos de integração na estimativa da disponibilidade de graus-dia para essa cultura, na região de Botucatu, Estado de São Paulo. Os modelos de integração, com discretização de 6 em 6 h, apresentaram resultados satisfatórios na estimativa de graus-dia. As metodologias tradicionais apresentaram desempenhos satisfatórios quanto à estimativa de grausdia com base na curva de temperatura horária para cada dia e para os agrupamentos de três, sete, 15 e 30 dias. Pelo método numérico de integração, a região de Botucatu, Estado de São Paulo, apresentou disponibilidade térmica anual média de 1.070,6 GD para a cultura da cana-de-açúcar.

An analytical and numerical study of plane change maneuvers using aerodynamic force

The goal of the present work is to analyze space missions that use the terrestrial atmosphere to accomplish orbital maneuvers that involve a plane change. A set of analytical solutions is presented for the variation of the orbital elements due to a single passage through the atmosphere, assuming that the interval the spacecraft travels through the atmosphere is not too large. The study considers both the lift influence on the spacecraft orbit as well as drag. The final equations are tested with numerical integration and can be considered in accordance with the numerical results whenever the perigee height is larger than a critical value. Next, a numerical study of the ratio between the velocity increment required to correct the semimajor axis decay due to the atmospheric passage and the velocity variation required to obtain the change in the inclination is also presented. This analysis can be used to decide if a maneuver passing through the atmosphere can decrease the fuel consumption of the mission and, in the cases where this technique can be used, if a multiple passage is more efficient than a single passage.

Artificial satellites orbits in 2:1 resonance: GPS constellation

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

GPS Satellites Orbits: Resonance

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The Orbital Dynamics of Synchronous Satellites: Irregular Motions in the 2:1 Resonance

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

Effective potential of two coupled binary matter wave bright solitons with spatially modulated nonlinearity

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

Correlation times in stochastic equations with delayed feedback and multiplicative noise

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Approximations for the generalized temperature integral: a method based on quadrature rules

The generalized temperature integral I(m, x) appears in non-isothermal kinetic analysis when the frequency factor depends on the temperature. A procedure based on Gaussian quadrature to obtain analytical approximations for the integral I(m, x) was proposed. The results showed good agreement between the obtained approximation values and those obtained by numerical integration. Unless other approximations found in literature, the methodology presented in this paper can be easily generalized in order to obtain approximations with the maximum of accurate.

GCHF basis sets and their application in the electronic structure study of PrMnO3

The approach called generator coordinate Hartree-Fock (GCHF) method is used in the selection of Gaussian basis set [25s18p for O ((3)p), 31s21p14d for Mn (S-6), and 33s22p16d9f for Pr ((4)J)] for atoms. The role of the weight functions in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are contracted to (25s18p/9s5p), (31s21p14d/9s6p4d), and (33s22pl6d9f118sl2p5d3f) by segmented contraction scheme of Dunning and they are utilized in calculations of Restricted-Open-HF (ROHF) Total and Orbital energies of the (MnO+1)-Mn-3 and (PrO+1)-Pr-1 fragments, to evaluate their quality in molecular studies. The addition of one d polarization function in the contracted (9s5p) basis set for O(P-3) atom and their application with the contracted (9s6p4d), (18s21p5d3f) basis sets for Mn (S-6) and Pr-Pr ((4)j) atoms lead to the electronic structure study of PrMnO3. The dipole moment, the total energy, and total atomic charges properties were calculated and were carried out at ROHF level with the [PrMnO3](2) fragment. The calculated values show that PrMnO3 does not present piezoelectric properties. (C) 2004 Elsevier B.V. All rights reserved.

Ab initio study of high tridymite by the formalism generator coordinate Hartree-Fock

The Generator Coordinate Hartree-Fock (GCHF) Method is applied to generate extended 14s 8p and 17s 11p Gaussian basis sets for the atoms O and Si, respectively. The role of the weight functions in the assessment of the numerical integration range of the GCHF is shown. The Gaussian basis sets are contracted to [6s4p] O atom and [8s5p] Si atom by the Dunning's segmented contraction scheme. To evaluate the quality of our contracted [6s4p] and [8s5p] bases in molecular calculations we accomplish calculations of total and orbital energies in the Hartree-Fock-Roothaan method for O-2 and SiO molecules. We compare the results obtained with the our (14s 8p) and (17s 11p) bases sets with the of 6-311G basis and with values from the literature. The addition of one d polarization function in the silicon basis and its utilization with the basis for oxygen leads to the calculation of electronic properties and IR Spectrum of high tridymite in space group D-3d. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Gaussian basis sets to the theoretical study of the electronic structure of perovskite (LaMnO3)

The Generator Coordinate Hartree-Fock (GCHF) method is applied to generate extended (20s14p), (30s19p13d), and (31s23p18d) Gaussian basis sets for the 0, Mn, and La atoms, respectively. The role of the weight functions (WFs) in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are then contracted to [5s3p] and [11s6p6d] for 0 and Mn atoms, respectively, and [17s11p7d] for La atom by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaan (HFR) method for (MnO1+)-Mn-5 and (LaO1+)-La-1 fragments. The results obtained with the contracted basis sets are compared with values obtained with the extended basis sets. The addition of one d polarization function in the contracted basis set for 0 atom and its utilization with the contracted basis sets for Mn and La atoms leads to the calculations of dipole moment and total atomic charges of perovskite (LaMnO3). The calculations were performed at the HFR level with the crystal [LaMnO3](2) fragment in space group C-2v the values of dipole moment, total energy, and total atomic charges showed that it is reasonable to believe that LaMnO3 presents behaviour of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

Determination of interfacial co-ion concentration in ionic micelles by chemical trapping: Halide concentration at the interface of sodium dodecyl sulfate micelles

Products from the spontaneous reaction of a long-chain arenediazonium salt, 2,6-dimethyl-4-hexadecylbenzenediazonium tetrafluoroborate(16-ArN2BF4), in aqueous micellar solutions of sodium dodecyl sulfate (SDS)? are used to estimate the local concentration of chloride and bromide ions at the micellar surface. The arenediazonium ion, 16-ArN2+, which is totally bound to the SDS micelle, reacts by rate-determining loss of N-2 to give an aryl cation that traps available nucleophiles, i,e., H2O, Cl-, and Br-, to give stable phenol, 16-ArOH, and halobenzene products, 16-ArCl and 16-ArBr, respectively. Product yields, determined by HPLC, are related to local concentrations using calibration curves obtained from independent standards. The local concentrations determined by this method are consistent with co-ion concentrations calculated, using a cell model, by numerical integration of the Poisson-Boltzmann equation (PBE) taking into account salt-induced micellar growth. The salt dependence of the intel facial concentrations of Cl- and Br- are identical. indicating no specific interactions in the interfacial co-ion compartment. PBE calculations predict that, in micellar SDS, increasing the concentration of a particular halide salt (NaX) at constant concentration of another halide (NaY) should result in an increase in the local concentrations of both co-ions. Using this chemical-trapping method, this prediction was demonstrated experimentally.

Simple generalization of Wisdom's perturbative method

A simple generalization of Wisdom's perturbative method, as originally proposed by Wisdom (1985), is obtained. Any number of resonant cosines can be handled and the method can also accommodate more involved disturbing functions. Averaged trajectories are easily obtained by drawing level curves of the action. Here, the method is first tested for simple models of 3:1 and 2:1 resonant problems. Comparisons with numerical integration and surface-section curves show very good agreements.

Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.