Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

Predictor-Corrector Primal-Dual Interior Point Method for Solving Economic Dispatch Problems: A Postoptimization Analysis

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

Short Lyapunov time: a method for identifying confined chaos

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

NUMERICAL SIMULATION of CONVECTION-DIFFUSION PROBLEMS BY THE CONTROL-VOLUME-BASED FINITE-ELEMENT METHOD

This work presents a numerical study of the tri-dimensional convection-diffusion equation by the control-volume-based on finite-element method using quadratic hexahedral elements. Considering that the equation governing this problem in its main variable may represent several properties, including temperature, turbulent kinetic energy, viscous dissipation rate of the turbulent kinetic energy, specific dissipation rate of the turbulent kinetic energy, or even the concentration of a contaminant in a given medium, among others, the wide applicability of this problem is thus evidenced. Three cases of temperature distributions will be studied specifically in this work, in addition to one case of pollutant dispersion upon analysis of the concentration of a contaminant in a fixed flow point. Some comparisons will be carried out against works found in the open literature, while others will be done according to each phenomenon characteristics.

A SIMPLE NUMERICAL-METHOD FOR HYDROSTATIC INCOMPRESSIBLE MODELS WITH RIGID LIDS

A simple and easily implemented method is developed to keep the vertical velocity equal to zero at the bottom and top of hydrostatic incompressible numerical models. The pressure is computed at the top by correcting its value given in the previous time step so that the vertical integral of the horizontal divergence is zero at each column. Numerical experiments that exhibit small time variations of pressure at the top are able to simplify the algorithm and save computer time. Numerical simulations illustrate the method effectiveness for a horizontal deformation-induced frontogenesis.

A combined wavelet-element free Galerkin method for numerical calculations of electromagnetic fields

A combined wavelet-element free Galerkin (EFG) method is proposed for solving electromagnetic EM) field problems. The bridging scales are used to preserve the consistency and linear independence properties of the entire bases. A detailed description of the development of the discrete model and its numerical implementations is given to facilitate the reader to. understand the proposed algorithm. A numerical example to validate the proposed method is also reported.

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.

A new implementation of population based incremental learning method for optimizations in electromagnetics

To enhance the global search ability of population based incremental learning (PBIL) methods, it is proposed that multiple probability vectors are to be included on available PBIL algorithms. The strategy for updating those probability vectors and the negative learning and mutation operators are thus re-defined correspondingly. Moreover, to strike the best tradeoff between exploration and exploitation searches, an adaptive updating strategy for the learning rate is designed. Numerical examples are reported to demonstrate the pros and cons of the newly implemented algorithm.

Evaluation of in vitro gas production of roughages and their mixtures using the curves subtraction method

The in vitro gas production of four single roughages and their paired combinations (1:1 on dry matter basis) were evaluated. Two roughage samples (100 mg) per treatment were fermented with ruminal fluid during a 48 h incubation period. Total 48 h gas volumes of fermentation dry matter (DM), neutral detergent fiber (NDF) and soluble compounds in neutral detergent (NDS) were for sugarcane = 16.8, 11.2, 6.9 mL; sugarcane + corn silage = 20.1, 12.6, 9.1 mL; sugarcane + 60-day elephantgrass = 16.5, 17.6 mL; sugarcane + 180-day elephantgrass = 13.8, 8.2, 5.9 mL; corn silage = 18.8, 16.8, 4.7 mL; corn silage + 60-day elephantgrass = 16.3, 15.4, 2.4 mL; corn silage + 180-day elephantgrass = 16.1, 11.8, 4.2 mL; 60-day elephantgrass = 16.9, 19.0 mL and 180-day elephantgrass = fermented 10.7, 12.2 mL, respectively. The NDS gas production was not possible to estimate for sugarcane + 60-day elephantgrass, 60-day elephantgrass and 180-day elephantgrass. The present data shows that the curves subtraction method can be an option to evaluate the contribution of the soluble fractions in roughages to digestion kinetics. However, this method underestimates the NDS gas contribution when roughages are low in crude protein and soluble carbohydrates. It is advisable to directly apply the two-compartmental mathematical model to the digestion curves for roughage DM, when determining the NDS gas volume and the digestion rate. This method is more straightforward and accurate when compared to the curve subtraction method. Non-structural carbohydrates combined with fiber and protein promoted a positive associative effect in sugarcane + corn silage (50:50) mixture. Therefore, it can be concluded that the soluble fraction of roughages greatly contributes to gas production. (C) 2004 Elsevier B.V. All rights reserved.

A new implementation of Population Based Incremental Learning method for optimization studies in electromagnetics

To enhance the global search ability of Population Based Incremental Learning (PBIL) methods, It Is proposed that multiple probability vectors are to be Included on available PBIL algorithms. As a result, the strategy for updating those probability vectors and the negative learning and mutation operators are redefined as reported. Numerical examples are reported to demonstrate the pros and cons of the newly Implemented algorithm. ©2006 IEEE.

Evaluation of a photographic method to measure dental angulation

Objetivo: analisar a confiabilidade e reprodutibilidade de um método simplificado para análise da angulação dentária que faz uso de fotografias digitalizadas de modelos de gesso. Métodos: foram realizadas fotografias digitalizadas e padronizadas de modelos de gesso, posteriormente transportadas para um programa gráfico de leitura de ângulos, para a obtenção das medidas. Tais procedimentos foram repetidos para avaliação do erro do método casual e para a análise da reprodutibilidade por meio da Correlação Intraclasse. A amostra constituiu-se de 12 indivíduos com dentição permanente completa e não tratados ortodonticamente, sendo seis do sexo masculino e seis do feminino. As análises foram feitas bilateralmente, gerando 24 medidas. Resultados: o erro casual mostrou uma variação de 0,77 a 2,55º para a angulação dos dentes. A análise estatística revelou que o método apresenta uma excelente reprodutibilidade (r = 0,65 - 0,91; p < 0,0001) para todos os dentes, exceto para os pré-molares superiores, mas ainda assim estatisticamente significativa (p < 0,001). Conclusão: o método proposto apresenta confiabilidade suficiente para justificar seu uso no desenvolvimento de pesquisas científicas, bem como na prática clínica.

3D stability orbits close to 433 Eros using an effective polyhedral model method

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Complex modal analysis of a slender vertical rotor by finite elements method

In this paper, natural frequencies were analyzed (axial, torsional and flexural) and frequency response of a vertical rotor with a hard disk at the edge through the classical modal and complex analysis. The equation that rules the movement was obtained through the Lagrangian formulation. The model considered the effects of bending, torsion and axial deformation of the shaft, besides the gravitational and gyroscopic effects. The finite element method was used to discretize the structure into hollow cylindrical elements with 12 degrees of freedom. Mass, stiffness and gyroscopic matrices were explained consistently. The classical modal analysis, usually applied to stationary structures, does not consider an important characteristic of rotating machinery which are the methods of forward and backward whirl. Initially, through the traditional modal analysis, axial and torsional natural frequencies were obtained in a static shaft, since they do not suffer the influence of gyroscopic effects. Later research was performed by complex modal analysis. This type of tool, based on the use of complex coordinates to describe the dynamic behavior of rotating shaft, allows the decomposition of the system in two submodes, backward and forward. Thus, it is possible to clearly visualize that the orbit and direction of the precessional motion around the line of the rotating shaft is not deformed. A finite element program was developed using MATLAB (TM) and numerical simulations were performed to validate this model. Natural frequencies and directional frequency forced response (dFRF) were obtained using the complex modal analysis for a simple vertical rotor and also for a typical drill string used in the construction of oil wells.

Interior-exterior point method with global convergence strategy for solving the reactive optimal power flow problem

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Interferometric measurements of nanometric displacements in a Piezoelectric Flextensional Actuator by using the new J1...J5 method

In this work, nanometric displacement amplitudes of a Piezoelectric Flextensional Actuator (PFA) designed using the topology optimization technique and operating in its linear range are measured by using a homodyne Michelson interferometer. A new improved version of the J1...J4 method for optical phase measurements, named J1...J5 method, is presented, which is of easier implementation than the original one. This is a passive phase detection scheme, unaffected by signal fading, source instabilities and changes in visibility. Experimental results using this improvement were compared with those obtained by using the J1... J4, J1...J6(pos) and J1...J 6(neg) methods, concluding that the dynamic range is increased while maintaining the sensitivity. Analysis based on the 1/f voltage noise and random fading show the new method is more stable to phase drift than all those methods. © 2012 IEEE.