Fractional polynomial response surface models

Second-order polynomial models have been used extensively to approximate the relationship between a response variable and several continuous factors. However, sometimes polynomial models do not adequately describe the important features of the response surface. This article describes the use of fractional polynomial models. It is shown how the models can be fitted, an appropriate model selected, and inference conducted. Polynomial and fractional polynomial models are fitted to two published datasets, illustrating that sometimes the fractional polynomial can give as good a fit to the data and much more plausible behavior between the design points than the polynomial model. © 2005 American Statistical Association and the International Biometric Society.

Hyaluronidase from the venom of the social wasp Polybia paulista (Hymenoptera, Vespidae): Cloning, structural modeling, purification, and immunological analysis

In this study, we describe the cDNA cloning, sequencing, and 3-D structure of the allergen hyaluronidase from Polybia paulista venom (Pp-Hyal). Using a proteomic approach, the native form of Pp-Hyal was purified to homogeneity and used to produce a Pp-specific polyclonal antibody. The results revealed that Pp-Hyal can be classified as a glycosyl hydrolase and that the full-length Pp-Hyal cDNA (1315 bp; GI: 302201582) is similar (80-90%) to hyaluronidase from the venoms of endemic Northern wasp species. The isolated mature protein is comprised of 338 amino acids, with a theoretical pI of 8.77 and a molecular mass of 39,648.8 Da versus a pI of 8.13 and 43,277.0 Da indicated by MS. The Pp-Hyal 3D-structural model revealed a central core (α/β)7 barrel, two sulfide bonds (Cys 19-308 and Cys 185-197), and three putative glycosylation sites (Asn79, Asn187, and Asn325), two of which are also found in the rVes v 2 protein. Based on the model, residues Ser299, Asp107, and Glu109 interact with the substrate and potential epitopes (five conformational and seven linear) located at surface-exposed regions of the structure. Purified native Pp-Hyal showed high similarity (97%) with hyaluronidase from Polistes annularis venom (Q9U6V9). Immunoblotting analysis confirmed the specificity of the Pp-Hyal-specific antibody as it recognized the Pp-Hyal protein in both the purified fraction and P. paulista crude venom. No reaction was observed with the venoms of Apis mellifera, Solenopsis invicta, Agelaia pallipes pallipes, and Polistes lanio lanio, with the exception of immune cross-reactivity with venoms of the genus Polybia (sericea and ignobilis). Our results demonstrate cross-reactivity only between wasp venoms from the genus Polybia. The absence of cross-reactivity between the venoms of wasps and bees observed here is important because it allows identification of the insect responsible for sensitization, or at least of the phylogenetically closest insect, in order to facilitate effective immunotherapy in allergic patients. © 2013 Elsevier Ltd.

Control of Limit Cycle Oscillation in a Three Degrees of Freedom Airfoil Section Using Fuzzy Takagi-Sugeno Modeling

This work presents a strategy to control nonlinear responses of aeroelastic systems with control surface freeplay. The proposed methodology is developed for the three degrees of freedom typical section airfoil considering aerodynamic forces from Theodorsen's theory. The mathematical model is written in the state space representation using rational function approximation to write the aerodynamic forces in time domain. The control system is designed using the fuzzy Takagi-Sugeno modeling to compute a feedback control gain. It useds Lyapunov's stability function and linear matrix inequalities (LMIs) to solve a convex optimization problem. Time simulations with different initial conditions are performed using a modified Runge-Kutta algorithm to compare the system with and without control forces. It is shown that this approach can compute linear control gain able to stabilize aeroelastic systems with discontinuous nonlinearities.

Inverse Modeling of GPS Multipath for Snow Depth Estimation-Part I: Formulation and Simulations

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

Surface coverage in wireless sensor networks based on delaunay tetrahedralization

In this work is presented a new method for sensor deployment on 3D surfaces. The method was structured on different steps. The first one aimed discretizes the relief of interest with Delaunay algorithm. The tetrahedra and relative values (spatial coordinates of each vertex and faces) were input to construction of 3D Voronoi diagram. Each circumcenter was calculated as a candidate position for a sensor node: the corresponding circular coverage area was calculated based on a radius r. The r value can be adjusted to simulate different kinds of sensors. The Dijkstra algorithm and a selection method were applied to eliminate candidate positions with overlapped coverage areas or beyond of surface of interest. Performance evaluations measures were defined using coverage area and communication as criteria. The results were relevant, once the mean coverage rate achieved on three different surfaces were among 91% and 100%.

Coupling remote sensing bio-optical and three-dimensional hydrodynamic modeling to study the phytoplankton dynamics in a tropical hydroelectric reservoir

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

Computer graphics applied to anatomy: a study of two bio-cad modeling methods on finite element analysis of human edentulous hemi-mandible

Modeling is a step to perform a finite element analysis. Different methods of model construction are reported in literature, as the Bio-CAD modeling. The purpose of this study was to perform a model evaluation and application using two methods of Bio-CAD modeling from human edentulous hemi-mandible on the finite element analysis. From CT scans of dried human skull was reconstructed a stereolithographic model. Two methods of modeling were performed: STL conversion approach (Model 1) associated to STL simplification and reverse engineering approach (Model 2). For finite element analysis was used the action of lateral pterygoid muscle as loading condition to assess total displacement (D), equivalent von-Mises stress (VM) and maximum principal stress (MP). Two models presented differences on the geometry regarding surface number (1834 (model 1); 282 (model 2)). Were observed differences in finite element mesh regarding element number (30428 nodes/16683 elements (model 1); 15801 nodes/8410 elements (model 2). D, VM and MP stress areas presented similar distribution in two models. The values were different regarding maximum and minimum values of D (ranging 0-0.511 mm (model 1) and 0-0.544 mm (model 2), VM stress (6.36E-04-11.4 MPa (model 1) and 2.15E-04-14.7 MPa (model 2) and MP stress (-1.43-9.14 MPa (model 1) and -1.2-11.6 MPa (model 2). From two methods of Bio-CAD modeling, the reverse engineering presented better anatomical representation compared to the STL conversion approach. The models presented differences in the finite element mesh, total displacement and stress distribution.

Influence of summertime mesoscale convective systems on the heat balance and surface mixed layer dynamics of a large Amazonian hydroelectric reservoir

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

Insights into the morphology of the Serra da Cangalha impact structure from geophysical modeling

Forward modeling is commonly applied to gravity field data of impact structures to determine the main gravity anomaly sources. In this context, we have developed 2.5-D gravity models of the Serra da Cangalha impact structure for the purpose of investigating geological bodies/structures underneath the crater. Interpretation of the models was supported by ground magnetic data acquired along profiles, as well as by high resolution aeromagnetic data. Ground magnetic data reveal the presence of short-wavelength anomalies probably related to shallow magnetic sources that could have been emplaced during the cratering process. Aeromagnetic data show that the basement underneath the crater occurs at an average depth of about 1.9 km, whereas in the region beneath the central uplift it is raised to 0.51 km below the current surface. These depths are also supported by 2.5-D gravity models showing a gentle relief for the basement beneath the central uplift area. Geophysical data were used to provide further constraints for numeral modeling of crater formation that provided important information on the structural modification that affected the rocks underneath the crater, as well as on shock-induced modifications of target rocks. The results showed that the morphology is consistent with the current observations of the crater and that Serra da Cangalha was formed by a meteorite of approximately 1.4 km diameter striking at 12 km s-1.

Mathematical modeling and analysis of the particle interactions in the direct filtration using the colloidal theory

This work used the colloidal theory to describe forces and energy interactions of colloidal complexes in the water and those formed during filtration run in direct filtration. Many interactions of particle energy profiles between colloidal surfaces for three geometries are presented here in: spherical, plate and cylindrical; and four surface interactions arrangements: two cylinders, two spheres, two plates and a sphere and a plate. Two different situations were analyzed, before and after electrostatic destabilization by action of the alum sulfate as coagulant in water studies samples prepared with kaolin. In the case were used mathematical modeling by extended DLVO theory (from the names: Derjarguin-Landau-Verwey-Overbeek) or XDLVO, which include traditional approach of the electric double layer (EDL), surfaces attraction forces or London-van der Waals (LvdW), esteric forces and hydrophobic forces, additionally considering another forces in colloidal system, like molecular repulsion or Born Repulsion and Acid-Base (AB) chemical function forces from Lewis.

Post-breakup tectonics in southeast Brazil from thermochronological data and combined inverse-forward thermal history modeling

The continental margin of southeast Brazil is elevated. Onshore Tertiary basins and Late Cretaceous/Paleogene intrusions are good evidence for post breakup tectono-magmatic activity. To constrain the impact of post-rift reactivation on the geological history of the area, we carried out a new thermochronological study. Apatite fission track ages range from 60.7 +/- 1.9 Ma to 129.3 +/- 4.3 Ma, mean track lengths from 11.41 +/- 0.23 mu m to 14.31 +/- 0.24 mu m and a subset of the (U-Th)/He ages range from 45.1 +/- 1.5 to 122.4 +/- 2.5 Ma. Results of inverse thermal history modeling generally support the conclusions from an earlier study for a Late Cretaceous phase of cooling. Around the onshore Taubate Basin, for a limited number of samples, the first detectable period of cooling occurred during the Early Tertiary. The inferred thermal histories for many samples also imply subsequent reheating followed by Neogene cooling. Given the uncertainty of the inversion results, we did deterministic forward modeling to assess the range of possibilities of this Tertiary part of the thermal history. The evidence for reheating seems to be robust around the Taubate Basin, but elsewhere the data cannot discriminate between this and a less complex thermal history. However, forward modeling results and geological information support the conclusion that the whole area underwent cooling during the Neogene. The synchronicity of the cooling phases with Andean tectonics and those in NE Brazil leads us to assume a plate-wide compressional stress that reactivated inherited structures. The present-day topographic relief of the margin reflects a contribution from post-breakup reactivation and uplift.

Conformation analysis of a surface loop that controls active site access in the GH11 xylanase A from Bacillus subtilis

Xylanases (EC 3.2.1.8 endo-1,4-glycosyl hydrolase) catalyze the hydrolysis of xylan, an abundant hemicellulose of plant cell walls. Access to the catalytic site of GH11 xylanases is regulated by movement of a short beta-hairpin, the so-called thumb region, which can adopt open or closed conformations. A crystallographic study has shown that the D11F/R122D mutant of the GH11 xylanase A from Bacillus subtilis (BsXA) displays a stable "open" conformation, and here we report a molecular dynamics simulation study comparing this mutant with the native enzyme over a range of temperatures. The mutant open conformation was stable at 300 and 328 K, however it showed a transition to the closed state at 338 K. Analysis of dihedral angles identified thumb region residues Y113 and T123 as key hinge points which determine the open-closed transition at 338 K. Although the D11F/R122D mutations result in a reduction in local inter-intramolecular hydrogen bonding, the global energies of the open and closed conformations in the native enzyme are equivalent, suggesting that the two conformations are equally accessible. These results indicate that the thumb region shows a broader degree of energetically permissible conformations which regulate the access to the active site region. The R122D mutation contributes to the stability of the open conformation, but is not essential for thumb dynamics, i.e., the wild type enzyme can also adapt to the open conformation.

Dynamics of the cold surface layer of polythermal Storglaciären, Sweden

Polythermal glaciers, i.e. glaciers with a combination of ice at and below the freezing point, are widespread in arctic and subarctic environments. The polythermal structure has major implications for glacier hydrology, ice flow and glacial erosion. However, the interplay of factors governing its spatial and temporal variations such as net mass balance, ice advection and water content in the ice is poorly investigated and as yet not fully understood. This study deals with a thorough investigation of the polythermal regime on Storglaciären, northern Sweden, a small valley glacier with a cold surface layer in the ablation area. Extensive field work was performed including mapping of the cold surface layer using ground-penetrating radar, ice temperature measurements, mass balance and ice velocity measurements. Analyses of these data combined with numerical modelling were used specifically to investigate the spatial and temporal variability of the cold surface layer, the spatial distribution of the water content just below the cold surface layer transition, the effect of radar frequency on the detection of the surface layer, and the sensitivity of the cold surface layer to changes in forcing. A comparison between direct temperature measurements in boreholes and ground-penetrating surveys shows that the radar-inferred cold-temperate transition depth is within ±1 m from the melting point of ice at frequencies above ~300 MHz. At frequencies below ~155 MHz, the accuracy degrades because of reduced scattering efficiency that occurs when the scatterers become much smaller compared to the wavelength. The mapped spatial pattern of the englacial cold-temperate transition boundary is complex. This pattern reflects the observed spatial variation in net loss of ice at the surface by ablation and vertical advection of ice, which is suggested to provide the predominant forcing of the cold surface layer thickness pattern. This is further supported by thermomechanical modeling of the cold surface layer, which indicates high sensitivity of the cold surface layer thickness to changes in vertical advection rates. The water content is the least investigated quantity that is relevant for the thermal regime of glaciers, but also the most difficult to assess. Spatial variability of absolute water content in the temperate ice immediately below the cold surface layer on Storglaciären was determined by combining relative estimates of water content from ground-penetrating radar data with absolute determination from temperature measurements and the thermal boundary condition at the freezing front. These measurements indicate large-scale spatial variability in the water content, which seems to arise from variations in entrapment of water at the firn-ice transition. However, this variability cannot alone explain the spatial pattern in the thermal regime on Storglaciären. Repeated surveys of the cold surface layer show a 22% average thinning of the cold surface layer on Storglaciären between 1989 and 2001. Transient thermomechanical modeling results suggest that the cold surface layer adapts to new equilibrium conditions in only a few decades after a perturbation in the forcing is introduced. An increased winter air temperature since mid-1980s seems to be the cause of the observed thinning of the cold surface layer. Over the last decades, mass balance measurements indicate that the glacier has been close to a steady state. The quasi-steady state situation is also reflected in the vertical advection, which shows no significant changes during the last decades. Increased winter temperatures at the ice surface would result in a slow-down of the formation of cold ice at the base of the cold surface layer and lead to a larger imbalance between net loss of ice at the surface and freezing of temperate ice at the cold-temperate transition.

Centrifuge Modeling of Sandy Slopes

Slope failure occurs in many areas throughout the world and it becomes an important problem when it interferes with human activity, in which disasters provoke loss of life and property damage. In this research we investigate the slope failure through the centrifuge modeling, where a reduced-scale model, N times smaller than the full-scale (prototype), is used whereas the acceleration is increased by N times (compared with the gravity acceleration) to preserve the stress and the strain behavior. The aims of this research “Centrifuge modeling of sandy slopes” are in extreme synthesis: 1) test the reliability of the centrifuge modeling as a tool to investigate the behavior of a sandy slope failure; 2) understand how the failure mechanism is affected by changing the slope angle and obtain useful information for the design. In order to achieve this scope we arranged the work as follows: Chapter one: centrifuge modeling of slope failure. In this chapter we provide a general view about the context in which we are working on. Basically we explain what is a slope failure, how it happens and which are the tools available to investigate this phenomenon. Afterwards we introduce the technology used to study this topic, that is the geotechnical centrifuge. Chapter two: testing apparatus. In the first section of this chapter we describe all the procedures and facilities used to perform a test in the centrifuge. Then we explain the characteristics of the soil (Nevada sand), like the dry unit weight, water content, relative density, and its strength parameters (c,φ), which have been calculated in laboratory through the triaxial test. Chapter three: centrifuge tests. In this part of the document are presented all the results from the tests done in centrifuge. When we talk about results we refer to the acceleration at failure for each model tested and its failure surface. In our case study we tested models with the same soil and geometric characteristics but different angles. The angles tested in this research were: 60°, 75° and 90°. Chapter four: slope stability analysis. We introduce the features and the concept of the software: ReSSA (2.0). This software allows us to calculate the theoretical failure surfaces of the prototypes. Then we show in this section the comparisons between the experimental failure surfaces of the prototype, traced in the laboratory, and the one calculated by the software. Chapter five: conclusion. The conclusion of the research presents the results obtained in relation to the two main aims, mentioned above.

The meccano method for isogeometric solid modeling and applications

[EN]We present a new method to construct a trivariate T-spline representation of complex solids for the application of isogeometric analysis. We take a genus-zero solid as a basis of our study, but at the end of the work we explain the way to generalize the results to any genus solids. The proposed technique only demands a surface triangulation of the solid as input data. The key of this method lies in obtaining a volumetric parameterization between the solid and the parametric domain, the unitary cube. To do that, an adaptive tetrahedral mesh of the parametric domain is isomorphically transformed onto the solid by applying a mesh untangling and smoothing procedure...