Characteristics Method Using Cubic Spline Interpolation for Advection Diffusion Equation
|
|
- Rosemary Russell
- 6 years ago
- Views:
Transcription
1 Characteristics Method Using Cubic Spline Interpolation for Advection Diffusion Equation Tung-Lin Tsai 1 ; Jinn-Chuang Yang, M.ASCE 2 ; and Liang-Hsiung Huang, A.M.ASCE 3 Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. Abstract: The characteristics method by using the cubic-spline interpolation is comparable to the Holly Preissmann scheme in solving the advection portion of the advection diffusion equation. In order to conduct a cubic-spline interpolation, an additional constraint must be specified at each endpoint. In general, four types of endpoint constraints are available, i.e., the first derivative, second derivative, quadratic, and not-a-knot constraints. The goal of this paper is to examine each type of endpoint constraints. Two hypothetical cases are used to conduct the investigation. Among the four types of constraints examined herein, the not-a-knot constraint and the first derivative constraint with high-order finite difference approximation yield the best results. However, as far as accuracy and simple implementation are concerned the not-a-knot constraint should be the best choice in solving the advection diffusion equation DOI: /ASCE :6580 CE Database subject headings: Advection; Diffusion; Interpolation. Introduction The advection diffusion equation is one of the governing equations in solving mass transport in river, lakes, oceans, and groundwater. A variety of numerical methods are available to solve the advection diffusion equation. However, none of them is considered satisfactory. Yet the primary difficulty arises from the combined hyperbolic advection and parabolic diffusion nature. The split-operator approach in which the advection and diffusion processes are separately computed using different numerical schemes has been pursued by many numerical modelers. In the splitoperator approach, since the diffusion process can be accurately computed by several numerical methods, such as the Crank Nicholson central difference scheme, the Crank Nicholson Galerkin finite element method, and some others. Hence, the accuracy of solving the advection diffusion equation will mainly depend on the computed results of the advection process. The characteristics-based Holly Preissmann scheme Holly and Preissmann 1977 with a fourth-order Hermitic interpolation between grid points is considered one of the best methods because of less numerical oscillation and damping in modeling the advection process. The Holly Preissmann scheme was successfully applied to a number of hydraulic problems Holly and Usseglio-Polatera 1984; Toda and Holly 1987; Holly and Rahuel 1990; Yang and Hsu It was further examined by Glass and Rodi In spite of all the successes, the Holly Preissmann 1 Research Assistant Professor, Natural Hazard Mitigation Research Center, National Chiao Tung Univ., Hsinchu, Taiwan 30010, R.O.C. 2 Professor, Dept. of Civ. Engrg., National Chiao Tung Univ., Hsinchu, Taiwan 30010, R.O.C. 3 Professor, Dept. of Civ. Engrg., National Taiwan Univ., Taipei, Taiwan 10617, R.O.C. Note. Discussion open until November 1, Separate discussions must be submitted for individual papers. To extend the closing date by one month, a written request must be filed with the ASCE Managing Editor. The manuscript for this technical note was submitted for review and possible publication on November 12, 2002; approved on October 30, This technical note is part of the Journal of Hydraulic Engineering, Vol. 130, No. 6, June 1, ASCE, ISSN / 2004/ /$ scheme achieves its high accuracy at the expense of lengthy computation time for tackling spatial derivatives with auxiliary equations. To obviate the need of solving the first spatial derivatives of the auxiliary equation, an alternative with the use of cubic-spline interpolation function was proposed by Schohl and Holly 1991, Karpik and Crokett 1997, and Stefanovic and Stefan The cubic-spline interpolation with the accuracy comparable to the Hermitic interpolation may achieve some reduction in computational time and memory. The cubic-spline interpolation function must pass through each given data location or node and be continuous in its first and second derivatives at interior nodes. Thus, it is necessary to specify two additional endpoint i.e., upstream and downstream boundaries constraints. In general, there exist four different types of endpoint constraints, i.e., first derivative, second derivative, quadratic, and not-a-knot endpoint constraints e.g., Gerald and Wheatley 1994; Knott 1999; Kvasov The natural cubicspline interpolation with neglect of second derivative at the endpoints is most frequently used e.g., Schohl and Holly 1991; Karpik and Crokett 1997; Stefanovic and Stefan However, in all their cases, the computational domains are so large that the influences of boundary are irrelevant, and therefore any type of endpoint conditions may obtain the same result. Further, the natural cubic-spline interpolation may yield the end cubics to approach linearity at their extremities and flatten the interpolating curve too much due to neglect of second derivative at endpoints. The main goal of this paper is to examine which type of endpoint constraints in solving the advection diffusion equation is the best, especially when the influences of boundary are not negligible. Numerical Framework of Advection Diffusion Simulation The one-dimensional advection diffusion equation can be expressed as U t D 2 2 (1) 580 / JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004
2 Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. Fig. 1. Schematic grid diagram for cubic-spine interpolation where the scalar function (x,t) may represent, for example, temperature or concentration at position x and time t with flow velocity U and diffusion coefficient D. By applying the splitoperator approach Yanenko 1968; Sobey 1984; Holly and Ussglio-Polatera 1984; Lin and Falconer 1997, Eq. 1 may be approximated with combination of two kinds of transport processes: a pure advection process described by t U 0 (2) and a pure diffusion process written as follows: t D 2 2 (3) The advantage of the split-operator approach is that it makes the choice of the most appropriate scheme for each transport process, whereas the disadvantage is the violation of the fact that advection and diffusion occur simultaneously. Solution of Advection Process The characteristics method with interpolating approximation of solution is widely used for solving the advection process. Eq. 2 can be written in terms of total derivative, i.e., D s /Dt, as D s 0 (4) Dt along dx U (5) dt Integration of the above equations yields h f (6) along x h f UtCr (7) where and t represent the grid size and time step, respectively see Fig. 1. Cr is the Courant number. The schematic diagram of the characteristic trajectory is shown in Fig. 1. h is the unknown concentration of grid point h at time level n, which is to be solved. f is the concentration of grid point f at time level n 1, in which concentrations of all grid points are known. Since, in general, the foot of the trajectory, x f, does not coincide with grid points, one must employ some form of interpolation to obtain the required concentration, f. Thus, the accuracy of the solution is significantly related to the choice of the interpolation function. One may develop a cubic-spline interpolation function for evaluating h corresponding to all the known concentrations at time level n1, that is, n1 i, i0,1,...,m shown in Fig. 1. In the cubic-spline interpolation, the second derivative is a continuous piecewise linear function and can be expressed as x n1 i1 x i S 1 S i1, x i i,x i1 (8) i where i i1 i, i0,...,m1; x i is the coordinate of grid point; n1 (x) denotes the function of second derivative of concentration; and S i and S i1 represent the second derivative at grid points i and i1, respectively. Integrating Eq. 8 twice and substituting the values of concentrations at grid points i and i1 yields the expression for the cubic function n1 (x) on i,x i1 as follows: n1 i1 3 i 3 S i S 6 i1 i 6 i x i1 i 2 n1 i i S i 6 2 n1 i i1 S i1 x i (9) 6 i The second derivative of grid points, S i, in Eq. 9 can be found by applying the continuity of the first derivative at interior nodes as follows: i1 S i1 2 i1 i S i i S i1 6 i1 n1 n1 i i n1 n1 i1 i i1, i1,...,m1 (10) The system of Eq. 10 is underdetermined as it contains only M1 equations for finding M1 unknowns. To complete this system two additional constraints i.e., two endpoint constraints for S 0 and S M are required. Four types of endpoint constraints Gerald and Wheatley 1994; Knott 1999; Kvasov 2000 may be used to produce the various types of cubic-spline interpolations. They are stated below. 1. First derivative endpoint constraints, 2 0 S 0 1 S n1 n (11a) and 2 M1 S M M2 S M1 6 M n1 n1 M M1 M1 (11b) 2. Second derivative endpoint constraints, S and S M 2 M 2 (12) 3. Quadratic endpoint constraints, S 0 S 1 and S M S M1 (13) 4. Not-a-knot endpoint constraints, and S S 1 0 S 2 1 S M M2 M1 S M1 M1 S M2 (14) M2 In Eqs. 11 and 12, the unknown first and second derivatives could be obtained by finite difference approximation corresponding to grid points near the boundary. However, the frequently JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004 / 581
3 Table 1. Root Mean Square Error of Various Endpoint Constraints in Pure Advection Test *0.02 *0.025 * Natural Quadratic Not-a-knot Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. First derivative st order nd order rd order th order th order Second derivative 1st order nd order rd order th order Note: *grid size and Courant numbercr0.3, 0.6, and 0.9. used natural cubic-spline interpolation simply takes S 0 S M 0 in Eq. 12 and neglects the second derivative at endpoints. Thus, it makes the end cubics approach linearity at their extremities. On the other hand, by substituting Eq. 13 into Eq. 9, one can clearly find that this endpoint constraint is equivalent to assuming that the end cubics approach to quadratic curves at their extremities. The not-a-knot endpoint constraint as shown in Eq. 14 represents the continuity of the third derivative at the nodes x 1 and x M1 de Boor In other words, two cubic segments that join at the node x 1 are adjacent parts of the same cubic curve. The identical result as that at node x 1 is also yielded at node x M1. Solution of Diffusion Step The second step in the split-operator approach is the diffusion transport of the advected concentration field as shown in Eq. 3. The computation of the diffusion process can be accurately executed by using a variety of finite-difference and finite-element numerical schemes. The Crank Nicholson second-order central difference scheme is used here. The discretized equation 3 for a uniform grid size can be written as n1 i a D t 2 2 n1 i1 2 n1 i n1 i1 n i1 2 n n i i1 (15) where a is a value of concentration obtained from the advection process at i point. This would form a tridiagonal system of algebraic equations, which can be solved by the efficient Thomas algorithm. Evaluation of Different Endpoint Constraints Calculation of Transport of Sine Function To investigate the computational performances of cubic-spine interpolation with different endpoint constraints, a nondimensional advection diffusion equation with uniform flow velocity, U, and diffusion coefficient, D, is considered as follows: * * t* * 1 2 * P * 2 (16) where */ 0, x*/l, t*ut/l are nondimensional concentration, coordinate, and time; 0 and L are characteristic concentration and length, respectively; and PUL/D represents Péclet number. Under the initial condition and the boundary conditions **,t*0sin2* (17) **0,t*sin2t*exp4 2 2 t*/p (18) **1,t*sin 21t*exp4 2 2 t*/p (19) the exact solution to 16 is **,t*sin 2*t*exp4 2 2 t/p (20) where represents the periodical number of sine function in the computational domain. In this study, 5 is used. First, the pure advection equation i.e., P ) is considered. Numerical results in terms of root mean square rms error with different Courant numbers and uniform grid of 0.02, 0.025, and 0.01 for 200 time steps by the natural, quadratic, and not-a-knot constraints as well as the first and second derivative constraints with different order finite difference approximation are, respectively, given in Table 1. Fig. 2 shows the simulated results of various constraints with Courant number of 0.3 and uniform grid of It is observed that the natural constraint induces larger rms error and numerical diffusion than the other constraints, excluding the first derivative constraint with first order finite difference approximation. This may be because the natural cubic-spline interpolation flattens the interpolating curve too much at the ends due to the assumption of neglecting second order derivative at endpoints, which makes the end cubics approach linearity at their extremities. The simulated results of not-a-knot constraint are better than those of the quadratic constraint. It is clearly found that the results of the first and second derivative endpoint constraints depend on the derivative approximations, especially for the former constraint. However, the more 582 / JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004
4 Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. Fig. 2. Computational results of various endpoint constraints for pure advection equation (*0.025, Cr0.3) accurate derivative approximations may not yield better results. The first and second derivative endpoint constraints seem to have better results when the fourth and second order finite difference approximations are applied, respectively. This may be due to the inconsistency of accuracy between the derivative approximations and the cubic-spline interpolation. For the first derivative endpoint constraint, the first order finite difference approximation produces larger numerical damping than the natural constraint. The second order finite difference approximation has worse simulated results than the quadratic constraint. The simulated results of the third and higher order finite difference approximations are comparable to those of the not-a-knot constraint. On the other hand, the second derivative constraint with finite difference approximation has similar simulated results in comparison with the quadratic constraint. In addition, one can find from Table 1 that the accuracy of simulated results for all the endpoint constraints increases significantly when the grid size decreases from to For the advection diffusion equation, the computational results of different endpoints constraints for 100 time steps and Table 2. RMS Error of Various Endpoint Constraints for Advection- Diffusion Equation (*0.025, Cr0.6) Peclet number 1/P /P0.002 Natural Quadratic Not-a-knot First derivative 1st order nd order rd order th order th order Second derivative 1st order nd order rd order th order Cr0.6 are shown in Table 2 for 1/P and 1/P0.002, respectively. For large Péclet number i.e., the advection dominates the transport process the performances of simulated results by all the constraints are similar to those of pure advection equation mentioned above. However, for small Péclet number, all the constraints yield comparable results that are better than those of large Péclet number. In this example, one can find that the natural constraint induces large numerical damping and rms error among those constraints considered. The simulated results of the quadratic and second derivative constraints with finite difference approximation are comparable. The not-a-knot constraint and the first derivative constraint with high-order finite difference approximation seem to have the best results in comparison with the other constraints. Since the results of the first derivative constraint are strongly related to the derivative approximation, as far as accuracy and simple implementation are concerned one may conclude that the not-a-knot constraint should be the best choice for solving the advection diffusion equation by characteristics method with cubic-spline interpolation. In the following example, further examination of the natural constraint and not-a-knot constraint for the transport of Gaussian distribution are conducted. Calculation of Transport of Gaussian Distribution With the following initial condition, exp exp (21) the solution of 1 is 1Ut2,t exp 2 2Dt 2 2 2Dt exp 2Ut Dt (22) where standard deviation of Gaussian distribution; and x 1 and x 2 represent the locations of peaks of Gaussian distribution. JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004 / 583
5 the natural condition. However, for the large diffusion coefficient, the computational results by the two constraints are almost identical to the exact solution. So, from the overall examination, the not-a-knot constraint should be suggested to replace the frequently used natural constraint as an alternative for solving the advection diffusion equation by characteristics method with cubic-spline interpolation. Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. Fig. 3. Computational results of Gaussian distribution for natural and not-a-knot endpoint constraints (U0.7 m/s,100 m, t 100 s): a D0 m 2 /s, b D0.35 m 2 /s, c D3.5 m 2 /s In this problem, the computational domain is x0,x b, x b is the location of downstream boundary where the domain of simulation is long enough so that the boundary effect can be ignored. A grid size of 100 m, a time step of 100 s, and a standard deviation of 120 m are used to conduct the simulation. The Gaussian distribution is advected for 50,000 s under a constant flow velocity U0.7 m/s. The two peak positions of this Gaussian distribution are at x m and x m, respectively. Thus, one can know that the Gaussian distribution is out of computational domain in the beginning of simulation. With the transport process, the Gaussian distribution passes through the upstream boundary (x0) into the computational domain. Numerical results by the natural and not-a-knot constraints are shown in Fig. 3 for D0, 0.35, and 3.5 m 2 /s, respectively. Fig. 3 demonstrates that for the advection dominated transport the natural constraint induces larger numerical damping than the not-aknot constraint. In spite of a little numerical oscillation, the nota-knot constraint has better simulated results in comparison with 584 / JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004 Conclusions The cubic-spline interpolation function not only needs to pass through each given data node but also needs to be continuous in its first and second derivatives at interior nodes, an additional constraint must therefore be specified at each endpoint. In general, the first derivative, second derivative, quadratic, and not-aknot constraints are four different types of endpoint constraints. Two hypothetical cases for transport of sine function and Gaussian distribution are used to investigate each type of constraint. The results show that the frequently used natural constraint induces larger numerical damping and rms error than the other constraints. The quadratic constraint and the second derivative constraint with finite difference approximation have comparable results. It is also found that the results of the first derivative constraint strongly depend on the derivative approximation. The nota-knot constraint and the first derivative constraint with higher order finite difference approximation seem to have the best results in comparison with the other constraints. In conclusion, as far as accuracy and simple implementation are concerned the not-a-knot constraint should be the best choice in solving the advection diffusion equation. Notation The following symbols are used in this technical note: Cr Courant number; D diffusion coefficient; L characteristic length; P Péclet number; S second derivative; U flow velocity component; computational grid interval; t time increment; standard deviation; concentration; 0 characteristic concentration; * nondimensional concentration; and periodical number of sine function in the computational domain. Subscripts i x-directional computational point index. Superscripts n time step index. References DeBoor, C A practical guide to spline. Springer, New York. Gerald, C. F., and Wheatley, P. O Applied numerical analysis. Addison-Wesley, New York.
6 Downloaded from ascelibrary.org by National Chiao Tung University on 04/30/14. Copyright ASCE. For personal use only; all rights reserved. Glass, J., and Rodi, W A higher order numerical scheme for scalar transport. Comput. Methods Appl. Mech. Eng., 31, Holly, F. M., Jr., and Preissmann, A Accurate calculation of transport in two dimensions. J. Hydraul. Div., Am. Soc. Civ. Eng., 10311, Holly, F. M., Jr., and Usseglio-Polatera, J Dispersion simulation in two-dimensional tidal flow. J. Hydraul. Eng., 1107, Holly, F. M., Jr., and Rahuel, J. L New numerical/physical framework for mobile-bed modeling, part I-Numerical and physical principles. J. Hydraul. Res., 284, Karpik, S. R., and Crockett, S. R Semi-Lagrangian algorithm for two-dimensional advection-diffusion equation on curvilinear coordinate meshes. J. Hydraul. Eng., 1235, Knott, G. D Interpolating cubic splines. Birkhauser, Boston. Kvasov, B. I Shape-preserving spline approximation. Scientific Publishing, Singapore. Lin, B., and Falconer, R. A Tidal flow and transport modeling using ULTIMATE QUICKEST scheme. J. Hydraul. Eng., 1234, Schohl, G. A., and Holly, F. M., Jr Cubic-spline interpolation in Lagrangian advection computation. J. Hydraul. Eng., 1172, Sobey, R. J Numerical alternatives in transient stream response. J. Hydraul. Eng., 1106, Stefanovic, D. L., and Stefan, H. G Accurate two-dimensional simulation of advective-diffusive-reactive transport. J. Hydraul. Eng., 1279, Toda, K., and Holly, F. M., Jr Hybrid numerical method for linear advection-diffusion. Microsoftw. Eng., 34, Yang, J. C., and Hsu, E. L On the use of the reach-back characteristics method for calculation of dispersion. Int. J. Numer. Methods Fluids, 12, Yaneko, N. N The method of fractional steps: The solution of problems of mathematical physics in several variables, Springer, New York. JOURNAL OF HYDRAULIC ENGINEERING ASCE / JUNE 2004 / 585
AN ACCURATE COMPUTATION FOR RAPIDLY VARIED FLOW IN AN OPEN CHANNEL
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 14, 361-374 (1992) AN ACCURATE COMPUTATION FOR RAPIDLY VARIED FLOW IN AN OPEN CHANNEL JINN-CHUANG YANG AND KUN-NAN CHEN Department of Civil Engineering,
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More information99 International Journal of Engineering, Science and Mathematics
Journal Homepage: Applications of cubic splines in the numerical solution of polynomials Najmuddin Ahmad 1 and Khan Farah Deeba 2 Department of Mathematics Integral University Lucknow Abstract: In this
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 13: The Lecture deals with:
The Lecture deals with: Some more Suggestions for Improvement of Discretization Schemes Some Non-Trivial Problems with Discretized Equations file:///d /chitra/nptel_phase2/mechanical/cfd/lecture13/13_1.htm[6/20/2012
More informationChapter 6. Petrov-Galerkin Formulations for Advection Diffusion Equation
Chapter 6 Petrov-Galerkin Formulations for Advection Diffusion Equation In this chapter we ll demonstrate the difficulties that arise when GFEM is used for advection (convection) dominated problems. Several
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationVideo 11.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 11.1 Vijay Kumar 1 Smooth three dimensional trajectories START INT. POSITION INT. POSITION GOAL Applications Trajectory generation in robotics Planning trajectories for quad rotors 2 Motion Planning
More informationCorrected/Updated References
K. Kashiyama, H. Ito, M. Behr and T. Tezduyar, "Massively Parallel Finite Element Strategies for Large-Scale Computation of Shallow Water Flows and Contaminant Transport", Extended Abstracts of the Second
More informationLecture 2.2 Cubic Splines
Lecture. Cubic Splines Cubic Spline The equation for a single parametric cubic spline segment is given by 4 i t Bit t t t i (..) where t and t are the parameter values at the beginning and end of the segment.
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More informationNUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD
18th Engineering Mechanics Division Conference (EMD007) NUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD Abstract S. Fu University of Texas at Austin, Austin,
More informationNumerical Simulation of Flow around a Spur Dike with Free Surface Flow in Fixed Flat Bed. Mukesh Raj Kafle
TUTA/IOE/PCU Journal of the Institute of Engineering, Vol. 9, No. 1, pp. 107 114 TUTA/IOE/PCU All rights reserved. Printed in Nepal Fax: 977-1-5525830 Numerical Simulation of Flow around a Spur Dike with
More informationAdvective and conservative semi-lagrangian schemes on uniform and non-uniform grids
Advective and conservative semi-lagrangian schemes on uniform and non-uniform grids M. Mehrenberger Université de Strasbourg and Max-Planck Institut für Plasmaphysik 5 September 2013 M. Mehrenberger (UDS
More information8 Piecewise Polynomial Interpolation
Applied Math Notes by R. J. LeVeque 8 Piecewise Polynomial Interpolation 8. Pitfalls of high order interpolation Suppose we know the value of a function at several points on an interval and we wish to
More informationHandout 4 - Interpolation Examples
Handout 4 - Interpolation Examples Middle East Technical University Example 1: Obtaining the n th Degree Newton s Interpolating Polynomial Passing through (n+1) Data Points Obtain the 4 th degree Newton
More informationA Comparison of Some Numerical Methods for the Advection-Diffusion Equation
Res Lett Inf Math Sci, 26, Vol1, pp49-62 Available online at http://iimsmasseyacnz/research/letters/ 49 A Comparison of Some Numerical Methods for the Advection-Diffusion Equation M Thongmoon 1 & R McKibbin
More informationAsymptotic Error Analysis
Asymptotic Error Analysis Brian Wetton Mathematics Department, UBC www.math.ubc.ca/ wetton PIMS YRC, June 3, 2014 Outline Overview Some History Romberg Integration Cubic Splines - Periodic Case More History:
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationA Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
Scientific Papers, University of Latvia, 2010. Vol. 756 Computer Science and Information Technologies 207 220 P. A Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
More informationInterpolation by Spline Functions
Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves
More informationCombos: Advection-Diffusion and operator splitting
96 Chapter 7 Combos: Advection-Diffusion and operator splitting Up to this point we have treated advection and diffusion as two separate processes and considered appropriate numerical techniques for solving
More informationlecture 8 Groundwater Modelling -1
The Islamic University of Gaza Faculty of Engineering Civil Engineering Department Water Resources Msc. Groundwater Hydrology- ENGC 6301 lecture 8 Groundwater Modelling -1 Instructor: Dr. Yunes Mogheir
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationCurve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur
Curve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur Email: jrkumar@iitk.ac.in Curve representation 1. Wireframe models There are three types
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationFour equations are necessary to evaluate these coefficients. Eqn
1.2 Splines 11 A spline function is a piecewise defined function with certain smoothness conditions [Cheney]. A wide variety of functions is potentially possible; polynomial functions are almost exclusively
More informationJ. Vira, M. Sofiev SILAM winter school, February 2013, FMI
Numerical aspects of the advection-diffusion equation J. Vira, M. Sofiev SILAM winter school, February 2013, FMI Outline Intro Some common requirements for numerical transport schemes Lagrangian approach
More informationTutorial Four Discretization Part 1
Discretization Part 1 4 th edition, Jan. 2018 This offering is not approved or endorsed by ESI Group, ESI-OpenCFD or the OpenFOAM Foundation, the producer of the OpenFOAM software and owner of the OpenFOAM
More informationACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION
Scientific Research of the Institute of Mathematics and Computer Science ACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION Ewa Węgrzyn-Skrzypczak, Tomasz Skrzypczak Institute of Mathematics, Czestochowa
More informationNatural Quartic Spline
Natural Quartic Spline Rafael E Banchs INTRODUCTION This report describes the natural quartic spline algorithm developed for the enhanced solution of the Time Harmonic Field Electric Logging problem As
More informationAbstract. Introduction. Numerical Techniques for Coextrusion Simulation
COMPARISON OF MESH PARTITIONING TECHNIQUE WITH LEVEL-SET METHOD FOR COEXTRUSION SIMULATION Mahesh Gupta 1, 2 1. Michigan Technological University, Houghton, MI 49931 2. Plastic Flow, LLC, Hancock, MI 49930
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 5:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture5/5_1.htm 1 of 1 6/20/2012 12:22 PM The Lecture deals with: Explicit and Implicit Methods file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture5/5_2.htm
More informationRational Bezier Curves
Rational Bezier Curves Use of homogeneous coordinates Rational spline curve: define a curve in one higher dimension space, project it down on the homogenizing variable Mathematical formulation: n P(u)
More informationA Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete Data
Applied Mathematical Sciences, Vol. 1, 16, no. 7, 331-343 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/1.1988/ams.16.5177 A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete
More informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More informationADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM. Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s
ADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s. 14-24 ABSTRACT In recent years considerable interest has been shown in the development
More informationMathematical Model and Surface Deviation of Cylindrical Gears With Curvilinear Shaped Teeth Cut by a Hob Cutter
Jui-Tang Tseng Graduate Student Chung-Biau Tsay Professor Department of Mechanical Engineering, National Chiao Tung University, Hsinchu, Taiwan 3000, Republic of China Mathematical Model Surface Deviation
More information08 - Designing Approximating Curves
08 - Designing Approximating Curves Acknowledgement: Olga Sorkine-Hornung, Alexander Sorkine-Hornung, Ilya Baran Last time Interpolating curves Monomials Lagrange Hermite Different control types Polynomials
More informationThe Application of Spline Functions and Bézier Curves to AGV Path Planning
IEEE ISIE 2005, June 20-23, 2005, Dubrovnik, Croatia The Application of Spline Functions and Bézier Curves to AGV Path Planning K. Petrinec, Z. Kova i University of Zagreb / Faculty of Electrical Engineering
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationEvaluating the knot vector to synthesize the cam motion using NURBS
Evaluating the knot vector to synthesize the cam motion using NURBS T.T.N. Nguyen 1, S. Kurtenbach 2, M. Hüsing 3 and B. Corves 4 Department of Mechanism Theory and Dynamics of Machines, RWTH Aachen University,
More informationNonparametric regression using kernel and spline methods
Nonparametric regression using kernel and spline methods Jean D. Opsomer F. Jay Breidt March 3, 016 1 The statistical model When applying nonparametric regression methods, the researcher is interested
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 Tobler s Law All places are related, but nearby places are related more than distant places Corollary:
More informationTYPE 529: RADIANT SLAB
TYPE 529: RADIANT SLAB General Description This component is intended to model a radiant floor-heating slab, embedded in soil, and containing a number of fluid filled pipes. The heat transfer within the
More informationThe Spalart Allmaras turbulence model
The Spalart Allmaras turbulence model The main equation The Spallart Allmaras turbulence model is a one equation model designed especially for aerospace applications; it solves a modelled transport equation
More information15.10 Curve Interpolation using Uniform Cubic B-Spline Curves. CS Dept, UK
1 An analysis of the problem: To get the curve constructed, how many knots are needed? Consider the following case: So, to interpolate (n +1) data points, one needs (n +7) knots,, for a uniform cubic B-spline
More informationSpatial Interpolation - Geostatistics 4/3/2018
Spatial Interpolation - Geostatistics 4/3/201 (Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Distance between pairs of points Lag Mean Tobler s Law All places are related, but nearby places
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More information1\C 1 I)J1mptMIll 'betl~flbejlli
1\C 1 )J1mptMll 'betl~flbejll l~~1hlbd ada]t6gh -or~«ejf. '~illrlf J~..6 '~~~ll!~ 4iJ~ "Mf:i',nl.Nqr2l' ~':l:mj}.i~:tv t.~l '\h Dr. N.Homsup, Abstract n this paper, two high-order FDTD schemes are developed
More informationAPPM/MATH Problem Set 4 Solutions
APPM/MATH 465 Problem Set 4 Solutions This assignment is due by 4pm on Wednesday, October 16th. You may either turn it in to me in class on Monday or in the box outside my office door (ECOT 35). Minimal
More informationAdvanced Numerical Methods for Numerical Weather Prediction
Advanced Numerical Methods for Numerical Weather Prediction Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943-5502 phone: (831) 656-4882 fax: (831) 656-4769 e-mail: giraldo@nrlmry.navy.mil
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We develop weighted essentially
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationResearch Article Data Visualization Using Rational Trigonometric Spline
Applied Mathematics Volume Article ID 97 pages http://dx.doi.org/.//97 Research Article Data Visualization Using Rational Trigonometric Spline Uzma Bashir and Jamaludin Md. Ali School of Mathematical Sciences
More informationConsider functions such that then satisfies these properties: So is represented by the cubic polynomials on on and on.
1 of 9 3/1/2006 2:28 PM ne previo Next: Trigonometric Interpolation Up: Spline Interpolation Previous: Piecewise Linear Case Cubic Splines A piece-wise technique which is very popular. Recall the philosophy
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More informationDeficient Quartic Spline Interpolation
International Journal of Computational Science and Mathematics. ISSN 0974-3189 Volume 3, Number 2 (2011), pp. 227-236 International Research Publication House http://www.irphouse.com Deficient Quartic
More informationNUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES
NUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES Ben R. Hodges 1 (A.M., ASCE) and Sarah Kelly Delavan 2 ABSTRACT Analysis of numerical diffusion and dissipation rates in a hydrostatic
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More information(Refer Slide Time: 00:02:24 min)
CAD / CAM Prof. Dr. P. V. Madhusudhan Rao Department of Mechanical Engineering Indian Institute of Technology, Delhi Lecture No. # 9 Parametric Surfaces II So these days, we are discussing the subject
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationSparse Grids. Lab 1. Discretization
Lab 1 Sparse Grids Lab Objective: Sparse Grids are an important tool when dealing with highdimensional problems. Computers operate in discrete space, not in continuous space. It is important to choose
More informationPositivity Preserving Interpolation of Positive Data by Rational Quadratic Trigonometric Spline
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 10, Issue 2 Ver. IV (Mar-Apr. 2014), PP 42-47 Positivity Preserving Interpolation of Positive Data by Rational Quadratic
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationFebruary 2017 (1/20) 2 Piecewise Polynomial Interpolation 2.2 (Natural) Cubic Splines. MA378/531 Numerical Analysis II ( NA2 )
f f f f f (/2).9.8.7.6.5.4.3.2. S Knots.7.6.5.4.3.2. 5 5.2.8.6.4.2 S Knots.2 5 5.9.8.7.6.5.4.3.2..9.8.7.6.5.4.3.2. S Knots 5 5 S Knots 5 5 5 5.35.3.25.2.5..5 5 5.6.5.4.3.2. 5 5 4 x 3 3.5 3 2.5 2.5.5 5
More informationC. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract
Advanced Computational Methods and Experiments in Heat Transfer XIII 15 A numerical study of heat diffusion using the Lagrangian particle SPH method and the Eulerian Finite-Volume method: analysis of convergence,
More informationSimulation of Time-Lag Permeation Experiments Using Finite Differences
Proceedings of the nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, pril 30 May 1 015 Paper No. 181 Simulation of Time-Lag Permeation Experiments Using Finite
More informationComputational Physics PHYS 420
Computational Physics PHYS 420 Dr Richard H. Cyburt Assistant Professor of Physics My office: 402c in the Science Building My phone: (304) 384-6006 My email: rcyburt@concord.edu My webpage: www.concord.edu/rcyburt
More informationRiver inundation modelling for risk analysis
River inundation modelling for risk analysis L. H. C. Chua, F. Merting & K. P. Holz Institute for Bauinformatik, Brandenburg Technical University, Germany Abstract This paper presents the results of an
More informationInvestigation of Shape Parameter for Exponential Weight Function in Moving Least Squares Method
2 (2017) 17-24 Progress in Energy and Environment Journal homepage: http://www.akademiabaru.com/progee.html ISSN: 2600-7662 Investigation of Shape Parameter for Exponential Weight Function in Moving Least
More informationConsiderations about level-set methods: Accuracy, interpolation and reinitialization
Considerations about level-set methods: Accuracy, interpolation and reinitialization Gustavo C. Buscaglia Coll: E. Dari, R. Ausas, L. Ruspini Instituto de Ciências Matemáticas e de Computação Universidade
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationApproximations for Linear Tenth-order Boundary Value Problems through Polynomial and Non-polynomial Cubic Spline Techniques
Proceedings of the Pakistan Academy of Sciences 52 4): 389 396 2015) Copyright Pakistan Academy of Sciences ISSN: 0377-2969 print), 2306-1448 online) Pakistan Academy of Sciences Research Article Approximations
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationEffectiveness of Element Free Galerkin Method over FEM
Effectiveness of Element Free Galerkin Method over FEM Remya C R 1, Suji P 2 1 M Tech Student, Dept. of Civil Engineering, Sri Vellappaly Natesan College of Engineering, Pallickal P O, Mavelikara, Kerala,
More informationDigital Image Processing. Prof. P.K. Biswas. Department of Electronics & Electrical Communication Engineering
Digital Image Processing Prof. P.K. Biswas Department of Electronics & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Image Segmentation - III Lecture - 31 Hello, welcome
More informationEnvironmental Modelling & Software
Environmental Modelling & Software 26 (2011) 1298e1308 Contents lists available at ScienceDirect Environmental Modelling & Software journal homepage: www.elsevier.com/locate/envsoft Implementation of quadratic
More information1 Trajectories. Class Notes, Trajectory Planning, COMS4733. Figure 1: Robot control system.
Class Notes, Trajectory Planning, COMS4733 Figure 1: Robot control system. 1 Trajectories Trajectories are characterized by a path which is a space curve of the end effector. We can parameterize this curve
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationComputer Project 3. AA Computational Fluid Dyanmics University of Washington. Mishaal Aleem March 17, 2015
Computer Project 3 AA 543 - Computational Fluid Dyanmics University of Washington Mishaal Aleem March 17, 2015 Contents Introduction........................................... 1 3.1 Grid Generator.......................................
More informationA NUMERICAL SOLUTION OF ONE DIMENSIONAL HEAT EQUATION USING CUBIC B-SPLINE BASIS FUNCTIONS
A NUMERICAL SOLUTION OF ONE DIMENSIONAL HEAT EQUATION USING CUBIC B-SPLINE BASIS FUNCTIONS 1 SHARANJEET DHAWAN, 2 SHEO KUMAR 1 Department of mathematics 2 Dr. B.R.Ambedkar National Institute of Technology,
More informationAdvances in Water Resources
Advances in Water Resources 53 (213) 198 27 Contents lists available at SciVerse ScienceDirect Advances in Water Resources journal homepage: www.elsevier.com/locate/advwatres An efficient time-domain approach
More informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationUsing the Discrete Ordinates Radiation Model
Tutorial 6. Using the Discrete Ordinates Radiation Model Introduction This tutorial illustrates the set up and solution of flow and thermal modelling of a headlamp. The discrete ordinates (DO) radiation
More informationIncompressible Viscous Flow Simulations Using the Petrov-Galerkin Finite Element Method
Copyright c 2007 ICCES ICCES, vol.4, no.1, pp.11-18, 2007 Incompressible Viscous Flow Simulations Using the Petrov-Galerkin Finite Element Method Kazuhiko Kakuda 1, Tomohiro Aiso 1 and Shinichiro Miura
More information