Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

Global minimization using an Augmented Lagrangian method with variable lower-level constraints

A novel global optimization method based on an Augmented Lagrangian framework is introduced for continuous constrained nonlinear optimization problems. At each outer iteration k the method requires the epsilon(k)-global minimization of the Augmented Lagrangian with simple constraints, where epsilon(k) -> epsilon. Global convergence to an epsilon-global minimizer of the original problem is proved. The subproblems are solved using the alpha BB method. Numerical experiments are presented.

Determination of sorbate and benzoate in beverage samples by capillary electrophoresis - Optimization of the method with inspection of ionic mobilities

The aim of this study was to develop a fast capillary electrophoresis method for the determination of benzoate and sorbate ions in commercial beverages. In the method development the pH and constituents of the background electrolyte were selected using the effective mobility versus pH curves. As the high resolution obtained experimentally for sorbate and benzoate in the studies presented in the literature is not in agreement with that expected from the ionic mobility values published, a procedure to determine these values was carried out. The salicylate ion was used as the internal standard. The background electrolyte was composed of 25 mmol L(-1) tris(hydroxymethyl)aminomethane and 12.5 mmol L(-1) 2-hydroxyisobutyric acid, atpH 8.1.Separation was conducted in a fused-silica capillary(32 cm total length and 8.5 cm effective length, 50 mu m I.D.), with short-end injection configuration and direct UV detection at 200 nm for benzoate and salicylate and 254 nm for sorbate ions. The run time was only 28 s. A few figures of merit of the proposed method include: good linearity (R(2) > 0.999), limit of detection of 0.9 and 0.3 mg L(-1) for benzoate and sorbate, respectively, inter-day precision better than 2.7% (n =9) and recovery in the range 97.9-105%. Beverage samples were prepared by simple dilution with deionized water (1:11, v/v). Concentrations in the range of 197-401 mg L(-1) for benzoate and 28-144 mg L(-1) for sorbate were found in soft drinks and tea. (c) 2008 Elsevier B.V. All rights reserved.

MULTIVARIATE STATISTICAL METHOD IN OPTIMIZATION OF PROTEIN EXTRACTION FROM SOYBEAN FLOUR WITH DAIRY WHEY

The aim objective of this project was to evaluate the protein extraction of soybean flour in dairy whey, by the multivariate statistical method with 2(3) experiments. Influence of three variables were considered: temperature, pH and percentage of sodium chloride against the process specific variable ( percentage of protein extraction). It was observed that, during the protein extraction against time and temperature, the treatments at 80 degrees C for 2h presented great values of total protein (5.99%). The increasing for the percentage of protein extraction was major according to the heating time. Therefore, the maximum point from the function that represents the protein extraction was analysed by factorial experiment 2(3). By the results, it was noted that all the variables were important to extraction. After the statistical analyses, was observed that the parameters as pH, temperature, and percentage of sodium chloride, did not sufficient for the extraction process, since did not possible to obtain the inflection point from mathematical function, however, by the other hand, the mathematical model was significant, as well as, predictive.

Reduced gradient method combined with augmented Lagrangian and barrier for the optimal power flow problem

A new approach for solving the optimal power flow (OPF) problem is established by combining the reduced gradient method and the augmented Lagrangian method with barriers and exploring specific characteristics of the relations between the variables of the OPF problem. Computer simulations on IEEE 14-bus and IEEE 30-bus test systems illustrate the method. (c) 2007 Elsevier Inc. All rights reserved.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Masticatory afficiency in denture wearers with bilateral balanced occlusion and canine guidance

The maintenance of masticatory function is especially important for patients wearing complete dentures due to their limitations. Thus, the bilateral balanced occlusal concept is used to achieve greater masticatory efficiency. However, a critical review of the literature reveals that there is not sufficient scientific evidence to support bilateral balanced occlusion as the most appropriate occlusal concept in complete dentures. Therefore, the aim of this study was to evaluate the masticatory efficiency in complete dentures wearers with bilateral balanced occlusion and canine guidance. A double-blinded controlled crossover clinical trial was conducted. The sample was composed by 24 edentulous patients who wore sets of complete dentures with both occlusal concepts during equal periods of 3 months. Objective data were collected through the masticatory efficiency test performed by the colorimetric method with the beads, in which capsules of a synthetic material enclosing fuchsine-containing granules were used. Subjective data were recorded by patient's ratings of their chewing function. No significant statistical difference was found for masticatory efficiency (p=0.095) between the two occlusal concepts studied. The results suggest that bilateral balanced occlusion does not improve the masticatory efficiency in complete denture wearers.

Mineralogical characterization of a highly-weathered soil by the Rietveld Method

The mineralogical characterization through mineral quantification of Brazilian soils by X-ray diffraction data using the Rietveld Method is not common. A mineralogical quantification of an Acric Ferralsol from the Ponta Grossa region, state of Paraná, Brazil, was carried out using this Method with X-Ray Diffraction data to verify if this method was suitable for mineral quantification of a highly-weathered soil. The A, AB and B3 horizons were fractioned to separate the different particle sizes: clay, silt, fine sand (by Stokes Law) and coarse sand fractions (by sieving), with the procedure free of chemical treatments. X-ray Fluorescence, Inductively Coupled Plasma Atomic Emission Spectrometry, Infrared Spectroscopy and Mössbauer Spectroscopy were used in order to assist the mineral identification and quantification. The Rietveld Method enabled the quantification of the present minerals. In a general way, the quantitative mineralogical characterization by the Rietveld Method revealed that quartz, gibbsite, rutile, hematite, goethite, kaolinite and halloysite were present in the clay and silt fractions of all horizons. The silt fractions of the deeper horizons were different from the more superficial ones due to the presence of large amounts of quartz. The fine and the coarse sand fractions are constituted mainly by quartz. Therefore, a mineralogical quantification of the finer fraction (clay and silt) by the Rietveld Method was successful.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

Electron collisions with the CH(2)O-H(2)O complex

We report cross sections for elastic collisions of low-energy electrons with the CH(2)O-H(2)O complex. We employed the Schwinger multichannel method with pseudopotentials in the static-exchange and in the static-exchange-polarization approximations for energies from 0.1 to 20 eV. We considered four different hydrogen-bonded structures for the complex that were generated by classical Monte Carlo simulations. Our aim is to investigate the effect of the water molecule on the pi* shape resonance of formaldehyde. Previous studies reported a pi* shape resonance for CH(2)O at around 1 eV. The resonance positions of the complexes appear at lower energies in all cases due to the mutual polarization between the two molecules. This indicates that the presence of water may favor dissociation by electron impact and may lead to an important effect on strand breaking in wet DNA by electron impact.

Gravitational wave recoil in Robinson-Trautman spacetimes

We consider the gravitational recoil due to nonreflection-symmetric gravitational wave emission in the context of axisymmetric Robinson-Trautman spacetimes. We show that regular initial data evolve generically into a final configuration corresponding to a Schwarzschild black hole moving with constant speed. For the case of (reflection-)symmetric initial configurations, the mass of the remnant black hole and the total energy radiated away are completely determined by the initial data, allowing us to obtain analytical expressions for some recent numerical results that have appeared in the literature. Moreover, by using the Galerkin spectral method to analyze the nonlinear regime of the Robinson-Trautman equations, we show that the recoil velocity can be estimated with good accuracy from some asymmetry measures (namely the first odd moments) of the initial data. The extension for the nonaxisymmetric case and the implications of our results for realistic situations involving head-on collision of two black holes are also discussed.

Simulation of sugarcane residue decomposition and aboveground growth

Due to the worldwide increase in demand for biofuels, the area cultivated with sugarcane is expected to increase. For environmental and economic reasons, an increasing proportion of the areas are being harvested without burning, leaving the residues on the soil surface. This periodical input of residues affects soil physical, chemical and biological properties, as well as plant growth and nutrition. Modeling can be a useful tool in the study of the complex interactions between the climate, residue quality, and the biological factors controlling plant growth and residue decomposition. The approach taken in this work was to parameterize the CENTURY model for the sugarcane crop, to simulate the temporal dynamics of aboveground phytomass and litter decomposition, and to validate the model through field experiment data. When studying aboveground growth, burned and unburned harvest systems were compared, as well as the effect of mineral fertilizer and organic residue applications. The simulations were performed with data from experiments with different durations, from 12 months to 60 years, in Goiana, TimbaA(0)ba and Pradpolis, Brazil; Harwood, Mackay and Tully, Australia; and Mount Edgecombe, South Africa. The differentiation of two pools in the litter, with different decomposition rates, was found to be a relevant factor in the simulations made. Originally, the model had a basically unlimited layer of mulch directly available for decomposition, 5,000 g m(-2). Through a parameter optimization process, the thickness of the mulch layer closer to the soil, more vulnerable to decomposition, was set as 110 g m(-2). By changing the layer of mulch at any given time available for decomposition, the sugarcane residues decomposition simulations where close to measured values (R (2) = 0.93), contributing to making the CENTURY model a tool for the study of sugarcane litter decomposition patterns. The CENTURY model accurately simulated aboveground carbon stalk values (R (2) = 0.76), considering burned and unburned harvest systems, plots with and without nitrogen fertilizer and organic amendment applications, in different climates and soil conditions.

Mechanical properties and cytotoxicity of 3Y-TZP bioceramics reinforced with Al(2)O(3) particles

The influence of Al(2)O(3) addition and sintering parameters on the mechanical properties and cytotoxicity of tetragonal ZrO(2)-3 mol% Y(2)O(3) ceramics was evaluated. Samples containing 0, 10, 20 and 30 wt.% of Al(2)O(3) particles were prepared by cold uniaxial pressing (80 MPa) and sintered in air at 1500, 1550 and 1600 degrees C for 120 min. The effects of the sintering conditions on the microstructure were analyzed by X-ray diffraction analysis and scanning electron microscopy. Hardness and fracture toughness were determined by the Vickers indentation method and the mechanical resistance by four-point bending tests. As a preliminary biological evaluation, ""in vitro"" cytotoxicity tests were realized to determine the cytotoxic level of the ZrO(2)-Al(2)O(3) composites, using the neutral red uptake method with NCTC clones L929 from the American Type Culture Collection (ATCC) bank. Fully dense ceramic materials were obtained with a hardness ranging between 1340 HV and 1585 HV, depending on the amount of Al(2)O(3) in the ZrO(2) matrix. On the other hand, no significant influence of the Al(2)O(3) addition on fracture toughness was observed, exhibiting values near 8 MPa m(1/2) for all compositions and sintering conditions studied. The non-cytotoxic behavior, the elevated fracture toughness, the good bending strength (sigma(f) = 690 MPa) and the elevated Weibull`s modulus (m = 11) exhibited by the material, show that these ceramic composites are highly suitable biomaterials for dental implant applications. (C) 2008 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Floquet stability analysis of the flow around an oscillating cylinder

This investigative work is concerned with the flow around a circular cylinder submitted to forced transverse oscillations. The goal is to investigate how the transition to turbulence is initiated in the wake for cases with different Reynolds numbers (Re) and displacement amplitudes (A). For each Re the motion frequency is kept constant, close to the Strouhal number of the flow around a fixed cylinder at the same Re. Stability analysis of two-dimensional periodic flows around a forced-oscillating cylinder is carried out with respect to three-dimensional infinitesimal perturbations. The procedure consists of performing a Floquet type analysis of time-periodic base flows, computed using the spectral/hp element method. With the results of the Floquet calculations, considerations regarding the stability of the system are drawn, and the form of the instability at its onset is obtained. The critical Reynolds number is observed to change with the amplitude of oscillation. With respect to instabilities, unstable modes with the same symmetry as mode A of a fixed cylinder are observed; however, they present different wavelengths. Also, the instabilities observed for the oscillating cylinder are distinctively stronger in the braid shear layers. Other unstable modes similar to mode B are found. Quasi-periodic modes are observed in the 2S wake, and subharmonic mode occurrences are reported in P + S wakes. (C) 2009 Elsevier Ltd. All rights reserved.