Proximal methods for nonlinear programming: double regularization and inexact subproblems

This paper describes the first phase of a project attempting to construct an efficient general-purpose nonlinear optimizer using an augmented Lagrangian outer loop with a relative error criterion, and an inner loop employing a state-of-the art conjugate gradient solver. The outer loop can also employ double regularized proximal kernels, a fairly recent theoretical development that leads to fully smooth subproblems. We first enhance the existing theory to show that our approach is globally convergent in both the primal and dual spaces when applied to convex problems. We then present an extensive computational evaluation using the CUTE test set, showing that some aspects of our approach are promising, but some are not. These conclusions in turn lead to additional computational experiments suggesting where to next focus our theoretical and computational efforts.

On Approximate KKT Condition and its Extension to Continuous Variational Inequalities

In this work, we introduce a necessary sequential Approximate-Karush-Kuhn-Tucker (AKKT) condition for a point to be a solution of a continuous variational inequality, and we prove its relation with the Approximate Gradient Projection condition (AGP) of Garciga-Otero and Svaiter. We also prove that a slight variation of the AKKT condition is sufficient for a convex problem, either for variational inequalities or optimization. Sequential necessary conditions are more suitable to iterative methods than usual punctual conditions relying on constraint qualifications. The AKKT property holds at a solution independently of the fulfillment of a constraint qualification, but when a weak one holds, we can guarantee the validity of the KKT conditions.

A Lagrangian Identification of the Main Sources of Moisture Affecting Northeastern Brazil during Its Pre-Rainy and Rainy Seasons

This work examines the sources of moisture affecting the semi-arid Brazilian Northeast (NEB) during its pre-rainy and rainy season (JFMAM) through a Lagrangian diagnosis method. The FLEXPART model identifies the humidity contributions to the moisture budget over a region through the continuous computation of changes in the specific humidity along back or forward trajectories up to 10 days period. The numerical experiments were done for the period that spans between 2000 and 2004 and results were aggregated on a monthly basis. Results show that besides a minor local recycling component, the vast majority of moisture reaching NEB area is originated in the south Atlantic basin and that the nearby wet Amazon basin bears almost no impact. Moreover, although the maximum precipitation in the ""Poligono das Secas'' region (PS) occurs in March and the maximum precipitation associated with air parcels emanating from the South Atlantic towards PS is observed along January to March, the highest moisture contribution from this oceanic region occurs slightly later (April). A dynamical analysis suggests that the maximum precipitation observed in the PS sector does not coincide with the maximum moisture supply probably due to the combined effect of the Walker and Hadley cells in inhibiting the rising motions over the region in the months following April.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

A Lagrangian identification of major sources of moisture over Central Brazil and La Plata Basin

[1] This work examines the main sources of moisture over Central Brazil and La Plata Basin during the year through a new Lagrangian diagnosis method which identifies the humidity contributions to the moisture budget over a region. This methodology computes budgets of evaporation minus precipitation by calculating changes in the specific humidity along back-trajectories for the previous 10 d. The origin of all air masses residing over each region was tracked during a period of 5 years (2000-2004). These regions were selected because they coincide with two centers of action of a known dipole precipitation variability mode observed in different temporal scales (from intra seasonal up to inter decadal timescales) and are related to the climatic variability of the South American Monsoon System. The results suggested the importance of the tropical south Atlantic as a moisture source for Central Brazil, and of recycling for La Plata basin. It seems that the Tropical South Atlantic plays an important role as a moisture source for Central Brazil and La Plata basin along the year, particularly during the austral summer. The north Atlantic is also an additional source for both regions during the austral summer.

A Lagrangian relaxation approach to a coupled lot-sizing and cutting stock problem

Industrial production processes involving both lot-sizing and cutting stock problems are common in many industrial settings. However, they are usually treated in a separate way, which could lead to costly production plans. In this paper, a coupled mathematical model is formulated and a heuristic method based on Lagrangian relaxation is proposed. Computational results prove its effectiveness. (C) 2009 Elsevier B.V. All rights reserved.

Automatic camera control in virtual environments augmented using multiple sparse videos

Automated virtual camera control has been widely used in animation and interactive virtual environments. We have developed a multiple sparse camera based free view video system prototype that allows users to control the position and orientation of a virtual camera, enabling the observation of a real scene in three dimensions (3D) from any desired viewpoint. Automatic camera control can be activated to follow selected objects by the user. Our method combines a simple geometric model of the scene composed of planes (virtual environment), augmented with visual information from the cameras and pre-computed tracking information of moving targets to generate novel perspective corrected 3D views of the virtual camera and moving objects. To achieve real-time rendering performance, view-dependent textured mapped billboards are used to render the moving objects at their correct locations and foreground masks are used to remove the moving objects from the projected video streams. The current prototype runs on a PC with a common graphics card and can generate virtual 2D views from three cameras of resolution 768 x 576 with several moving objects at about 11 fps. (C)2011 Elsevier Ltd. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

In vitro evaluation of an alternative method to bond molar tubes

Despite the advances in bonding materials, many clinicians today still prefer to place bands on molar teeth. Molar bonding procedures need improvement to be widely accepted clinically. OBJECTIVE: The purpose of this study was to evaluate the shear bond strength when an additional adhesive layer was applied on the occlusal tooth/tube interface to provide reinforcement to molar tubes. MATERIAL AND METHODS: Sixty third molars were selected and allocated to the 3 groups: group 1 received a conventional direct bond followed by the application of an additional layer of adhesive on the occlusal tooth/tube interface, group 2 received a conventional direct bond, and group 3 received a conventional direct bond and an additional cure time of 10 s. The specimens were debonded in a universal testing machine. The results were analyzed statistically by ANOVA and Tukey's test (α=0.05). RESULTS: Group 1 had a significantly higher (p<0.05) shear bond strength compared to groups 2 and 3. No difference was detected between groups 2 and 3 (p>0.05). CONCLUSIONS: The present in vitro findings indicate that the application of an additional layer of adhesive on the tooth/tube interface increased the shear bond strength of the bonded molar tubes.

Synthesis and biocompatibility of an experimental glass ionomer cement prepared by a non-hydrolytic sol-gel method

The aims of this study were to demonstrate the synthesis of an experimental glass ionomer cement (GIC) by the non-hydrolytic sol-gel method and to evaluate its biocompatibility in comparison to a conventional glass ionomer cement (Vidrion R). Four polyethylene tubes containing the tested cements were implanted in the dorsal region of 15 rats, as follows: GI - experimental GIC and GII - conventional GIC. The external tube walls was considered the control group (CG). The rats were sacrificed 7, 21 and 42 days after implant placement for histopathological analysis. A four-point (I-IV) scoring system was used to graduate the inflammatory reaction. Regarding the experimental GIC sintherization, thermogravimetric and x-ray diffraction analysis demonstrated vitreous material formation at 110oC by the sol-gel method. For biocompatibility test, results showed a moderate chronic inflammatory reaction for GI (III), severe for GII (IV) and mild for CG (II) at 7 days. After 21 days, GI presented a mild reaction (II); GII, moderate (III) and CG, mild (II). At 42 days, GI showed a mild/absent inflammatory reaction (II to I), similar to GII (II to I). CG presented absence of chronic inflammatory reaction (I). It was concluded that the experimental GIC presented mild/absent tissue reaction after 42 days, being biocompatible when tested in the connective tissue of rats.

Method for determination of root curvature radius using cone-beam computed tomography images

This article describes and discusses a method to determine root curvature radius by using cone-beam computed tomography (CBCT). The severity of root canal curvature is essential to select instrument and instrumentation technique. The diagnosis and planning of root canal treatment have traditionally been made based on periapical radiography. However, the higher accuracy of CBCT images to identify anatomic and pathologic alterations compared to panoramic and periapical radiographs has been shown to reduce the incidence of false-negative results. In high-resolution images, the measurement of root curvature radius can be obtained by circumcenter. Based on 3 mathematical points determined with the working tools of Planimp® software, it is possible to calculate root curvature radius in both apical and coronal directions. The CBCT-aided method for determination of root curvature radius presented in this article is easy to perform, reproducible and allows a more reliable and predictable endodontic planning, which reflects directly on a more efficacious preparation of curved root canals.

Microleakage in conservative cavities varying the preparation method and surface treatment

OBJECTIVE: To assess microleakage in conservative class V cavities prepared with aluminum-oxide air abrasion or turbine and restored with self-etching or etch-and-rinse adhesive systems. Materials and Methods: Forty premolars were randomly assigned to 4 groups (I and II: air abrasion; III and IV: turbine) and class V cavities were prepared on the buccal surfaces. Conditioning approaches were: groups I/III - 37% phosphoric acid; groups II/IV - self-priming etchant (Tyrian-SPE). Cavities were restored with One Step Plus/Filtek Z250. After finishing, specimens were thermocycled, immersed in 50% silver nitrate, and serially sectioned. Microleakage at the occlusal and cervical interfaces was measured in mm and calculated by a software. Data were subjected to ANOVA and Tukey's test (α=0.05). RESULTS: Marginal seal provided by air abrasion was similar to high-speed handpiece, except for group I. There was SIGNIFICANT difference between enamel and dentin/cementum margins for to group I and II: air abrasion. The etch-and-rinse adhesive system promoted a better marginal seal. At enamel and dentin/cementum margins, the highest microleakage values were found in cavities treated with the self-etching adhesive system. At dentin/cementum margins, high-speed handpiece preparations associated with etch-and-rinse system provided the least dye penetration. CONCLUSION: Marginal seal of cavities prepared with aluminum-oxide air abrasion was different from that of conventionally prepared cavities, and the etch-and-rinse system promoted higher marginal seal at both enamel and dentin margins.

Effect of insertion method on knoop hardness of high viscous glass ionomer cements

The aim of this study was to assess the Knoop hardness of three high viscous glass ionomer cements: G1 - Ketac Molar; G2 - Ketac Molar Easymix (3M ESPE) and G3 - Magic Glass ART (Vigodent). As a parallel goal, three different methods for insertion of Ketac Molar Easymix were tested: G4 - conventional spatula; G5 - commercial syringe (Centrix) and G6 - low-cost syringe. Ten specimens of each group were prepared and the Knoop hardness was determined 5 times on each specimen with a HM-124 hardness machine (25 g/30 s dwell time) after 24 h, 1 and 2 weeks. During the entire test period, the specimens were stored in liquid paraffin at 37ºC. Significant differences were found between G3 and G1/G2 (two-way ANOVA and Tukey's post hoc test; p<0.01). There was no significant difference in the results among the multiple ways of insertion. The glass ionomer cement Magic Glass ART showed the lowest hardness, while the insertion technique had no significant influence on hardness.

A method to determine volatile contaminants in polyethylene terephthalate (PET) packages by HDC-GC-FID and its application to post-consumer materials

A simple and low cost method to determine volatile contaminants in post-consumer recycled PET flakes was developed and validated by Headspace Dynamic Concentration and Gas Chromatography-Flame Ionization Detection (HDC-GC-FID). The analytical parameters evaluated by using surrogates include: correlation coefficient, detection limit, quantification limit, accuracy, intra-assay precision, and inter-assay precision. In order to compare the efficiency of the proposed method to recognized automated techniques, post-consumer PET packaging samples collected in Brazil were used. GC-MS was used to confirm the identity of the substances identified in the PET packaging. Some of the identified contaminants were estimated in the post-consumer material at concentrations higher than 220 ng.g-1. The findings in this work corroborate data available in the scientific literature pointing out the suitability of the proposed analytical method.

An Eulerian Immersed Boundary Method for flow simulations over stationary and moving rigid bodies

The fluid flow over bodies with complex geometry has been the subject of research of many scientists and widely explored experimentally and numerically. The present study proposes an Eulerian Immersed Boundary Method for flows simulations over stationary or moving rigid bodies. The proposed method allows the use of Cartesians Meshes. Here, two-dimensional simulations of fluid flow over stationary and oscillating circular cylinders were used for verification and validation. Four different cases were explored: the flow over a stationary cylinder, the flow over a cylinder oscillating in the flow direction, the flow over a cylinder oscillating in the normal flow direction, and a cylinder with angular oscillation. The time integration was carried out by a classical 4th order Runge-Kutta scheme, with a time step of the same order of distance between two consecutive points in x direction. High-order compact finite difference schemes were used to calculate spatial derivatives. The drag and lift coefficients, the lock-in phenomenon and vorticity contour plots were used for the verification and validation of the proposed method. The extension of the current method allowing the study of a body with different geometry and three-dimensional simulations is straightforward. The results obtained show a good agreement with both numerical and experimental results, encouraging the use of the proposed method.