Mesh Configuration Research Articles

FEM-BEM Coupling for the MHD Pipe Flow in an Exterior Region

In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM-BEM coupling approach. The coupled partial differential equations defined for the interior medium are transformed into homogenous modified Helmholtz equations. For the exterior medium on an infinite region, the Laplace equation is considered for the exterior magnetic field. Unknowns in the equations are also related with the corresponding Dirichlet and Neumann type coupled boundary conditions. Unknown values of the magnetic field on the boundary and for the exterior region are obtained by using BEM, and the unknown velocity and magnetic field inside the pipe are obtained by using SUPG type stabilized FEM. Computations are carried for very high values of magnetic Reynolds numbers Rm1, Reynolds number Re and magnetic pressure Rh of the fluid. The results show that using stabilized method enables us to get stable and accurate numerical approximations consistent with the physical configuration of the problem over rough mesh which also results a cheap computational cost.

Network on Chip for 3D Mesh Structure with Enhanced Security Algorithm in HDL Environment

Network on chip (NOC) architecture is an approach to develop large and complex systems on a single chip. In this work 3D mesh topological structures has been implemented with enhanced TACIT Security in HDL environment. The architecture supports physical and architectural level design integration. Basic communication mechanism between resources is envisioned to be packet switched message passing through the switches. For the node identification in 3D networks 8 × 8 × 8 Mesh configuration is taken. A new Algorithm is proposed for secured data transmission among nodes. TACIT algorithm is used for data encryption and decryption and applicable for n bit block size and key. Design is implemented in Xilinx 14.2 VHDL software, and functional simulation was carried out in Modelsim 10.1 b, student edition. Hardware parameters such as size cost and timings are extracted from the design code. General Terms TACIT Security, Simulation, Synthesis, Integrated Circuits, Intellectual Property

An efficient meshing approach for the calculation of notch stresses

For the application of the notch stress approach, local stresses at the notch root have to be calculated. This demands a fine discretization of the finite element model in the notch leading for complex structures to large finite element models which are difficult to handle and take a high amount of calculation time. This paper contains the detailed investigation of the required finite element mesh in the notch area for application of the notch stress approach. Based on small-scale specimens, various meshing rules are applied and compared to each other. As the most important factor for the efficiency of a finite element mesh, the element shape in the notch root was identified. It should be chosen according to the stress gradients: The higher the stress gradients are, the smaller the element edge length should be modeled. The shape of the subsequent elements inside the volume has no major influence on the accuracy of the calculation. Recommendations for practical use are given for various element types, including commonly used tetrahedron and hexahedron element with linear and quadratic shape function. Expected numerical errors for the chosen mesh configuration are indicated.

Two-phase numerical study of the flow field formed in water pump sump: influence of air entrainment

In a pump sump it is imperative that the amount of non-homogenous flow and entrained air be kept to a minimum. Free air-core vortex occurring at a water-intake pipe is an important problem encountered in hydraulic engineering. These vortices reduce pump performances, may have large effects on the operating conditions and lead to increase plant operating costs.This work is an extended study starting from 2006 in LML and published by ISSA and al. in 2008, 2009 and 2010. Several cases of sump configuration have been numerically investigated using two specific commercial codes and based on the initial geometry proposed by Constantinescu and Patel. Fluent and Star CCM+ codes are used in the previous studies. The results, obtained with a structured mesh, were strongly dependant on main geometrical sump configuration such as the suction pipe position, the submergence of the suction pipe on one hand and the turbulence model on the other hand. Part of the results showed a good agreement with experimental investigations already published. Experiments, conducted in order to select best positions of the suction pipe of a water-intake sump, gave qualitative results concerning flow disturbances in the pump-intake related to sump geometries and position of the pump intake. The purpose of this paper is to reproduce the flow pattern of experiments and to confirm the geometrical parameter that influences the flow structure in such a pump. The numerical model solves the Reynolds averaged Navier-Stokes (RANS) equations and VOF multiphase model. STAR CCM+ with an adapted mesh configuration using hexahedral mesh with prism layer near walls was used. Attempts have been made to calculate two phase unsteady flow for stronger mass flow rates and stronger submergence with low water level in order to be able to capture air entrainment. The results allow the knowledge of some limits of numerical models, of mass flow rates and of submergences for air entrainment. In the validation of this numerical model, emphasis was placed on the prediction of the number, location, size and strength of the various types of vortices coming from the free surface. Contours of vorticity at free surface, air cores, isoline of pressure surface were particularly examined for some cases. Streamlines issued from the free surface and the volume of fraction of air allows visualizing the air entrainment.

Stretchable and highly sensitive graphene-on-polymer strain sensors

The use of nanomaterials for strain sensors has attracted attention due to their unique electromechanical properties. However, nanomaterials have yet to overcome many technological obstacles and thus are not yet the preferred material for strain sensors. In this work, we investigated graphene woven fabrics (GWFs) for strain sensing. Different than graphene films, GWFs undergo significant changes in their polycrystalline structures along with high-density crack formation and propagation mechanically deformed. The electrical resistance of GWFs increases exponentially with tensile strain with gauge factors of ~103 under 2~6% strains and ~106 under higher strains that are the highest thus far reported, due to its woven mesh configuration and fracture behavior, making it an ideal structure for sensing tensile deformation by changes in strain. The main mechanism is investigated, resulting in a theoretical model that predicts very well the observed behavior.

Randomized Partially-Minimal Routing: Near-Optimal Oblivious Routing for 3-D Mesh Networks

The increasing viability of 3-D silicon integration technology has opened new opportunities for chip architecture innovations. One direction is in the extension of 2-D mesh-based tiled chip-multiprocessor architectures into three dimensions. This paper focuses on efficient routing algorithms for such 3-D mesh networks. Existing routing algorithms suffer from either poor worst-case throughput (DOR, ROMM) or poor latency (VAL). Although the minimal routing algorithm O1TURN proposed in already achieves near-optimal worst-case throughput for 2-D mesh networks, the optimality result does not extend to higher dimensions. For 3-D and higher dimensional meshes, the worst-case throughput of O1TURN degrades tremendously. The main contribution of this paper is a new oblivious routing algorithm for 3-D mesh networks called randomized partially-minimal (RPM) routing. RPM provably achieves optimal worst-case throughput for 3-D meshes when the network radix is even and within a factor of 1/k2 of optimal worst-case throughput when is odd. Finally, whereas VAL achieves optimal worst-case throughput at a penalty factor of 2 in average latency over DOR, RPM achieves (near) optimal worst-case throughput with a much smaller factor of 1.33. For practical asymmetric 3-D mesh configurations where the number of device layers are fewer than the number of tiles along the edge of a layer, the average latency of RPM reduces to just a factor of 1.11 to 1.19 of DOR. Additionally, a variant of RPM called randomized minimal first (RMF) routing is proposed, which leverages the inherent load-balancing properties of the network traffic to further reduce packet latency without compromising throughput.

A variational multiscale method with subscales on the element boundaries for the Helmholtz equation

SUMMARYIn this paper, we apply the variational multiscale method with subgrid scales on the element boundaries to the problem of solving the Helmholtz equation with low‐order finite elements. The expression for the subscales is obtained by imposing the continuity of fluxes across the interelement boundaries. The stabilization parameter is determined by performing a dispersion analysis, yielding the optimal values for the different discretizations and finite element mesh configurations. The performance of the method is compared with that of the standard Galerkin method and the classical Galerkin least‐squares method with very satisfactory results. Some numerical examples illustrate the behavior of the method. Copyright © 2012 John Wiley & Sons, Ltd.

Comparison of different higher order finite element schemes for the simulation of Lamb waves

Structural Health Monitoring (SHM) applications call for both efficient and powerful numerical tools to predict the behavior of ultrasonic guided waves. When considering waves in thin-walled structures, so called Lamb waves, conventional linear or quadratic pure displacement finite elements soon reach their limits. The spatial as well as temporal discretisation, required to obtain good quality results has to be very fine. This results in enormous computational costs (computational time and memory storage requirements) when ultrasonic wave propagation problems are solved in the time domain. To resolve this issue several higher order finite element methods with polynomial degrees p>2 are proposed. The objective of the current article is to develop such higher order schemes and to verify their capabilities with respect to accuracy and numerical performance. To the best of the authors’ knowledge such comparison has not been reported in literature, yet. Specifically, spectral elements based on Lagarange polynomials (SEM), p-elements using the normalized integrals of the Legendre polynomials (p-FEM) and isogeometric elements utilizing non-uniform rational B-splines (NURBS, N-FEM) are discussed in this paper. By solving a two-dimensional benchmark problem, their advantages and drawbacks with respect to Lamb wave propagation are highlighted. The results of the convergence studies are then used to derive guidelines for estimating the optimal element size for a given finite element type and polynomial degree template. These findings serve the purpose to determine the optimal mesh configuration a priori and thus, save a considerable amount of computational effort. The proposed guideline is then tested on a three-dimensional structure with a conical hole showing an excellent agreement with the predicted behaviour.

Multifunctional graphene woven fabrics

Tailoring and assembling graphene into functional macrostructures with well-defined configuration are key for many promising applications. We report on a graphene-based woven fabric (GWF) prepared by interlacing two sets of graphene micron-ribbons where the ribbons pass each other essentially at right angles. By using a woven copper mesh as the template, the GWF grown from chemical vapour deposition retains the network configuration of the copper mesh. Embedded into polymer matrices, it has significant flexibility and strength gains compared with CVD grown graphene films. The GWFs display both good dimensional stability in both the warp and the weft directions and the combination of film transparency and conductivity could be optimized by tuning the ribbon packing density. The GWF creates a platform to integrate a large variety of applications, e.g., composites, strain sensors and solar cells, by taking advantages of the special structure and properties of graphene.

Experience with using the LICON methodology for predicting long term creep behaviour in materials

The LICON methodology is a newly developed approach for predicting the lifetime of materials under creep loading condition. The LICON method has gained attention for application to newly developed materials, for which there is little existing creep-rupture data, and for their welded structures. The LICON method predicts long term uniaxial creep strength using the results from several short duration creep crack incubation tests in conjunction with the outcome of a mechanical analysis (e.g. involving finite element analysis, at least initially for a new material and/or multi-axial specimen geometry). In this study, the sensitivity of results from the LICON method to input conditions for the associated finite element analysis is investigated. Two important considerations for creep finite element analyses are the mesh configuration and the type of creep constitutive model used. In order to evaluate the sensitivity of the outcome of the approach to finite element mesh condition, twelve different configurations in terms of size, geometry and interpolating formulation are examined. Moreover, to examine the sensitivity of the outcome of the approach to the type of creep equation employed, three constitutive models have been evaluated. Finally, the sensitivity of the LICON approach to the evaluated conditions have been discussed with reference to the information contained in a comprehensive creep-rupture database available for a 1%CrMoV steel at 550 °C.

Selectivity of commercial compared to larger mesh and square mesh trawl codends for four fish species in the Aegean Sea

Summary Differences in size selectivity of commercial (40 mm diamond mesh, 40D), larger mesh (48 mm diamond mesh, 48D), and square mesh codends (40 mm square mesh, 40S) for hake (Merluccius merluccius), greater forkbeard (Phycis blennoides), blackbelly rosefish (Helicolenus dactylopterus dactylopterus) and four-spot megrim (Lepidorhombus boscii) were investigated in the Aegean Sea. The study was conducted using the covered codend method. Data were analysed taking between-haul variations into account. Results showed that changing from a 40D to 48D codend significantly improved mean L50 values, with increases of about 22% for hake, 8% for greater forkbeard, 20% for blackbelly rosefish, and 11% for four-spot megrim (P = 0.000). A change from diamond to square mesh configuration with the same 40 mm netting also significantly increased the mean L50 values, with 45% for hake, 36% for greater forkbeard, and 25% for blackbelly rosefish (P = 0.000). For four-spot megrim the mean L50 value was about 10% lower, but this difference was statistically insignificant (P > 0.05). Despite the increases in L50 values, this study concludes that the selectivity of the 48D and 40S codends is still not sufficient to release fish smaller than length at first maturity for these four species.

A higher order finite element including transverse normal strain for linear elastic composite plates with general lamination configurations

This paper describes a higher-order global–local theory for thermal/mechanical response of moderately thick laminated composites with general lamination configurations. In-plane displacement fields are constructed by superimposing the third-order local displacement field to the global cubic displacement field. To eliminate layer-dependent variables, interlaminar shear stress compatibility conditions have been employed, so that the number of variables involved in the proposed model is independent of the number of layers of laminates. Imposing shear stress free condition at the top and the bottom surfaces, derivatives of transverse displacement are eliminated from the displacement field, so that C 0 interpolation functions are only required for the finite element implementation. To assess the proposed model, the quadratic six-node C 0 triangular element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate. Comparing to various existing laminated plate models, it is found that simple C 0 finite elements with non-zero normal strain could produce more accurate displacement and stresses for thick multilayer composite plates subjected to thermal and mechanical loads. Finally, it is remarked that the proposed model is quite robust, such that the finite element results are not sensitive to the mesh configuration and can rapidly converge to 3-D elasticity solutions using regular or irregular meshes.

An Experimental and Numerical Study of Ballistic Impacts on a Turbine Casing Material at Varying Temperatures

An experimental and numerical study of ballistic impacts on steel plates at various temperatures (700 °C, 400 °C and room temperature) has been carried out. The motivation for this work is the blade-off event that may occur inside a jet engine turbine. However, as a first attempt to understand this complex loading process, a somewhat simpler approach is carried out in the present work. The material used in this study is the FV535 martensitic stainless steel, which is one of the most commonly used materials for turbine casings. Based on material test data, a Modified Johnson-Cook (MJC) model was calibrated for numerical simulations using the LS-DYNA explicit finite element code. To check the mesh size sensitivity, 2D axisymmetric finite element models with three different mesh sizes and configurations were used for the various temperatures. Two fixed meshes with 64 and 128 elements over the 2 mm thick plate and one mesh with 32 elements over the thickness with adaptive remeshing were used in the simulations. Both the formation of adiabatic shear bands in the perforation process and the modeling of the thermal softening effects at high temperatures have been found crucial in order to achieve good results.

Performance characterization of nonlinear optimization methods for mesh quality improvement

We characterize the performance of gradient- and Hessian-based optimization methods for mesh quality improvement. In particular, we consider the steepest descent and Polack-Ribiere conjugate gradient methods which are gradient based. In the Hessian-based category, we consider the quasi-Newton, trust region, and feasible Newton methods. These techniques are used to improve the quality of a mesh by repositioning the vertices, where the overall mesh quality is measured by the sum of the squares of individual elements according to the aspect ratio metric. The effects of the desired degree of accuracy in the improved mesh, problem size, initial mesh configuration, and heterogeneity in element volume on the performance of the optimization solvers are characterized on a series of tetrahedral meshes.

Recommendations for grounding systems in lightning protection systems

This paper presents some practical recommendations for designing grounding systems as part of an integral protection system against lightning strikes. These recommendations are made taking into account the results of academic software whose development was based on hybrid electromagnetic method and the method of moments. This paper presents the results of impedance and transient voltage for triangle, wye, counterpoises and mesh configurations. Some recommendation are made concerning the use and characteristics of earth electrodes, for example effective counterpoise length, where to locate grounding rods, where to connect down conductors in a mesh and the potential difference between points on the same grounding electrodes. These recommendations guide a systems' designer to ensure greater benefit from grounding setups without wasting money.

Comparison between spring network models and continuum constitutive laws: Application to the large deformation of a capsule in shear flow

A capsule is a liquid drop enclosed by a solid, deformable membrane. To analyze the deformation of a capsule accurately, both the fluid mechanics of the internal and external fluids and the solid mechanics of the membrane must be solved precisely. Recently, many researchers have used discrete spring network models to express the membrane mechanics of capsules and biological cells. However, it is unclear whether such modeling is sufficiently accurate to solve for capsule deformation. This study examines the correlations between the mechanical properties of the discrete spring network model and continuum constitutive laws. We first compare uniaxial and isotropic deformations of a two-dimensional (2D) sheet, both analytically and numerically. The 2D sheet is discretized with four kinds of mesh to analyze the effect of the spring network configuration. We derive the relationships between the spring constant and continuum properties, such as the Young modulus, Poisson ratio, area dilation modulus, and shear modulus. It is found that the mechanical properties of spring networks are strongly dependent on the mesh configuration. We then calculate the deformation of a capsule under inflation and in a simple shear flow in the Stokes flow regime, using various membrane models. To achieve high accuracy in the flow calculation, a boundary-element method is used. Comparing the results between the different membrane models, we find that it is hard to express the area incompressibility observed in biological membranes using a simple spring network model.

Degradation, infection and heat effects on polypropylene mesh for pelvic implantation: what was known and when it was known.

Many properties of polypropylene mesh that are causative in producing the complications that our patients are experiencing were published in the literature prior to the marketing of most currently used mesh configurations and mesh kits. These factors were not sufficiently taken into account prior to the sale of these products for use in patients. This report indicates when this information was available to both mesh kit manufacturers and the Food and Drug Administration.

Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation

Haemodynamic factors, in particular wall shear stresses (WSSs) may have significant impact on growth and rupture of cerebral aneurysms. Without a means to measure WSS reliably in vivo, computational fluid dynamic (CFD) simulations are frequently employed to visualise and quantify blood flow from patient-specific computational models. With increasing interest in integrating these CFD simulations into pretreatment planning, a better understanding of the validity of the calculations in respect to computation parameters such as volume element type, mesh size and mesh composition is needed. In this study, CFD results for the two most common aneurysm types (saccular and terminal) are compared for polyhedral- vs. tetrahedral-based meshes and discussed regarding future clinical applications. For this purpose, a set of models were constructed for each aneurysm with spatially varying surface and volume mesh configurations (mesh size range: 5119–258, 481 volume elements). WSS distribution on the model wall and point-based velocity measurements were compared for each configuration model. Our results indicate a benefit of polyhedral meshes in respect to convergence speed and more homogeneous WSS patterns. Computational variations of WSS values and blood velocities are between 0.84 and 6.3% from the most simple mesh (tetrahedral elements only) and the most advanced mesh design investigated (polyhedral mesh with boundary layer).

Experimental and code based study of beam–column behaviour of equal leg single-angles

This study intends to establish the behaviour of equal leg single-angles under bi-axial bending and compression experimentally and analytically using nonlinear finite element analysis techniques. In the experimental system, 100×100×10 nominally sized equal leg single-angles were tested. The legs of angles were placed in a flange-down position, and their axes were restrained laterally as a simply supported beam. These simply supported beams were loaded with a concentrated vertical load applied at the middle of the length and an axial load parallel to the beam axis until failure was reached. The FEM model of the experimental system was generated with ABAQUS software. Suitable finite element mesh configurations, boundary conditions and nonlinear analysis solution techniques were determined to identify the system behaviour. The test and ABAQUS software results were interpreted using some code specifications and one design procedure.

The impact of mesh adaptivity on the gravity current front speed in a two-dimensional lock-exchange

Numerical simulations of the two-dimensional lock-exchange flow are used to evaluate the performance of adaptive meshes as implemented in the non-hydrostatic, finite-element model Fluidity-ICOM. The lock-exchange is a widely studied laboratory-scale set-up that produces two horizontally propagating gravity currents and incorporates key physical processes associated with gravity currents over many scales, including ocean overflows. The Froude number (non-dimensional front speed) is used to assess simulations performed on structured-fixed, unstructured-fixed and unstructured-adaptive meshes and different adaptive mesh configurations are compared. Fluidity-ICOM successfully captures the flow dynamics, including the development of Kelvin–Helmholtz billows. Mesh adapts are guided by a metric which is key to the ability of an adaptive mesh to represent the flow. The metric employed in Fluidity-ICOM is simple, based on the curvature of the solution fields and user-defined solution field weights. Good representation of the gravity current front region is essential to the quality of the solution and for the adaptive meshes this is achieved by reducing the horizontal velocity field weight near the boundaries. Adaptive meshes that are configured in this way are seen to perform as well as high-resolution fixed meshes whilst using at least one order of magnitude fewer nodes. The Froude numbers also compare well with previously published values determined from experimental, numerical and theoretical approaches. The substantial reduction in the number of nodes used by the adaptive meshes is particularly encouraging as it suggests that even greater gains may be achieved in three-dimensional simulations and larger-scale problems. Results show that successful use of the adaptive mesh approach employed requires a clear understanding of the physics of the system and the metric. These considerations will be vital to the effective application of adaptive mesh approaches in numerical modelling of more complex ocean flows.