Alternate treatments of jacobian singularities in polar coordinates within finite-difference schemes

Jacobian singularities of differential operators in curvilinear coordinates occur when the Jacobian determinant of the curvilinear-to-Cartesian mapping vanishes, thus leading to unbounded coefficients in partial differential equations. Within a finite-difference scheme, we treat the singularity at the pole of polar coordinates by setting up complementary equations. Such equations are obtained by either integral or smoothness conditions. They are assessed by application to analytically solvable steady-state heat-conduction problems.

Scales, Universality and Finite-Range Correction in Three-body Systems

The scale invariance manifested by the weakly-bound Efimov states implies that all the Efimov spectrum can be merged in a single scaling function. By considering this scaling function, the ratio between two consecutive energy levels, E3 (N+1) and E3 (N), can be obtained from a two-body low-energy observable (usually the scattering length a), given in units of the three-body energy level N. The zero-ranged scaling function is improved by incorporating finite range corrections in first order of r0/a (r0 is the potential effective range). The critical condition for three-identical bosons in s-wave, when the excited E3 (N+1) state disappears in the 2 + 1 threshold, is given by √E2/E3 (N) ≈ 0.38+0.12(r0/a). © 2012 Springer-Verlag.

Spin squeezing and entanglement via finite-dimensional discrete phase-space description

We show how mapping techniques inherent to N2-dimensional discrete phase spaces can be used to treat a wide family of spin systems which exhibits squeezing and entanglement effects. This algebraic framework is then applied to the modified Lipkin-Meshkov-Glick (LMG) model in order to obtain the time evolution of certain special parameters related to the Robertson- Schrödinger (RS) uncertainty principle and some particular proposals of entanglement measure based on collective angular-momentum generators. Our results reinforce the connection between both the squeezing and entanglement effects, as well as allow to investigate the basic role of spin correlations through the discrete representatives of quasiprobability distribution functions. Entropy functionals are also discussed in this context. The main sequence correlations → entanglement → squeezing of quantum effects embraces a new set of insights and interpretations in this framework, which represents an effective gain for future researches in different spin systems. © 2013 World Scientific Publishing Company.

Study of the ℤ3 symmetry in QCD at finite temperature and chemical potential using a worm algorithm

Traditional Monte Carlo simulations of QCD in the presence of a baryon chemical potential are plagued by the complex phase problem and new numerical approaches are necessary for studying the phase diagram of the theory. In this work we consider a ℤ3 Polyakov loop model for the deconfining phase transition in QCD and discuss how a flux representation of the model in terms of dimer and monomer variable solves the complex action problem. We present results of numerical simulations using a worm algorithm for the specific heat and two-point correlation function of Polyakov loops. Evidences of a first order deconfinement phase transition are discussed. © 2013 American Institute of Physics.

Comparison of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible - A computed tomography-based three-dimensional finite element analysis

A finite element analysis was used to compare the effect of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible. Four models of an human mandible were constructed. In the OR (O'ring) group, the mandible was restored with an overdenture retained by four unsplinted implants with O'ring attachment; in the BC (bar-clip) -C and BC groups, the mandibles were restored with overdentures retained by four splinted implants with bar-clip anchor associated or not with two distally placed cantilevers, respectively; in the FD (fixed denture) group, the mandible was restored with a fixed full-arch four-implant-supported prosthesis. Models were supported by the masticatory muscles and temporomandibular joints. A 100-N oblique load was applied on the left first molar. Von Mises (σvM), maximum (σmax) and minimum (σmin) principal stresses (in MPa) analyses were obtained. BC-C group exhibited the highest stress values (σvM=398.8, σmax=580.5 and σmin=-455.2) while FD group showed the lowest one (σvM=128.9, σmax=185.9 and σmin=-172.1). Within overdenture groups, the use of unsplinted implants reduced the stress level in the implant/prosthetic components (59.4% for σvM, 66.2% for σmax and 57.7% for σmin versus BC-C group) and supporting tissues (maximum stress reduction of 72% and 79.5% for σmax, and 15.7% and 85.7% for σmin on the cortical and trabecular bones, respectively). Cortical bone exhibited greater stress concentration than the trabecular bone for all groups. The use of fixed implant dentures and removable dentures retained by unsplinted implants to rehabilitate edentulous mandible reduced the stresses in the periimplant bone tissue, mucosa and implant/prosthetic components. © 2013 Elsevier Ltd.

Free tools and strategies for the generation of 3D finite element meshes: Modeling of the cardiac structures

The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context. © 2013 E. Pavarino et al.

Reduction of infinite dimensional systems to finite dimensions: Compact convergence approach

We consider parameter dependent semilinear evolution problems for which, at the limit value of the parameter, the problem is finite dimensional. We introduce an abstract functional analytic framework that applies to many problems in the existing literature for which the study of asymptotic dynamics can be reduced to finite dimensions via the invariant manifolds technique. Some practical models are considered to show wide applicability of the theory. © 2013 Society for Industrial and Applied Mathematics.

A finite element approach for predicting the ultimate rotation capacity of RC beams

This paper presents a numerical approach to model the complex failure mechanisms that define the ultimate rotational capacity of reinforced concrete beams. The behavior in tension and compression is described by a constitutive damage model derived from a combination of two specific damage models [1]. The nonlinear behavior of the compressed region is treated by the compressive damage model based on the Drucker-Prager criterion written in terms of the effective stresses. The tensile damage model employs a failure criterion based on the strain energy associated with the positive part the effective stress tensor. This model is used to describe the behavior of very thin bands of strain localization, which are embedded in finite elements to represent multiple cracks that occur in the tensioned region [2]. The softening law establishes dissipation energy compatible with the fracture energy of the concrete. The reinforcing steel bars are modeled by truss elements with elastic-perfect plastic behavior. It is shown that the resulting approach is able to predict the different stages of the collapse mechanism of beams with distinct sizes and reinforcement ratios. The tensile damage model and the finite element embedded crack approach are able to describe the stiffness reduction due to concrete cracking in the tensile zone. The truss elements are able to reproduce the effects of steel yielding and, finally, the compressive damage model is able to describe the non-linear behavior of the compressive zone until the complete collapse of the beam due to crushing of concrete. The proposed approach is able to predict well the plastic rotation capacity of tested beams [3], including size-scale effects.

Influence of the implant diameter with different sizes of hexagon: Analysis by 3-dimensional finite element method

The aim of this study was to evaluate the stress distribution in implants of regular platforms and of wide diameter with different sizes of hexagon by the 3-dimensional finite element method. We used simulated 3-dimensional models with the aid of Solidworks 2006 and Rhinoceros 4.0 software for the design of the implant and abutment and the InVesalius software for the design of the bone. Each model represented a block of bone from the mandibular molar region with an implant 10 mm in length and different diameters. Model A was an implant 3.75 mm/regular hexagon, model B was an implant 5.00 mm/regular hexagon, and model C was an implant 5.00 mm/ expanded hexagon. A load of 200 N was applied in the axial, lateral, and oblique directions. At implant, applying the load (axial, lateral, and oblique), the 3 models presented stress concentration at the threads in the cervical and middle regions, and the stress was higher for model A. At the abutment, models A and B showed a similar stress distribution, concentrated at the cervical and middle third; model C showed the highest stresses. On the cortical bone, the stress was concentrated at the cervical region for the 3 models and was higher for model A. In the trabecular bone, the stresses were less intense and concentrated around the implant body, and were more intense for model A. Among the models of wide diameter (models B and C), model B (implant 5.00 mm/regular hexagon) was more favorable with regard to distribution of stresses. Model A (implant 3.75 mm/regular hexagon) showed the largest areas and the most intense stress, and model B (implant 5.00 mm/regular hexagon) showed a more favorable stress distribution. The highest stresses were observed in the application of lateral load.

Three-dimensional finite element analysis of stress distribution in retention screws of different crown-implant ratios

The retaining screw of the implant-supported dental prosthesis is the weakest point of the crown/implant system. Furthermore, crown height is another important factor that may increase the lever arm. Therefore, the aim of this study was to assess the stress distribution in implant prosthetic screws with different heights of the clinical crown of the prosthesis using the method of three-dimensional finite element analysis. Three models were created with implants (3.75 mm × 10 mm) and crowns (heights of 10, 12.5 and 15 mm). The results were visualised by means of von Mises stress maps that increased the crown heights. The screw structure exhibited higher levels of stresses in the oblique load. The oblique loading resulted in higher stress concentration when compared with the axial loading. It is concluded that the increase of the crown was damaging to the stress distribution on the screw, mainly in oblique loading. © 2013 Taylor & Francis.

Luminescence quantum efficiency at 1.5 μm of Er3+-doped tellurite glass determined by thermal lens spectroscopy

Erbium doped tellurite glasses (TeO2 + Li2O + TiO2) were prepared by conventional melt-quenching method to study the influence of the Er3+ concentration on the luminescence quantum efficiency (η) at 1.5 μm. Absorption and luminescence data were used to characterize the samples, and the η parameter was measured using the well-known thermal lens spectroscopy. For low Er3+ concentration, the measured values are around 76%, and the concentration behavior of η shows Er-Er and Er-OH- interactions, which agreed with the measured lifetime values. © 2013 Elsevier B.V. All rights reserved.

Stress analysis in bone tissue around single implants with different diameters and veneering materials: A 3-D finite element study

The aimof this study was to evaluate the stress distribution on bone tissue with a single prosthesis supported by implants of large and conventional diameter and presenting different veneering materials using the 3-D finite elementmethod. Sixteenmodels were fabricated to reproduce a bone block with implants, using two diameters (3.75 × 10 mmand 5.00 × 10 mm), four different veneering materials (composite resin, acrylic resin, porcelain, and NiCr crown), and two loads (axial (200 N) and oblique (100 N)). For data analysis, the maximum principal stress and vonMises criterion were used. For the axial load, the cortical bone in allmodels did not exhibit significant differences, and the trabecular bone presented higher tensile stresswith reduced implant diameter. For the oblique load, the cortical bone presented a significant increase in tensile stress on the same side as the loading for smaller implant diameters. The trabecular bone showed a similar but more discreet trend. There was no difference in bone tissue with different veneering materials. The veneering material did not influence the stress distribution in the supporting tissues of single implant-supported prostheses. The large-diameter implants improved the transference of occlusal loads to bone tissue and decreased stress mainly under oblique loads.Oblique loading was more detrimental to distribution stresses than axial loading. © 2013 Elsevier B.V. All rights reserved.

Aproximações hiperbólicas dos tempos de trânsito de raios paraxiais a um raio central refletido e difratado

A simulação de uma seção sísmica de afastamento nulo (AN) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhor a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda, recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Também, seguindo este novo conceito de imageamento sísmico está surgindo um método para simular seções com afastamento comum (AC) a partir de dados de cobertura múltipla usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento finito. Esta nova aproximação dos tempos de trânsito depende de cinco atributos cinemáticos. Neste trabalho, a partir da aproximação dos tempos de trânsito paraxiais em relação a um raio central com afastamento finito, derivamos uma nova equação do tempo de trânsito usando a condição de um ponto difrator em profundidade, reduzindo a equação original para quatro parâmetros. Para ambas aproximações (reflexão e difração), mostramos a superfície de empilhamento SRC com afastamento finito (SRC-AF). Considerando um modelo sintético, realizamos um estudo comparativo das aproximações dos tempos de trânsito para as quatro configurações sísmicas (fonte comum (FC), receptor comum (RC), ponto-médio-comum (PMC) e afastamento comum (AC)). Para analisar o comportamento do operador SRC-AF, quando este é perturbado, discutimos sua sensibilidade em relação a cada um dos cinco parâmetros (K₁, K₂, K₃, β_S e β_G). Esta análise de sensibilidade é realizada em duas formas: Sensibilidade através da primeira derivada e Sensibilidade no Empilhamento SRC-AF. Após realizar a análise de sensibilidade utilizamos uma nova condição, K₂ = 0 e assim, obtemos uma nova aproximação, agora dependente de três parâmetros. Usando essas aproximações hiperbólicas (em função de cinco, quatro e três parâmetros), propomos um algoritmo para a simulação de seções AC a partir de dados de cobertura múltipla. Finalmente, é apresentado um estudo da zona de Fresnel, com o objetivo de determinar a delimitação da abertura da superfície de empilhamento SRC-AF.

Migração usando equação da onda unidirecional de abertura arbitrária

Neste trabalho avaliamos uma classe de operadores de continuação de campos de onda, baseados em equações unidirecionais e com aplicação direta à migração sísmica. O método de representação de equações de onda unidirecionais, desenvolvido neste trabalho, é válido para abertura angular arbitrária, baseia-se no conceito de rigidez de um semiespaço, na transformação Dirichlet-Neumann e em sua discretização por elementos finitos. O método de construção dos operadores de continuação requer a introdução de variáveis auxiliares cujo número cresce em função da maior abertura angular desejada para o operador. Efetuamos a implementação no domínio do espaço e da frequência o que permite sua imediata paralelização. Baseados em experimentos numéricos, que avaliam a relação de dispersão e a resposta ao impulso do operador, propomos prescrições que permitem especificar o número de variáveis auxiliares e o passo de continuação para o operador de migração. A aplicação do algoritmo nos dados do modelo de domo salino da SEG-EAGE demonstra a capacidade do algoritmo em migrar refletores com forte mergulho em meios com forte variação lateral de velocidade.

Use of synthetic aperture radar for recognition of Coastal Geomorphological Features, land-use assessment and shoreline changes in Bragança coast, Pará, Northern Brazil

Imagens de radar de abertura sintética (SAR) vem sendo bem mais utilizadas do que antes nas aplicações de geociências em regiões tropicais úmidas. Nesta investigação, uma imagem RADARSAT-1, na banda C, polarização HH adquirida em 1998 foi usada para o mapeamento costeiro e avaliação da cobertura da terra na área de Bragança, norte do Brasil. Imagem do radar aerotransportado GEMS-1000, na banda X, polarização HH, adquirida em 1972 durante o projeto RADAM foi também utilizada para avaliar as variações costeiras ocorridas nas últimas três décadas. A pesquisa tem confirmado a utilidade da imagem RADARSAT-1 para o mapeamento geomorfológico e avaliação da cobertura da terra, particularmente em costas de manguezal de macromaré. Além disso, um novo método para estimar as variações da linha de costa baseado na superposição de vetores extraídos de diferentes imagens SAR, com alta acurácia geométrica, tem mostrado que a planície costeira de Bragança tem estado sujeita a severa erosão responsável pelo recuo de aproximadamente 32 km² e acreção de 20 km², resultando em uma perda de área de manguezal de aproximadamente 12 km². Como perspectiva de aplicação, dados SAR orbitais e aerotransportados provaram ser uma importante fonte de informação tanto para o mapeamento geomorfológico, quando para o monitoramento de modificações costeiras em ambientes tropicais úmidos.