The application of interpolating MLS approximations to the analysis of MHD flows

The element-free Galerkin method (EFGM) is a very attractive technique for solutions of partial differential equations, since it makes use of nodal point configurations which do not require a mesh. Therefore, it differs from FEM-like approaches by avoiding the need of meshing, a very demanding task for complicated geometry problems. However, the imposition of boundary conditions is not straightforward, since the EFGM is based on moving-least-squares (MLS) approximations which are not necessarily interpolants. This feature requires, for instance, the introduction of modified functionals with additional unknown parameters such as Lagrange multipliers, a serious drawback which leads to poor conditionings of the matrix equations. In this paper, an interpolatory formulation for MLS approximants is presented: it allows the direct introduction of boundary conditions, reducing the processing time and improving the condition numbers. The formulation is applied to the study of two-dimensional magnetohydrodynamic flow problems, and the computed results confirm the accuracy and correctness of the proposed formulation. (C) 2002 Elsevier B.V. All rights reserved.

Polymerization shrinkage: effects of boundary conditions and filling technique of resin composite restorations

Objective: To evaluate the linear polymerization shrinkage (LPS) and the effect of polymerization shrinkage of a resin composite and resin-dentin bond strength under different boundary conditions and filling techniques.Methods: Two cavities (4 x 4 x 2 MM) were prepared in bovine incisors (n = 30). The teeth were divided into three groups, according to boundary conditions: In group TE, the total-etch technique was used. In group EE, only enamel was conditioned, and in group NE, none of the watts of the cavities were conditioned. A two-step adhesive system was applied to all cavities. The resin composite was inserted in one (B) or three increments (1), and tight-cured with 600 mW/cm(2) (80 s). The LPS (%) was measured in the top-bottom direction, by placing a probe in contact with resin composite during curing. Enamel and total mean gap widths were measured (400 x) in three slices obtained after sectioning the restorations. Then, the slices were sectioned again, either to obtain sticks from the adhesive interface from the bottom of the cavity or to obtain resin composite sticks (0.8 mm(2)) to be tested for tensile strength (Kratos machine, 0.5 mm/min). The data was subjected to a two-way repeated measures ANOVA and Tukey's test for comparison of the means (alpha = 0.05).Results: the highest percentage of LPS was found for the TE when bulk fitted, and the lowest percentage of LPS was found in the Hand NE when incrementally fitted. The resin dentin bond strength was higher and the total mean gap width was tower for TE group; no significant effect was detected for the main factor fitting techniques. No difference was detected for the tensile strength of resin composite among the experimental groups.Conclusions: the filling technique is not able to minimize effects of the polymerization shrinkage, and bonding to the cavity watts is necessary to assure reduced mean gap width and high bond strength values. (C) 2004 Elsevier Ltd. All rights reserved.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

Temperature dependence of fractal dimension of grain boundary region in SnO2 based ceramics

Fractal dimensions of grain boundary region in doped SnO2 ceramics were determined based on previously derived fractal model. This model considers fractal dimension as a measure of homogeneity of distribution of charge carriers. Application of the derived fractal model enables calculation of fractal dimension using results of impedance spectroscopy. The model was verified by experimentally determined temperature dependence of the fractal dimension of SnO2 ceramics. Obtained results confirm that the non-Debye response of the grain boundary region is connected with distribution of defects and consequently with a homogeneity of a distribution of the charge carriers. Also, it was found that C-T-1 function has maximum at temperature at which the change in dominant type of defects takes place. This effect could be considered as a third-order transition.

Exact boundary control for the wave equation in a polyhedral time-dependent domain

We establish exact boundary controllability for the wave equation in a polyhedral domain where a part of the boundary moves slowly with constant speed in a small interval of time. The control on the moving part of the boundary is given by the conormal derivative associated with the wave operator while in the fixed part the control is of Neuman type. For initial state H-1 x L-2 we obtain controls in L-2. (C) 1999 Elsevier B.V. Ltd. All rights reserved.

Linear analysis of building floor structures by a BEM formulation based on Reissner's theory

In this work, the plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis is extended to the analysis of zoned plates in order to model a building floor structure. In the proposed formulation each sub-region defines a beam or a slab and depending on the way the sub-regions are represented, one can have two different types of analysis. In the simple bending problem all sub-regions are defined by their middle surface. on the other hand, for the coupled stretching-bending problem all sub-regions are referred to a chosen reference surface, therefore eccentricity effects are taken into account. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. The bending and stretching values defined on the interfaces are approximated along the beam width, reducing therefore the number of degrees of freedom. Then, in the proposed model the set of equations is written in terms of the problem values on the beam axis and on the external boundary without beams. Finally some numerical examples are presented to show the accuracy of the proposed model.

Ferroelectric state analysis in grain boundary of Na0.85Li0.15NbO3 ceramic

The electric and dielectric properties of the grain boundary of Na0.85Li0.15NbO3 lead-free ferroelectric-semiconductor perovskite were investigated. The impedance spectroscopy was carried out as a function of a thermal cycle. The sodium lithium niobate was synthesized by a chemical route based on the evaporation method. Dense ceramic, relative density of 97%, was prepared at 1423 K for 2 h in air atmosphere. ac measurements were carried out in the frequency range of 5 Hz-13 MHz and from 673 to 1023 K. Theoretical adjust of the impedance data was performed to deriving the electric parameters of the grain boundary. The electric conductivity follows the Arrhenius law, with activation energy values equal to 1.55 and 1.54 eV for heating and cooling cycle, respectively. The nonferroelectric state of the grain boundary and its correlation with symmetry are discussed in the temperature domain. (C) 2003 American Institute of Physics.

Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.

ANALYSIS OF CURVED EDGE PLATES BY THE BOUNDARY-ELEMENT METHOD

The aim of this work is to present a formulation of the boundary element method to analyse elastic and isotropic plates with curved boundaries. In this study the plate boundary is approximated, along each element, by a second degree polynomial relation or by a circular arch, in order to better represent the real boundary. The numerical integration is performed by the self-adaptive coordinate transformation proposed by Telles. The effective shear forces are approximated by concentrated reactions applied at the boundary element nodes, according to the alternative formulation introduced by Paiva. Some examples are presented to demonstrate the better accuracy obtained with the proposed elements.

INFLUENCE of THE EVAPORATOR POSITION IN STATIC DOMESTIC REFRIGERATORS

The temperature and velocity distributions of the air inside the cabinet of domestic refrigerators affect the quality of food products. If the consumer knows the location of warm and cold zones in the refrigerator, the products can be placed in the right zone. In addition, the knowledge of the thickness of thermal and hydrodynamic boundary layers near the evaporator and the other walls is also important. If the product is too close to the evaporator wall, freezing can occur, and if it is too close to warm walls, the products can be deteriorated. The aim of the present work is to develop a steady state computational fluid dynamics (CFD) model for domestic refrigerators working on natural convection regime. The Finite Volume Methodology is chosen as numerical procedure for discretizing the governing equations. The SIMPLE-Semi-Implicit Method for Pressure-Linked Equations algorithm applied to a staggered mesh was used for solving the pressure-velocity coupling problem. The Power-Law scheme is employed as interpolation function for the convective-diffusive terms, and the TDMA-Tri-Diagonal Matrix Algorithm is used to solve the systems of algebraic equations. The model is applied to a commercial static refrigerator, where the cabinet is considered an empty three-dimensional rectangular cavity with one drawer at the bottom of the cabinet, but without shelves. In order to analyze the velocity and temperature fields of the air flow inside the cabinet the evaporator temperature, Te, was varied from -20 degrees C to 0 degrees C, and nine different evaporator positions are evaluated for evaporator temperature of -15 degrees C. The cooling capacity of the evaporator for the steady state regime is also computed for each case. One can conclude that the vertical positioning of the evaporator inside the cabinet plays an important role on the temperature distribution inside the cabinet.

Non-linear analysis of reinforced concrete plates by the boundary element method and continuum damage mechanics

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

Fine structure of the boundary tissue of the seminiferous tubuli of the vampire bat (Desmodus rotundus, G.; Microchiroptera).

Structurally the boundary tissue of the vampire bat seminiferous tubuli showed 2 to 5 layers of connective tissue in which elongated contractile cells and lamellar and/or fibrillar collagen were noticed. This boundary tissue forms the seminiferous tubular lamina propria. Its structure was more complex around the seminiferous tubuli near the Capsula testicularis than between the adjacent and contiguous tubuli into the testicular lobuli. The whole ultrastructural organization of the seminiferous lamina propria was described here.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.