On asymptotic behavior at infinity and the finite section method for integral equations on the half-line

e consider integral equations on the half-line of the form and the finite section approximation to x obtained by replacing the infinite limit of integration by the finite limit β. We establish conditions under which, if the finite section method is stable for the original integral equation (i.e. exists and is uniformly bounded in the space of bounded continuous functions for all sufficiently large β), then it is stable also for a perturbed equation in which the kernel k is replaced by k + h. The class of perturbations allowed includes all compact and some non-compact perturbations of the integral operator. Using this result we study the stability and convergence of the finite section method in the space of continuous functions x for which ()()()=−∫∞dttxt,sk)s(x0()syβxβx()sxsp+1 is bounded. With the additional assumption that ()(tskt,sk−≤ where ()()(),qsomefor,sassOskandRLkq11>+∞→=∈− we show that the finite-section method is stable in the weighted space for ,qp≤≤0 provided it is stable on the space of bounded continuous functions. With these results we establish error bounds in weighted spaces for x - xβ and precise information on the asymptotic behaviour at infinity of x. We consider in particular the case when the integral operator is a perturbation of a Wiener-Hopf operator and illustrate this case with a Wiener-Hopf integral equation arising in acoustics.

Gradient recovery in adaptive finite-element methods for parabolic problems

We derive energy-norm a posteriori error bounds, using gradient recovery (ZZ) estimators to control the spatial error, for fully discrete schemes for the linear heat equation. This appears to be the �rst completely rigorous derivation of ZZ estimators for fully discrete schemes for evolution problems, without any restrictive assumption on the timestep size. An essential tool for the analysis is the elliptic reconstruction technique.Our theoretical results are backed with extensive numerical experimentation aimed at (a) testing the practical sharpness and asymptotic behaviour of the error estimator against the error, and (b) deriving an adaptive method based on our estimators. An extra novelty provided is an implementation of a coarsening error "preindicator", with a complete implementation guide in ALBERTA in the appendix.

A frequency-independent boundary element method for scattering by two-dimensional screens and apertures

We propose and analyse a hybrid numerical–asymptotic hp boundary element method (BEM) for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high-frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom N increases, and that to achieve any desired accuracy it is sufficient to increase N in proportion to the square of the logarithm of the frequency as the frequency increases (standard BEMs require N to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary ‘breakwater’ problem of propagation through an aperture in an infinite sound-hard screen.

High Royds: an integrated, analytical approach for mapping the unmarked burials of a pauper cemetery

Applying geophysical techniques to detect and map the physical extent of individual unmarked graves proves difficult in many cases. The success of individual geophysical techniques for detecting unmarked graves may be due to a poor understanding of the nature of the graves themselves, the context in which they lie in, and temporal changes to the burial state. Given the unpredictability of these variables, it is surprising that grave prospection is often undertaken using only a single method. This paper presents a multi-methodological survey strategy for detecting unmarked burials and utilises an analytical approach for visualising and evaluating survey results.

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.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

Financial time series analysis : Chaos and neurodynamics approach

This work aims at combining the Chaos theory postulates and Artificial Neural Networks classification and predictive capability, in the field of financial time series prediction. Chaos theory, provides valuable qualitative and quantitative tools to decide on the predictability of a chaotic system. Quantitative measurements based on Chaos theory, are used, to decide a-priori whether a time series, or a portion of a time series is predictable, while Chaos theory based qualitative tools are used to provide further observations and analysis on the predictability, in cases where measurements provide negative answers. Phase space reconstruction is achieved by time delay embedding resulting in multiple embedded vectors. The cognitive approach suggested, is inspired by the capability of some chartists to predict the direction of an index by looking at the price time series. Thus, in this work, the calculation of the embedding dimension and the separation, in Takens‘ embedding theorem for phase space reconstruction, is not limited to False Nearest Neighbor, Differential Entropy or other specific method, rather, this work is interested in all embedding dimensions and separations that are regarded as different ways of looking at a time series by different chartists, based on their expectations. Prior to the prediction, the embedded vectors of the phase space are classified with Fuzzy-ART, then, for each class a back propagation Neural Network is trained to predict the last element of each vector, whereas all previous elements of a vector are used as features.

Development of a new analytical approach based on cellulose membrane and chelator for differentiation of labile and inert metal species in aquatic systems

A new procedure was developed in this study, based on a system equipped with a cellulose membrane and a tetraethylenepentamine hexaacetate chelator (MD-TEPHA) for in situ characterization of the lability of metal species in aquatic systems. To this end, the DM-TEPHA system was prepared by adding TEPHA chelator to cellulose bags pre-purified with 1.0 mol L-1 of HCl and NaOH solutions. After the MD-TEPHA system was sealed, it was examined in the laboratory to evaluate the influence of complexation time (0-24 h), pH (3.0, 4.0, 5.0, 6.0 and 7.0), metal ions (Cu, Cd, Fe, Mn and Ni) and concentration of organic matter (15, 30 and 60 mg L-1) on the relative lability of metal species by TEPHA chelator. The results showed that Fe and Cu metals were complexed more slowly by TEPHA chelator in the MD-TEPHA system than were Cd, Ni and Mn in all pH used. It was also found that the pH strongly influences the process of metal complexation by the MD-TEPHA system. At all the pH levels, Cd, Mn and Ni showed greater complexation with TEPHA chelator (recovery of about 95-75%) than did Cu and Fe metals. Time also affects the lability of metal species complexed by aquatic humic substances (AHS); while Cd, Ni and Mn showed a faster kinetics, reaching equilibrium after about 100 min, and Cu and Fe approached equilibrium after 400 min. Increasing the AHS concentration decreases the lability of metal species by shifting the equilibrium to AHS-metal complexes. Our results indicate that the system under study offers an interesting alternative that can be applied to in situ experiments for differentiation of labile and inert metal species in aquatic systems. (c) 2006 Elsevier B.V. All rights reserved.

Otimização de forma aplicando B-splines sob critério integral de tensões

This work proposes a computational methodology to solve problems of optimization in structural design. The application develops, implements and integrates methods for structural analysis, geometric modeling, design sensitivity analysis and optimization. So, the optimum design problem is particularized for plane stress case, with the objective to minimize the structural mass subject to a stress criterion. Notice that, these constraints must be evaluated at a series of discrete points, whose distribution should be dense enough in order to minimize the chance of any significant constraint violation between specified points. Therefore, the local stress constraints are transformed into a global stress measure reducing the computational cost in deriving the optimal shape design. The problem is approximated by Finite Element Method using Lagrangian triangular elements with six nodes, and use a automatic mesh generation with a mesh quality criterion of geometric element. The geometric modeling, i.e., the contour is defined by parametric curves of type B-splines, these curves hold suitable characteristics to implement the Shape Optimization Method, that uses the key points like design variables to determine the solution of minimum problem. A reliable tool for design sensitivity analysis is a prerequisite for performing interactive structural design, synthesis and optimization. General expressions for design sensitivity analysis are derived with respect to key points of B-splines. The method of design sensitivity analysis used is the adjoin approach and the analytical method. The formulation of the optimization problem applies the Augmented Lagrangian Method, which convert an optimization problem constrained problem in an unconstrained. The solution of the Augmented Lagrangian function is achieved by determining the analysis of sensitivity. Therefore, the optimization problem reduces to the solution of a sequence of problems with lateral limits constraints, which is solved by the Memoryless Quasi-Newton Method It is demonstrated by several examples that this new approach of analytical design sensitivity analysis of integrated shape design optimization with a global stress criterion purpose is computationally efficient

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

Finite Element Modeling for Development and Optimization of a Bone Plate for Mandibular Fracture in Dogs

This study aimed to develop a plate to treat fractures of the mandibular body in dogs and to validate the project using finite elements and biomechanical essays. Mandible prototypes were produced with 10 oblique ventrorostral fractures (favorable) and 10 oblique ventrocaudal fractures (unfavorable). Three groups were established for each fracture type. Osteosynthesis with a pure titanium plate of double-arch geometry and blocked monocortical screws offree angulanon were used. The mechanical resistance of the prototype with unfavorable fracture was lower than that of the fcworable fracture. In both fractures, the deflection increased and the relative stiffness decreased proportionally to the diminishing screw number The finite element analysis validated this plate study, since the maximum tension concentration observed on the plate was lower than the resistance limit tension admitted by the titanium. In conclusion, the double-arch geometry plate fixed with blocked monocortical screws has sufficient resistance to stabilize oblique,fractures, without compromising mandibular dental or neurovascular structures. J Vet Dent 24 (7); 212 - 221, 2010

Finite Element Analysis to Compare Complete Denture and Implant-Retained Overdentures With Different Attachment Systems

This finite element analysis compared stress distribution on complete dentures and implant-retained overdentures with different attachment systems. Four models of edentulous mandible were constructed: group A (control), complete denture; group B, overdenture retained by 2 splinted implants with bar-clip system; group C, overdenture retained by 2 unsplinted implants with o'ring system; and group D, overdenture retained by 2 splinted implants with bar-clip and 2 distally placed o'ring system. Evaluation was performed on Ansys software, with 100-N vertical load applied on central incisive teeth. The lowest maximum general stress value (in megapascal) was observed in group A (64.305) followed by groups C (119.006), D (258.650), and B (349.873). The same trend occurred it) supporting tissues with the highest stress value for cortical bone. Unsplinted implants associated with the o'ring attachment system showed the lowest maximum stress values among all overdenture groups. Furthermore, o'ring system also improved stress distribution when associated with bar-clip system.

Three-Dimensional Finite Element Analysis of Vertical and Angular Misfit in Implant-Supported Fixed Prostheses

Purpose: Three-dimensional finite element analysis was used to evaluate the effect of vertical and angular misfit in three-piece implant-supported screw-retained fixed prostheses on the biomechanical response in the peri-implant bone, implants, and prosthetic components. Materials and Methods: Four three-dimensional models were fabricated to represent a right posterior mandibular section with one implant in the region of the second premolar (2PM) and another in the region of the second molar (2M). The implants were splinted by a three-piece implant-supported metal-ceramic prosthesis and differed according to the type of misfit, as represented by four different models: Control = prosthesis with complete fit to the implants; UAM (unilateral angular misfit) = prosthesis presenting unilateral angular misfit of 100 pm in the mesial region of the 2M; UVM (unilateral vertical misfit) = prosthesis presenting unilateral vertical misfit of 100 pm in the mesial region of the 2M; and TVM (total vertical misfit) = prosthesis presenting total vertical misfit of 100 pm in the platform of the framework in the 2M. A vertical load of 400 N was distributed and applied on 12 centric points by the software Ansys, ie, a vertical load of 150 N was applied to each molar in the prosthesis and a vertical load of 100 N was applied at the 2PM. Results: The stress values and distribution in peri-implant bone tissue were similar for all groups. The models with misfit exhibited different distribution patterns and increased stress magnitude in comparison to the control. The highest stress values in group UAM were observed in the implant body and retention screw. The groups UVM and TVM exhibited high stress values in the platform of the framework and the implant hexagon, respectively. Conclusions: The three types of misfit influenced the magnitude and distribution of stresses. The influence of misfit on peri-implant bone tissue was modest. Each type of misfit increased the stress values in different regions of the system. INT J ORAL MAXILLOFAC IMPLANTS 2011;26:788-796

Biomechanical evaluation of internal and external hexagon platform switched implant-abutment connections: An in vitro laboratory and three-dimensional finite element analysis

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