The background of error estimation and adaptivity in finite element computations
|
|
- Beatrice Thompson
- 5 years ago
- Views:
Transcription
1 Comput. Methods Appl. Mech. Engrg. 195 (2006) The background of error estimation and adaptivity in finite element computations Olgierd C. Zienkiewicz 1 School of Engineering, University of Wales, Swansea, UK Received 29 June 2004; accepted 10 July 2004 Abstract The paper summarizes the present state of the art concerning error estimation and adaptive remeshing. Recent work shows that methods based on recovery result on most accurate error estimators. Ó 2005 Elsevier B.V. All rights reserved. Keywords: Error estimation; Adaptive remeshing 1. Introduction This paper is being written to open the special issue of this journal concerned with reporting the major achievements of the conference called ADMOS held in Göteborg, Sweden, from September 29 to October 1, It is appropriate when considering the present state of the art of any subject, to be informed about its background and origins. Error estimation in numerical computations is obviously as old as the numerical computations themselves. The very first paper by Richardson [1] published in 1910 is the first paper for practical computations utilizing finite differences which presents a general procedure by which errors can be estimated. The original Richardson method starts from the fact that errors in the solution of any finite difference scheme usually depend of some power of the mesh size used (or today the size of finite elements). The order of the error is thus usually given and depends as well on the solution method adopted. In the Richardson procedure a set of meshes, which nest within each other, are used to provide the solution, then with a knowledge of the convergence rate the errors can be estimated. 1 Unesco Professor, CIMNE, UPC, Barcelona, Spain /$ - see front matter Ó 2005 Elsevier B.V. All rights reserved. doi: /j.cma
2 208 O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) Unfortunately, the procedures used by Richardson and his early followers are not easily adaptable to the computer age. At the present stage the time and resources available for resolution of the problem on a finer mesh do not exist as generally all the funds available for computation have been used in the first solution. What is therefore required is to find the error in the existing solution which certainly cannot be repeated on a mesh of half the size. 2. Residual and recovery based methods 2.1. Error estimation The work, originally introduced by Babuška and Rheinboldt [2,3], to estimate errors, considers local residuals of the numerical solution. By investigating the residuals occurring in a patch of elements or even in a single element it becomes possible to estimate the errors which arise locally, usually in such norm as the energy norm. The original approaches for error estimation and indeed for correcting these by adaptive refinement started by simply looking at the elements with largest error and dividing these to achieve some acceptable accuracy. Since this early work many players entered the field and today the procedures available for error estimation are essentially reduced to two kinds. These are, first, the residual error estimators continuing the original suggestions of Babuška but the emphasis is now on use of self-equilibrating patches. Here the work of Ainsworth and Oden [4,5] and others is important. Such residual approaches have gained universal mathematical aplomb and for some years formed the basis of acceptable error estimators. However, in the second approach, an alternative was possible and this was suggested by myself and my collaborator Dr. J.Z. Zhu. This is the use of a recovery process to obtain more accurate representation of the unknowns. In the original process [6], introduced in 1987, simple averaging and the so-called L 2 projection recovery of Brauchli and Oden [7] were used to estimate errors. Although these methods were quite successful, better recovery procedures were sought. It was only in 1992that a very accurate recovery process based on superconvergence was achieved [8 10]. This so-called superconvergent patch recovery (SPR) method led to very accurate recovered values which generally converged at a higher rate than the original solution and which gave a solid basis for error estimation. In the first uses of the SPR the standard process was to find, in a patch corresponding to the solution surrounding each node or element, a set of points at which the accuracy was higher such as the superconvergence points and then to pass a least squares approximation of the function of the stresses or gradients (as in the case of thermal problems) through such values. With the knowledge of these more accurate stresses or gradients of the solution, the energy norm could easily be estimated and indeed so could the local values of stresses and gradients. Today the methodology is widely used, industrially, following the same practice but a variant of it is sometimes desirable. In that variant the superconvergent values are taken not as the stresses or gradients of the solution, but on the values of the original functions themselves. It is known for instance that superconvergence will always occur in these values for triangular elements at the vertex nodes, whereas in the same triangular elements such superconvergence for the stresses or gradients does not exist. Therefore for triangular elements it is preferable to include in the patch superconvergence of the displacements or the basic function itself and then to find its gradients by establishing a suitable polynomial of higher order. This procedure of obtaining the recovered values is currently widely used and many papers by Wiberg (who is a co-editor of this issue of the Journal) have contributed very much to this recovery methodology [11 15]. Although both the recovery and residual based procedures have proven their worth it was necessary to establish which of the particular methods proved most accurate. To test these procedures Babuška together
3 with a team of researchers led by Professor Strouboulis at the University of Texas at Austin, produced a general methodology. This methodology is essentially based on what I like to call the Babǔska patch test in which arbitrary fields of finite elements, assumed to be repeatable in space, are subjected to the next higher order basic polynomial than contained in the original shape functions. Obviously if the finite elements used include quadratics any solution which is quadratic will be exactly returned through the numerical computation and again equally the main errors will be concentrated only in the cubic terms. Thus it is necessary, in a repeatable sense, to impose only cubic polynomials on the patch and test whether errors arise and, if so, their magnitude. By doing this it is easy to establish whether the procedures capture errors exactly or simply show the errors for inexact solutions. Here Babuška proposed to value each error estimator with a robustness index and after very many tests concluded that the robustness index was definitively optimal for the recovery methods, with the residual methodology giving much poorer results. This work is reported in a series of papers [16 20] and in the book available in 2001 [21] Adaptivity O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) Once the error estimator process has been achieved it is natural to seek methods by which the design can be improved. This process is not as simple as first implied by the work of the pioneers who simply considered a subdivision of elements in which error indicators were too large. There should be a logical way in which adaptivity can proceed. This is something to which attention was only sparingly given in the early days. It seems to me that we now have a method of obtaining a solution which gives certain accuracy with such norms as energy. The whole process today would be summarized as follows. First the errors in each element would be estimated, if for instance energy norm is sought, and the total energy norm error for the whole problem will be established. Once the decision is taken that the total error permissible should be less than a certain value, it is a simple matter to search the design field for a new solution in which the total error satisfies this requirement. The simplest process is to try to achieve an equal error distribution between all elements as with such equal error distribution the results are most economical [6,22 25]. Considering the total permissible error, the error can be divided between the existing number of elements used and it is then found whether more elements need to be used. Generally this will be the case but at the same time, in some other areas of the problem we may be over generous with subdivision. For this reason it is best in our experience to redesign the mesh completely and use the results of such a calculation to gain a new subdivision. Of course this is simple, if such requirements as those placed by mathematicians where the total energy norm has to be limited to some value, is used. It is much more difficult to provide answers for other quantities of interest such as local stresses, groups of particular forces, etc. How to get this is obviously something that occupies minds of mesh designers and I should not discuss it further, except to mention that in many cases we are not concerned by overall errors or errors in such norms as energy but wish to look at another kind of norm and possibly to a local value of some quantity, forgetting all the others. Much progress is today being achieved in solution of problems of goal oriented adaptivity or estimation in the error on quantities of interest. And this by itself will form a subject of some papers at this conference [26 29] and I anticipate new ideas will be presented. 3. Bounds on quantities of interest Although we have shown that excellent estimators of errors exist today, many researchers are striving to know that these estimators are not only close but that they are bounded. The strain energy was one of the first quantities in which bounds could be established. Here the classical work of Fraeijs de Veubeke in the mid-1960s is of vital importance [30,31].
4 210 O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) It was quickly realized by Fraeijs de Veubeke that the standard (displacement) procedures from which structural analysis usually started would provide a lower bound of the strain energy contained in the structure and thus the solution always underestimated the value of strain energy. He therefore sought procedures which could solve the same structural problem by concentrating on so-called complementary energy this would allow the establishment of solutions in which the strain energy would always be overestimated. This process proved very difficult as equilibrating solutions have to be established at all stages. I had the opportunity of collaborating with Fraeijs de Veubeke and suggested a possible way, useful for many two dimensional problems, by the use of stress functions and slab analogy [32]. Nevertheless the methodology did not succeed as a practical way of providing the bounds of strain energy in an actual analysis. Much later when the residual method was being applied to determine error in structural analysis it was realized that once again opportunity would arise for establishing bounds. The residuals are nothing else but a measure by which the numerical solution fails to satisfy the differential equations of the problem. By using a local solution, generally based on single elements, the error can be estimated locally and the total error be obtained by combining these estimates for all elements. The completely independent solution for the displacement, stresses, etc. established by the residuals provides the measure of the error. This solution can be carried out in a number of ways and here a departure from the original, say classical displacement, method can be made. As many local problems are solved, it is assumed the total error is a combination of local (patch) solutions. It seems that one of the first to extend the concept of establishing upper bounds is Kelly in 1984 [33] and in his subsequent work. He endeavored to obtain solutions for the residual placed as a load with equilibrating methodologies. Two similar alternative approaches, though stemming from completely different origins, were proposed by Ladevèze [34,35] and almost simultaneously by Bank and Weiser [36]. These were later adopted by Ainsworth and Oden forming a large part to their book [37]. The adaptive methodologies so far produced are very effective when the basic quantity of interest is a simple one, such as strain energy or the energy norm. However, when the quantity of interest is more localized and if for instance it is the displacement at some part of the structure, rather than an overall measure of stresses as it is in the case of energy, different procedures arise and pure examination of energy errors does not suffice or it is not very selective in showing how to correct the answers. For this reason much effort has been given in recent years to discussing the possibilities of bounds and optimizing the manner in which such localized adaptive goals of analysis can be solved. Much of the present work in this field concentrates on such methodologies and indeed many papers presented later in this special issue deal with these problems. The first to give attention to the possible extension of norms was a series of papers presented by Babuška and Miller in 1984 [38 40]. These papers laid the foundations for much of the work started some 10 years later and which today occupies much of our interest. Here it appears that the first full extension of the methodology is due to Peraire who published in 1997 with many papers following [41 46]. By extending the ideas introduced in the Babuška and Miller papers, Peraire et al. show that whatever the quantity of interest is, it is always possible to establish an adjoint problem which can be solved on the same mesh of the original problem but now with different loads for dealing with the accuracy. Such adjoint problems may be called differently and names such as extraction problems and dual problems are also used. Although many people are now entering the field and the methodology has been followed by Oden and Prudhomme in a series of papers, it appears that only Peraire so far has extended the approach to non self-adjoint problems such as fluid dynamics [43,45,46]. Here again, the extension of Babuška and Miller goes much further and I believe that a number of papers included in this issue are following this line of attack. Clearly, much remains to be done and the range of possibilities to appear is quite large. The equilibrated methods have always provided upper bounds for such quantities as strain energy, etc. Similar position of the bounding occurs if we look at the energy in the adjoint approaches. However, some interest now goes back to placing satisfactorily the lower bounds and thus bracket the solutions. Of course,
5 O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) lower bounds of zero values have been used and this is rather defeatist because they do not provide a tight specification. It is my belief that more precise bounds are required and here the work of Díez, who himself is a co-editor of this issue, and Huerta seems to lead the way [47 49]. Although the bounding of quantities of interest involve quite sophisticated mathematics, it is possible that simpler approaches will evolve in general engineering practice, thus leading possibly to a parallel with the main error estimators which have been so satisfactorily introduced for industry in general. There is apparent need for such error estimators in fluid dynamics and in many other problems. We look forward to some interesting work being presented for this in the future. 4. Closure Clearly the concept of bounds and local solutions can be extended to other quantities of interest beyond strain energy, such as may arise for instance in determining a local displacement, a maximum stress, a fracture mechanics quantity, etc. Indeed in extending the procedures to non self-adjoint problems in which symmetry does not exist the concept of residual is paramount, however, it does not need to be the only way of establishing the kind of error bound. As clearly the residuals of the original solution need not be calculated to obtain the best answers, already it has been shown that much better estimates of error and accuracy are obtained by the process of recovery. And it seems to me we should endeavour to combine the achievements of a good recovery process with the possibilities of bounding by local solutions. It is clear to me that recovered solution values could be used in any procedure to find a set of residuals corresponding to the recovered solution and not to the original numerical one and, thus, to establish bounds on the recovered solution which has been shown to be very precise. Here, incidentally, an energy norm bound obtaining local solutions now would be more efficacious than starting from the so-called equilibrated-residual basis. The reader will have observed that in many papers included in this symposium the question of bounds and of extending the possibilities beyond the simple strain energy is possible. Some of the work concerning fluid mechanics should be of special interest here as clearly in such problems a bound no longer exists in the energy norm way. References [1] L.F. Richardson, The approximate arithmetical solution by finite differences of physical problems, Trans. Roy. Soc. (London) A 210 (1910) [2] I. Babuška, C. Rheinboldt, A-posteriori error estimates for the finite element method, Int. J. Numer. Methods Engrg. 12(1978) [3] I. Babuška, C. Rheinboldt, Adaptive approaches and reliability estimates in finite element analysis, Comput. Methods Appl. Mech. Engrg (1979) [4] M. Ainsworth, J.T. Oden, A procedure for a posteriori error estimation for h p finite element methods, Comput. Methods Appl. Mech. Engrg. 101 (1992) [5] M. Ainsworth, J.T. Oden, A unified approach to a posteriori error estimation using element residual methods, Numerische Mathematik 65 (1993) [6] O.C. Zienkiewicz, J.Z. Zhu, A simple error estimator and adaptive procedure for practical engineering analysis, Int. J. Numer. Methods Engrg. 24 (1987) [7] H.J. Brauchli, J.T. Oden, On the calculation of consistent stress distributions in finite element applications, Int. J. Numer. Methods Engrg. 3 (1971) [8] O.C. Zienkiewicz, J.Z. Zhu, The superconvergent patch recovery (SPR) and adaptive finite element refinement, Comput. Methods Appl. Mech. Engrg. 101 (1992) [9] O.C. Zienkiewicz, J.Z. Zhu, The superconvergent patch recovery (SPR) and a posteriori error estimates. Part 1: The recovery technique, Int. J. Numer. Methods Engrg. 33 (1992)
6 212 O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) [10] O.C. Zienkiewicz, J.Z. Zhu, The superconvergent patch recovery (SPR) and a posteriori error estimates. Part 2: Error estimates and adaptivity, Int. J. Numer. Methods Engrg. 33 (1992) [11] N.-E. Wiberg, L.F. Zeng, X.D. Li, Error estimation and adaptivity in elastodynamics, Comput. Methods Appl. Mech. Engrg. 101 (1992) [12] N.-E. Wiberg, F. Abdulwahab, Patch recovery based on superconvergent derivatives and equilibrium, Int. J. Numer. Methods Engrg. 36 (1993) [13] N.-E. Wiberg, F. Abdulwahab, S. Ziukas, Enhanced superconvergent patch recovery incorporating equilibrium and boundary conditions, Int. J. Numer. Methods Engrg. 37 (1994) [14] N.-E. Wiberg, X.D. Li, Superconvergent patch recovery of finite element solutions and a posteriori l 2 norm error estimate, Commun. Numer. Methods Engrg. 10 (1994) [15] N.-E. Wiberg, F. Abdulwahab, X.D. Li, Error estimation and adaptive procedures based on superconvergent patch recovery (SPR), Arch. Computat. Methods Engrg. 4 (3) (1997) [16] I. Babuška, T. Strouboulis, C.S. Upadhyay, A model study of the quality of a posteriori error estimators for linear elliptic problems. Error estimation in the interior of patchwise uniform grids of triangles, Comput. Methods Appl. Mech. Engrg. 114 (1994) [17] I. Babuška, T. Strouboulis, C.S. Upadhyay, S.K. Gangaraj, K. Copps, Validation of a posteriori error estimators by numerical approach, Int. J. Numer. Methods Engrg. 37 (1994) [18] I. Babuška, T. Strouboulis, C.S. Upadhyay, S.K. Gangaraj, K. Copps, An objective criterion for assessing the reliability of a posteriori error estimators in finite element computations, USACM Bull. 7 (1994) [19] I. Babuška, T. Strouboulis, C.S. Upadhyay, K.S. Gangaraj, A posteriori estimation and adaptive control of the pollution error in the h-version of the finite element method, Int. J. Numer. Methods Engrg. 38 (1995) [20] I. Babuška, T. Strouboulis, K.S. Gangaraj, C.S. Upadhyay, Pollution error in the h-version of the finite element method and the local quality of the recovered derivatives, Comput. Methods Appl. Mech. Engrg. 140 (1995) [21] I. Babuška, T. Strouboulis, The Finite Element Method and Its Reliability, Clarendon Press, Oxford, [22] L.Y. Li, P. Bettess, J.W. Bull, T. Bond, I. Applegarth, Theoretical formulations for adaptive finite element formulations, Commun. Numer. Methods Engrg. 11 (1995) [23] L.Y. Li, P. Bettess, Notes on mesh optimal criteria in adaptive finite element computations, Commun. Numer. Methods Engrg. 11 (1995) [24] E. Oñate, G. Bugeda, A study of mesh optimality criteria in adaptive finite element analysis, Engrg. Computat. 10 (1993) [25] P. Díez, A. Huerta, A unified approach to remeshing strategies for finite element h-adaptivity, Comput. Methods Appl. Mech. Engrg. 176 (1999) [26] A. Huerta, N. Parés, P. Díez, Flux free residual a posteriori error estimators yielding simultaneous upper and lower bounds of the error, in: N.E. Wiberg, P. Díez, (Eds.), Adaptive Modeling and Simulation, vol. 1, CIMNE, October 2003, pp [27] E. Florentin, L. Gallimard, P. Ladevèze, J.P. Pelle, P. Rougeot, Local constitutive relation error for mises stresses in 3d elasticity, in: N.E. Wiberg, P. Díez, (Eds.), Adaptive Modeling and Simulation, vol. 1, CIMNE, October 2003, pp [28] J. Peraire, Computing bounds for functional outputs of exact solutions of partial differential equations, in: N.E. Wiberg, P. Díez, (Eds.), Adaptive Modeling and Simulation, vol. 1, CIMNE, October 2003, pp [29] N. Parés, J. Bonet, A. Huerta, J. Peraire, Guaranteed bounds for linear outputs of interest in elasticity, in: N.E. Wiberg, P. Díez, (Eds.), Adaptive Modeling and Simulation, vol. 1, CIMNE, October 2003, pp [30] B. Fraeijs de Veubeke, Displacement and equilibrium models in finite element method, in: O.C. Zienkiewicz, G.S. Holister (Eds.), Stress Analysis, John Wiley and Sons, Chichester, 1965, pp , Chapter 9. [31] B. Fraeijs de Veubeke, Bending and stretching of plates. Special models for upper and lower bounds, in: Proc. First Conf. Matrix Methods in Structural Mechanics, vol. AFFDL-TR-66-80, Wright Patterson Air force Base, Ohio, October 1966, pp [32] B. Fraeijs de Veubeke, O.C. Zienkiewicz, Strain energy bounds in finite element analysis by slab analogy, J. Strain Anal. 2 (1967) [33] D.W. Kelly, The self-equilibration of residuals and complementary a posteriori error estimates in the finite element method, Int. J. Numer. Methods Engrg. 20 (1984) [34] P. Ladevèze, D. Leguillon, Error estimate procedure in the finite element method and applications, SIAM J. Numer. Anal. 20 (3) (1983) [35] P. Ladevèze, G. Coffignal, J.P. Pelle, Accuracy of elastoplastic and dynamic analysis, in: I. Babuška, O.C. Zienkiewicz, J. Gago, E.R. de A. Oliviera, (Eds.), Accuracy Estimates and Adaptive Refinements in Finite Element Computations, 1986, Chapter 11. [36] R.E. Bank, A. Weiser, Some a posteriori error estimators for elliptic partial differential equations, Math. Computat. 44 (1985) [37] M. Ainsworth, J.T. Oden, A Posteriori Error Estimation in Finite Element Analysis, John Wiley and Sons, New York, [38] I. Babuška, A. Miller, The post-processing approach in the finite element method. Part 1: Calculation of displacements, stresses and other higher derivatives of displacements, Int. J. Numer. Methods Engrg. 20 (1984)
7 O.C. Zienkiewicz / Comput. Methods Appl. Mech. Engrg. 195 (2006) [39] I. Babuška, A. Miller, The post-processing approach in the finite element method. Part 2: The calculation of stress intensity factors, Int. J. Numer. Methods Engrg. 20 (1984) [40] I. Babuška, A. Miller, The post-processing approach in the finite element method. Part 3: A posteriori error estimates and adaptive mesh selection, Int. J. Numer. Methods Engrg. 20 (1984) [41] M. Paraschivoiu, J. Peraire, A.T. Patera, A posteriori finite element bounds for linear-functional outputs of elliptic partial differential equations, Comput. Methods Appl. Mech. Engrg. 150(1 4): Symposium on Advances in Computational Mechanics, vol. 2(Austin, TX, 1997). [42] J. Peraire, A.T. Patera, Asymptotic a posteriori finite element bounds for the outputs of noncoercive problems: the Helmholtz and Burgers equations, Comput. Methods Appl. Mech. Engrg. 171 (1 2) (1999) [43] L. Machiels, J. Peraire, A.T. Patera, A posteriori finite-element output bounds for the incompressible Navier Stokes equations: application to a natural convection problem, J. Computat. Phys. 172 (2) (2001) [44] Y. Maday, A.T. Patera, J. Peraire, A general formulation for a posteriori bounds for output functionals of partial differential equations; application to the eigenvalue problem, C. R. Acad. Sci. Paris Sér. I Math. 328 (9) (1999) [45] A.T. Patera, J. Peraire, A general Lagrangian formulation for the computation of a posteriori finite element bounds, in: Error Estimation and Adaptive Discretization Methods in Computational Fluid Dynamics, Lecture Notes Computer Science Engineering, vol. 25, Springer, Berlin, 2003, pp [46] L. Machiels, J. Peraire, A.T. Patera, Output bound approximations for partial differential equations; application to the incompressible Navier Stokes equations, in: Industrial and Environmental Applications of Direct and Large-eddy Simulation (Istanbul 1998), Lecture Notes in Physics, vol. 529, Springer, Berlin, 1999, pp [47] P. Díez, J.J. Egozcue, A. Huerta, A posteriori error estimation for standard finite element analysis, Comput. Methods Appl. Mech. Engrg. 163 (1998) [48] A. Huerta, P. Díez, Error estimation including pollution assessment for nonlinear finite element analysis, Comput. Methods Appl. Mech. Engrg. 181 (2000) [49] P. Díez, N. Parés, A. Huerta, Recovering lower bounds of the error postprocessing implicit residual a posteriori error estimates, Int. J. Numer. Methods Engrg. 56 (2003)
Department of Computing and Software
Department of Computing and Software Faculty of Engineering McMaster University Assessment of two a posteriori error estimators for steady-state flow problems by A. H. ElSheikh, S.E. Chidiac and W. S.
More informationClassic Reprints Series
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 52:287 342 (DOI: 10.1002/nme.339) Classic Reprints Series Displacement and equilibrium models in the nite element
More informationVisualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium
Visualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium Abstract The aim of this paper is to show how to use
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationSHAPE SENSITIVITY ANALYSIS INCLUDING QUALITY CONTROL WITH CARTESIAN FINITE ELEMENT MESHES
VI International Conference on Adaptive Modeling and Simulation ADMOS 2013 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) SHAPE SENSITIVITY ANALYSIS INCLUDING QUALITY CONTROL WITH CARTESIAN
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More information3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis
Adv. Theor. Appl. Mech., Vol. 1, 2008, no. 3, 131-139 3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis Yunhua Luo Department of Mechanical & Manufacturing Engineering, University of
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationOn a nested refinement of anisotropic tetrahedral grids under Hessian metrics
On a nested refinement of anisotropic tetrahedral grids under Hessian metrics Shangyou Zhang Abstract Anisotropic grids, having drastically different grid sizes in different directions, are efficient and
More informationROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL ELEMENTS
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng, 15, 493±499 (1999) ROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationAS 5850 Finite Element Analysis
AS 5850 Finite Element Analysis INTRODUCTION Prof. IIT Madras Goal of engineering computations Perform analysis and design of physical systems and processes subjected to imposed conditions (or loads) and
More informationOn the Matlab S-FEM code for 2D problems using T3 and Q4 elements
On the Matlab S-FEM code for 2D problems using T3 and Q4 elements These codes were developed by Liu, Nguyen and workers. The detailed theoretical background, formulation and implementation procedure are
More informationEffective adaptation of hexahedral mesh using local refinement and error estimation
Key Engineering Materials Vols. 243-244 (2003) pp. 27-32 online at http://www.scientific.net (2003) Trans Tech Publications, Switzerland Online Citation available & since 2003/07/15 Copyright (to be inserted
More informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationLETTERS TO THE EDITOR
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, VOL. 7, 135-141 (1983) LETTERS TO THE EDITOR NUMERICAL PREDICTION OF COLLAPSE LOADS USING FINITE ELEMENT METHODS by S. W. Sloan
More informationAdaptive Refinement of Quadrilateral Finite Element Meshes Based on MSC.Nastran Error Measures
Adaptive Refinement of Quadrilateral Finite Element Meshes Based on MSC.Nastran Error Measures Mark E. Botkin GM R&D Center Warren, MI 48090-9055 Rolf Wentorf and B. Kaan Karamete Rensselaer Polytechnic
More informationcourse outline basic principles of numerical analysis, intro FEM
idealization, equilibrium, solutions, interpretation of results types of numerical engineering problems continuous vs discrete systems direct stiffness approach differential & variational formulation introduction
More informationApplication of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution
Application of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution P. BURDA a,, J. NOVOTNÝ b,, J. ŠÍSTE a, a Department of Mathematics Czech University of Technology
More informationIsogeometric Analysis of Fluid-Structure Interaction
Isogeometric Analysis of Fluid-Structure Interaction Y. Bazilevs, V.M. Calo, T.J.R. Hughes Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA e-mail: {bazily,victor,hughes}@ices.utexas.edu
More informationINTERIOR POINT METHOD BASED CONTACT ALGORITHM FOR STRUCTURAL ANALYSIS OF ELECTRONIC DEVICE MODELS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationA STUDY ON ENERGY NORM AS AN ALTERNATE APPROACH TO MESH CONVERGENCE ON CIRCUMFERENTIAL DISK DOVETAIL OF AIRCRAFT ENGINE COMPRESSOR
International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue 13, December 2018, pp. 1293 1301, Article ID: IJMET_09_13_131 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=13
More informationControl Volume Finite Difference On Adaptive Meshes
Control Volume Finite Difference On Adaptive Meshes Sanjay Kumar Khattri, Gunnar E. Fladmark, Helge K. Dahle Department of Mathematics, University Bergen, Norway. sanjay@mi.uib.no Summary. In this work
More informationA new 8-node quadrilateral spline finite element
Journal of Computational and Applied Mathematics 195 (2006) 54 65 www.elsevier.com/locate/cam A new 8-node quadrilateral spline finite element Chong-Jun Li, Ren-Hong Wang Institute of Mathematical Sciences,
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM. Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s
ADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s. 14-24 ABSTRACT In recent years considerable interest has been shown in the development
More informationFOURTH ORDER COMPACT FORMULATION OF STEADY NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS
International Journal of Mechanical Engineering and Technology (IJMET Volume 9 Issue 10 October 2018 pp. 179 189 Article ID: IJMET_09_10_11 Available online at http://www.iaeme.com/ijmet/issues.asp?jtypeijmet&vtype9&itype10
More informationSmooth finite elements
Smooth finite elements seamless handling of incompressibility, distorted and polygonal meshes; links with equilibrium methods Stéphane Bordas * Nguyen-Xuan Hung ** Nguyen-Dang Hung *** * University of
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationREINER HORST. License or copyright restrictions may apply to redistribution; see
MATHEMATICS OF COMPUTATION Volume 66, Number 218, April 1997, Pages 691 698 S 0025-5718(97)00809-0 ON GENERALIZED BISECTION OF n SIMPLICES REINER HORST Abstract. A generalized procedure of bisection of
More informationAdaptive finite element refinement
15 Adaptive finite element refinement 1 5.1 Introduction In the previous chapter we have discussed at some length various methods of recovery by which the finite element solution results could be made
More informationTICAM - Texas Institute for Computational and Applied Mathematics. The University of Texas at Austin. Taylor Hall Abstract
A New Cloud-Based hp Finite Element Method J. T. Oden, C. A. M. Duarte y and O. C. Zienkiewicz z TICAM - Texas Institute for Computational and Applied Mathematics The University of Texas at Austin Taylor
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationSolution for Euler Equations Lagrangian and Eulerian Descriptions
Solution for Euler Equations Lagrangian and Eulerian Descriptions Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com Abstract We find an exact solution for the system of Euler equations, supposing
More informationModeling Skills Thermal Analysis J.E. Akin, Rice University
Introduction Modeling Skills Thermal Analysis J.E. Akin, Rice University Most finite element analysis tasks involve utilizing commercial software, for which you do not have the source code. Thus, you need
More informationmethod Dietmar Gallistl
Snapshots of modern mathematics from Oberwolfach 13/2016 T h e a d a p t i ve f i n i t e e l e m e n t method Dietmar Gallistl Computer simulations of many physical phenomena rely on approximations by
More information3D Finite Element Software for Cracks. Version 3.2. Benchmarks and Validation
3D Finite Element Software for Cracks Version 3.2 Benchmarks and Validation October 217 1965 57 th Court North, Suite 1 Boulder, CO 831 Main: (33) 415-1475 www.questintegrity.com http://www.questintegrity.com/software-products/feacrack
More informationSDC. Engineering Analysis with COSMOSWorks. Paul M. Kurowski Ph.D., P.Eng. SolidWorks 2003 / COSMOSWorks 2003
Engineering Analysis with COSMOSWorks SolidWorks 2003 / COSMOSWorks 2003 Paul M. Kurowski Ph.D., P.Eng. SDC PUBLICATIONS Design Generator, Inc. Schroff Development Corporation www.schroff.com www.schroff-europe.com
More informationCOMPONENT-BASED REDUCED BASIS SIMULATIONS
COMPONENT-BASED REDUCED BASIS SIMULATIONS 2016 AKSELOS S.A. EPFL Innovation Park 1015 Lausanne, Switzerland www.akselos.com info@akselos.com community.akselos.com Abstract Akselos enables extremely fast
More informationCURVILINEAR MESH GENERATION IN 3D
CURVILINEAR MESH GENERATION IN 3D Saikat Dey, Robert M. O'Bara 2 and Mark S. Shephard 2 SFA Inc. / Naval Research Laboratory, Largo, MD., U.S.A., dey@cosmic.nrl.navy.mil 2 Scientific Computation Research
More informationStudy of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket
RESEARCH ARTICLE OPEN ACCESS Study of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket Gowtham K L*, Shivashankar R. Srivatsa** *(Department of Mechanical Engineering, B. M.
More informationADVANCES IN ADAPTIVE METHODS IN COMPUTATIONAL FLUID MECHANICS. J. Tinsley Oden. Abstract
ADVANCES N ADAPTVE METHODS N COMPUTATONAL FLUD MECHANCS J. Tinsley Oden Texas nstitute for Computational Mechanics The University of Texas at Austin Austin, Texas 78712 Abstract Recent developments in
More informationAdaptive Nonlinear Finite Element Analyses of High Strength Steel Cladding Systems
Missouri University of Science and Technology Scholars' Mine International Specialty Conference on Cold- Formed Steel Structures (2002) - 16th International Specialty Conference on Cold-Formed Steel Structures
More informationRELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS
International Conference on Economic Engineering and Manufacturing Systems Braşov, 26 27 November 2009 RELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS Galina TODOROVA, Valentin
More informationOPTIMIZATION, OPTIMAL DESIGN AND DE NOVO PROGRAMMING: DISCUSSION NOTES
OPTIMIZATION, OPTIMAL DESIGN AND DE NOVO PROGRAMMING: DISCUSSION NOTES MILAN ZELENY Introduction Fordham University, New York, USA mzeleny@fordham.edu Many older texts, with titles like Globally Optimal
More information2 neighbor interaction is necessary; thus, simplifying implementation on a parallel computer. Gradient jumps may be shared between elements with a com
A POSTERIORI FINITE ELEMENT ERROR ESTIMATION FOR DIFFUSION PROBLEM SLIMANE ADJERID y, BELKACEM BELGUENDOUZ z, AND JOSEPH E. FLAHERTY y Abstract. Adjerid et al. [2] and Yu [9, 2] show that a posteriori
More informationSpeed-up of Parallel Processing of Divisible Loads on k-dimensional Meshes and Tori
The Computer Journal, 46(6, c British Computer Society 2003; all rights reserved Speed-up of Parallel Processing of Divisible Loads on k-dimensional Meshes Tori KEQIN LI Department of Computer Science,
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More informationDomain Decomposition and hp-adaptive Finite Elements
Domain Decomposition and hp-adaptive Finite Elements Randolph E. Bank 1 and Hieu Nguyen 1 Department of Mathematics, University of California, San Diego, La Jolla, CA 9093-011, USA, rbank@ucsd.edu. Department
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More informationIntroduction to the Finite Element Method (3)
Introduction to the Finite Element Method (3) Petr Kabele Czech Technical University in Prague Faculty of Civil Engineering Czech Republic petr.kabele@fsv.cvut.cz people.fsv.cvut.cz/~pkabele 1 Outline
More informationTopology Optimization of Two Linear Elastic Bodies in Unilateral Contact
2 nd International Conference on Engineering Optimization September 6-9, 2010, Lisbon, Portugal Topology Optimization of Two Linear Elastic Bodies in Unilateral Contact Niclas Strömberg Department of Mechanical
More informationLocking-Free Smoothed Finite Element Method with Tetrahedral/Triangular Mesh Rezoning in Severely Large Deformation Problems
Locking-Free Smoothed Finite Element Method with Tetrahedral/Triangular Mesh Rezoning in Severely Large Deformation Problems Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology (Japan) P. 1 P. 1 Motivation
More informationFirst Order Analysis for Automotive Body Structure Design Using Excel
Special Issue First Order Analysis 1 Research Report First Order Analysis for Automotive Body Structure Design Using Excel Hidekazu Nishigaki CAE numerically estimates the performance of automobiles and
More informationChallenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure
Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment
More informationInvestigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models. David M. Trujillo 1. December 2005
Investigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models by David M. Trujillo 1 December 2005 1 Consultant, TRUCOMP, Fountain Valley, California trucomp@earthlink.net Abstract
More informationAdaptive numerical methods
METRO MEtallurgical TRaining On-line Adaptive numerical methods Arkadiusz Nagórka CzUT Education and Culture Introduction Common steps of finite element computations consists of preprocessing - definition
More informationPrecise FEM solution of corner singularity using adjusted mesh applied to 2D flow
Precise FEM solution of corner singularity using adjusted mesh applied to 2D flow Jakub Šístek, Pavel Burda, Jaroslav Novotný Department of echnical Mathematics, Czech echnical University in Prague, Faculty
More informationCommutative filters for LES on unstructured meshes
Center for Turbulence Research Annual Research Briefs 1999 389 Commutative filters for LES on unstructured meshes By Alison L. Marsden AND Oleg V. Vasilyev 1 Motivation and objectives Application of large
More informationConjectures concerning the geometry of 2-point Centroidal Voronoi Tessellations
Conjectures concerning the geometry of 2-point Centroidal Voronoi Tessellations Emma Twersky May 2017 Abstract This paper is an exploration into centroidal Voronoi tessellations, or CVTs. A centroidal
More informationIt has been widely accepted that the finite element
Advancements of Extrusion Simulation in DEFORM-3D By G. Li, J. Yang, J.Y. Oh, M. Foster, and W. Wu, Scientific Forming Tech. Corp.; and P. Tsai and W. Chang, MIRDC Introduction It has been widely accepted
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationTruss structural configuration optimization using the linear extended interior penalty function method
ANZIAM J. 46 (E) pp.c1311 C1326, 2006 C1311 Truss structural configuration optimization using the linear extended interior penalty function method Wahyu Kuntjoro Jamaluddin Mahmud (Received 25 October
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationCURVES OF CONSTANT WIDTH AND THEIR SHADOWS. Have you ever wondered why a manhole cover is in the shape of a circle? This
CURVES OF CONSTANT WIDTH AND THEIR SHADOWS LUCIE PACIOTTI Abstract. In this paper we will investigate curves of constant width and the shadows that they cast. We will compute shadow functions for the circle,
More informationPressure Correction Scheme for Incompressible Fluid Flow
AALTO UNIVERSITY School of Chemical Technology CHEM-E7160 Fluid Flow in Process Units Pressure Correction Scheme for Incompressible Fluid Flow Ong Chin Kai 620503 Lee De Ming Benedict 620448 Page 1 Abstract
More informationLOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements
LOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements Sanjay Patel, Claus Hoff, Mark Gwillim The MacNeal-Schwendler Corporation Abstract In large finite element models
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationA New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes
A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationLARGE DEFORMATION ANALYSIS IN GEOMECHANICS USING ADAPTIVE FINITE ELEMENT METHODS
LARGE DEFORMATION ANALYSIS IN GEOMECHANICS USING ADAPTIVE FINITE ELEMENT METHODS By Mina Kardani B.Sc., Software Engineering A Thesis submitted for the Degree of Doctor of Philosophy at The University
More informationAbstract. 1 Where We've Been
Structural optimization: where we've been and where we're going G. N. Vanderplaats Vanderplaats Research & Development, Inc., 1767 S. 8thStreet, Swfz 2% Co/oWo Sprmgj, CO 80906 C/&4 EMail: vanderplaats@vma.com
More informationAPPROXIMATION AND COMPARISON OF ORDINARY DIFFERENTIAL EQUATION USING BY NEW ITERATION METHODS
Review of the Air Force Academy No.1 (33)/2017 APPROXIMATION AND COMPARISON OF ORDINARY DIFFERENTIAL EQUATION USING BY NEW ITERATION METHODS Necdet BİLDİK *, Yasemin BAKIR ** * Celal Bayar University,
More informationCase C3.1: Turbulent Flow over a Multi-Element MDA Airfoil
Case C3.1: Turbulent Flow over a Multi-Element MDA Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
More informationTHE FINITE ELEMENT METHOD AND ITS APPLICATIONS IN ENGINEERING USING ANSYS Course contents:
THE FINITE ELEMENT METHOD AND ITS APPLICATIONS IN ENGINEERING USING ANSYS Course contents: The finite element method (FEM) has become a staple for predicting and simulating the physical behavior of complex
More informationAN IMPROVED METHOD TO MODEL SEMI-ELLIPTICAL SURFACE CRACKS USING ELEMENT MISMATCH IN ABAQUS
AN IMPROVED METHOD TO MODEL SEMI-ELLIPTICAL SURFACE CRACKS USING ELEMENT MISMATCH IN ABAQUS R. H. A. Latiff and F. Yusof School of Mechanical Engineering, UniversitiSains, Malaysia E-Mail: mefeizal@usm.my
More informationFatigue Crack Growth Simulation using S-version FEM
Copyright c 2008 ICCES ICCES, vol.8, no.2, pp.67-72 Fatigue Crack Growth Simulation using S-version FEM M. Kikuchi 1,Y.Wada 2, A. Utsunomiya 3 and Y. Li 4 Summary Fatigue crack growth under mixed mode
More informationCollocation and optimization initialization
Boundary Elements and Other Mesh Reduction Methods XXXVII 55 Collocation and optimization initialization E. J. Kansa 1 & L. Ling 2 1 Convergent Solutions, USA 2 Hong Kong Baptist University, Hong Kong
More informationThe Cell Method: Quadratic Interpolation with Tetrahedra for 3D Scalar Fields
Copyright 2013 Tech Science Press CMES, vol.94, no.4, pp.279-300, 2013 The Cell Method: Quadratic Interpolation with Tetrahedra for 3D Scalar Fields Martino Pani 1, Fulvia Taddei 1 Abstract: The Cell Method
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More informationMATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC
Advanced Mathematics and Mechanics Applications Using MATLAB Third Edition Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hulman Institute of Technology David Halpern University of Alabama
More informationCOMPONENT-BASED REDUCED BASIS SIMULATIONS
COMPONENT-BASED REDUCED BASIS SIMULATIONS www.akselos.com info@akselos.com White Paper Abstract Akselos enables extremely fast simulations of large and complex parametrized systems. The two main ingredients
More informationA MEASURE OF THE AREA OF A HOMOGENEOUS RANDOM SURFACE IN SPACE*
404 NOTES [Vol. XII, No. 4 where fi is a constant, where a(t), b(t), c{t) are arbitrary functions of t, and where Y = (y b) cos Qt + (z c) sin tit, Z = (y b) sin Qi + (z c) cos 01. Comparison of (3) with
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationParallel hp-finite Element Simulations of 3D Resistivity Logging Instruments
Parallel hp-finite Element Simulations of 3D Resistivity Logging Instruments M. Paszyński 1,3, D. Pardo 1,2, L. Demkowicz 1, C. Torres-Verdin 2 1 Institute for Computational Engineering and Sciences 2
More information