GAMS MINLP Solver comparisons and some improvements to the AlphaECP algorithm

I doktorsavhandlingen undersöks förmågan att lösa hos ett antal lösare för optimeringsproblem och ett antal svårigheter med att göra en rättvis lösarjämförelse avslöjas. Dessutom framläggs några förbättringar som utförts på en av lösarna som heter GAMS/AlphaECP. Optimering innebär, i det här sammanhanget, att finna den bästa möjliga lösningen på ett problem. Den undersökta klassen av problem kan karaktäriseras som svårlöst och förekommer inom ett flertal industriområden. Målet har varit att undersöka om det finns en lösare som är universellt snabbare och hittar lösningar med högre kvalitet än någon av de andra lösarna. Det kommersiella optimeringssystemet GAMS (General Algebraic Modeling System) och omfattande problembibliotek har använts för att jämföra lösare. Förbättringarna som presenterats har utförts på GAMS/AlphaECP lösaren som baserar sig på skärplansmetoden Extended Cutting Plane (ECP). ECP-metoden har utvecklats främst av professor Tapio Westerlund på Anläggnings- och systemteknik vid Åbo Akademi.

Influence of data acquisition geometry on soybean spectral response simulated by the prosail model

View angle and directional effects significantly affect reflectance and vegetation indices, especially when daily images collected by large field-of-view (FOV) sensors like the Moderate Resolution Imaging Spectroradiometer (MODIS) are used. In this study, the PROSAIL radiative transfer model was chosen to evaluate the impact of the geometry of data acquisition on soybean reflectance and two vegetation indices (Normalized Difference Vegetation Index - NDVI and Enhanced Vegetation Index -EVI) by varying biochemical and biophysical parameters of the crop. Input values for PROSAIL simulation were based on the literature and were adjusted by the comparison between simulated and real satellite soybean spectra acquired by the MODIS/Terra and hyperspectral Hyperion/Earth Observing-One (EO-1). Results showed that the influence of the view angle and view direction on reflectance was stronger with decreasing leaf area index (LAI) and chlorophyll concentration. Because of the greater dependence on the near-infrared reflectance, the EVI was much more sensitive to viewing geometry than NDVI presenting larger values in the backscattering direction. The contrary was observed for NDVI in the forward scattering direction. In relation to the LAI, NDVI was much more isotropic for closed soybean canopies than for incomplete canopies and a contrary behavior was verified for EVI.

Software for calculation of reservoir active capacity with the sequent-peak algorithm

It is presented a software developed with Delphi programming language to compute the reservoir's annual regulated active storage, based on the sequent-peak algorithm. Mathematical models used for that purpose generally require extended hydrological series. Usually, the analysis of those series is performed with spreadsheets or graphical representations. Based on that, it was developed a software for calculation of reservoir active capacity. An example calculation is shown by 30-years (from 1977 to 2009) monthly mean flow historical data, from Corrente River, located at São Francisco River Basin, Brazil. As an additional tool, an interface was developed to manage water resources, helping to manipulate data and to point out information that it would be of interest to the user. Moreover, with that interface irrigation districts where water consumption is higher can be analyzed as a function of specific seasonal water demands situations. From a practical application, it is possible to conclude that the program provides the calculation originally proposed. It was designed to keep information organized and retrievable at any time, and to show simulation on seasonal water demands throughout the year, contributing with the elements of study concerning reservoir projects. This program, with its functionality, is an important tool for decision making in the water resources management.

Effect of a simulated heat wave in thermal and aerial environment broiler-rearing environment

ABSTRACT Global warming increases the occurrence of events such as extreme heat waves. Research on thermal and air conditions affecting broiler-rearing environment are important to evaluate the animal welfare under extreme heat aiming mitigation measures. This study aimed at evaluating the effect of a simulated heat wave, in a climatic chamber, on the thermal and air environment of 42-day-old broilers. One hundred and sixty broilers were housed and reared for 42 days in a climatic chamber; the animals were divided into eight pens. Heat wave simulation was performed on the 42nd day, the period of great impact and data sampling. The analyzed variables were room and litter temperatures, relative humidity, concentrations of oxygen, carbon monoxide and ammonia at each pen. These variables were assessed each two hours, starting at 8 am, simulating a day heating up to 4 pm, when it is reached the maximum temperature. By the results, we concluded that increasing room temperatures promoted a proportional raise in litter temperatures, contributing to ammonia volatilization. In addition, oxygen concentrations decreased with increasing temperatures; and the carbon monoxide was only observed at temperatures above 27.0 °C, relative humidity higher than 88.4% and litter temperatures superior to 30.3 °C.

Testing Platform Development for Large Scale Solar Integration

Currently, a high penetration level of Distributed Generations (DGs) has been observed in the Danish distribution systems, and even more DGs are foreseen to be present in the upcoming years. How to utilize them for maintaining the security of the power supply under the emergency situations, has been of great interest for study. This master project is intended to develop a control architecture for studying purposes of distribution systems with large scale integration of solar power. As part of the EcoGrid EU Smart Grid project, it focuses on the system modelling and simulation of a Danish representative LV network located in Bornholm island. Regarding the control architecture, two types of reactive control techniques are implemented and compare. In addition, a network voltage control based on a tap changer transformer is tested. The optimized results after applying a genetic algorithm to five typical Danish domestic loads are lower power losses and voltage deviation using Q(U) control, specially with large consumptions. Finally, a communication and information exchange system is developed with the objective of regulating the reactive power and thereby, the network voltage remotely and real-time. Validation test of the simulated parameters are performed as well.

Algorithm to determine the intersection curves between bezier surfaces by the solution of multivariable polynomial system and the differential marching method

The determination of the intersection curve between Bézier Surfaces may be seen as the composition of two separated problems: determining initial points and tracing the intersection curve from these points. The Bézier Surface is represented by a parametric function (polynomial with two variables) that maps a point in the tridimensional space from the bidimensional parametric space. In this article, it is proposed an algorithm to determine the initial points of the intersection curve of Bézier Surfaces, based on the solution of polynomial systems with the Projected Polyhedral Method, followed by a method for tracing the intersection curves (Marching Method with differential equations). In order to allow the use of the Projected Polyhedral Method, the equations of the system must be represented in terms of the Bernstein basis, and towards this goal it is proposed a robust and reliable algorithm to exactly transform a multivariable polynomial in terms of power basis to a polynomial written in terms of Bernstein basis .

Numerical Simulation of the Cavitation in the Hydrodynamic Lubrication of Journal Bearings: a Parallel Algorithm

In this paper we present an algorithm for the numerical simulation of the cavitation in the hydrodynamic lubrication of journal bearings. Despite the fact that this physical process is usually modelled as a free boundary problem, we adopted the equivalent variational inequality formulation. We propose a two-level iterative algorithm, where the outer iteration is associated to the penalty method, used to transform the variational inequality into a variational equation, and the inner iteration is associated to the conjugate gradient method, used to solve the linear system generated by applying the finite element method to the variational equation. This inner part was implemented using the element by element strategy, which is easily parallelized. We analyse the behavior of two physical parameters and discuss some numerical results. Also, we analyse some results related to the performance of a parallel implementation of the algorithm.

Neutron tomography using projection data obtained by Monte Carlo simulation for nondestructive evaluation

This work present the application of a computer package for generating of projection data for neutron computerized tomography, and in second part, discusses an application of neutron tomography, using the projection data obtained by Monte Carlo technique, for the detection and localization of light materials such as those containing hydrogen, concealed by heavy materials such as iron and lead. For tomographic reconstructions of the samples simulated use was made of only six equal projection angles distributed between 0º and 180º, with reconstruction making use of an algorithm (ARIEM), based on the principle of maximum entropy. With the neutron tomography it was possible to detect and locate polyethylene and water hidden by lead and iron (with 1cm-thick). Thus, it is demonstrated that thermal neutrons tomography is a viable test method which can provide important interior information about test components, so, extremely useful in routine industrial applications.

Streamside management zone (SMZ) efficiency in herbicide retention from simulated surface flow

Plot-scale overland flow experiments were conducted to evaluate the efficiency of streamside management zones (SMZs) in retaining herbicides in runoff generated from silvicultural activities. Herbicide retention was evaluated for five different slopes (2, 5, 10, 15, and 20%), two cover conditions (undisturbed O horizon and raked surface), and two periods under contrasting soil moisture conditions (summer dry and winter wet season) and correlated to O horizon and site conditions. Picloram (highly soluble in water) and atrazine (moderately sorbed into soil particles) at concentrations in the range of 55 and 35 µg L-1 and kaolin clay (approximately 5 g L-1) were mixed with 13.000 liters of water and dispersed over the top of 5 x 10 m forested plots. Surface flow was collected 2, 4, 6, and 10 m below the disperser to evaluate the changes in concentration as it moved through the O horizon and surface soil horizon-mixing zone. Results showed that, on average, a 10 m long forested SMZ removed around 25% of the initial concentration of atrazine and was generally ineffective in reducing the more soluble picloram. Retention of picloram was only 6% of the applied quantity. Percentages of mass reduction by infiltration were 36% for atrazine and 20% for picloram. Stronger relationships existed between O horizon depth and atrazine retention than in any other measured variable, suggesting that better solid-solution contact associated with flow through deeper O horizons is more important than either velocity or soil moisture as a determinant of sorption.

A comparative analysis of preprocessing techniques of cardiac event series for the study of heart rhythm variability using simulated signals

In the present study, using noise-free simulated signals, we performed a comparative examination of several preprocessing techniques that are used to transform the cardiac event series in a regularly sampled time series, appropriate for spectral analysis of heart rhythm variability (HRV). First, a group of noise-free simulated point event series, which represents a time series of heartbeats, was generated by an integral pulse frequency modulation model. In order to evaluate the performance of the preprocessing methods, the differences between the spectra of the preprocessed simulated signals and the true spectrum (spectrum of the model input modulating signals) were surveyed by visual analysis and by contrasting merit indices. It is desired that estimated spectra match the true spectrum as close as possible, showing a minimum of harmonic components and other artifacts. The merit indices proposed to quantify these mismatches were the leakage rate, defined as a measure of leakage components (located outside some narrow windows centered at frequencies of model input modulating signals) with respect to the whole spectral components, and the numbers of leakage components with amplitudes greater than 1%, 5% and 10% of the total spectral components. Our data, obtained from a noise-free simulation, indicate that the utilization of heart rate values instead of heart period values in the derivation of signals representative of heart rhythm results in more accurate spectra. Furthermore, our data support the efficiency of the widely used preprocessing technique based on the convolution of inverse interval function values with a rectangular window, and suggest the preprocessing technique based on a cubic polynomial interpolation of inverse interval function values and succeeding spectral analysis as another efficient and fast method for the analysis of HRV signals

Effects of tryptophan depletion on anxiety induced by simulated public speaking

Several lines of evidence point to the participation of serotonin (5HT) in anxiety. Its specific role, however, remains obscure. The objective of the present study was to evaluate the effect of reducing 5HT-neurotransmission through an acute tryptophan depletion on anxiety induced by a simulated public speaking (SPS) test. Two groups of 14-15 subjects were submitted to a 24-h diet with a low or normal content of tryptophan and received an amino acid mixture without (TRY-) or with (TRY+) tryptophan under double-blind conditions. Five hours later they were submitted to the SPS test. The state-trait anxiety inventory (STAI) and the visual analogue mood scale (VAMS) were used to measure subjective anxiety. Both scales showed that SPS induced a significant increase in anxiety. Although no overall difference between groups was found, there was a trend (P = 0.078) to an interaction of group x gender x phases of the SPS, and a separate analysis of each gender showed an increase in anxiety measured by the STAI in females of the TRY- group. The results for the female TRY- group also suggested a greater arousing effect of the SPS test. In conclusion, the tryptophan depletion procedure employed in the present study did not induce a significant general change in subjective anxiety, but tended to induce anxiety in females. This suggests a greater sensitivity of the 5HT system to the effects of the procedure in this gender.

Yam (Discorea sp) drying with different cuts and temperatures: experimental and simulated results

The yam (Discorea sp) is a tuber rich in carbohydrates, vitamins and mineral salts, besides several components that serve as raw material for medicines. It grows well in tropical and subtropical climates and develops well in zones with an annual pluvial precipitation of around 1300mm, and with cultural treatments, its productivity can exceed 30t/ha. When harvested, the tubers possess about 70% of moisture, and are merchandised "in natura", in the atmospheric temperature, which can cause its fast deterioration. The present work studied the drying of the yam in the form of slices of 1.0 and 2.5cm thickness, as well as in the form of fillets with 1.0 x 1.0 x 5.0cm, with the drying air varying from 40 to 70°C. The equating of the process was accomplished, allowing to simulate the drying as a function of the conditions of the drying air and of the initial and final moisture of the product. Also investigated was the expense of energy as function of the air temperature. The drying in the form of fillets, with the air in a temperature range between 45 and 50°C, was shown to be the most viable process when combining both the quality of the product and the expense of energy.

An artificial neural network model for prediction of quality characteristics of apples during convective dehydration

In this study, the effects of hot-air drying conditions on color, water holding capacity, and total phenolic content of dried apple were investigated using artificial neural network as an intelligent modeling system. After that, a genetic algorithm was used to optimize the drying conditions. Apples were dried at different temperatures (40, 60, and 80 °C) and at three air flow-rates (0.5, 1, and 1.5 m/s). Applying the leave-one-out cross validation methodology, simulated and experimental data were in good agreement presenting an error < 2.4 %. Quality index optimal values were found at 62.9 °C and 1.0 m/s using genetic algorithm.

Extracting Ramachandran torsional angle distributions from 2D NMR data using Tikhonov regularization /

Solid state nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for studying structural and dynamical properties of disordered and partially ordered materials, such as glasses, polymers, liquid crystals, and biological materials. In particular, twodimensional( 2D) NMR methods such as ^^C-^^C correlation spectroscopy under the magicangle- spinning (MAS) conditions have been used to measure structural constraints on the secondary structure of proteins and polypeptides. Amyloid fibrils implicated in a broad class of diseases such as Alzheimer's are known to contain a particular repeating structural motif, called a /5-sheet. However, the details of such structures are poorly understood, primarily because the structural constraints extracted from the 2D NMR data in the form of the so-called Ramachandran (backbone torsion) angle distributions, g{^,'4)), are strongly model-dependent. Inverse theory methods are used to extract Ramachandran angle distributions from a set of 2D MAS and constant-time double-quantum-filtered dipolar recoupling (CTDQFD) data. This is a vastly underdetermined problem, and the stability of the inverse mapping is problematic. Tikhonov regularization is a well-known method of improving the stability of the inverse; in this work it is extended to use a new regularization functional based on the Laplacian rather than on the norm of the function itself. In this way, one makes use of the inherently two-dimensional nature of the underlying Ramachandran maps. In addition, a modification of the existing numerical procedure is performed, as appropriate for an underdetermined inverse problem. Stability of the algorithm with respect to the signal-to-noise (S/N) ratio is examined using a simulated data set. The results show excellent convergence to the true angle distribution function g{(j),ii) for the S/N ratio above 100.

Surface effect ferromagnetism in pure and reduced strontium titanate

A room temperature ferromagnetic hysteresis is observed in single crystal strontium titanate substrates as purchased from several manufacturers. It was found that polishing all sides of the substrates removed this observed hysteresis, suggesting that the origin of the ferromagnetic behavior resides on the surface of the substrates. X-ray diffraction and energy dispersive x-ray spectra were measured however they were unable to detect any impurity phases. In similar semiconducting oxides it was previously suggested that ferromagnetism could originate in oxygen vacancies or from disorder within the single crystal. To this end substrates were annealed in both air and vacuum in a range of temperatures (600°C to 1100°G) to both create bulk oxygen vacancies and to heal surface damage. Annealing in vacuum was found to create a measureable number of oxygen vacancies however their creation could not be correlated to the ferromagnetic signal of the substrate. Annealing in air was found to effect the remnant moment of the substrate as well as the width of the x-ray diffraction peaks on the unpolished face, weakly suggesting a relation between surface based disorder and ferromagnetism. Argon ion bombardment was employed to create a layer of surface disorder in the polished crystal, however it was not found to induce ferromagnetism. It was found that acid etching was sufficient to remove the ferromagnetism from as purchased samples and similarly simulated handling with stainless steel tweezers was sufficient to re-create the ferromagnetism. It is suggested that the origin of this ferromagnetism in SrTi03 is surface contaminants (mainly iron).