Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

A method for the calculation of the adsorbed phase volume and pseudo-saturation pressure from adsorption isotherm data on activated carbon

We propose a new method for evaluating the adsorbed phase volume during physisorption of several gases on activated carbon specimens. We treat the adsorbed phase as another equilibrium phase which satisfies the Gibbs equation and hence assume that the law of rectilinear diameters is applicable. Since invariably the bulk gas phase densities are known along measured isotherms, the constants of the adsorbed phase volume can be regressed from the experimental data. We take the Dubinin-Astakhov isotherm as the model for verifying our hypothesis since it is one of the few equations that accounts for adsorbed phase volume changes. In addition, the pseudo-saturation pressure in the supercritical region is calculated by letting the index of the temperature term in Dubinin's equation to be temperature dependent. Based on over 50 combinations of activated carbons and adsorbates (nitrogen, oxygen, argon, carbon dioxide, hydrocarbons and halocarbon refrigerants) it is observed that the proposed changes fit experimental data quite well.

Structural vibration control using resistively shunted piezoceramics

Application of piezoceramic materials in actuation and sensing of vibration is of current interest. Potential and more popular applications of piezoceramics are probably in the field of active vibration control. However, the objective of this work is to investigate the effect of shunted piezoceramics as passive vibration control devices when bonded to a host structure. Resistive shunting of a piezoceramic bonded to a cantilevered duralumin beam has been investigated. The piezoceramic is connected in parallel to an electrical network comprising of resistors and inductors. The piezoceramic is a capacitor that stores and discharges electrical energy that is transformed from the mechanical motion of the structure to which it is bonded. A resistor across the piezoceramic would be termed as a resistively shunted piezoceramic. Similarly, an inductor across the piezoceramic is termed as a resonantly shunted piezoceramic. In this study, the effect of resistive shunting on the nature of damping enhancement to the host structure has been investigated. Analytical studies are presented along with experimental results.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.

Influence of matrix characteristics on fracture toughness of high volume fraction Al2O3/Al–AlN composites

The role of matrix microstructure on the fracture of Al-alloy composites with 60 vol% alumina particulates was studied. The matrix composition and microstructure were systematically varied by changing the infiltration temperature and heat treatment. Characterization was carried out by a combination of metallography, hardness measurements, and fracture studies conducted on compact tension specimens to study the fracture toughness and crack growth in the composites. The composites showed a rise in crack resistance with crack extension (R curves) due to bridges of intact matrix ligaments formed in the crack wake. The steady-state or plateau toughness reached upon stable crack growth was observed to be more sensitive to the process temperature rather than to the heat treatment. Fracture in the composites was predominantly by particle fracture, extensive deformation, and void nucleation in the matrix. Void nucleation occurred in the matrix in the as-solutionized and peak-aged conditions and preferentially near the interface in the underaged and overaged conditions. Micromechanical models based on crack bridging by intact ductile ligaments were modified by a plastic constraint factor from estimates of the plastic zone formed under indentations, and are shown to be adequate in predicting the steady-state toughness of the composite.

Fracture and R-curves in high volume fraction Al2O3/Al composites

Fracture toughness and fracture mechanisms in Al2O3/Al composites are described. The unique flexibility offered by pressureless infiltration of molten Al alloys into porous alumina preforms was utilized to investigate the effect of microstructural scale and matrix properties on the fracture toughness and the shape of the crack resistance curves (R-curves). The results indicate that the observed increment in toughness is due to crack bridging by intact matrix ligaments behind the crack tip. The deformation behavior of the matrix, which is shown to be dependent on the microstructural constraints, is the key parameter that influences both the steady-state toughness and the shape of the R-curves. Previously proposed models based on crack bridging by intact ductile particles in a ceramic matrix have been modified by the inclusion of an experimentally determined plastic constraint factor (P) that determines the deformation of the ductile phase and are shown to be adequate in predicting the toughness increment in the composites. Micromechanical models to predict the crack tip profile and the bridge lengths (L) correlate well with the observed behavior and indicate that the composites can be classified as (i) short-range toughened and (ii) long-range toughened on the basis of their microstructural characteristics.

Dynamics of pulled desorption with effects of excluded-volume interaction: The p-Laplacian diffusion equation and its exact solution

We analyze the dynamics of desorption of a polymer molecule which is pulled at one of its ends with force f, trying to desorb it. We assume a monomer to desorb when the pulling force on it exceeds a critical value f(c). We formulate an equation for the average position of the n-th monomer, which takes into account excluded-volume interaction through the blob-picture of a polymer under external constraints. The approach leads to a diffusion equation with a p-Laplacian for the propagation of the stretching along the chain. This has to be solved subject to a moving boundary condition. Interestingly, within this approach, the problem can be solved exactly in the trumpet, stem-flower and stem regimes. In the trumpet regime, we get tau = tau(0)n(d)(2), where n(d) is the number of monomers that have desorbed at the time tau. tau(0) is known only numerically, but for f close to f(c), it is found to be tau(0) similar to f(c)/(f(2/3) - f(c)(2/3)) If one used simple Rouse dynamics, this result would change to tau similar to f(c)n(d)(2)/(f - f(c)). In the other regimes too, one can find exact solution, and interestingly, in all regimes tau similar to n(d)(2). Copyright (C) EPLA, 2011

Efficient Algorithms for Combinatorial Auctions with Volume Discounts Arising in Web Service Composition

Web services are now a key ingredient of software services offered by software enterprises. Many standardized web services are now available as commodity offerings from web service providers. An important problem for a web service requester is the web service composition problem which involves selecting the right mix of web service offerings to execute an end-to-end business process. Web service offerings are now available in bundled form as composite web services and more recently, volume discounts are also on offer, based on the number of executions of web services requested. In this paper, we develop efficient algorithms for the web service composition problem in the presence of composite web service offerings and volume discounts. We model this problem as a combinatorial auction with volume discounts. We first develop efficient polynomial time algorithms when the end-to-end service involves a linear workflow of web services. Next we develop efficient polynomial time algorithms when the end-to-end service involves a tree workflow of web services.

Volume Change Behaviour of Soil Influenced with Sulfuric Acid

The paper brings out the role of calcium carbonate (CaCO3) on the volume change behaviour of natural black cotton soil with 1N sulfuric acid (H2SO4) as pore fluid. Natural black cotton soil contained predominantly montmorillonite [Ca0.2(Al,Mg)2Si4 O10 (OH)2 .4H2O] along with other minerals such as amesite [(Mg Fe)2 Al (Si Al)2 O5 (OH)4], kalsilite [KAlSiO4] and quartz [SiO2]. The calcitic soil, reacted with H2SO4 during consolidation testing, showed the presence of the new mineral yavapaiite [K Fe(SO4)2]. Consequently, the carbonate soil treated with 1N H2SO4 led to higher swell at seating load and more compression upon loading than the soil with no carbonate. The swelling increased with increase in the amount of carbonate present in the soil.

Volume Change behaviour in Calcitic Soil influenced with Sulphuric Acid using Artificial Neural Networks

This paper elucidates the methodology of applying artificial neural network model (ANNM) to predict the percent swell of calcitic soil in sulphuric acid solutions, a complex phenomenon involving many parameters. Swell data required for modelling is experimentally obtained using conventional oedometer tests under nominal surcharge. The phases in ANN include optimal design of architecture, operation and training of architecture. The designed optimal neural model (3-5-1) is a fully connected three layer feed forward network with symmetric sigmoid activation function and trained by the back propagation algorithm to minimize a quadratic error criterion.The used model requires parameters such as duration of interaction, calcite mineral content and acid concentration for prediction of swell. The observed strong correlation coefficient (R2 = 0.9979) between the values determined by the experiment and predicted using the developed model demonstrates that the network can provide answers to complex problems in geotechnical engineering.

Closed form solutions for element equilibrium and flexibility matrices of eight node rectangular plate bending element using integrated force method

Closed form solutions for equilibrium and flexibility matrices of the Mindlin-Reissner theory based eight-node rectangular plate bending element (MRP8) using integrated Force Method (IFM) are presented in this paper. Though these closed form solutions of equilibrium and flexibility matrices are applicable to plate bending problems with square/rectangular boundaries, they reduce the computational time significantly and give more exact solutions. Presented closed form solutions are validated by solving large number of standard square/rectangular plate bending benchmark problems for deflections and moments and the results are compared with those of similar displacement-based eight-node quadrilateral plate bending elements available in the literature. The results are also compared with the exact solutions.

Transient electric field distribution within the working volume of an NEMP simulator

Nuclear electro-magnetic pulse (NEMP) simulators which are used in the simulation of transient electromagnetic fields due to a high altitude nuclear detonation are generally excited with a double exponential high voltage pulse. This results in a current distribution on the wires of the simulator and hence a transient electric field in the working volume of the simulator where the test object is kept. It is found that for the simulator under study, the current distribution is non-uniform and so is the field distribution along the width of the simulator in the working volume. To make the current distribution uniform, several methods have been suggested and the results of these methods are analyzed and suitable conclusions are arrived at from those results.