Numerical simulation of powder mixed electric discharge machining (PMEDM) using finite element method

In the present paper, an axisymmetric two-dimensional model for powder mixed electric discharge machining (PMEDM) has been developed using the finite element method (FEM). The model utilizes the several important aspects such as temperature-sensitive material properties, shape and size of heat source (Gaussian heat distribution), percentage distribution of heat among tool, workpiece and dielectric fluid, pulse on/off time, material ejection efficiency and phase change (enthalpy) etc. to predict the thermal behaviour and material removal mechanism in PMEDM process. The developed model first calculates the temperature distribution in the workpiece material using ANSYS (version 5.4) software and then material removal rate (MRR) is estimated from the temperature profiles. The effect of various process parameters on temperature distributions along the radius and depth of the workpiece has been reported. Finally, the model has been validated by comparing the theoretical MRR with the experimental one obtained from a newly designed experimental setup developed in the laboratory.     

New adaptive mixed finite element method (AMFEM)

The need to develop reliable and efficient adaptive algorithms using mixed finite element methods arises from various applications in fluid dynamics and computational continuum mechanics. In order to save degrees of freedom, not all but just some selected set of finite element domains are refined and hence the fundamental question of convergence requires a new mathematical argument as well as the question of optimality. We will present a new adaptive algorithm for mixed finite element methods to solve the model Poisson problem, for which optimal convergence can be proved. The a posteriori error control of mixed finite element methods dates back to Alonso (1996) Error estimators for a mixed method. and Carstensen (1997) A posteriori error estimate for the mixed finite element method. The error reduction and convergence for adaptive mixed finite element methods has already been proven by Carstensen and Hoppe (2006) Error Reduction and Convergence for an Adaptive Mixed Finite Element Method, Convergence analysis of an adaptive nonconforming finite element methods. Recently, Chen, Holst and Xu (2008) Convergence and Optimality of Adaptive Mixed Finite Element Methods. presented convergence and optimality for adaptive mixed finite element methods following arguments of Rob Stevenson for the conforming finite element method. Their algorithm reduces oscillations, before applying and a standard adaptive algorithm based on usual error estimation. The proposed algorithm does this in a natural way, by switching between the reduction of either the estimated error or oscillations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A finite element method for nearly incompressible elasticity problems

A finite element method is considered for dealing with nearly incompressible material. In the case of large deformations the nonlinear character of the volumetric contribution has to be taken into account. The proposed mixed method avoids volumetric locking also in this case and is robust for (with being the well-known Lamé constant). Error estimates for the -norm are crucial in the control of the nonlinear terms.

     

A symmetric mixed finite element method for nearly incompressible elasticity based on biorthogonal systems

We present a symmetric version of the nonsymmetric mixed finite element method presented in (Lamichhane, ANZIAM J 50 (2008), C324–C338) for nearly incompressible elasticity. The displacement–pressure formulation of linear elasticity is discretized using a Petrov–Galerkin discretization for the pressure equation in (Lamichhane, ANZIAM J 50 (2008), C324–C338) leading to a non‐symmetric saddle point problem. A new three‐field formulation is introduced to obtain a symmetric saddle point problem which allows us to use a biorthogonal system. Working with a biorthogonal system, we can statically condense out all auxiliary variables from the saddle point problem arriving at a symmetric and positive‐definite system based only on the displacement. We also derive a residual based error estimator for the mixed formulation of the problem. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

A mixed virtual element method for nearly incompressible linear elasticity equations

Over the last few years deep artificial neural networks (ANNs) have very successfully been used in numerical simulations for a wide variety of computational problems including computer vision, image classification, speech recognition, natural language processing, as well as computational advertisement. In addition, it has recently been proposed to approximate solutions of high-dimensional partial differential equations (PDEs) by means of stochastic learning problems involving deep ANNs. There are now also a few rigorous mathematical results in the scientific literature which provide error estimates for such deep learning based approximation methods for PDEs. All of these articles provide spatial error estimates for ANN approximations for PDEs but do not provide error estimates for the entire space-time error for the considered ANN approximations. It is the subject of the main result of this article to provide space-time error estimates for deep ANN approximations of Euler approximations of certain perturbed differential equations. Our proof of this result is based (i) on a certain ANN calculus and (ii) on ANN approximation results for products of the form \([0,T]\times \mathbb {R}^{d}\ni (t,x){\kern -.5pt}\mapsto {\kern -.5pt} tx{\kern -.5pt}\in {\kern -.5pt} \mathbb {R}^{d}\) where \(T{\kern -.5pt}\in {\kern -.5pt} (0,\infty )\), \(d{\kern -.5pt}\in {\kern -.5pt} \mathbb {N}\), which we both develop within this article.

     

An optimally convergent adaptive mixed finite element method

We prove convergence and optimal complexity of an adaptive mixed finite element algorithm, based on the lowest-order Raviart–Thomas finite element space. In each step of the algorithm, the local refinement is either performed using simple edge residuals or a data oscillation term, depending on an adaptive marking strategy. The inexact solution of the discrete system is controlled by an adaptive stopping criterion related to the estimator.     

Numerical analysis of incompressible wormhole propagation with mass-preserving characteristic mixed finite element procedure

In this article, we construct a new combined characteristic mixed finite element procedure to simulate the incompressible wormhole propagation. In this procedure, we use the classical mixed finite element method to solve the pressure equation and a modified mass-preserving characteristic finite element method for the solute transport equation, and solve the porosity function straightly by the given concentration. This combined method not only keeps mass balance globally but also preserves maximum principle for the porosity. We considered the corresponding convergence and derive the optimal L²-norm error estimate. Finally, we present some numerical examples to confirm theoretical analysis.

     

An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method

An approximation scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using two methods. Standard mixed finite element is used for the Darcy velocity equation. A characteristics-mixed finite element method is presented for the concentration equation. Characteristic approximation is applied to handle the convection part of the concentration equation, and a lowest-order mixed finite element spatial approximation is adopted to deal with the diffusion part. Thus, the scalar unknown concentration and the diffusive flux can be approximated simultaneously. In order to derive the optimal L²$L^2$-norm error estimates, a post-processing step is included in the approximation to the scalar unknown concentration. This scheme conserves mass globally; in fact, on the discrete level, fluid is transported along the approximate characteristics. Numerical experiments are presented finally to validate the theoretical analysis.     

On using meshes of mixed element types in adaptive finite element analysis

Four different automatic mesh generators capable of generating either triangular meshes or hybrid meshes of mixed element types have been used in the mesh generation process. The performance of these mesh generators were tested by applying them to the adaptive finite element refinement procedure. It is found that by carefully controlling the quality and grading of the quadrilateral elements, an increase in efficiency over pure triangular meshes can be achieved. Furthermore, if linear elements are employed, an optimal hybrid mesh can be obtained most economically by a combined use of the mesh coring technique suggested by Lo and Lau and a selective removal of diagonals from the triangular element mesh. On the other hand, if quadratic elements are used, it is preferable to generate a pure triangular mesh first, and then obtain a hybrid mesh by merging of triangles.     

Numerical simulation of transversely isotropic material using adaptive FEM

Lightweight construction plays an important role in the global task to save energy. A common approach to reduce the weight of structures is the use of composite materials like fibre reinforced polymers (FRP). FRP can be characterised by transversely isotropic material behaviour, a special case of anisotropy. In our arcticle we present the necessary efforts to include such material behaviour into an existing adaptive FEM code. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new family of stable elements for nearly incompressible elasticity based on a mixed Petrov-Galerkin finite element formulation

Summary Adding to the classical Hellinger Reissner formulation another residual form of the equilibrium equation, a new Petrov-Galerkin finite element method is derived. It fits within the framework of a mixed finite element method and is proved to be stable for rather general combinations of stress and displacement interpolations, including equal-order discontinuous stress and continuous displacement interpolations which are unstable within the Galerkin approach. Error estimates are presented using the Babuka-Brezzi theory and numerical results confirm these estimates as well as the good accuracy and stability of the method.Dedicated to Professor Ivo Babuka on the occasion of his sixtieth birthdayPrepared for the conference on: The Impact of Mathematical Analysis on the Solution of Engineering Problems. University of Maryland, September 1986.     

Convergence of an adaptive mixed finite element method for convection-diffusion-reaction equations

We prove the convergence of an adaptive mixed finite element method (AMFEM) for (nonsymmetric) convection-diffusion-reaction equations. The convergence result holds for the cases where convection or reaction is not present in convection- or reaction-dominated problems. A novel technique of analysis is developed by using the superconvergence of the scalar displacement variable instead of the quasi-orthogonality for the stress and displacement variables, and without marking the oscillation dependent on discrete solutions and data. We show that AMFEM is a contraction of the error of the stress and displacement variables plus some quantity. Numerical experiments confirm the theoretical results.     

An adaptive least-squares mixed finite element method for nonlinear parabolic problems

A least-squares mixed finite element method for nonlinear parabolic problems is investigated in terms of computational efficiency. An a posteriori error estimator, which is needed in an adaptive refinement algorithm, was composed with the least-squares functional, and a posteriori errors were effectively estimated.     

Error reduction and convergence for an adaptive mixed finite element method

An adaptive mixed finite element method (AMFEM) is designed to guarantee an error reduction, also known as saturation property: after each refinement step, the error for the fine mesh is strictly smaller than the error for the coarse mesh up to oscillation terms. This error reduction property is established here for the Raviart-Thomas finite element method with a reduction factor uniformly for the norm of the flux errors. Our result allows for linear convergence of a proper adaptive mixed finite element algorithm with respect to the number of refinement levels. The adaptive algorithm surprisingly does not require any particular mesh design, unlike the conforming finite element method. The new arguments are a discrete local efficiency and a quasi-orthogonality estimate. The proof does not rely on duality or on regularity.

     

A robust boundary element method for nearly incompressible linear elasticity

Summary. For a Dirichlet boundary value problem in linear elasticity we consider a boundary element method which is robust for nearly incompressible materials. Based on the spectral properties of the single layer potential for the Stokes problem we introduce an orthogonal splitting of the trial space. The resulting variational problem is then well conditioned and can be discretized by using standard boundary element methods. Mathematics Subject Classification (1991):65N38     

A stress-velocity least-squares mixed finite element formulation for incompressible elastodynamics

In the framework of solving elastodynamic problems using a least-squares mixed finite element method (LSFEM) the implementation of a stress-velocity formulation for small strains is introduced and discussed in the present contribution. The element formulation is based on a first-order div – grad system, with the balance equation of momentum and the constitutive law as the governing equations. Application of the L₂-norm to the two residuals leads to a functional depending on stresses and velocities. Different time discretization schemes are considered, a scalar weighting is introduced and chosen in dependency of the different time discretization methods. In a numerical example the influence of the time integration method, the chosen time step width and the related weighing factor are investigated for a two-dimensional problem. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical analysis of a mixed finite element method for a flow-transport problem

We consider a generic flow-transport system of partial differential equations, which has wide application in the waste disposal industry. The approximation of this system, using a finite element method for the brine, radionuclides, and heat combined with a mixed finite element method for the pressure and velocity, is analyzed. Optimal order error estimates in H¹ and L² are derived. The error analysis is given with no restriction on the diffusion tensor. That is, we have included the effects of molecular diffusion and dispersion. © 1996 John Wiley & Sons, Inc.     

A Galerkin/least-square finite element formulation for nearly incompressible elasticity/stokes flow

A Galerkin/least-square finite element formulation (GLS) is used to study mixed displacement-pressure formulation of nearly incompressible elasticity. In order to fully incorporate the effect of the residual-based stabilized term to the weak form, the second derivatives of shape functions were also derived and accounted, which can accurately discretize the residual term and improve the GLS method as well as the Petrov–Galerkin method. The numerical studies show that improved stabilized method can effectively remove volumetric locking problem for incompressible elasticity and stabilize the pressure field for stokes flow. When apply GLS to study material nonlinearity, the derivative of tangent modulus at the integration point will be required. Both advantage and disadvantage of using GLS method for nearly incompressible elasticity/stokes flow were demonstrated.     

2D simulation of fluid-structure interaction using finite element method

This paper deals with pressure-based finite element analysis of fluid–structure systems considering the coupled fluid and structural dynamics. The present method uses two-dimensional fluid elements and structural line elements for the numerical simulation of the problem. The equations of motion of the fluid, considered inviscid and compressible, are expressed in terms of the pressure variable alone. The solution of the coupled system is accomplished by solving the two systems separately with the interaction effects at the fluid–solid interface enforced by an iterative scheme. Non-divergent pressure and displacement are obtained simultaneously through iterations. The Galerkin weighted residual method-based FE formulation and the iterative solution procedure are explained in detail followed by some numerical examples. Numerical results are compared with the existing solutions to validate the code for sloshing with fluid–structure coupling.