Linking Correlations Spanning Adjacent Applicability Domains
|
|
- Barrie Harvey
- 5 years ago
- Views:
Transcription
1 Ian David Lockhart Bogle and Michael Fairweather (Editors), Proceedings of the 22nd European Symposium on Computer Aided Process Engineering, June 2012, London Elsevier B.V. All rights reserved. Linking Correlations Spanning Adjacent Applicability Domains Tareg M. Alsoudani and I. D. L. Bogle Centre for Process Systems Engineering, University College London, Dept. of Chemical Engineering, Torrington Place, London, WC1E 7JE, U.K. Abstract Discontinuities between distinct regions, described by different equation sets, cause difficulties for solvers. We present a new algorithm that eliminates integrator discontinuities through two steps. First, it determines the optimum switch point between two empirical correlations spanning two or more domains. The optimum switch point is determined by searching for a jump point that minimizes a discontinuity between results obtained through calculating a parameter using two or more different methods at an intersection or proimity of their validity domains. Second, it links the two adjacent discontinuous domains with an interpolating polynomial. When applied, this approach should eliminate the need for conventional integrators to discretize and then link discontinuities through the conventional approach that generates interpolating polynomials based on state variables. The new approach tackles each discontinuity at the level of the constitutive equations rather than at the level of state variables, hence eliminating errors associated with polynomials generated at state variable level. Keywords: discontinuity, domains overlap. 1. Introduction and Literature Survey Literature is rich in empirical correlations, spanning different domains, to calculate certain parameters. Each of these correlations spans a certain validity domain with a specified accuracy. Theoretically, parameters calculated through different correlations should have eact or very close values at the intersection points of the correlations' applicability domains. However, this is not the case with most empirical formula; resulting in large discontinuities when solving spatio-temporal problems. A typical eample, illustrating the discussion, would be the calculation of heat transfer coefficient at the intersections between laminar and turbulent regimes. A process can be thought of as a comple system that is described by, mostly, continuous mathematical functions (algebraic or differential). The solution of these mathematical functions, usually through integration, brings an insight into the behaviour of the process under study. However, the continuity of these mathematical functions is sometimes broken by internal or eternal influences. A rapid phase change or flow reversal are eamples of an internally generated discontinuity in a ODE/DEA system whereas switching a pump on or off can be considered as an eternal influence that raises a discontinuity. Regardless of the form or source of discontinuity, it needs to be resolved either before starting to integrate the ODE/DAE system (if possible) or whenever it is encountered during the evolution of the integration process. Eternal influences affecting model continuity are usually resolved through the use of the modelling language, especially modern modelling languages. Once a DAE/ODE solver encounters an eternal
2 2 T. Alsoudani and I.D.L. Bogle discontinuity, it reinitializes state variables based on the new conditions post discontinuity. Eamples of such cases include changes in process boundary conditions, changes in process network through flow rerouting or pump switching, relief valve triggering, etc. The literature in handling internal discontinuities can be classified into two types. The first type (Type I) handles discontinuities generically through detection and resolution algorithms that are usually embedded within the solver. The solver usually detects the discontinuity through observing a jump in one or more of the state variables although the origin of the discontinuity might not necessarily be the state variable. The solver, then, tries to bridge between the two discontinuous states either through discretizing the discontinuous domain or through generating an interpolating polynomial in the discontinuous variable as illustrated in the figure below (Borst, 2008). For a comprehensive review of state event detection algorithms, the reader may refer to (Park and Barton, 1996), (Mao and Petzold, 2002) and (Archibald, et al, 2008). (Javey, 1988) reports three methods for resolving discontinuities once detected: changing equations to those post discontinuity conditions, halving the step in an effort to find a resolution at a finer grid, and reinitializing the ODE/DAE system to conditions post discontinuity. The second type (Type II) of discontinuity handling techniques deals with detection and resolution of discontinuities through structuring the ODE/DAE system in a way that minimizes, and hopefully eliminates, discontinuities. The literature is etremely sparse on this type of discontinuity detection and resolution. The first encountered attempt into resolving such discontinuity was by Brackbill et al. (1992). Borst (2008) emphasized that the use of regulating functions derived from the physics of the problem (Type II) will better eliminate discontinuities than the sole use of discretization techniques. 2. Methodology In this work, we try to resolve type II internal discontinuities arising from the use of logical epressions such as (1a) by devising discontinuity detection and resolution algorithms. Such epressions usually arise when simulating adjacent domains defined by different empirical correlations. The modeler usually starts with an epression of the form (1a.) and then transforms it into form (1b.) during model construction. Assuming a < b, c < d and b c in epression (1b.) will satisfy (the last condition enforces overlapping or touching domains), g c g b to maintain continuity in dimension. We are then free to select any g that satisfies g [ b, c]. However, an arbitrary choice will f ( ) f 1 = f 2 ( ), [ a,b] ( ), [ c,d] 1a. 1b. If (<g) then (Domain I) f() = () Elseif ( g) then (Domain II) f() = ()
3 Linking Correlations Spanning Adjacent Applicability Domains 3 probably result in selecting g points that increase or create discontinuities as illustrated in (2a) and (2b), respectively. Nevertheless, this freedom in selecting g facilitates transforming the search for the best g into an optimization problem of the form: min. () () s.t. [ b, c] The solution of this optimization problem will result in g that minimizes the discontinuity between non-intersecting and (2a); or eliminates the discontinuity between intersecting and (2b). a c g ' g b d 2a. 2b. For multiple functions f i (i=1 n) spanning multiple adjacent domains [a i, b i ] and satisfying a i b i-1 (i=2 n), there eist n-1 intersecting g i points that need to be located by applying the minimization algorithm at each intersection interval. Once g i points are identified results can be reflected directly in the logical epression (1b). We can then rely on type I discontinuity handlers to resolve any discontinuity in a state variable resulting from (2a) type of discontinuity in the state variable s respective constitutive equation. However, a better solution would be to eliminate the discontinuity at the constitutive equation level. This is achievable through linking and by an interpolating polynomial at the location of the discontinuity. Any interpolating polynomial would satisfy zero order smoothness. However, 2 nd order smoothness (usually required in integration routines) is only achieved through a few of the readily available interpolating functions such as cubic splines and cubic hermite interpolating polynomials. In subsequent discussion, we will demonstrate the concept using hermite interpolating polynomials although any 2 nd order interpolating polynomial would fit. Three points are usually enough to construct an interpolating polynomial. However, to ensure better closure of the interpolating polynomial with functions and, a fourpoint interpolating polynomial is selected as illustrated in (3a.) and (3b.). The points are a c g g ' b d f ( ) Rela tive to { p=0.0 p =0.1 p =0.2 p =0.3 p =0.4 } p =0.4 p =0.5 p =0.3 p =0.2 p =0.1 p =0.0 f ( ) p=0.0 p=0.1 p=0.2 p=0.3 p=0.4 p=0.5 Rela tive to 3a. Three point interpolation 3b. Four point interpolation
4 4 T. Alsoudani and I.D.L. Bogle If (<g-1.5h) then (Domain I) f() = () elseif (g-1.5h g-1.5h) f() = Interpolating polynomial Elseif (>g) then (Domain II) f() = () 4. equally separated by a distance h and equally distributed between the right and left sides of the g point location. Two of the points lie eactly on the respective functions and and located at distances g-1.5g and g+1.5h, respectively. The other two points are located at distances g-0.5h and g+0.5h. To ensure proper closure of the interpolating polynomial with and, interpolating points values of the latter locations will correspond to the values of their respective functions after adding p (g)- (g) to the lower valued function and subtracting p (g)- (g) from the higher valued one, where p ranges from 0 to 0.5. The p factor is a tuning parameter that can be used to increase smoothness at the link points between the interpolating polynomial and the respective functions as illustrated in (3a.) and (3b.). The logical epression in (1b) then converts to (4) and the algorithm terminates. The 1D detection algorithm can be epanded to cover 2D functions f(,y) through treating each dimension separately. However, the 2D discontinuity resolution algorithm differs because the independent variables and y move freely in a 2D space based on evolution of simulation run. This simulation-dependent movement makes it hard to determine where the discontinuity would be hit during a simulation run and hence it is hard to fi the coordinates of the interpolating polynomial. Thus, we recommend tracking the movement trajectory of the independent variables and y through a vector v i at each simulation integration step and etrapolating the trajectory to the discontinuity plane in order to continuously update the coordinates of the interpolating polynomial points as in (5a). P o v 1 P1 v 2 P 2 v n f (, y) P n B v p A P n 1 P n 2 C y 5a. 5b. 3. An Eample To test the theory, we used a simplified model of the paraffins isomerization reactor that was part of the patent published by Minkkinen (1993). The model is constructed using the gproms modeling language (PSE,2011). We tested the effect of transition from Laminar to Turbulent flow regimes on the wall heat transfer coefficient. For Laminar
5 Linking Correlations Spanning Adjacent Applicability Domains 5 flow, we used the simplified constant heat-flu equation of Nu d = For turbulent flow we used the Gnielinski correlation (Keith, 2000): / where: ^ 2 A plot of the Nu d versus Re and Pr is illustrated in (5b.). Since minimum Re for Gnielinski correlation is 2300, this value was taken as the separation value between Laminar and Turbulent flow regimes. We started the simulation at a feed velocity corresponding to a Reynolds number o030. Then, during the course of the simulation, we ramped the velocity to cross the 2300 laminar boundary. Unlike other simulation packages that detect discontinuities only at the state variable level, gproms is smart enough to detect discontinuities at the constitutive equation level. However, gproms still performs reinitialization once a discontinuity is detected. This reinitialization is a form of discretization that we discussed earlier (Type I). To apply the two dimensional solution we proposed, we needed to link gproms to a C++ code through gproms Foreign Object Interface. gproms passes the required parameters to compute Nu d and the C++ code determines whether it needs to calculate Nu d using Laminar, Turbulent, or the interpolating polynomial equation. Once the simulation is run with the C++ code, gproms previously reported discontinuity disappeared and hence gproms reinitialization was not performed. The evolution of the simulation vector is superimposed on the Nu d mesh. Simulation results also demonstrate a noticeable decline in overall simulation computations due to elimination of model reinitialization; although it is not conclusive. To conclude, we were able to demonstrate an algorithm that detects and resolves discontinuities resulting from the use of correlations describing behaviors of adjacent domains. The solution to the 1D detection and resolution problem is unique for a fied set of correlations irrespective of the model. Thus the algorithm can be eecuted prior to starting the simulation run. For the 2D case, the solution is dependent on the evolution of the independent variables trajectory. We think the 2D solution can be further etended to cover higher dimension correlations. Also this approach is best embedded inside the simulator since it allows a simulation compiler to automatically generate the code containing the interpolating polynomial from that written by the modeler. References 1. Archibald, R., A. Gelb and J. Yoon, "Determining the Locations and Discontinuities in the Derivatives of Functions", Applied Numerical Mathematics, vol. 58,, pp , Borst, R., "Challenges in computational materials science: Multiple scales, multi-physics and evolving discontinuities", Computational Materials Science, vol. 43,, pp. 1-15, Brackbill, J.U., D.B. Kothe and C. Zemach, "A Continuum Method for Modeling Surface Tension", Journal of Computational Physics, vol. 100,, pp , gproms Modelling Language, Copyright Process Systems Enterprise Ltd. 5. Kreith, Frank, CRC Handbook of Thermal Engineering, CRC Press, Mao, G. and L.R. Petzold, "Efficient Integration over Discontinuities for Differential- Algebraic Systems", Computers and Mathematics with Applics, vol. 43,, pp , Minkkinen, A., L. Mank, and S. Jullian, "Process for the Isomerization of of C5/C6 Normal Paraffins with Recycling of Normal Paraffins", U. S. Patent 5,233,120 (1993). 8. Park, T. and P.I. Barton, "State Event Location in Differential-Algebraic Models", ACM Transactions on Modeling & Computer Simulation, vol. 6, issue 2, pp , 1996.
Numerical Analysis Fall. Numerical Differentiation
Numerical Analysis 5 Fall Numerical Differentiation Differentiation The mathematical definition of a derivative begins with a difference approimation: and as is allowed to approach zero, the difference
More informationBoundary/Contour Fitted Grid Generation for Effective Visualizations in a Digital Library of Mathematical Functions
Boundary/Contour Fitted Grid Generation for Effective Visualizations in a Digital Library of Mathematical Functions Bonita Saunders Qiming Wang National Institute of Standards and Technology Bureau Drive
More informationA novel method for identification of critical points in flow sheet synthesis under uncertainty
Ian David Lockhart Bogle and Michael Fairweather (Editors), Proceedings of the nd European Symposium on Computer Aided Process Engineering, 17-0 June 01, London. 01 Elsevier B.V. All rights reserved A
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 informationKochanek-Bartels Cubic Splines (TCB Splines)
Kochanek-Bartels Cubic Splines (TCB Splines) David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More information3.5 Indeterminate Forms. Consider the function given by the rule f(x) = Whenever direct substitution into lim x c f(x) yields a 0 0
3.5 Indeterminate Forms Introduction; Consider the function given by the rule f() = ; its graph is shown below. 0 0 situations Clearly, = 3. Note, however, that if one merely tried to plug in 0 for when
More informationModeling with Uncertainty Interval Computations Using Fuzzy Sets
Modeling with Uncertainty Interval Computations Using Fuzzy Sets J. Honda, R. Tankelevich Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, CO, U.S.A. Abstract A new method
More informationComputational Fluid Dynamics (CFD) Simulation in Air Duct Channels Using STAR CCM+
Available onlinewww.ejaet.com European Journal of Advances in Engineering and Technology, 2017,4 (3): 216-220 Research Article ISSN: 2394-658X Computational Fluid Dynamics (CFD) Simulation in Air Duct
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 informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More 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 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 informationSUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS
SUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS M. Hofer and H. Pottmann Institute of Geometry Vienna University of Technology, Vienna, Austria hofer@geometrie.tuwien.ac.at, pottmann@geometrie.tuwien.ac.at
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationSurface Approximation and Interpolation via Matrix SVD
Surface Approximation and Interpolation via Matrix SVD Andrew E. Long and Clifford A. Long Andy Long (longa@nku.edu) is at Northern Kentucky University. He received his Ph.D. in applied mathematics at
More informationFinding a Best Fit Plane to Non-coplanar Point-cloud Data Using Non Linear and Linear Equations
AIJSTPME (013) 6(): 17-3 Finding a Best Fit Plane to Non-coplanar Point-cloud Data Using Non Linear and Linear Equations Mulay A. Production Engineering Department, College of Engineering, Pune, India
More informationA nodal based evolutionary structural optimisation algorithm
Computer Aided Optimum Design in Engineering IX 55 A dal based evolutionary structural optimisation algorithm Y.-M. Chen 1, A. J. Keane 2 & C. Hsiao 1 1 ational Space Program Office (SPO), Taiwan 2 Computational
More informationKeyword: Quadratic Bézier Curve, Bisection Algorithm, Biarc, Biarc Method, Hausdorff Distances, Tolerance Band.
Department of Computer Science Approximation Methods for Quadratic Bézier Curve, by Circular Arcs within a Tolerance Band Seminar aus Informatik Univ.-Prof. Dr. Wolfgang Pree Seyed Amir Hossein Siahposhha
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 informationPlotting technologies for the study of functions of two real variables
Plotting technologies for the study of functions of two real variables David Zeitoun 1 and Thierry Dana-Picard 2 1 Department of Mathematics, Orot College of Education, Rehovot, Israel, ed.technologie@gmail.com
More informationCFD Analysis on Heat Transfer Through Different Extended Surfaces
CFD Analysis on Heat Transfer Through Different Extended Surfaces Ravindra Kondaguli 1 1 Department of Mechanical Engineering BLDECET Vijayapur Abstract: The present work includes CFD analysis and comparison
More informationChapter 3 Path Optimization
Chapter 3 Path Optimization Background information on optimization is discussed in this chapter, along with the inequality constraints that are used for the problem. Additionally, the MATLAB program for
More informationFUZZY INFERENCE SYSTEMS
CHAPTER-IV FUZZY INFERENCE SYSTEMS Fuzzy inference is the process of formulating the mapping from a given input to an output using fuzzy logic. The mapping then provides a basis from which decisions can
More information3D object comparison with geometric guides for Interactive Evolutionary CAD
Loughborough University Institutional Repository 3D object comparison with geometric guides for Interactive Evolutionary CAD This item was submitted to Loughborough University's Institutional Repository
More informationA SCATTERED DATA INTERPOLANT FOR THE SOLUTION OF THREE-DIMENSIONAL PDES
European Conference on Computational Fluid Dnamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périau (Eds) c TU Delft, The Netherlands, 26 A SCATTERED DATA INTERPOLANT FOR THE SOLUTION OF THREE-DIMENSIONAL
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationThe method of rationalizing
Roberto s Notes on Differential Calculus Chapter : Resolving indeterminate forms Section The method of rationalizing What you need to know already: The concept of it and the factor-and-cancel method of
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationGrid and Mesh Generation. Introduction to its Concepts and Methods
Grid and Mesh Generation Introduction to its Concepts and Methods Elements in a CFD software sstem Introduction What is a grid? The arrangement of the discrete points throughout the flow field is simpl
More informationThe method of rationalizing
Roberto s Notes on Differential Calculus Chapter : Resolving indeterminate forms Section The method of rationalizing What you need to know already: The concept of it and the factor-and-cancel method of
More informationGlobal Solution of Mixed-Integer Dynamic Optimization Problems
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 25 Elsevier Science B.V. All rights reserved. Global Solution of Mixed-Integer Dynamic Optimization
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationAnalyzing and Segmenting Finger Gestures in Meaningful Phases
2014 11th International Conference on Computer Graphics, Imaging and Visualization Analyzing and Segmenting Finger Gestures in Meaningful Phases Christos Mousas Paul Newbury Dept. of Informatics University
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationProgressive Surface Modeling Based On 3D Motion Sketch
Progressive Surface Modeling Based On 3D Motion Sketch SHENGFENG IN, and DAVID K WRIGHT School of Engineering and Design Brunel University Uxbridge, Middlesex UB8 3PH UK Abstract: - This paper presents
More informationPreparing for AS Level Further Mathematics
Preparing for AS Level Further Mathematics Algebraic skills are incredibly important in the study of further mathematics at AS and A level. You should therefore make sure you are confident with all of
More informationLecture 1 GENERAL INTRODUCTION: HISTORICAL BACKGROUND AND SPECTRUM OF APPLICATIONS
Lecture 1 GENERAL INTRODUCTION: HISTORICAL BACKGROUND AND SPECTRUM OF APPLICATIONS 1.1 INTRODUCTION Analysis of physical problems in any area of engineering and science involves a multipronged approach:
More informationPROTECTION AGAINST MODELING AND SIMULATION UNCERTAINTIES IN DESIGN OPTIMIZATION NSF GRANT DMI
PROTECTION AGAINST MODELING AND SIMULATION UNCERTAINTIES IN DESIGN OPTIMIZATION NSF GRANT DMI-9979711 Bernard Grossman, William H. Mason, Layne T. Watson, Serhat Hosder, and Hongman Kim Virginia Polytechnic
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 informationC1M0 Introduction to Maple Assignment Format C1M1 C1M1 Midn John Doe Section 1234 Beginning Maple Syntax any
CM0 Introduction to Maple Our discussion will focus on Maple 6, which was developed by Waterloo Maple Inc. in Waterloo, Ontario, Canada. Quoting from the Maple 6 Learning Guide, Maple is a Symbolic Computation
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationSection 1.1: Four Ways to Represent a Function
Section.: Four Ways to Represent a Function. The Definition of a Function Functions are one of the most basic tools in mathematics, so we start by considering the definition of a function and all related
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationTowards An Integrated Classification of Exceptions
Towards An Integrated Classification of Annelise Janse van Rensburg a Karen Renaud b University of South Africa a aisg@lantic.net b renaukv@unisa.ac.za Abstract : refer to the situations that are not modelled
More informationNumerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers
Numerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers R. Jendrny, H. Damanik, O. Mierka, S. Turek Institute of Applied Mathematics (LS III), TU
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationMODELLING-Choosing a model
MODELLING-Choosing a model Model categories When using a model to help in the design process, it is important that the right type of model is used. Using the wrong type of model can waste computing power
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 informationDIMENSIONAL SYNTHESIS OF SPATIAL RR ROBOTS
DIMENSIONAL SYNTHESIS OF SPATIAL RR ROBOTS ALBA PEREZ Robotics and Automation Laboratory University of California, Irvine Irvine, CA 9697 email: maperez@uci.edu AND J. MICHAEL MCCARTHY Department of Mechanical
More informationSolve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method
Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method P.B. Choksi 1 and A.K. Pathak 2 1 Research Scholar, Rai University,Ahemdabad. Email:pinalchoksey@gmail.com 2 H.O.D.
More informationThe Eikonal Equation
The Eikonal Equation Numerical efficiency versus computational compleity Shu-Ren Hysing III Institute of Applied Mathematics LSIII University of Dortmund Level set - methodology By embedding an interface
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 informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationNAG Library Chapter Introduction. D04 Numerical Differentiation
D4 Numerical Differentiation Introduction D4 NAG Library Chapter Introduction D4 Numerical Differentiation Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Description of the
More information3D Modeling Parametric Curves & Surfaces
3D Modeling Parametric Curves & Surfaces Shandong University Spring 2012 3D Object Representations Raw data Point cloud Range image Polygon soup Solids Voxels BSP tree CSG Sweep Surfaces Mesh Subdivision
More informationSolution Recording and Playback: Vortex Shedding
STAR-CCM+ User Guide 6663 Solution Recording and Playback: Vortex Shedding This tutorial demonstrates how to use the solution recording and playback module for capturing the results of transient phenomena.
More informationConnected Minimal Acceleration Trigonometric Curves
Connected Minimal Acceleration Trigonometric Curves Tony Barrera Barrera Kristiansen AB Anders Hast Creative Media Lab, University of Gävle Ewert Bengtsson Centre for Image Analysis Uppsala University
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationLimits and Derivatives (Review of Math 249 or 251)
Chapter 3 Limits and Derivatives (Review of Math 249 or 251) 3.1 Overview This is the first of two chapters reviewing material from calculus; its and derivatives are discussed in this chapter, and integrals
More informationLecture 3.2 Methods for Structured Mesh Generation
Lecture 3.2 Methods for Structured Mesh Generation 1 There are several methods to develop the structured meshes: Algebraic methods, Interpolation methods, and methods based on solving partial differential
More informationIntroduction to and calibration of a conceptual LUTI model based on neural networks
Urban Transport 591 Introduction to and calibration of a conceptual LUTI model based on neural networks F. Tillema & M. F. A. M. van Maarseveen Centre for transport studies, Civil Engineering, University
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 informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationcomputational Fluid Dynamics - Prof. V. Esfahanian
Three boards categories: Experimental Theoretical Computational Crucial to know all three: Each has their advantages and disadvantages. Require validation and verification. School of Mechanical Engineering
More informationDECONFLICTION AND SURFACE GENERATION FROM BATHYMETRY DATA USING LR B- SPLINES
DECONFLICTION AND SURFACE GENERATION FROM BATHYMETRY DATA USING LR B- SPLINES IQMULUS WORKSHOP BERGEN, SEPTEMBER 21, 2016 Vibeke Skytt, SINTEF Jennifer Herbert, HR Wallingford The research leading to these
More informationGetting Started with montaj GridKnit
Getting Started with montaj GridKnit Today s explorers use regional and local geophysical compilations of magnetic, radiometric and apparent resistivity data to perform comprehensive geologic interpretations
More informationData Mining Approaches to Characterize Batch Process Operations
Data Mining Approaches to Characterize Batch Process Operations Rodolfo V. Tona V., Antonio Espuña and Luis Puigjaner * Universitat Politècnica de Catalunya, Chemical Engineering Department. Diagonal 647,
More informationAshwin Shridhar et al. Int. Journal of Engineering Research and Applications ISSN : , Vol. 5, Issue 6, ( Part - 5) June 2015, pp.
RESEARCH ARTICLE OPEN ACCESS Conjugate Heat transfer Analysis of helical fins with airfoil crosssection and its comparison with existing circular fin design for air cooled engines employing constant rectangular
More informationLarge Data Analysis via Interpolation of Functions: Interpolating Polynomials vs Artificial Neural Networks
American Journal of Intelligent Systems 2018, 8(1): 6-11 DOI: 10.5923/j.ajis.20180801.02 Large Data Analysis via Interpolation of Functions: Interpolating Polynomials vs Artificial Neural Networks Rohit
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 informationMaximizing an interpolating quadratic
Week 11: Monday, Apr 9 Maximizing an interpolating quadratic Suppose that a function f is evaluated on a reasonably fine, uniform mesh {x i } n i=0 with spacing h = x i+1 x i. How can we find any local
More informationSOLIDWORKS Flow Simulation 2015
An Introduction to SOLIDWORKS Flow Simulation 2015 John E. Matsson, Ph.D. SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following websites
More informationComputational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy. Jordan A.
Computational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy Jordan A. Denev Abstract: The present paper is a first one from a series of papers
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 informationComputational Simulation of the Wind-force on Metal Meshes
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Computational Simulation of the Wind-force on Metal Meshes Ahmad Sharifian & David R. Buttsworth Faculty
More informationCHAPTER 6 Parametric Spline Curves
CHAPTER 6 Parametric Spline Curves When we introduced splines in Chapter 1 we focused on spline curves, or more precisely, vector valued spline functions. In Chapters 2 and 4 we then established the basic
More informationRobots are built to accomplish complex and difficult tasks that require highly non-linear motions.
Path and Trajectory specification Robots are built to accomplish complex and difficult tasks that require highly non-linear motions. Specifying the desired motion to achieve a specified goal is often a
More informationFinite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2
Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,
More informationRemark. Jacobs University Visualization and Computer Graphics Lab : ESM4A - Numerical Methods 331
Remark Reconsidering the motivating example, we observe that the derivatives are typically not given by the problem specification. However, they can be estimated in a pre-processing step. A good estimate
More informationDynamic Memory Usage Optimization using ILP
Dynamic emory Usage Optimization using IP A. Allam and J. Ramanuam 2 Electrical Engineering Dept., Assiut University, Egypt 2 Electrical and Computer Engineering Dept, ouisiana State University, USA atef
More informationUntil now we have worked with flat entities such as lines and flat polygons. Fit well with graphics hardware Mathematically simple
Curves and surfaces Escaping Flatland Until now we have worked with flat entities such as lines and flat polygons Fit well with graphics hardware Mathematically simple But the world is not composed of
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
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 informationJoe Wingbermuehle, (A paper written under the guidance of Prof. Raj Jain)
1 of 11 5/4/2011 4:49 PM Joe Wingbermuehle, wingbej@wustl.edu (A paper written under the guidance of Prof. Raj Jain) Download The Auto-Pipe system allows one to evaluate various resource mappings and topologies
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More informationOperators to calculate the derivative of digital signals
9 th IMEKO TC 4 Symposium and 7 th IWADC Workshop July 8-9, 3, Barcelona, Spain Operators to calculate the derivative of digital signals Lluís Ferrer-Arnau, Juan Mon-Gonzalez, Vicenç Parisi-Baradad Departament
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 informationImportant Properties of B-spline Basis Functions
Important Properties of B-spline Basis Functions P2.1 N i,p (u) = 0 if u is outside the interval [u i, u i+p+1 ) (local support property). For example, note that N 1,3 is a combination of N 1,0, N 2,0,
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More 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 informationModule 7 VIDEO CODING AND MOTION ESTIMATION
Module 7 VIDEO CODING AND MOTION ESTIMATION Lesson 20 Basic Building Blocks & Temporal Redundancy Instructional Objectives At the end of this lesson, the students should be able to: 1. Name at least five
More informationResearch Article A Family of Even-Point Ternary Approximating Schemes
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 97, pages doi:.5//97 Research Article A Family of Even-Point Ternary Approximating Schemes Abdul Ghaffar and Ghulam
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
More informationExact Optimal Solution of Fuzzy Critical Path Problems
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Vol. 6, Issue (June 0) pp. 5 67 (Previously, Vol. 6, Issue, pp. 99 008) Applications and Applied Mathematics: An International Journal
More informationEfficient computation of source magnetic scalar potential
Adv. Radio Sci., 4, 59 63, 2006 Author(s) 2006. This work is licensed under a Creative Commons License. Advances in Radio Science Efficient computation of source magnetic scalar potential W. Hafla, A.
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 informationArxiv. Golden interpolation. Ying He a, Jincai Chang b, A B S T R A C T
Ariv (2018) Contents lists available at Ariv Ariv ariv homepage: https://arivorg/ Golden interpolation Interpolation is a classic problem In the interpolation problem, it is required to construct a function
More information