DATA STRUCTURES FOR THE NONLINEAR SUBGRID STABILIZATION METHOD
|
|
- Wesley Newman
- 5 years ago
- Views:
Transcription
1 DATA STRUCTURES FOR TE NONLINEAR SUBGRID STABILIZATION METOD Isaac P. Santos Lucia Catabriga Departamento de Informática Universidade Federal do Espírito Santo - UFES Av. Fernando Ferrari, 54 - Vitória, ES, Brazil CEP: {ipsantos,luciac}@inf.ufes.br Regina C. Almeida Departamento de Mecânica Computacional Laboratório Nacional de Computação Científica - LNCC/MCT Av. Getúlio Vargas, Petrópolis, RJ, Brazil CEP: rcca@lncc.br Abstract. In this work we present a new implementation strateg for the Nonlinear Subgrid Stabilization (NSGS) method (Santos and Almeida, 27, 28). The NSGS method is a methodolog for approimating transport problems, based on a multiscale decomposition of the approimation space. The two-scale viscosit model is obtained b adding to the Galerkin formulation a nonlinear operator acting onl on the unresolved mesh scales. The resulting method is stable and does not require an tune-up parameter. owever, it requires the solution of a sstem twice as large as that related to the final resolved scale resolution. To avoid this drawback, we propose here a new strateg to approimate the effects of the unresolved scales on the resolved scales combined with an element-based and edge-based data structures. The resulting methodolog ields superior performance, which is demonstrated b a variet of numerical eperiments covering the regimes of dominant advection and dominant absorption. Kewords: Nonlinear subgrid viscosit, Edge-based data structure, Krlov iterative methods
2 . INTRODUCTION The numerical solution of transport equations using the Galerkin formulation ma ehibit global spurious oscillations when the diffusive term is smaller than the advective or reactive ones. More accurate and stable results can be obtained using stabilized or multiscale methods. In a general sense, each of these alternative methodologies are based on adding some tpe of linear stabilization to the standard Galerkin formulation and usuall depends on the definition of one or more tuning parameters. The multiscale approach consists of decomposing the approimation space into a resolved coarse scale and an unresolved subgrid scale spaces such that the weak form of the problem is split into two sub-problems: one for the coarse scales and one for the subgrid scales. From this point of view, the attempts to recover stabilit of the resolved scale solution ma be interpreted as a wa of improving the simulation b considering the effects of the smallest scales on the larger ones Brezzi and Marini (22). The Nonlinear Subgrid Stabilization (NSGS) method is a step towards the development of two-scale methods whose stabilit and convergence properties do not rel on tune-up parameters. Its first version was presented in (Santos and Almeida, 27) for advectiondiffusion problems and was based on adding a nonlinear artificial viscosit term onl to the subgrid scales. The amount of the subgrid viscosit depends on assuming an additive decomposition of the velocit field and solving the resolved scale equation at the element level, ielding a self adaptive method. A linear subgrid stabilization was kept when the resolved scale solution was accurate enough. The NSGS method was etended in (Santos and Almeida, 28) to advection-diffusionreaction problems and some improvements were introduced in smooth regions so as to recover the Galerkin formulation for regular resolved scale solutions. To improve the convergence of the linearization process, an alternative iterative algorithm to solve the nonlinear problem was also proposed. The resulting self adapted methodolog has ver good stabilit and convergence properties for a wide range of transport problems, but presents the drawback of requiring the solution of a sstem twice as large as that associated with the resolved scale resolution. To overcome this drawback, in this work we present a strateg for building the NSGS method using element-based (EBE) and edge-based (EDS) implementations (Catabriga et al., 998; Catabriga and Coutinho, 22). The present strateg locall approimates the effects of the subgrid scales on the resolved scales b performing a static condensation of the subgrid degrees of freedom of the associated macro element. Solution efficienc is improved b using an edge based data structure and the final algebraic sstem is solved b using the generalized minimal residual method (GMRES) (Saad, 995). The outline of this paper is as follows. We briefl describe the Nonlinear Subgrid Stabilization (NSGS) method in Section 2. Section 3 presents the proposed strateg of incorporating the subgrid effects into the resolved scale and Section 4 presents the related element-based (EBE) and edge-based (EDS) data structures. Numerical eamples are conducted in Section 5 and Section 6 concludes this paper b pointing out the main achievements, difficulties and perspectives. 2. TE NONLINEAR SUBGRID STABILIZATION METOD Consider the stead scalar advection-diffusion-reaction problem whose solution satisfies ǫ u + β u + σu = f in Ω; () u = on Ω, (2) where Ω R d (d = 2) is an open bounded domain with a Lipschitz boundar Ω and unit outward normal n, β is the velocit field, σ is the reaction coefficient, < ǫ is the
3 (constant) diffusion coefficient, and f is the source term. For simplicit, we onl consider homogeneous Dirichlet boundar conditions and it is assumed that β W, (Ω), σ L (Ω) and f L 2 (Ω). The weak formulation of the problem ()-(2) reads: find u (Ω) such that where B(u, v) = (f, v) v (Ω), (3) B(u, v) = ǫ( u, v) + (β u, v) + (σu, v) (4) and (, ) stands for the usual inner product in L 2 (Ω). It is also assumed that σ and that there eists a constant σ such that σ 2 β σ >. (5) The La-Milgram Lemma implies that the problem (3) has an unique solution. To define the discrete model, we consider a triangular partition T = {T } of the domain Ω where stands for the mesh parameter. From each triangle T T, four triangles are created b connecting the midpoints of the edges. We set h = /2 and denote b T h = {T h } the resulting finer triangulation. A two-level piecewise linear finite element approimation is defined b introducing the following two spaces: X = {u (Ω) u T P (T), T T }; (6) X h = {u h (Ω) u h Th P (T h ), T h T h }. (7) Moreover, we introduce an additional discrete space Xh decomposition holds: X h = X X h, X h, such that the following (8) where X is the resolved (coarse) scale space whereas Xh is the subgrid (fine) scale space. We define also the projection operator P : X h X so that Xh = (I P )X h, where I is the identit operator. Then, given u h X h and u X such that u h and u coincide in the coarse scale nodes, the space decomposition (8) implies that u h = (I P )u h = u h u, u h X h. The couple (X, X h ) will be referred to as the two-level P setting (see Fig. ) Guermond (999). T X nodes X h nodes Figure : Schematic representation of the two-level P setting
4 The Nonlinear Subgrid Stabilization Method (NSGS) (see Santos and Almeida (27, 28) for details) is given b Find u h X h such that B(u h, v h ) + T h T h D(u, u h, v h ) = (f, v h), v h X h, (9) where the non-linear operator D(u, u h, v h ) in each element T h T h is given b with D(u, u h, v h T ) = ξ(u ) u h v h dω, () h ξ(u ) = 2 µ(h) β h and µ(h) stands for the subgrid characteristic lehgth. For the two-level P setting, we use µ(h) = h/2. The subgrid velocit field β h is given b β h = { R(u ) u 2 u, if u ;, otherwise, () where R(u ) = ǫ u + β u + σu f. (2) Setting ξ(u ) = c b h, with c b > constant, we obtain the Linear Subgrid Stabilization (SGS) method (Guermond, 999, 2; Laton, 22). When u =, the formulation (9) recovers the standard Galerkin method. The NSGS method is solved using an iterative procedure for which the SGS solution is the initial guess (with c b = ) and ξ(u ) (or β h ) is delaed one iteration (see Algorithm ). The iterative NSGS process is defined b: Given u n h, we find u n h satisfing B(u n h, v h) + ) u;n h vh dω = (f, v h), v h X h, (3) T h T h T h where ) = c b,t h µ(h) is locall adjusted depending on the residual of the resolved scale solution. The iterative NSGS procedure is given b (Algorithm ), where u i and β;i h stand for the approimate resolved scale solution and the subgrid velocit field, respectivel, at iteration i; maiter is the maimum number of iterations and tol is the prescribed tolerance. Thus, the desired resolved solution u is obtained when the convergence of (Algorithm ) is attained. At each iteration, a linear sstem of the form A = b, with A nonsmmetric, has to be solved in order to determine u i (step 2). The sstem dimension is equal to the number of degrees freedom of the mesh T h. As mentioned before, this is the method major drawback and ma be overcome b appling the strateg proposed in the net section. 3. APPROXIMATING TE SUBGRID EFFECTS In this section we propose a new implementation strateg for the NSGS method, in which the non-linear formulation (9) is solved for the coarse mesh T. Let ũ be the approimation of the resolved scale solution u. ũ is obtained b locall approimating the effects of each
5 Algorithm Require: This algorithm receive as input u - calculated using SGS with c b = Ensure: The output is the approimate solution u : for (each element T h T h ) do 2: c T h c b 3: end for 4: i 5: repeat 6: i i + 7: for (each element T h T h ) do 8: determine β ;i 9: c i T h 2 [ h c i T h + 2 β ;i h : c b,th c i T h : end for 2: determine ( u i ) 3: until i maiter or ma u i ;j u i ;j tol j=,...,dof ] X h v 3 X v 3 T 4 h v 6 v 5 P v +v 3 2 T 4 h v 2 +v 3 2 T h T 3 h v v 2 v 4 T 2 h T h T 3 h v v 2 v +v 2 2 Figure 2: Macro element T T and its corresponding micro elements T h,t 2 h,t 3 h and T 4 h T h: relationship among the degrees of freedom of T and T h, due to the projection operator P. T 2 h
6 subgrid degree of freedom belonging to T T. To show how this approimation is performed, let us first notice that both linear and linearized operators in (3) can be built b summing up the contributions of each macro element T T. Each 6 6 local matri results from assembling the contributions of its corresponding micro elements Th i T, i =, 2, 3, 4 (see Fig. 2). Moreover, the degrees of freedom of the micro elements are connected to the resolved scale nodes of the related macro element due to the linearized operator D(u n ; u;n h, vh ), which prevents the use of an edge based structure. This can be clearl seen b rewriting it in the following wa: ) u;n h vh dω = ) (un h un ) (v h v ) dω T T T T T T = where I = III = T T T T ) un h v hdω; II = T ) un v hdω and IV = ( I II III + IV T ) un h v dω; ) un v dω. ), (4) The term I is similar to the diffusive term coming from the bilinear form B(u n h, v h). In the other three terms (II, III, IV ) the T resolved degrees of freedom are coupled with the T h degrees of freedom. For each Th i T, i =, 2, 3, 4, the former are obtained b projecting the latter into the macro element T using the projection operator P (Fig. 2). After assembling the contributions of each micro element T h T, the local problem associated with each T T ma be defined as K T U T = F T. (5) The local vector of unknowns U T = {u T h;, ut h;2, ut h;3, ut h;4, ut h;5, ut h;6 }t ma also be written as U T = {u T ;, ut ;2, ut ;3, ut h;4, ut h;5, ut h;6 }t, since u h and u coincide in the coarse scale nodes. The local sstem (5) can be partitioned in the following wa [ ][ ] [ ] A B U F = (6) C D U 2 F 2 where U = {u T ;, ut ;2, ut ;3 }t and U 2 = {u T h;4, ut h;5, ut h;6 }t. This local sstem shows quite well the coupling of the resolved scale and the subgrid nodes. Although we are interested in the resolved scale solution, such coupling should not be disregard. Besides, such local connection spreads to the macro element neighbors due to the support of the interpolation functions (twolevel P setting). This propert prevents the use of an EBE and EDS data structure. owever, we ma approimate the problem b disregarding the global subgrid connectivit and performing a static condensation of the unknowns U 2 at each macro element level. Thus, assuming that the matri D is nonsingular, the condensed local problem becomes K T U T = F T (7) where K T = (A BD C) is the 3 3 local matri which approimates the macro element matri of T T. The respective local vector is F T = (F BD F 2 ). After assembling the contributions of all T T, the following new global linear sstem is obtained KŨ = F, (8)
7 where K = nel A T= KT ; F = nel A T= FT ; Ũ = nel A T= UT. (9) In this epression, nel is the total number of macro elements of the coarse mesh T. The linear sstem (8) is then solved for the resolved scale unknown vector Ũ = {ũ ;,..., ũ ;Neq } t, where Neq is the total number of resolved scale nodes. Such approach ields a remarkable decrease of the computational cost, which can also be improved b using an EBE or EDS data structure. The resulting methodolog can be seen as an approimation of the original NSGS method. Although lacking a full mathematical eplanation, the numerical results are promising, as will be discussed in the following. 4. ELEMENT-BASED AND EDGE-BASED STRUCTURES Element-based and Edge-based structures have been etensivel used in finite element implementations, resulting in considerable improvements comparing to standard implementations. The success of this solution strateg requires an efficient implementation of matri-vector products and the choice of a suitable preconditioner. Generall, the edgebased data structure reduces processing time and requires around one half of the storage area to hold the coefficient matri when compared to an enhanced element-based implementation. We perform an implementation of the NSGS method using element-based and edge-based implementations. The conventional finite element data structure associates to each triangle e is its connectivit, that is, the mesh nodes I, J and K. In the edge-based data structure each edge s is associated to the adjacent elements e and f, thus to the nodes I, J, K and L, as shown in Figure 3. J K e s f L I Figure 3: Elements adjacent to edge s, formed b nodes I e J. Each element matri can be disassembled into its contributions to three edges, s, s + and s + 2, with connectivities IJ, JK and KI, that is, = + +, (2) } {{ } element e } {{ } edge s } {{ } edge s + } {{ } edge s + 2
8 where and are matri coefficients defined from (3). Thus, all the contributions belonging to edge s will be present in the adjacent elements e and f. The resulting edge matri is the sum of the corresponding sub-element matrices containing all the contributions to nodes I and J, that is, [ ] [ ] [ ] = +. (2) }{{} edge s }{{} element e }{{} element f Considering a conventional elementwise description of a given finite element mesh, the topological informations are manipulated, generating a new edge-based mesh description. Thus, the assembled global matri given in equation (9) ma be written now as, K = nedges A s= Ks where nedges is the total number of edges of the macro mesh T. The edge matri K s is obtained from the contributions of all the element matrices K T that share the edge s. In the element b element (EBE) implementation strateg, the coefficients of the global matri are stored in each macro element matri as defined b (9). The global matri K is stored in a compact form of size nel 9. On the other hand, in the edge-based (EDS) strateg the coefficients of the global matri, defined b (22), are also stored in a compact form of size nedge NUMERICAL RESULTS In this section we evaluate the numerical performance of the proposed Condensed NSGS on some academic test cases. In all of them the computational domain Ω = (, ) (, ) is discretized b using triangular meshes with 2 2, 4 4 and 8 8 cells, where each cell is subdivided in two triangles (macro elements). Besides comparing the resolved solution in the lights of the condensed and the original NSGS, we investigate the efficienc of both EBE and EDS implementations of the new methodolog. The tolerance of the GMRES algorithm is 7. It is used 5 vectors in the Krlov basis in eamples and 2 and vectors in the eample 3. The convergence of the NSGS method is attained for a prescribed tolerance (tol = 2 to the eamples and 2, and tol = 4 to the eample 3) or maiter = Eample : Advection-diffusion problem This eample simulates a two-dimensional advection dominated advection-diffusion problem with ǫ = 2, β = (, ) and f = σ =. The Dirichlet boundar conditions are given b {,.3 ; u(, ) = u(, ) = u(, ) = and u(, ) =, >.3. These conditions ield a solution with an internal laer in the direction of the velocit field starting at (.3,.) and an eponential eternal laer at =. Figure 4 shows the standard NSGS solution and the EBE/EDS-based condensed NSGS solution using the coarsest mesh (2 2). Thus, the mesh T is composed b 8 macro elements with 36 unknowns, while the mesh T h is formed b 32 elements with 52 unknowns. Note that the EBE/EDS-based condensed NSGS solution, obtained on T, is quite accurate. (22)
9 Table presents the computational performance of the element-based (EBE) and edgebased (EDS) data structures, respectivel, for the condensed NSGS method. In Tab., Neq is the number of unknowns related to coarse mesh and the CPU times are reported. The edgebased approach is about 6% faster than the element-based one. For T with 2 2 and 4 4, the CPU times for the standard NSGS (resolved for T h, direct solver) is 96.6 and seconds, respectivel. Thus, either using EBE or EDS, remarkable computational improvements are obtained with the condensed NSGS strateg (a) standard NSGS (b) EBE/EDS-based condensed NSGS Figure 4: Eample - NSGS method solutions. Table : Eample : Computational costs - EBE and EDS structures Mesh(T ) Neq Time (s) EBE EDS Eample 2: Advection-diffusion problem with source term This eample simulates a two-dimensional advection dominated advection-diffusion problem with ǫ = 9, β = (, ), σ = and a constant source term f =. We set u = on all boundar so that the eact solution is a 45 o slope, possessing parabolic laers at =, = and an eponential laer at =. The standard NSGS solution and the edge-based condensed NSGS solution using the mesh with 2 2 cells are shown in Fig. 5. The standard NSGS method presents some oscillations in the neighborhood of the parabolic laers ( = and = ). The condensed NSGS solution presents a diffusive behavior in the these laers and a small overshoot in =. owever, there is no oscillation at the eternals laers. Table 2 presents the computational performance for the condensed NSGS method considering both data structures. The edge-based approach is about 6% faster than the elementbased one, as in the eample. For T with 2 2 and 4 4, the CPU times for the standard NSGS (resolved for T h, direct solver) is and seconds, respectivel.
10 .2.2 (a) standard NSGS (b) EBE/EDS-based condensed NSGS Figure 5: Eample 2 - NSGS method solutions. Table 2: Eample 2: Computational costs: EBE and EDS structures Mesh(T ) Neq Time (s) EBE EDS Eample 3: Diffusion - reaction problem In this eample, we consider the singular perturbed case where ǫ = 2, σ = and f =.5, with the following boundar conditions: u(, ) = u(, ) = and u(, ) = u(, ) =. The Galerkin solution obtained on macro mesh 2 2 is shown in Fig. 6(a), which presents some locallized oscillations in the neighborhood of the eternal laers. Fig. 6(b) and 6(c) show the corresponding standard and condensed NSGS solution, respectivel. Although some oscillations still remain at eternal laers, the edge-based condensed NSGS solution is more accurate than the NSGS standard one. The computational performance for the condensed EBE and EDS based NSGS implementation are presented in Tab. 3. For this eample the element-based data structure is more effective than the edge-based, that is, the EBE CPU time for all meshes is smaller than the EDS CPU time. Usuall, EDS is more effective than EBE due to the matri-vector product float-point operations. owever, as the number of the GMRES iterations is ver small in the present eample, the solution required fewer number of matri-vector products floatpoint operations. Thus, the computational costs are mainl concentrated on the generation of the edge-based mesh. Table 3: Eample 3: Computational costs: EBE and EDS structures Mesh(T ) Neq() Time (s) EBE EDS
11 (a) Galerkin (b) standard NSGS (c) EBE/EDS-based condensed NSGS Figure 6: Eample 3 - Galerkin and NSGS methods solutions.
12 6. CONCLUSIONS We propose a new implementation strateg for the Nonlinear Subgrid Stabilization (NSGS) method. The NSGS method is a free parameter two-level subgrid method with good stabilit and convergence properties, owing to a subgrid viscosit that acts onl on the unresolved scales. It requires the solution of linear sstems associated with a mesh with characteristic length h = /2 to obtain a solution whose resolution is. The Condensed NSGS method proposed here introduces a procedure to overcome this drawback. This is performed b locall approimating the effects of the unresolved scale, ielded b a static condensation of the not null subgrid degrees of freedom at the macro element level. Then, the Condensed NSGS method can be seen as an approimation of the original NSGS method. Although lacking a full mathematical eplanation, such simple procedure results in a remarkable computational gain. Moreover, some numerical eperiments indicate that it ields etra regularization. Since the resolved scale solution is obtained directl for the macro finite element mesh, we incorporate the EBE and EDS data structures, ielding additional computational efficienc. Contrasting other techniques that use tune-up parameters to improve accurac, the present free Condensed EBE/EDS NSGS method combines simplicit with computational efficienc. Some important issues associated with the convergence of the nonlinear procedure, the mathematical and numerical analsis, application to nonlinear and transient problems remain to be investigated in forthcoming works. Acknowledgements The first author would like to thank the fellowship provided b the Brazilian Government, through the Agenc CNPq/Brazil Pos-Doc Junior (PDJ) 5649/27-3. The second and third authors acknowledge the financial support from the Agenc CNPq/Brazil, contract/grant numbers 6265/26-5 and 3573/27-, respectivel. REFERENCES Brezzi, F. & Marini, L. D., 22. Augmented spaces, two-level methods and stabilizing subgrids. International Journal for Numerical Methods in Fluids, vol. 4, pp Catabriga, L. & Coutinho, A., 22. Implicit SUPG solution of Euler equations using edgebased data structures. Computer Methods in Applied Mechanics and Engineering, vol. 9, pp Catabriga, L., Martins, M. A. D., Coutinho, A. L. G. A., & Alves, J. L. D., 998. Clustered edge-b-edge preconditioners for non-smmetric finite element equations. In Onate, E. & Idelsohn, S. R., eds, Computational Mechanics: New Trends and Applications, pp. 4, Buenos Aires, Argentina. 4th World Congress on Computational Mechanics, CD-ROM. Guermond, J.-L., 999. Stabilization of galerkin approimations of transport equation b subgrid modeling. Mathematical Modelling and Numerical Analsis, vol. 33, pp Guermond, J.-L., 2. Subgrid stabilization of galerkin approimations of linear monotone operators. IMA Journal of Numerical Analsis, vol. 2, pp Laton, W. J., 22. A connection between subgrid scale edd viscosit and mied methods. Appl. Math. and Comput., vol. 33, pp Saad, Y., 995. Iterative methods for sparse linear sstems. PWS Publishing Compan.
13 Santos, I. P. & Almeida, R. C., 27. A nonlinear subgrid method for advection-diffusion problems. Computer Methods in Applied Mechanics and Engineering, vol. 96, pp Santos, I. P. & Almeida, R. C., 28. A free parameter subgrid scale model for transport problems. Computer Methods in Applied Mechanics and Engineering.
Krylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
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 informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
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 informationIX CONGRESSO BRASILEIRO DE ENGENHARIA E CIÊNCIAS TÉRMICAS. 9th BRAZILIAN CONGRESS OF THERMAL ENGINEERING AND SCIENCES
IX CONGRESSO BRASILEIRO DE ENGENHARIA E CIÊNCIAS TÉRMICAS 9th BRAZILIAN CONGRESS OF THERMAL ENGINEERING AND SCIENCES Paper CIT02-0184 A NUMERICAL METHODOLOGY TO SOLVE THERMAL POLLUTION PRO- BLEMS R. S.
More informationImproved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations
American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More informationAN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures AB = BA = I,
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 7 MATRICES II Inverse of a matri Sstems of linear equations Solution of sets of linear equations elimination methods 4
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationEvaluating the LCD algorithm for solving linear systems of equations arising from implicit SUPG formulation of compressible
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2003; 00: 22 [Version: 2002/09/8 v2.02] Evaluating the LCD algorithm for solving linear systems of equations arising
More informationCOMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER PARABOLIC EQUATIONS. Lili Ju. Xinfeng Liu.
DISCRETE AND CONTINUOUS doi:13934/dcdsb214191667 DYNAMICAL SYSTEMS SERIES B Volume 19, Number 6, August 214 pp 1667 1687 COMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER
More informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More informationRANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More information2.2 SEPARABLE VARIABLES
44 CHAPTER FIRST-ORDER DIFFERENTIAL EQUATIONS 6 Consider the autonomous DE 6 Use our ideas from Problem 5 to find intervals on the -ais for which solution curves are concave up and intervals for which
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 informationA high order Padé ADI method for unsteady convection-diffusion equations
Center for Turbulence Research Annual Research Briefs 2005 85 A high order Padé ADI method for unstead convection-diffusion equations B D. You 1. Motivation and objectives The unstead convection-diffusion
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
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 informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationFRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG. Lehrstuhl für Informatik 10 (Systemsimulation)
FRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG INSTITUT FÜR INFORMATIK (MATHEMATISCHE MASCHINEN UND DATENVERARBEITUNG) Lehrstuhl für Informatik 1 (Sstemsimulation) Efficient hierarchical grid coarsening
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationDomain decomposition on different levels of the Jacobi-Davidson method
hapter 5 Domain decomposition on different levels of the Jacobi-Davidson method Abstract Most computational work of Jacobi-Davidson [46], an iterative method suitable for computing solutions of large dimensional
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationSEPARABLE EQUATIONS 2.2
46 CHAPTER FIRST-ORDER DIFFERENTIAL EQUATIONS 4. Chemical Reactions When certain kinds of chemicals are combined, the rate at which the new compound is formed is modeled b the autonomous differential equation
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 information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationParameterized Joint Densities with Gaussian and Gaussian Mixture Marginals
Parameterized Joint Densities with Gaussian and Gaussian Miture Marginals Feli Sawo, Dietrich Brunn, and Uwe D. Hanebeck Intelligent Sensor-Actuator-Sstems Laborator Institute of Computer Science and Engineering
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationThe approximation of piecewise linear membership functions and Łukasiewicz operators
Fuzz Sets and Sstems 54 25 275 286 www.elsevier.com/locate/fss The approimation of piecewise linear membership functions and Łukasiewicz operators József Dombi, Zsolt Gera Universit of Szeged, Institute
More informationThree-Dimensional Explicit Parallel Finite Element Analysis of Functionally Graded Solids under Impact Loading. Ganesh Anandakumar and Jeong-Ho Kim
Three-Dimensional Eplicit Parallel Finite Element Analsis of Functionall Graded Solids under Impact Loading Ganesh Anandaumar and Jeong-Ho Kim Department of Civil and Environmental Engineering, Universit
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationAn elementary model for the validation of flamelet approximations in non-premixed turbulent combustion
Combust. Theor Modelling 4 () 89. Printed in the UK PII: S364-783()975- An elementar model for the validation of flamelet approimations in non-premied turbulent combustion A Bourliou and A J Majda Département
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationDISCONTINUOUS GALERKIN DISCRETIZATION WITH EMBEDDED BOUNDARY CONDITIONS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol.3(23), No., pp.35 58 c 23 Editorial Board of the Journal Computational Methods in Applied Mathematics Joint Co. Ltd. DISCONTINUOUS GALERKIN DISCRETIZATION
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationDemonstrate solution methods for systems of linear equations. Show that a system of equations can be represented in matrix-vector form.
Chapter Linear lgebra Objective Demonstrate solution methods for sstems of linear equations. Show that a sstem of equations can be represented in matri-vector form. 4 Flowrates in kmol/hr Figure.: Two
More informationRELATIONS AND FUNCTIONS through
RELATIONS AND FUNCTIONS 11.1.2 through 11.1. Relations and Functions establish a correspondence between the input values (usuall ) and the output values (usuall ) according to the particular relation or
More informationA Padé approximation to the scalar wavefield extrapolator for inhomogeneous media
A Padé approimation A Padé approimation to the scalar wavefield etrapolator for inhomogeneous media Yanpeng Mi, Zhengsheng Yao, and Gary F. Margrave ABSTRACT A seismic wavefield at depth z can be obtained
More informationUNCORRECTED SAMPLE PAGES. 3Quadratics. Chapter 3. Objectives
Chapter 3 3Quadratics Objectives To recognise and sketch the graphs of quadratic polnomials. To find the ke features of the graph of a quadratic polnomial: ais intercepts, turning point and ais of smmetr.
More informationA Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 A Meshless Simulations for D Nonlinear Reaction-diffusion Brusselator Sstem Ahmad Shirzadi, Vladimir Sladek, Jan Sladek 3 Abstract: This paper
More informationLESSON #42 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART 2 COMMON CORE ALGEBRA II
LESSON #4 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART COMMON CORE ALGEBRA II You will recall from unit 1 that in order to find the inverse of a function, ou must switch and and solve for. Also,
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
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 informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More informationBifurcations of the Controlled Escape Equation
Bifurcations of the Controlled Escape Equation Tobias Gaer Institut für Mathematik, Universität Augsburg 86135 Augsburg, German gaer@math.uni-augsburg.de Abstract In this paper we present numerical methods
More informationChapter 14 Truss Analysis Using the Stiffness Method
Chapter 14 Truss Analsis Using the Stiffness Method Structural Mechanics 2 ept of Arch Eng, Ajou Univ Outline undamentals of the stiffness method Member stiffness matri isplacement and force transformation
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
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 informationDeterminants. We said in Section 3.3 that a 2 2 matrix a b. Determinant of an n n Matrix
3.6 Determinants We said in Section 3.3 that a 2 2 matri a b is invertible if and onl if its c d erminant, ad bc, is nonzero, and we saw the erminant used in the formula for the inverse of a 2 2 matri.
More informationExact Differential Equations. The general solution of the equation is f x, y C. If f has continuous second partials, then M y 2 f
APPENDIX C Additional Topics in Differential Equations APPENDIX C. Eact First-Order Equations Eact Differential Equations Integrating Factors Eact Differential Equations In Chapter 6, ou studied applications
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More information6 = 1 2. The right endpoints of the subintervals are then 2 5, 3, 7 2, 4, 2 9, 5, while the left endpoints are 2, 5 2, 3, 7 2, 4, 9 2.
5 THE ITEGRAL 5. Approimating and Computing Area Preliminar Questions. What are the right and left endpoints if [, 5] is divided into si subintervals? If the interval [, 5] is divided into si subintervals,
More informationREORDERING EFFECTS ON PRECONDITIONED KRYLOV METHODS IN AMR SOLUTIONS OF FLOW AND TRANSPORT
Proc. of the XXVII Iberian Latin-American Congress on Computational Methods in Engineering CILAMCE 2006 Brazilian Assoc. for Comp. Mechanics & Latin American Assoc. of Comp. Methods in Engineering Belém,
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationFATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES FATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS L. Samuelsson * * Volvo Aero Corporation, Sweden Kewords:
More informationc 2016 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 38, No. 2, pp. A765 A788 c 206 Societ for Industrial and Applied Mathematics ROOTS OF BIVARIATE POLYNOMIAL SYSTEMS VIA DETERMINANTAL REPRESENTATIONS BOR PLESTENJAK AND MICHIEL
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 information5.6. Differential equations
5.6. Differential equations The relationship between cause and effect in phsical phenomena can often be formulated using differential equations which describe how a phsical measure () and its derivative
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationElliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations
5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have
More informationAdditional Topics in Differential Equations
0537_cop6.qd 0/8/08 :6 PM Page 3 6 Additional Topics in Differential Equations In Chapter 6, ou studied differential equations. In this chapter, ou will learn additional techniques for solving differential
More informationA DLM/FD Method for Simulating Particle Motion in Viscoelastic Fluid
A Method for Simulating Particle Motion in Viscoelastic Fluid Tsorng-Wha Pan and Roland Glowinski Department of Mathematics Universit of Houston Houston, Teas Ma 3, 6 International Workshop on Fluid-Structure
More informationA Hybrid Numerical Method for High Contrast. Conductivity Problems. Liliana Borcea, and George C. Papanicolaou y. June 10, 1997.
A Hbrid Numerical Method for High Contrast Conductivit Problems Liliana Borcea, and George C. Papanicolaou June, 997 Abstract We introduce and test etensivel a numerical method for computing ecientl current
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More informationChapter 9 Vocabulary Check
9 CHAPTER 9 Eponential and Logarithmic Functions Find the inverse function of each one-to-one function. See Section 9.. 67. f = + 68. f = - CONCEPT EXTENSIONS The formula = 0 e kt gives the population
More informationAffine transformations
Reading Optional reading: Affine transformations Brian Curless CSE 557 Autumn 207 Angel and Shreiner: 3., 3.7-3. Marschner and Shirle: 2.3, 2.4.-2.4.4, 6..-6..4, 6.2., 6.3 Further reading: Angel, the rest
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More information