A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

Structure and stress of Re(11%21); chiral terraces at a racemic surface

The surface structure and morphology of the clean Re(11%21) surface has been investigated through combined low energy electron diffraction intensity analysis of data taken at multiple angles of incidence, scanning tunneling microscopy, and first-principles density functional calculations. The results show how this globally racemic surface terminates in two chirally distinct terraces, which show largescale out-of-plane atomic relaxations and in-plane lateral movement of the uppermost atoms. We further identify and discuss the initial stages of step bunching upon adsorption of oxygen that leads ultimately to the large-scale faceting of the surface. Finally, we present calculations of surface stress and the response to applied surface strain, which suggest routes to the exertion of control over the expression of chirality at the surface.

Reducing noise associated with the Monte Carlo Independent Column Approximation for weather forecasting models

The Monte Carlo Independent Column Approximation (McICA) is a flexible method for representing subgrid-scale cloud inhomogeneity in radiative transfer schemes. It does, however, introduce conditional random errors but these have been shown to have little effect on climate simulations, where spatial and temporal scales of interest are large enough for effects of noise to be averaged out. This article considers the effect of McICA noise on a numerical weather prediction (NWP) model, where the time and spatial scales of interest are much closer to those at which the errors manifest themselves; this, as we show, means that noise is more significant. We suggest methods for efficiently reducing the magnitude of McICA noise and test these methods in a global NWP version of the UK Met Office Unified Model (MetUM). The resultant errors are put into context by comparison with errors due to the widely used assumption of maximum-random-overlap of plane-parallel homogeneous cloud. For a simple implementation of the McICA scheme, forecasts of near-surface temperature are found to be worse than those obtained using the plane-parallel, maximum-random-overlap representation of clouds. However, by applying the methods suggested in this article, we can reduce noise enough to give forecasts of near-surface temperature that are an improvement on the plane-parallel maximum-random-overlap forecasts. We conclude that the McICA scheme can be used to improve the representation of clouds in NWP models, with the provision that the associated noise is sufficiently small.

Stress distribution in orthodontic mini-implants

The objective of this study was to evaluate the stress distribution in the resin in contact with the spirals of cylindrical and conical mini-implants, when submitted to lateral load and insertion torsion. A photoelastic model was fabricated using transparent gelatin to simulate the alveolar bone. The model was observed with a plane polariscope and photographically recorded before and after activation of the two screws with a lateral force and torsion. The lateral force application caused bending moments on both mini-implants, with the uprising of fringes or isochromatics, characteristics of stresses, along the threads of the mini-implants and in the apex. When the torsion was exerted in the mini-implants, a great concentration of stress upraised close to the apex. The conclusion was that, comparing conical with cylindrical mini-implants under lateral load, the stresses were similar on the traction sides. The differences appear (1) on the apex, where the cylindrical mini-implant showed a greater concentration of stress, and (2) along the spirals, in the compression side, where the conical mini-implant showed a greater concentration of stress. The greater part of the stress produced by both mini-implants, after torsion load in insertion, were concentrated on the apex. With the cylindrical mini-implant, the greater concentration of tension was close to the apex, while with the conical one, the stresses were distributed along a greater amount of apical threads.

Study of stress distribution with finite element method in a root with prefabricated threaded split shaft post - Part I

A finite element analysis was carried out to study the role of prefabricated threaded split shaft post (Flexi-Post) on dentinal stress in pulpless tooth. Three dimensional plane strain model of mesio-distal section of a human maxillary central incisor without restoration was analysed with the MSC/NASTRAN (MacNeal/ Schwendler) general purpose finite analysis program was executed on a microcomputer. The model as discretized into 48.954 axisymmetric finite elements defined by 10.355 nodes. Each element was assigned unique elastic properties to represent the materials modeled. Homogeneity, isotropy and linear elasticity were assume for all material. A simulation of static load of 100N was applied to the incisal edge of the post; vertical. Maximal principal stresses and von Mises equivalent stress were calculated. Using the element analysis model employed in this study, the following can be concluded concerning threaded split shaft post (Flexi-Post): Maximum principal stresses in dentin were located at cervical place and at the post apex. The apical threads of the post not redirecting stresses away from the root.

Inclusão de pilares numa formulação do método dos elementos de contorno para análise linear de pavimentos de edifícios

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

Therapeutic concentration of morphine reduces oxidative stress in glioma cell line

Morphine is a potent analgesic opioid used extensively for pain treatment. During the last decade, global consumption grew more than 4-fold. However, molecular mechanisms elicited by morphine are not totally understood. Thus, a growing literature indicates that there are additional actions to the analgesic effect. Previous studies about morphine and oxidative stress are controversial and used concentrations outside the range of clinical practice. Therefore, in this study, we hypothesized that a therapeutic concentration of morphine (1 μM) would show a protective effect in a traditional model of oxidative stress. We exposed the C6 glioma cell line to hydrogen peroxide (H₂O₂) and/or morphine for 24 h and evaluated cell viability, lipid peroxidation, and levels of sulfhydryl groups (an indicator of the redox state of the cell). Morphine did not prevent the decrease in cell viability provoked by H₂O₂) but partially prevented lipid peroxidation caused by 0.0025% H₂O₂) (a concentration allowing more than 90% cell viability). Interestingly, this opioid did not alter the increased levels of sulfhydryl groups produced by exposure to 0.0025% H₂O₂), opening the possibility that alternative molecular mechanisms (a direct scavenging activity or the inhibition of NAPDH oxidase) may explain the protective effect registered in the lipid peroxidation assay. Our results demonstrate, for the first time, that morphine in usual analgesic doses may contribute to minimizing oxidative stress in cells of glial origin. This study supports the importance of employing concentrations similar to those used in clinical practice for a better approximation between experimental models and the clinical setting.

Impact of Quantity of Resin, C-factor, and Geometry on Resin Composite Polymerization Shrinkage Stress in Class V Restorations

Objective: This study evaluated the effect of quantity of resin composite, C-factor, and geometry in Class V restorations on shrinkage stress after bulk fill insertion of resin using two-dimensional finite element analysis.Methods: An image of a buccolingual longitudinal plane in the middle of an upper first premolar and supporting tissues was used for modeling 10 groups: cylindrical cavity, erosion, and abfraction lesions with the same C-factor (1.57), a second cylindrical cavity and abfraction lesion with the same quantity of resin (QR) as the erosion lesion, and then all repeated with a bevel on the occlusal cavosurface angle. The 10 groups were imported into Ansys 13.0 for two-dimensional finite element analysis. The mesh was built with 30,000 triangle and square elements of 0.1 mm in length for all the models. All materials were considered isotropic, homogeneous, elastic, and linear, and the resin composite shrinkage was simulated by thermal analogy. The maximum principal (MPS) and von Mises stresses (VMS) were analyzed for comparing the behavior of the groups.Results: Different values of angles for the cavosurface margin in enamel and dentin were obtained for all groups and the higher the angle, the lower the stress concentration. When the groups with the same C-factor and QR were compared, the erosion shape cavity showed the highest MPS and VMS values, and abfraction shape, the lowest. A cavosurface bevel decreased the stress values on the occlusal margin. The geometry factor overcame the effects of C-factor and QR in some situations.Conclusion: Within the limitations of the current methodology, it is possible to conclude that the combination of all variables studied influences the stress, but the geometry is the most important factor to be considered by the operator.

Frontal Plane Biomechanics in Males and Females with and without Patellofemoral Pain

NAKAGAWA, T. H., E. T. U. MORIYA, C. D. MACIEL, and F. V. SERRAO. Frontal Plane Biomechanics in Males and Females with and without Patellofemoral Pain. Med. Sci. Sports &ere., Vol. 44, No. 9, pp. 1747-1755, 2012. Purpose: The study's purpose was to compare trunk, pelvis, hip, and knee frontal plane biomechanics in males and females with and without patellofemoral pain syndrome (PFPS) during stepping. Methods: Eighty recreational athletes were equally divided into four groups: female PFPS, female controls, male PFPS, and male controls. Trunk, pelvis, hip, and knee frontal plane kinematics and activation of the gluteus medius were evaluated at 15 degrees, 30 degrees, 45 degrees, and 60 degrees of knee flexion during the downward and upward phases of the stepping task. Isometric hip abductor torque was also evaluated. Results: Females showed increased hip adduction and knee abduction at all knee flexion angles, greater ipsilateral trunk lean and contralateral pelvic drop from 60 degrees of knee flexion till the end of the stepping task (P = 0.027-0.001), diminished hip abductor torque (P < 0.001), and increased gluteus medius activation than males (P = 0.008-0.001). PFPS subjects presented increased knee abduction at all the angles evaluated; greater trunk, pelvis, and hip motion from 45 of knee flexion of the downward phase till the end of the maneuver; and diminished gluteus medius activation at 60 degrees of knee flexion, compared with controls (P = 0.034-0.001). Females with PFPS showed lower hip abductor torque compared with the other groups. Conclusions: Females presented with altered frontal plane biomechanics that may predispose them to knee injury. PFPS subjects showed frontal plane biomechanics that could increase the lateral patellofemoral joint stress at all the angles evaluated and could increase even more from 45 degrees of knee flexion in the downward phase untill the end of the maneuver. Hip abductor strengthening and motor control training should be considered when treating females with PFPS.

Comparison of active and passive forces of the pelvic floor muscles in women with and without stress urinary incontinence

Background: The reduction of the pelvic floor muscles (PFM) strength is a major cause of stress urinary incontinence (SUI). Objective: To compare active and passive forces, and vaginal cavity aperture in continent and stress urinary incontinent women. Method: The study included a total of thirty-two women, sixteen continent women (group 1 - G1) and sixteen women with SUI (group 2 - G2). To evaluate PFM passive and active forces in anteroposterior (sagittal plane) and left-right directions (frontal plane) a stainless steel specular dynamometer was used. Results: The anteroposterior active strength for the continent women (mean +/- standard deviation) (0.3 +/- 0.2 N) was greater compared to the values found in the evaluation of incontinent women (0.1 +/- 0.1 N). The left-right active strength (G1=0.43 +/- 0.1 N; G2=0.40 +/- 0.1 N), the passive force (G1=1.1 +/- 0.2 N; G2=1.1 +/- 0.3 N) and the vaginal cavity aperture (G1=21 +/- 3 mm; G2=24 +/- 4 mm) did not differ between groups 1 and 2. Conclusion: The function evaluation of PFM showed that women with SUI had a lower anteroposterior active strength compared to continent women.

Effect of sample pre-cracking method and notch geometry in plane strain fracture toughness tests as applied to a PMMA resin

The influence of test method factors (notch shape, square or angular, and pre-cracking method, by tapping onto or pressing a razor blade) on the results obtained in plane strain fracture toughness test according to standard ASTM D5045 using SENB specimens made of a commercial PMMA resin were investigated. Results were analyzed quantitatively by comparing the obtained K-IC values and qualitatively by observing their effect on the Moire fringes observed using photoelasticity, showing that, at 95% significance level, the K-IC values are affected by the pre-cracking method, with the most conservative value being obtained when natural pre-cracks were introduced by tapping onto a razor blade (K-IC = 1.15 +/- 0.11 MPa.m(0.5)). This correlates with a perturbation in the stress field close to the pre-crack tip observed in the photoelasticity test sample when it was introduced by pressing the razor blade. Surprisingly, notch geometry only slightly affects the results. (C) 2012 Elsevier Ltd. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Issues on the mean stress effect in fretting fatigue of a 7050-T7451 Al alloy posed by new experimental data

The aims of this work are: (i) to produce new experimental data for fretting fatigue considering the presence of a mean bulk stress and (ii) to assess two design methodologies against failure by fretting fatigue. Tests on a cylinder–flat contact configuration were conducted using a fretting apparatus mounted on a servo-hydraulic machine. The material used for both the pads and fatigue specimen was an aeronautical 7050-T7451 Al alloy. The experimental program was designed with all relevant parameters, apart from the mean bulk load (always applied before the contact loads), kept constant. The mean bulk stress varied from compressive to tensile values while maintaining a high peak pressure in order to encourage crack initiation. Two methodologies against fretting fatigue are proposed and confronted against the experimental data. The non-local stress-based methodology considers the evaluation of a critical plane fatigue criterion at the center of a process zone located beneath the contacting surfaces. The results showed that it correctly predicts crack initiation, but was not capable to provide successful prediction of the integrity of the specimens. Alternatively, we considered a crack arrest criterion which has the potential to provide a more complete description about the integrity of the specimens.

Quantum effects and dynamics in hydrogen-bonded systems: a first-principles approach to spectroscopic experiments

Computer simulations have become an important tool in physics. Especially systems in the solid state have been investigated extensively with the help of modern computational methods. This thesis focuses on the simulation of hydrogen-bonded systems, using quantum chemical methods combined with molecular dynamics (MD) simulations. MD simulations are carried out for investigating the energetics and structure of a system under conditions that include physical parameters such as temperature and pressure. Ab initio quantum chemical methods have proven to be capable of predicting spectroscopic quantities. The combination of these two features still represents a methodological challenge. Furthermore, conventional MD simulations consider the nuclei as classical particles. Not only motional effects, but also the quantum nature of the nuclei are expected to influence the properties of a molecular system. This work aims at a more realistic description of properties that are accessible via NMR experiments. With the help of the path integral formalism the quantum nature of the nuclei has been incorporated and its influence on the NMR parameters explored. The effect on both the NMR chemical shift and the Nuclear Quadrupole Coupling Constants (NQCC) is presented for intra- and intermolecular hydrogen bonds. The second part of this thesis presents the computation of electric field gradients within the Gaussian and Augmented Plane Waves (GAPW) framework, that allows for all-electron calculations in periodic systems. This recent development improves the accuracy of many calculations compared to the pseudopotential approximation, which treats the core electrons as part of an effective potential. In combination with MD simulations of water, the NMR longitudinal relaxation times for 17O and 2H have been obtained. The results show a considerable agreement with the experiment. Finally, an implementation of the calculation of the stress tensor into the quantum chemical program suite CP2K is presented. This enables MD simulations under constant pressure conditions, which is demonstrated with a series of liquid water simulations, that sheds light on the influence of the exchange-correlation functional used on the density of the simulated liquid.

High resolution study of local stress inside alumina - micro mechanical analysis using laser scanning confocal microscope

The aim of this work is to measure the stress inside a hard micro object under extreme compression. To measure the internal stress, we compressed ruby spheres (a-Al2O3: Cr3+, 150 µm diameter) between two sapphire plates. Ruby fluorescence spectrum shifts to longer wavelengths under compression and can be related to the internal stress by a conversion coefficient. A confocal laser scanning microscope was used to excite and collect fluorescence at desired local spots inside the ruby sphere with spatial resolution of about 1 µm3. Under static external loads, the stress distribution within the center plane of the ruby sphere was measured directly for the first time. The result agreed to Hertz’s law. The stress across the contact area showed a hemispherical profile. The measured contact radius was in accord with the calculation by Hertz’s equation. Stress-load curves showed spike-like decrease after entering non-elastic phase, indicating the formation and coalescence of microcracks, which led to relaxing of stress. In the vicinity of the contact area luminescence spectra with multiple peaks were observed. This indicated the presence of domains of different stress, which were mechanically decoupled. Repeated loading cycles were applied to study the fatigue of ruby at the contact region. Progressive fatigue was observed when the load exceeded 1 N. As long as the load did not exceed 2 N stress-load curves were still continuous and could be described by Hertz’s law with a reduced Young’s modulus. Once the load exceeded 2 N, periodical spike-like decreases of the stress could be observed, implying a “memory effect” under repeated loading cycles. Vibration loading with higher frequencies was applied by a piezo. Redistributions of intensity on the fluorescence spectra were observed and it was attributed to the repopulation of the micro domains of different elasticity. Two stages of under vibration loading were suggested. In the first stage continuous damage carried on until certain limit, by which the second stage, e.g. breakage, followed in a discontinuous manner.