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.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

Numerical analysis of different methods to calculate residual stresses in thin films based on instrumented indentation data

In this work, different methods to estimate the value of thin film residual stresses using instrumented indentation data were analyzed. This study considered procedures proposed in the literature, as well as a modification on one of these methods and a new approach based on the effect of residual stress on the value of hardness calculated via the Oliver and Pharr method. The analysis of these methods was centered on an axisymmetric two-dimensional finite element model, which was developed to simulate instrumented indentation testing of thin ceramic films deposited onto hard steel substrates. Simulations were conducted varying the level of film residual stress, film strain hardening exponent, film yield strength, and film Poisson's ratio. Different ratios of maximum penetration depth h(max) over film thickness t were also considered, including h/t = 0.04, for which the contribution of the substrate in the mechanical response of the system is not significant. Residual stresses were then calculated following the procedures mentioned above and compared with the values used as input in the numerical simulations. In general, results indicate the difference that each method provides with respect to the input values depends on the conditions studied. The method by Suresh and Giannakopoulos consistently overestimated the values when stresses were compressive. The method provided by Wang et al. has shown less dependence on h/t than the others.

Solvability near the characteristic set for a special class of complex vector fields

This work deals with the solvability near the characteristic set Sigma = {0} x S-1 of operators of the form L = partial derivative/partial derivative t+(x(n) a(x)+ ix(m) b(x))partial derivative/partial derivative x, b not equivalent to 0 and a(0) not equal 0, defined on Omega(epsilon) = (-epsilon, epsilon) x S-1, epsilon > 0, where a and b are real-valued smooth functions in (-epsilon, epsilon) and m >= 2n. It is shown that given f belonging to a subspace of finite codimension of C-infinity (Omega(epsilon)) there is a solution u is an element of L-infinity of the equation Lu = f in a neighborhood of Sigma; moreover, the L-infinity regularity is sharp.

Partially Observed Markov Random Fields Are Variable Neighborhood Random Fields

The present paper has two goals. First to present a natural example of a new class of random fields which are the variable neighborhood random fields. The example we consider is a partially observed nearest neighbor binary Markov random field. The second goal is to establish sufficient conditions ensuring that the variable neighborhoods are almost surely finite. We discuss the relationship between the almost sure finiteness of the interaction neighborhoods and the presence/absence of phase transition of the underlying Markov random field. In the case where the underlying random field has no phase transition we show that the finiteness of neighborhoods depends on a specific relation between the noise level and the minimum values of the one-point specification of the Markov random field. The case in which there is phase transition is addressed in the frame of the ferromagnetic Ising model. We prove that the existence of infinite interaction neighborhoods depends on the phase.

Carbon monoxide and related trace gases and aerosols over the Amazon Basin during the wet and dry seasons

We present the results of airborne measurements of carbon monoxide (CO) and aerosol particle number concentration (CN) made during the Balan double dagger o Atmosf,rico Regional de Carbono na Amazonia (BARCA) program. The primary goal of BARCA is to address the question of basin-scale sources and sinks of CO2 and other atmospheric carbon species, a central issue of the Large-scale Biosphere-Atmosphere (LBA) program. The experiment consisted of two aircraft campaigns during November-December 2008 (BARCA-A) and May-June 2009 (BARCA-B), which covered the altitude range from the surface up to about 4500 m, and spanned most of the Amazon Basin. Based on meteorological analysis and measurements of the tracer, SF6, we found that airmasses over the Amazon Basin during the late dry season (BARCA-A, November 2008) originated predominantly from the Southern Hemisphere, while during the late wet season (BARCA-B, May 2009) low-level airmasses were dominated by northern-hemispheric inflow and mid-tropospheric airmasses were of mixed origin. In BARCA-A we found strong influence of biomass burning emissions on the composition of the atmosphere over much of the Amazon Basin, with CO enhancements up to 300 ppb and CN concentrations approaching 10 000 cm(-3); the highest values were in the southern part of the Basin at altitudes of 1-3 km. The Delta CN/Delta CO ratios were diagnostic for biomass burning emissions, and were lower in aged than in fresh smoke. Fresh emissions indicated CO/CO2 and CN/CO emission ratios in good agreement with previous work, but our results also highlight the need to consider the residual smoldering combustion that takes place after the active flaming phase of deforestation fires. During the late wet season, in contrast, there was little evidence for a significant presence of biomass smoke. Low CN concentrations (300-500 cm(-3)) prevailed basinwide, and CO mixing ratios were enhanced by only similar to 10 ppb above the mixing line between Northern and Southern Hemisphere air. There was no detectable trend in CO with distance from the coast, but there was a small enhancement of CO in the boundary layer suggesting diffuse biogenic sources from photochemical degradation of biogenic volatile organic compounds or direct biological emission. Simulations of CO distributions during BARCA-A using a range of models yielded general agreement in spatial distribution and confirm the important contribution from biomass burning emissions, but the models evidence some systematic quantitative differences compared to observed CO concentrations. These mismatches appear to be related to problems with the accuracy of the global background fields, the role of vertical transport and biomass smoke injection height, the choice of model resolution, and reliability and temporal resolution of the emissions data base.

Local symmetries in gauge theories in a finite-dimensional setting

It is shown that the correct mathematical implementation of symmetry in the geometric formulation of classical field theory leads naturally beyond the concept of Lie groups and their actions on manifolds, out into the realm of Lie group bundles and, more generally, of Lie groupoids and their actions on fiber bundles. This applies not only to local symmetries, which lie at the heart of gauge theories, but is already true even for global symmetries when one allows for fields that are sections of bundles with (possibly) non-trivial topology or, even when these are topologically trivial, in the absence of a preferred trivialization. (C) 2012 Elsevier B.V. All rights reserved.

Gauge and integrable theories in loop spaces

We propose an integral formulation of the equations of motion of a large class of field theories which leads in a quite natural and direct way to the construction of conservation laws. The approach is based on generalized non-abelian Stokes theorems for p-form connections, and its appropriate mathematical language is that of loop spaces. The equations of motion are written as the equality of a hyper-volume ordered integral to a hyper-surface ordered integral on the border of that hyper-volume. The approach applies to integrable field theories in (1 + 1) dimensions, Chern-Simons theories in (2 + 1) dimensions, and non-abelian gauge theories in (2 + 1) and (3 + 1) dimensions. The results presented in this paper are relevant for the understanding of global properties of those theories. As a special byproduct we solve a long standing problem in (3 + 1)-dimensional Yang-Mills theory, namely the construction of conserved charges, valid for any solution, which are invariant under arbitrary gauge transformations. (C) 2012 Elsevier B.V. All rights reserved.

Quantum motion in superposition of Aharonov-Bohm with some additional electromagnetic fields

The structure of additional electromagnetic fields to the Aharonov-Bohm field, for which the Schrodinger, Klein-Gordon, and Dirac equations can be solved exactly are described and the corresponding exact solutions are found. It is demonstrated that aside from the known cases (a constant and uniform magnetic field that is parallel to the Aharonov-Bohm solenoid, a static spherically symmetrical electric field, and the field of a magnetic monopole), there are broad classes of additional fields. Among these new additional fields we have physically interesting electric fields acting during a finite time or localized in a restricted region of space. There are additional time-dependent uniform and isotropic electric fields that allow exact solutions of the Schrodinger equation. In the relativistic case there are additional electric fields propagating along the Aharonov-Bohm solenoid with arbitrary electric pulse shape. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4714352]

Application of the log-conformation tensor to three-dimensional time-dependent free surface flows

The numerical simulation of flows of highly elastic fluids has been the subject of intense research over the past decades with important industrial applications. Therefore, many efforts have been made to improve the convergence capabilities of the numerical methods employed to simulate viscoelastic fluid flows. An important contribution for the solution of the High-Weissenberg Number Problem has been presented by Fattal and Kupferman [J. Non-Newton. Fluid. Mech. 123 (2004) 281-285] who developed the matrix-logarithm of the conformation tensor technique, henceforth called log-conformation tensor. Its advantage is a better approximation of the large growth of the stress tensor that occur in some regions of the flow and it is doubly beneficial in that it ensures physically correct stress fields, allowing converged computations at high Weissenberg number flows. In this work we investigate the application of the log-conformation tensor to three-dimensional unsteady free surface flows. The log-conformation tensor formulation was applied to solve the Upper-Convected Maxwell (UCM) constitutive equation while the momentum equation was solved using a finite difference Marker-and-Cell type method. The resulting developed code is validated by comparing the log-conformation results with the analytic solution for fully developed pipe flows. To illustrate the stability of the log-conformation tensor approach in solving three-dimensional free surface flows, results from the simulation of the extrudate swell and jet buckling phenomena of UCM fluids at high Weissenberg numbers are presented. (C) 2012 Elsevier B.V. All rights reserved.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s-1 for 9.1 mu s. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.130.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

Vortices in the extended Skyrme-Faddeev model

We construct analytical and numerical vortex solutions for an extended Skyrme-Faddeev model in a (3 + 1) dimensional Minkowski space-time. The extension is obtained by adding to the Lagrangian a quartic term, which is the square of the kinetic term, and a potential which breaks the SO(3) symmetry down to SO(2). The construction makes use of an ansatz, invariant under the joint action of the internal SO(2) and three commuting U(1) subgroups of the Poincare group, and which reduces the equations of motion to an ordinary differential equation for a profile function depending on the distance to the x(3) axis. The vortices have finite energy per unit length, and have waves propagating along them with the speed of light. The analytical vortices are obtained for a special choice of potentials, and the numerical ones are constructed using the successive over relaxation method for more general potentials. The spectrum of solutions is analyzed in detail, especially its dependence upon special combinations of coupling constants.

delta-Superderivations of Semisimple Finite-Dimensional Jordan Superalgebras

In the paper, a complete description of the delta-derivations and the delta-superderivations of semisimple finite-dimensional Jordan superalgebras over an algebraically closed field of characteristic p not equal 2 is given. In particular, new examples of nontrivial (1/2)-derivations and odd (1/2)-superderivations are given that are not operators of right multiplication by an element of the superalgebra.

Causal amplitudes in the Schwinger model at finite temperature

We show, in the imaginary time formalism, that the temperature dependent parts of all the retarded (advanced) amplitudes vanish in the Schwinger model. We trace this behavior to the CPT invariance of the theory and give a physical interpretation of this result in terms of forward scattering amplitudes of on-shell thermal particles.