Stabilized and Coupled FEM/EFG Approximations for Fluid Problems
|
|
- Wesley Skinner
- 5 years ago
- Views:
Transcription
1 COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM4, Sept. 5-, 24, Beijing, China c 24 Tsinghua University Press & Springer-Verlag Stabilized and Coupled FEM/EFG Approximations for Fluid Problems T.P. Fries, H.G. Matthies Institute for Scientific Computing, Technische Universität Braunschweig Hans-Sommer-Str. 65, D 386 Braunschweig, Germany t.fries@tu-bs.de Abstract Coupled FEM/EFG approximations are employed for the solution of fluid problems, taking advantage of the individual features of meshfree and meshbased methods. Standard coupling approaches are modified in order to obtain shape functions that may be stabilized reliably. The coupled fluid solver can be applied straightforward to problems including moving and rotating obstacles. Key words: meshfree, meshless, coupling, fluid, Navier-Stokes INTRODUCTION Stabilized finite element methods have developed to be standard tools for the simulation of flow problems [5]. However, in the presence of rotating and moving obstacles, the maintenance of a conforming mesh may be almost impossible. Therefore, a comparably new and innovative class of methods will be used which enables the approximation of partial differential equations based on a set of nodes without the need for an additional mesh. However, these meshfree methods [, 7] are comparably time-consuming which limits their usefulness in practice. Throughout this paper the Element Free Galerkin (EFG) method [2] is considered as a popular member of meshfree methods. Coupling meshbased FEM approximations and meshfree EFG approximations enables one to use meshfree methods only in small parts of the domain, where a mesh is difficult to maintain, and standard meshbased methods in the rest of the domain, thereby combining the advantages of both methodologies. The coupling approaches of Huerta [] and Belytschko [3] are discussed. They are modified in order to obtain shape functions that are more suited for stabilization [8, 9]. The resulting formulation is successfully applied to a number of test cases. The plan of the paper is as follows: The next section gives an outline of the Element Free Galerkin (EFG) method [2], which employs shape functions built by the Moving Least Squares (MLS) concept [5]. Then, the approaches of Huerta [] and Belytschko [3] for coupling FEM and EFG approximations are discussed and modified to allow a stabilization of the coupled scheme. The following section compares the coupling approaches for the solution of the one-dimensional advection-diffusion equation, and a convergence test is carried out. Results from the coupled fluid solver are presented, and the advantage of this approach for the solution of complex flows with moving and rotating obstacles is displayed. All test cases show that the coupled approximations have the same order of convergence as pure FEM calculations, and that reliable and accurate solutions are obtained with the modified coupling approaches. OUTLINE OF THE ELEMENT FREE GALERKIN (EFG) METHOD Throughout this paper, we focus in particular on the Element Free Galerkin (EFG) method [2], however, most conclusions can be applied to other meshfree methods as well. The EFG in its original
2 version employs Moving Least Squares (MLS) approximations in a Bubnov-Galerkin formulation of a problem, i.e. test and shape functions of a weak form are identical. Also, a certain treatment of essential boundary conditions has been published as being part of the method [2], but this aspect is not considered here. Moving Least Squares Approximations The construction of meshfree approximations with the MLS methodology, introduced in [5], shall only be briefly outlined at this point, see [, 7] for a more detailed discussion. In the MLS, the domain is discretized by a set of particles I with corresponding dilatation parameters ρ i, which determine the supports Ω i of the resulting meshfree shape functions. One defines an approximation at an arbitrary fixed point x Ω of ũ (x, x) = p T (x) a (x), where p is a complete basis vector which in general consists of monomials. The length k of this vector depends on the dimension of the problem and the desired order of consistency of the resulting approximations. The vector a is the vector of unknown coefficients of the approximation at x. These unknowns are determined by a minimization of the following weighted error functional: J (a) = i I w (x x i ) [u (x i ) ũ (x i )] 2 = i I w (x x i ) [ u i p T (x i ) a (x) ] 2, which leads to a k k system of equations for a. The weighting function w (x x i ) is in general a bell-shaped function e.g. an approximation of the Gaussian function and is non-zero only in the corresponding support Ω i. It ensures the locality of the approximation. The resulting system of equation is w (x x i ) p (x i ) p T (x i ) a (x) i I = w (x x i ) p (x i ) u i, i I M (x) a (x) = B (x) u, hence a (x) = [M (x)] B (x) u. In words, the minimization of the error functional enables one to find a connection between the unknown coefficients of the approximation a with the nodal unknowns u. At this stage, the so-called moving procedure is applied [5], which in fact results into a replacement of x by x. Inserting a into the approximation gives ũ (x) = p T (x) [M (x)] B (x) u = p T (x) [M (x)] i I w (x x i ) p (x i ) u i, where the meshfree shape functions N i (x) = p T (x) [M (x)] w (x x i ) p (x i ) may directly be extracted. In the EFG method, these shape functions are as well taken as test functions in the weak form of a problem. EFG Approximations for Stabilized Weak Forms A straightforward usage of numerical methods based on the Bubnov-Galerkin principle may result in severe numerical problems and stabilization is required. Stabilized methods have developed to be standard tools in the numerical world [5]. They all share the common property of perturbing the test function of a weak form and multiply this modification with the residual of the differential equation under consideration, thereby maintaining the consistency of the formulation [8]. The perturbation of the test function leads to Petrov-Galerkin methods, i.e. test and trial functions are no longer identical. Stabilization for Advection-Dominated Problems It is a well-known fact that Bubnov-Galerkin formulations result for advection-dominated problems as they frequently arise in the context of fluid dynamics, in heavy oscillations, polluting the solution inside the whole domain, see e.g. [4]. These oscillations are apparent in classical meshbased methods, such as the Finite Element method, and occur analogously in meshfree methods. The Streamline-Upwind/Petrov-Galerkin (SUPG) [4] and Galerkin/Least-Squares (GLS) [4] stabilization are the most popular stabilization schemes to smooth oscillations in advection-dominated problems.
3 Ω el = Ω FEM Ω Ω FEM Ω Ω EFG FEM Γ EFG Γ Fig. Decomposition of the domain into Ω FEM, Ω EFG and Ω. Stabilization of Equal-order Interpolations Most desirable from a computational viewpoint are equal-order interpolations for the velocities and pressure. Then, every node has the same degrees of freedom, which are approximated by the same shape functions. The (most restrictive) governing stability criterion of equal-order interpolations for variational problems with constraints such as the incompressible Navier-Stokes equations is the Babuška-Brezzi condition [6]. Stabilization ideas often circumvent this condition (rather than satisfying it). It can be shown that the GLS stabilization may be used for the stabilization of equal-order interpolations as well [2]. However, SUPG stabilized formulations still require stabilization in order to enable equal-order interpolations. In practice, this may be done with the Pressure-Stabilizing/Petrov-Galerkin (PSPG) stabilization [3, 8]. Consequences Numerical studies for the solution of the incompressible Navier-Stokes equations with meshfree equal-order approximations show that SUPG/PSPG stabilization leads to less diffusive, i.e. slightly superior results than GLS [9]. It has also been shown that it is important for the success of the stabilization to employ meshfree shape functions with small dilatation parameters, i.e. small support sizes [8, 9]. A reliable stabilization of non-linear problems cannot be guaranteed with large supports of the meshfree shape functions. This fact is an important motivation for the modification of the existing coupling approaches in the following section. COUPLING EFG AND FEM The EFG method has been successfully applied to a large number of problems, where the maintenance of a conforming mesh may be almost impossible. However, the integration of the weak form is very time-consuming compared with standard meshbased methods such as the FEM. Therefore, EFG and FEM shall be coupled to employ EFG only in small areas, where a mesh is hardly to maintain and FEM in the rest of the domain. In the literature, there are most importantly two coupling approaches: the one introduced by Belytschko et al. in [3] and the other of Huerta et al. []. Preliminaries For a coupling of EFG and FEM, the domain Ω is decomposed into disjoint domains Ω el and Ω EFG, with the common boundary Γ EFG. The domain Ω el is discretized with standard quadrilateral finite elements. The union of all elements along Γ EFG is called the transition area Ω, so that Ω el may further be decomposed into the disjoint domains Ω FEM and Ω, connected by a boundary labeled Γ FEM ; clearly Ω FEM Ω EFG =. This situation can be depicted from Figure. Throughout this paper, consistency of first order is fulfilled by the set of meshbased, meshfree and coupled shape functions. This results in the ability of reproducing linear solutions exactly. Standard bilinear shape functions are chosen in the finite elements. The complete basis vector for the meshfree shape functions is p (x) = (, x) T in one dimension and p (x) = (, x, y) T in two dimensions. Coupling Approach of Huerta The coupling approach of Huerta [] considers the contributions of the meshbased FEM shape functions in the computation of the MLS shape functions by modified consistency conditions. The resulting, coupled set of shape functions is consistent up to the desired order.
4 original Huerta coupling modified Huerta coupling function values N * N EFG I FEM I * I EFG function values N * N EFG I FEM I * I EFG Ω FEM Γ FEM Ω * Γ EFG Ω EFG domain Ω Ω FEM Γ FEM Ω * Γ EFG Ω EFG domain Ω Fig. 2 Shape functions of the coupling approach of Huerta in the original and modified version. In the original approach, FEM nodes are placed in the standard way in the elements inside Ω FEM, however not in Ω. The corresponding meshbased shape functions of the FEM nodes remain unchanged, the coupling is considered only in the particles shape functions. EFG particles with corresponding supports Ω i may be arbitrarily distributed in Ω EFG and Ω. Then, the shape functions for the nodes and particles are computed as follows: FEM : N i = Ni FEM EFG : N i = N ( i EFG coupled : N i = p T (x) ) j I FEM j (x) p T (x j ) i I FEM i I EFG [M (x)] w (x x i ) p (x i ) i I, with I FEM = { i xi Ω } { } { FEM, I EFG = i Ω i Ω EFG and I } = i Ω i Ω el. In words, I EFG is the set of particles whose supports are fully inside Ω EFG, and I is the set of particles that have supports overlapping with elements. Ni FEM are the standard bilinear finite element shape functions, and Ni EFG are the standard MLS functions, defined in the previous section. Figure 2 shows the sets I FEM, I and I EFG and displays the resulting shape functions of this approach in a section of a one-dimensional domain with a regular node/particle distribution around the transition area Ω. Modification Instead of keeping the FEM shape function unchanged inside the transition area as in the original approach, one may additionally place EFG particles at the FEM node positions along { Γ FEM and superimpose the two shape functions at these nodes. That is, I FEM reduces to i xi Ω FEM \ Γ FEM}, and for the nodes along Γ FEM we define coupled : N i = ( p T (x) ) j I FEM j (x) p T (x j ) [M (x)] w (x x i ) p (x i ) + Ni FEM i : x i Γ FEM. Figure 2 shows the resulting shape functions of this approach. The important advantage of this modification is that smaller dilatation parameters are possible (although in this figure ρ i = 2.9 x has been taken for both approaches). For example, in the original approach and a regular distribution of particles and nodes in one dimension, one finds that for the regularity of the matrix M (x), dilatation parameters of ρ i > 2. x are required []. In contrast, with the modified approach ρ i >. x are sufficient. This holds analogously in multi-dimensional domains. The advantage can most importantly be realized for the solution of stabilized weak forms of non-linear partial differential equations. Numerical studies show for stabilized flow computations that sufficient and reliable stabilization is often only possible for about ρ i.7 x in regular node/particle distributions [9]. For the role of small dilatation parameters in the solution of stabilized problems see also [8].
5 original Belytschko coupling modified Belytschko coupling function values N * N EFG I FEM I * I EFG function values N * N EFG I FEM I * I EFG Ω FEM Γ FEM Ω * Γ EFG Ω EFG domain Ω Ω FEM Γ FEM Ω * Γ EFG Ω EFG domain Ω Fig. 3: Shape functions of the coupling approach of Belytschko in the original and modified version. Coupling Approach of Belytschko In the approach of Belytschko [3], meshfree and meshbased shape functions for the nodes i are computed and defined independently and coupled with help of a ramp function R (x). FEM nodes are placed in Ω el, in contrast to the coupling approach of Huerta, where they are placed in Ω FEM only. EFG particles with supports Ω i are distributed arbitrarily in Ω Ω EFG and may also be included inside Ω FEM. The latter will only affect the shape function inside Ω Ω EFG, i.e. particles with Ω i Ω FEM have no influence at all. Belytschko et al. define the ramp function as follows: R (x) =, x Ω FEM, x Ω EFG i I i (x), x Ω, I = { i xi Γ EFG }, i.e. it varies monotonically between and in Ω. The linear consistency of the resulting set of meshbased, meshfree and coupled shape functions is maintained [3]. The shape functions are defined as follows: FEM : N i = Ni FEM i I FEM EFG : N i = Ni EFG i I EFG coupled : N i = [ R (x)] Ni FEM + R (x) Ni EFG i I, { with I FEM = i x i Ω FEM, Ω { } { i Ω }, FEM I EFG = i Ω i Ω EFG and I = i Ω i Ω el }. The resulting shape functions are shown in Figure 3 and have been used with ρ i = 2.5 x in [3] for some test cases. One may note that EFG particles inside Ω FEM have an undesirable influence in Ω Ω EFG, leading to a numerically awkward form of the shape function (dashed line in Figure 3). Especially in the context of stabilization, it is neither clear from a mathematical viewpoint nor from an empirical viewpoint how to choose suitable stabilization parameters for these shape functions [8]. Modification In order to avoid these numerically awkward shape functions, this approach is modified slightly. The EFG particles are restricted to the area Ω Ω EFG. The resulting shape functions are then (for sufficiently small dilatation parameters) suitable to motivate the usage of standard stabilization parameters [8]. Consequently, a reliable stabilization is obtained. The resulting shape functions of the modified coupling approach of Belytschko is shown in the right part of Figure 3. NUMERICAL RESULTS One-dimensional Advection-Diffusion Equation The aim of this study is to show that the same order of convergence is obtained with all coupling approaches. The obtained convergence of order 2 in the L 2 -norm is equivalent to purely meshbased
6 n=2 n=35 n=49 L 2 error 3 convergence of the coupling schemes 2 Belyt., orig. Belyt., mod. Huerta, orig. Huerta, mod. pure FEM 4 Ω FEM Ω Ω EFG 5 2 node number Fig. 4 Convergence test of the coupling approaches. computations with bi-linear finite elements only. The test case is defined as follows. Approximate the stabilized one-dimensional advection-diffusion equation ( w + τc w ) ( c ũ ) x x K 2 ũ x f dω =, x (, ), 2 Ω with f = 2cπ cos (2πx) + 4Kπ 2 sin (2πx), c =, K = and boundary conditions u() = u() =. The exact solution is u (x) = sin (2πx). Throughout the convergence test, the ratio of the domains is kept constant at Ω F EM : Ω : Ω EF G = 6 : : 6, which may be seen in the left part of Figure 4. The right part shows the convergence results for the two original and modified coupling approaches of Huerta and Belytschko respectively. The rate of convergence remains the same than in pure FEM computations. The higher rates of convergence of a pure EFG computation cannot be reproduced, which is in agreement to [3] and []. Two-dimensional Incompressible Navier-Stokes Equations The SUPG/PSPG stabilized weak form of the instationary, incompressible Navier-Stokes equations is ( ) ui w i σ ij n j dγ = w i ϱ Γ Ω t + u u i w i j f i dω + σ ij dω + q u i dω x j Ω x j Ω x i ( u i + τ u k + ) [ ( ) q ui ϱ x k ϱ x i t + u u i j f i σ ] ij dω, x j x j ( u i x j with σ ij = pδ ij + µ + u j being the stress tensor. u i and p are approximations of the velocities and pressure, ϱ stands for the density, the external forces are f i = in the test cases presented herein. The Navier-Stokes equations are given in Arbitrary Lagrangian Eulerian (ALE) form, to consider possible mesh movements, such as in the rotor test case shown below. The first line of the NS equations is the Bubnov-Galerkin part, the second is the SUPG/PSPG stabilization, realized by a modification of the test function multiplied with the residual of the momentum equations. The boundary conditions are u i = g on the Dirichlet part and σ ij n j = h i on the Neumann part of the boundary. The stabilization parameters of the FEM shape functions are computed with τ = / ( ) t x i ) ( ) 2 ( ) 2 2 ui 4µ +, h el h 2 el as suggested in [7, 9]. This formula has also been applied to the meshfree and coupled shape functions, where the element length h el is replaced by the support length ρ. For details, see [8, 9].
7 a) b) c) u=, v= Ω Ω EFG Ω EFG Ω Ω EFG Ω Ω EFG.8.6 y vertical center velocity profile, Re=.4 reference coupled, 2x2 Ω coupled, 4x4 coupled, 6x6.2 pure FEM, 2x2 Ω FEM pure FEM, 4x4 pure FEM, 6x u Fig. 5: a) Driven cavity test case with domain decomposed into Ω FEM, Ω EFG and Ω, b) discretization with 4 4 node/particle distribution, c) convergence against reference solution. Driven Cavity The driven cavity test case is a standard benchmark for fluid problems. Reference solutions are given in []. A flow inside a quadratic domain Ω = (, ) (, ) with no-slip boundary conditions on the left, right and lower wall develops under a shear flow applied on the upper boundary until a stationary solution is reached. Figure 5 gives an outline of the problem and shows a discretization with 4 4 nodes/particles. The modified coupling approaches are applied, in order to enable the use of small dilatation parameters of ρ i =.3 x. For a Reynolds number of Re = convergence may not be reached for dilatation parameters ρ i 2. x, underlining the importance of the modified coupling versions. The results of both approaches are almost identical for this test case and therefore, only the result for Huerta coupling is shown. In Figure 5c) the convergence against the reference solution along the vertical center velocity profile is shown and coupled results are compared with pure FEM solutions. The Ω EFG holes are placed such that the center profile directly cuts through them. In between the EFG particles, linear interpolation has been applied for simplicity. One may see that coupling does not adversely affect the solution. Cylinder Flow The channel flow around a cylinder has been developed as a test case by Turek in [6]. The cylinder is placed slightly unsymmetrically in y-direction of the channel. For sufficiently high Reynolds numbers the well-known Kármán vortex street develops behind the cylinder. A quasi-stationary solution is obtained. Turek gives reference solutions for the lift and drag coefficients c L and c D of the cylinder. Figure 6 shows a sketch of the Kármán vortex street, together with the development of c D in time until a periodical solution is obtained. In the right part, the mesh around the cylinder is shown, together with the Ω EFG part, where particles are distributed. In the left part of Figure 7, the results for the drag coefficient obtained with the modified approach of Huerta and Belytschko and the pure FEM computation are compared ( t =.5). The horizontal lines show the limits, in which the exact value for the maximum of c D lies [6]. One may again see that the results are quite close together. The drag coefficient is slightly improved with the coupled approaches, the coupling approach of Huerta achieves somewhat better results than the approach of Belytschko. Both are slightly better than the pure FEM computation. Results for lift coefficient and pressure difference between the most left and right point on the cylinder surface are not shown here, because the drag coefficient turned out to be the most sensitive. For c L and p, the results between the coupled approaches and pure FEM are almost indistinguishable. The right part of Figure 7 shows the dependency of the drag coefficient on the time step t. The Strouhal number St = D v mt, with the diameter D =. of the cylinder, the average inflow from the left with v m = and the time T for 2 periods of c D (equals period of c L ), is displayed in the figure
8 a).4 y.2 Kármán vortex street c) b) x 3.3 drag coefficient 3.2 c D time [s] Fig. 6: Cylinder test case at Re =, a) the Kármán vortex street, b) development of the drag coefficient in time, c) discretization around the cylinder. drag coefficient drag coefficient c 3.22 D coupled Belytschko coupled Huerta pure FEM time [s] c 3.22 D Belyt., t=., St=.299 Belyt., t=.5, St=.2932 Belyt., t=.25, St=.2947 Huerta, t=., St=.299 Huerta, t=.5, St=.2932 Huerta, t=.25, St= time [s] Fig. 7: The left part compares the different modified coupling approaches and the pure FEM solution for t =.5, the right part shows the convergence in time of the coupled approaches. as well. A clear convergence against the reference Strouhal number of.295 St.35 may be seen and the amplitudes of the drag coefficient improve. Flow Around a Rotating Obstacle The previously described test cases verified the coupled fluid solver. However, they did not take advantage of the beneficial properties of this approach. The following test case replaces the cylinder of the previous test case by a rotating obstacle. Standard meshbased methods fail to give results due to the distortion of the mesh which must follow the rotation. However, this is no problem with the coupled fluid solver, where the rotating inner mesh and the stationary outer mesh are separated by a meshfree area. Figure 8 shows the discretization around the meshfree area and the resulting momentum around the rotor s center in dependence of the angle α of the inner mesh. Rather than giving all details of the test case settings, in here only the easiness with which the proposed coupled FEM/EFG approximations solve this complicated flow problem shall be shown. CONCLUSION In this paper, coupling approaches of meshfree and meshbased methods for the solution of flow problems are discussed. Meshfree methods are introduced considering the EFG method in particular. However, most conclusions can be applied to other meshfree methods as well. Standard coupling procedures are reviewed and modified such that the resulting shape functions are more suited for stabilization. This is a crucial aspect, because meshfree methods for stabilized, non-linear problems
9 moment around rotor center x 3.5 moment M angle α Fig. 8: Discretization around the rotating obstacle and resulting moment around the center of the rotor in dependence of the angle α. require certain attention [8, 9]. A convergence test for a one-dimensional test case shows that the coupled approaches achieve the same order of convergence than purely meshbased FEM calculations. The coupling approaches are then used for the solution of the two-dimensional incompressible Navier-Stokes equations. The coupled fluid solver is verified with standard test cases and is employed to solve a flow around a moving and rotating obstacle, showing the straightforward usability of this approach to complex flow problems. We conclude that coupled FEM/EFG approximations are a very promising tool for the simulation of complex flow problems. REFERENCES [] Belytschko, T.; Krongauz, Y.; Organ, D.; Fleming, M.; Krysl, P.: Meshless Methods: An Overview and Recent Developments. Comp. Methods Appl. Mech. Engrg., 39, 3 47, 996. [2] Belytschko, T.; Lu, Y.Y.; Gu, L.: Element-free Galerkin Methods. Internat. J. Numer. Methods Engrg., 37, , 994. [3] Belytschko, T.; Organ, D.; Krongauz, Y.: A Coupled Finite Element Element-free Galerkin Method. Comput. Mech., 7, 86 95, 995. [4] Brooks, A.N.; Hughes, T.J.R.: Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comp. Methods Appl. Mech. Engrg., 32, , 982. [5] Donea, J.; Huerta, A.: Finite Element Methods for Flow Problems. John Wiley & Sons, Chichester, 23. [6] Franca, L.P.; Hughes, T.J.R.: Two classes of mixed finite element methods. Comp. Methods Appl. Mech. Engrg., 69, 89 29, 988. [7] Fries, T.P.; Matthies, H.G.: Classification and Overview of Meshfree Methods. Informatikbericht-Nr. 23-3, Technical University Braunschweig, ( Brunswick, 23. [8] Fries, T.P.; Matthies, H.G.: A Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods. Informatikbericht-Nr. 24-, Technical University of Braunschweig, ( Brunswick, 24.
10 [9] Fries, T.P.; Matthies, H.G.: Meshfree Petrov-Galerkin Methods for the Incompressible Navier- Stokes Equations. In Meshfree Methods for Partial Differential Equations. (Griebel, M.; Schweitzer, M.A., Eds.), Springer Verlag, Berlin, 24 (to appear). [] Ghia, U.; Ghia, K.N.; Shin, C.T.: High-Re solutions for incompressible flow using the Navier- Stokes equations and a multi-grid method. J. Comput. Phys., 48, 387 4, 982. [] Huerta, A.; Fernández-Méndez, S.: Enrichment and Coupling of the Finite Element and Meshless Methods. Internat. J. Numer. Methods Engrg., 48, , 2. [2] Hughes, T.J.R.; Franca, L.P.: A new finite element formulation for computational fluid dynamics: VII. The Stokes problem with various well-posed boundary conditions: symmetric formulations that converge for all velocity/pressure spaces. Comp. Methods Appl. Mech. Engrg., 65, 85 96, 987. [3] Hughes, T.J.R.; Franca, L.P.; Balestra, M.: A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuška-Brezzi condition: a stable Petrov-Galerkin formulation of the Stokes problem accommodating equal-order interpolations. Comp. Methods Appl. Mech. Engrg., 59, 85 99, 986. [4] Hughes, T.J.R.; Franca, L.P.; Hulbert, G.M.: A new finite element formulation for computational fluid dynamics: VIII. The Galerkin/Least-squares method for advective-diffusive equations. Comp. Methods Appl. Mech. Engrg., 73, 73 89, 989. [5] Lancaster, P.; Salkauskas, K.: Surfaces Generated by Moving Least Squares Methods. Math. Comput., 37, 4 58, 98. [6] Schäfer, M.; Turek, S.: Benchmark Computations of Laminar Flow around a Cylinder. In Flow Simulation with High-Performance Computers II. (Hirschel, E.H., Ed.), Vieweg Verlag, Braunschweig, 996. [7] Shakib, F.; Hughes, T.J.R.; Johan, Z.: A new finite element formulation for computational fluid dynamics: X. The compressible Euler and Navier-Stokes equations. Comp. Methods Appl. Mech. Engrg., 89, 4 29, 99. [8] Tezduyar, T.E.; Mittal, S.; Ray, S.E.; Shih, R.: Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements. Comp. Methods Appl. Mech. Engrg., 95, , 992. [9] Tezduyar, T.E.; Osawa, Y.: Finite element stabilization parameters computed from element matrices and vectors. Comp. Methods Appl. Mech. Engrg., 9, 4 43, 2.
Meshfree Petrov-Galerkin Methods for the Incompressible Navier-Stokes Equations
Meshfree Petrov-Galerkin Methods for the Incompressible Navier-Stokes Equations Thomas-Peter Fries and Hermann Georg Matthies Institute of Scientific Computing, Technical University of Braunschweig, Brunswick,
More informationA Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods
ScientifiComputing A Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods Thomas-Peter Fries, Hermann G. Matthies Institute of Scientific Computing Technical University
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
More informationStabilized Formulations and Smagorinsky Turbulence Model for Incompressible Flows
WCCM V Fifth World Congress on Computational Mechanics July 7 2, 22, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Stabilized Formulations and Smagorinsky Turbulence Model for
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationSome remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes
Some remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes Stefano Micheletti Simona Perotto Marco Picasso Abstract In this paper we re-address the anisotropic recipe
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationStabilized finite elements for 3-D reactive flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2002/09/18 v1.01] Stabilized finite elements for 3-D reactive flows M. Braack, Th. Richter Institute
More informationThe Meshless Local Petrov-Galerkin (MLPG) Method for Solving Incompressible Navier-Stokes Equations
Copyright cfl 2001 Tech Science Press CMES, vol.2, no.2, pp.117-142, 2001 The Meshless Local Petrov-Galerkin (MLPG) Method for Solving Incompressible Navier-Stokes Equations H. Lin and S.N. Atluri 1 Abstract:
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationA variational multiscale stabilized formulation for the incompressible Navier Stokes equations
Comput Mech (2009) 44:145 160 DOI 10.1007/s00466-008-0362-3 ORIGINAL PAPER A variational multiscale stabilized formulation for the incompressible Navier Stokes equations Arif Masud Ramon Calderer Received:
More informationForchheimer law derived by homogenization of gas flow in turbomachines
Journal of Computational and Applied Mathematics 215 (2008) 467 476 www.elsevier.com/locate/cam Forchheimer law derived by homogenization of gas flow in turbomachines Gottfried Laschet ACCESS e.v., IntzestraYe
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6. Mathematik Preprint Nr. 228 A Variational Multiscale Method for Turbulent Flow Simulation with Adaptive Large Scale
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationSTABILIZED FEM SOLUTIONS OF MHD EQUATIONS AROUND A SOLID AND INSIDE A CONDUCTING MEDIUM
Available online: March 09, 2018 Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. Volume 68, Number 1, Pages 197 208 (2019) DOI: 10.1501/Commua1_0000000901 ISSN 1303 5991 http://communications.science.ankara.edu.tr/index.php?series=a1
More informationPoint interpolation method based on local residual formulation using radial basis functions
Structural Engineering and Mechanics, Vol. 14, No. 6 (2002) 713-732 713 Point interpolation method based on local residual formulation using radial basis functions G.R. Liu, L. Yan, J.G. Wang and Y.T.
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMulti-level hp-fem and the Finite Cell Method for the Navier-Stokes equations using a Variational Multiscale Formulation
Department of Civil, Geo and Environmental Engineering Chair for Computation in Engineering Prof. Dr. rer. nat. Ernst Rank Multi-level hp-fem and the Finite Cell Method for the Navier-Stokes equations
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationFinite calculus formulation for incompressible solids using linear triangles and tetrahedra
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 59:1473 1500 (DOI: 10.1002/nme.922) Finite calculus formulation for incompressible solids using linear triangles
More informationLocking in the incompressible limit: pseudo-divergence-free element free Galerkin
Locking in the incompressible limit: pseudodivergencefree element free Galerkin Yolanda Vidal Pierre Villon and Antonio Huerta Departament de Matemàtica Aplicada III, ES de Ingenieros de Caminos, Canales
More informationFinite element stabilization parameters computed from element matrices and vectors
Comput. Methods Appl. Mech. Engrg. 190 (2000) 411±430 www.elsevier.com/locate/cma Finite element stabilization parameters computed from element matrices and vectors Tayfun E. Tezduyar *, Yasuo Osawa Mechanical
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationConvectively Unstable Anti-Symmetric Waves in Flows Past Bluff Bodies
Copyright 29 Tech Science Press CMES, vol.53, no.2, pp.95-12, 29 Convectively Unstable Anti-Symmetric Waves in Flows Past Bluff Bodies Bhaskar Kumar 1 and Sanjay Mittal 1,2 Abstract: The steady flow past
More informationStabilization and shock-capturing parameters in SUPG formulation of compressible flows
Comput. Methods Appl. Mech. Engrg. 195 (006) 161 163 www.elsevier.com/locate/cma Stabilization and shock-capturing parameters in SUPG formulation of compressible flows Tayfun E. Tezduyar *, Masayoshi Senga
More informationAnalysis of Steady State Heat Conduction Problem Using EFGM
International Journal of Engineering and Management Research, Vol.-2, Issue-6, December 2012 ISSN No.: 2250-0758 Pages: 40-47 www.ijemr.net Analysis of Steady State Heat Conduction Problem Using EFGM Manpreet
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationElement diameter free stability parameters. for stabilized methods applied to uids
Element diameter free stability parameters for stabilized methods applied to uids by Leopoldo P. Franca Laboratorio Nacional de Computac~ao Cientica (LNCC/CNPq) Rua Lauro Muller 455 22290 Rio de Janeiro,
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationA FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS NUMBER FLUID FLOWS
Mathematical Models and Methods in Applied Sciences Vol. 12, No. 4 (2002) 525 539 c World Scientific Publishing Compan A FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationRESIDUAL-BASED TURBULENCE MODELS FOR INCOMPRESSIBLE FLOWS IN DOMAINS WITH MOVING BOUNDARIES RAMON CALDERER ELIAS DISSERTATION
RESIDUAL-BASED TURBULENCE MODELS FOR INCOMPRESSIBLE FLOWS IN DOMAINS WITH MOVING BOUNDARIES BY RAMON CALDERER ELIAS DISSERTATION Submitted in partial fulfillment of the requirements for the degree of Doctor
More informationA STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION
A STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION Eugenio Oñate International Centre for Numerical Methods in Engineering Universidad Politécnica
More informationA unified finite element formulation for compressible and incompressible flows using augmented conservation variables
. 2% mm d BB ELSEVIER Comput. Methods Appl. Mech. Engrg. 161 (1998) 229-243 Computer methods in applied mechanics and englneerlng A unified finite element formulation for compressible and incompressible
More informationSTABILIZED FINITE ELEMENT METHODS OF GLS TYPE FOR MAXWELL-B AND OLDROYD-B VISCOELASTIC FLUIDS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationSMOOTHED PARTICLE HYDRODYNAMICS METHOD IN MODELING OF STRUCTURAL ELEMENTS UNDER HIGH DYNAMIC LOADS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China SMOOTHE PARTICLE HYROYAMICS METHO I MOELIG OF STRUCTURAL ELEMETS UER HIGH YAMIC LOAS. Asprone *, F. Auricchio, A. Reali,
More informationSTABILIZED GALERKIN FINITE ELEMENT METHODS FOR CONVECTION DOMINATED AND INCOMPRESSIBLE FLOW PROBLEMS
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 29 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 1994 STABILIZED GALERIN FINITE ELEMENT METHODS FOR CONVECTION
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationCOMPUTATION OF MEAN DRAG FOR BLUFF BODY PROBLEMS USING ADAPTIVE DNS/LES
COMPUTATION OF MEAN DRAG FOR BLUFF BODY PROBLEMS USING ADAPTIVE DNS/LES JOHAN HOFFMAN Abstract. We compute the time average of the drag in two benchmark bluff body problems: a surface mounted cube at Reynolds
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationA hierarchical multiscale framework for problems with multiscale source terms
Available online at www.sciencedirect.com Comput. Methods Appl. Mech. Engrg. 197 (2008) 2692 2700 www.elsevier.com/locate/cma A hierarchical multiscale framework for problems with multiscale source terms
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationOptimal control in fluid mechanics by finite elements with symmetric stabilization
Computational Sciences Center Optimal control in fluid mechanics by finite elements with symmetric stabilization Malte Braack Mathematisches Seminar Christian-Albrechts-Universität zu Kiel VMS Worshop
More informationModélisation et discrétisation par éléments finis des discontinuités d interfaces pour les problèmes couplés
Modélisation et discrétisation par éléments finis des discontinuités d interfaces pour les problèmes couplés Application de XFEM pour les discontinuités mobiles Andreas Kölke Assistant Professor for Fluid-Structure
More informationUnsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation
Unsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation Mohamed S. Ebeida Carnegie Mellon University, Pittsburgh, PA 15213 Roger L. Davis and Roland W. Freund
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationA Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
Structural Engineering and Mechanics, Vol. 21, No. 3 (2005) 315-332 315 A Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationMinimal stabilization techniques for incompressible flows
Minimal stabilization tecniques for incompressible flows G. Lube 1, L. Röe 1 and T. Knopp 2 1 Numerical and Applied Matematics Georg-August-University of Göttingen D-37083 Göttingen, Germany 2 German Aerospace
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationSingularly Perturbed Partial Differential Equations
WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN 978-8-7378--9 MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics,
More informationA NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES
A NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES R. VERFÜRTH Abstract. In this note we look at constant-free a posteriori error estimates from a different perspective. We show that they can be interpreted
More informationOn spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P1 and Q1 finite elements
Volker John, Petr Knobloch On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P and Q finite elements Preprint No. MATH-knm-27/4 27.
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationLINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 445 461 c World Scientific Publishing Company LINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
More informationImproved discontinuity-capturing finite element techniques for reaction effects in turbulence computation
Comput. Mech. (26) 38: 356 364 DOI.7/s466-6-45- ORIGINAL PAPER A. Corsini F. Rispoli A. Santoriello T.E. Tezduyar Improved discontinuity-capturing finite element techniques for reaction effects in turbulence
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationReference values for drag and lift of a two-dimensional time-dependent ow around a cylinder
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2004; 44:777 788 (DOI: 10.1002/d.679) Reference values for drag and lift of a two-dimensional time-dependent ow around
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More information