NGsolve::Give me your element
|
|
- Garry Dorsey
- 6 years ago
- Views:
Transcription
1 NGsolve::Give me your element And let your eyes delight in my ways Jay Gopalakrishnan Portland State University Winter 2015 Download code (works only on version 6.0) for these notes from here. Give me your heart / And let your eyes delight in my ways. The Bible Jay Gopalakrishnan 1/18
2 Contents 1 Automatic differentiation 2 Shape functions 3 Finite elements and spaces 4 Orientation 5 Bicubic qualdrilateral element Jay Gopalakrishnan 2/18
3 Automatic differentiation Goal: Create variables that know how to (exactly) differentiate themselves. Idea: Differentiation obeys some rules (product rule, quotient rule etc) that we can implement by overloading operators, e.g., overload * to implement i (f g) = f ( i g) + g( i f ). Suppose an object representing x i (the ith coordinate) knows its value and the value of its derivatives, at any given point. Then, we can compute both the value and the value of derivatives of x i x j by overloading * as above. Jay Gopalakrishnan 3/18
4 A minimalist class for differentiation template<i n t D> c l a s s MyDiffVar {// My d i f f e r e n t i a b l e v a r i a b l e s // F i l e : d i f f e r e n t i a b l e s. hpp double Value ; double D e r i v a t i v e s [D ] ; p u b l i c : MyDiffVar ( ) { } ; MyDiffVar ( double x i, i n t i ) { // i th c o o r d i n a t e x i has Value = x i ; // grad = i th u n i t v e c t o r f o r ( auto & d : D e r i v a t i v e s ) d = 0. 0 ; D e r i v a t i v e s [ i ] = 1. 0 ; } double GetValue ( ) const { r e t u r n Value ; } double& SetValue ( ) { r e t u r n Value ; } double G e t D e r i v a t i v e ( i n t i ) const { r e t u r n D e r i v a t i v e s [ i ] ; } double& S e t D e r i v a t i v e ( i n t i ) { r e t u r n D e r i v a t i v e s [ i ] ; } } ; Jay Gopalakrishnan 4/18
5 Overload * for the differentiables Template implementation of implement i (f g) = f ( i g) + g( i f ): // implement p r o d u c t r u l e template<i n t D> MyDiffVar<D> operator ( const MyDiffVar<D> & f, const MyDiffVar<D> & g ) { } MyDiffVar<D> f g ; f g. S etvalue ( ) = f. GetValue ( ) g. GetValue ( ) ; f o r ( i n t i =0; i <D; i ++) f g. S e t D e r i v a t i v e ( i ) = f. GetValue ( ) g. G e t D e r i v a t i v e ( i ) + g. GetValue ( ) f. G e t D e r i v a t i v e ( i ) ; r e t u r n f g ; Quiz: Open the file and provide operators +,, and /. Jay Gopalakrishnan 5/18
6 Using your class #i n c l u d e d i f f e r e n t i a b l e s. hpp // F i l e d0. cpp using namespace s t d ; i n t main ( ) { MyDiffVar <2> x ( 0. 5, 0 ), y ( 2. 0, 1 ) ; } cout << x : << x << e n d l << y : << y << e n d l << x y : << x y << e n d l << x y y+y : << x y y+y << e n d l ; Using your simple class, you can now differentiate polynomial expressions built using x and y coordinates (or x i, i = 1,..., N, in N-dimensions). Jay Gopalakrishnan 6/18
7 Exercise! How would you modify differentiables.hpp so that you can also differentiate expressions like sin(xy)? Make sure your modified file compiles and runs correctly with this driver: #i n c l u d e d i f f e r e n t i a b l e s. hpp // F i l e d0x. cpp using namespace s t d ; i n t main ( ) { MyDiffVar <2> x ( 0. 5, 0 ), y ( 2. 0, 1 ) ; } cout << x : << x << e n d l << y : << y << e n d l << s i n ( x y )/ y : << s i n ( x y )/ y << e n d l ; Jay Gopalakrishnan 7/18
8 Netgen s AutoDiff class An implementation of these ideas is available in $NGSRC/netgen/libsrc/general/autodiff.hpp. Here is an example showing how to use it: #i n c l u d e <fem. hpp> // F i l e d1. cpp using namespace s t d ; i n t main ( ) { } AutoDiff <2> x ( 0. 5, 0 ), y ( 2. 0, 1 ) ; // x and y c o o r d s AutoDiff <2> b [ 3 ] = { x, y, 1 x y } ; // b a r y c e n t r i c c o o r d s cout << x : << x << e n d l << y : << y << e n d l << x y : << x y << e n d l << x y y+y : << x y y+y << e n d l << ( b0 b1 b2 1)/y : << ( b [ 0 ] b [ 1 ] b [ 2 ] 1)/ y << e n d l ; We will use AutoDiff and the following classes to program finite elements. Jay Gopalakrishnan 8/18
9 FlatVector, SliceVector, etc. #i n c l u d e <b l a. hpp> // F i l e : f l a t v e c. cpp using namespace s t d ; using namespace n g b l a ; i n t main ( ) { double mem [ ] = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 1 0 } ; F l a t V e c t o r <double> f 1 ( 2,mem) ; // A v e c t o r c l a s s t h a t s t e a l s F l a t V e c t o r <double> f 2 ( 2,mem+3); // memory from e l s e w h e r e. cout << f 1 : \ n << f 1 << e n d l ; // This p r i n t s 1, 2. cout << f 2 : \ n << f 2 << e n d l ; // This p r i n t s 5, 6. S l i c e V e c t o r <> s1 ( 4, 2,mem) ; // A l s o s t e a l s memory. cout << s1 : \ n << s1 << e n d l ; // This p r i n t s 1, 3, 5, 7. S l i c e V e c t o r <> s2 ( 5, 1,mem+4); // What i s t h i s? // : SliceVector class does not allocate or delete memory. Their constructors just create/copy pointers. Jay Gopalakrishnan 9/18
10 class ScalarFiniteElement template <i n t D> c l a s s S c a l a r F i n i t e E l e m e n t : p u b l i c F i n i t e E l e m e n t { v i r t u a l void CalcShape ( const I n t e g r a t i o n P o i n t & ip, S l i c e V e c t o r <> shape ) const = 0 ; v i r t u a l void CalcDShape ( const I n t e g r a t i o n P o i n t & ip, S l i c e M a t r i x <> dshape ) const = 0 ; //... } ; shape and dshape are cheap to pass by value as function arguments even when they contain many elements. Any derived finite element class must provide shape functions and their derivatives. Jay Gopalakrishnan 10/18
11 Visualizing finite element shape functions # FILE: shapes.pde geometry = square.in2d mesh = squaretrg.vol fespace v -type=h1ho -order=2 gridfunction u -fespace=v numproc shapetester nptest -gridfunction=u The numproc shapetester is an NGsolve tool to visualize global basis functions (called global shape functions) of an FESpace. Load this PDE file. Click Solve button before doing anything else. Look for a tiny window called Shape Tester that pops up. The number (0,1,... ) that you input in Shape Tester window determines which basis function will be set in gridfunction u. Got to Visual menu and pick gridfunction u to visualize. Jay Gopalakrishnan 11/18
12 Visualizing finite element shape functions # FILE: shapes.pde geometry = square.in2d mesh = squaretrg.vol fespace v -type=h1ho -order=2 gridfunction u -fespace=v numproc shapetester nptest -gridfunction=u 1, 2, etc. Global shape functions Jay Gopalakrishnan 11/18
13 Prepare to write your own finite element Study these files in the folder my_little_ngsolve: myelement.hpp, myelement.cpp, myhoelement.hpp, myhoelement.cpp All elements in a mesh are mapped from a fixed reference element. Pay particular attention to CalcShape(..) and CalcDshape(..). They give the values and derivatives of all local shape functions on the reference element. myfespace.hpp, myfespace.cpp, myhofespace.hpp, myhofespace.cpp Each global degree of freedom ( dof ) gives a global basis function and is associated to a geometrical object of the mesh (like a vertex, edge, or element). Pay particular attention to GetDofNrs(...), which return global dof-numbers connected to an element. Jay Gopalakrishnan 12/18
14 Homework Your assignment is to code the bicubic finite element Q 3 in bicubicelem.cpp. On the reference element, the unit square, this element consists of the space of functions Q 3 = span{x i y j : 0 i 3, 0 j 3}. Also code a bicubic finite element space (derived from FESpace), for any mesh of quadrilateral elements, in file bicubicspace.cpp Then, use your space to approximate the operator + I and solve a Neumann boundary value problem. Tabulate errors. The ensuing slides give you hints to complete this homework and suggest separating the work into smaller separate tasks. Jay Gopalakrishnan 13/18
15 Bicubic shape functions on unit square Task 1: In bicubicelem.cpp, provide shape functions. E.g., here is a valid basis set of shape functions (you may use others) : Vertex basis Edge basis φ 0 = (1 x)(1 y) φ 1 = (1 x)y φ 2 = x(1 y) φ 3 = xy φ 4 = (φ 0φ 1)φ 0 φ 5 = (φ 0φ 1)φ 1 φ 6 = (φ 1φ 3)φ 3 φ 7 = (φ 1φ 3)φ 1 : Interior basis φ 4 = (1 x)(1 y)xyφ 0 φ 4 = (1 x)(1 y)xyφ 1 :... Jay Gopalakrishnan 14/18
16 Bicubic shape functions on unit square Task 1: In bicubicelem.cpp, provide shape functions. Your basis expressions should go into the CalcShape member function. For the CalcDShape member function, you can use AutoDiff variables and the same expressions you need in CalcShape. Consider simplifying your code so that you only type the basis expressions once. Jay Gopalakrishnan 14/18
17 Orientation Task 2: In bicubicspace.cpp, write your finite element space. Remember to keep track of matching local and global orientation (go back and revise your bicubicelem.cpp if necessary). / What i s the l o c a l o r i e n t a t i o n? I s the o r d e r i n g o f v e r t i c e s and edges w i t h i n the r e f e r e n c e element v2 e1 v3 v3 e1 v2 o o o o e2 e3 e2 e3 o o o o or something v0 e0 v1, v0 e0 v1, e l s e? / What is the global orientation? NGsolve s mesh edges are directed/oriented. If edge shape functions from adjacent elements are not given in that orientation, then you may lose continuity! Jay Gopalakrishnan 15/18
18 Check your basis Task 3: Compile the code you wrote and make a shared library make libmyquad.so and check your basis functions on the three given quadrilateral mesh files. # FILE : bicubicshapes.pde geometry = square.in2d #mesh = squarequad1.vol.gz #mesh = squarequad2.vol.gz mesh = squarequad3.vol.gz shared = libmyquad define fespace v -type=myquadspace define gridfunction u -fespace=v numproc shapetester nptest -gridfunction=u Jay Gopalakrishnan 16/18
19 Solve a PDE Task 4: Using your finite element space, solve this boundary value problem: u + u = f on Ω u/ n = 0 on Ω. Hints: Do you know the variational formulation for this problem? You want to write a PDE file that mixes your finite element space with the NGSolve integrators. E.g., the NGSolve integrator laplace, can work with any finite element which provides CalcDShape, by dynamic polymorphism. Jay Gopalakrishnan 17/18
20 Solve a PDE Task 4: Using your finite element space, solve this boundary value problem: u + u = f on Ω This task includes these steps: u/ n = 0 on Ω. 1 Set f so that your exact solution is u = sin(πx) 2 sin(πy) 2. 2 Compute the L 2 (Ω) error (code this either in your own C++ numproc like we did before or find facilities to directly do it in the pde file). 3 Solve on mesh = squarequad3.vol.gz by loading your pde file and pressing the Solve button. Compute the L 2 (Ω)-error. Note it down. 4 Pressing the Solve button again to solve and compute the L 2 -error on a uniformly refined mesh. Note the L 2 -error. Repeat (until you can t). 5 What is the rate of convergence of L 2 -error with meshsize? Jay Gopalakrishnan 17/18
21 Project Student Team Project: Learn about the DPG method and download an implementation of it in GitHub. Your job is to extend it to quadrilateral elements. You will need to code a new finite element space that will serve as the test space for the DPG method. Details will be progressively made clear as you proceed with the project. Jay Gopalakrishnan 18/18
NGsolve::Come to the edge
NGsolve::Come to the edge Facet spaces and hybrid methods Jay Gopalakrishnan Portland State University Winter 2015 Download code (works only on version 6.0) for these notes from here. Come to the edge,
More informationNETGEN/NGSolve Manual
NETGEN/NGSolve Manual Appendix to the master thesis Numerische Lösungen elliptischer und parabolischer Differentialgleichungen in zwei und drei Dimensionen mit NETGEN/NGSolve Noam Arnold March 8, 2013
More informationTMA4220: Programming project - part 1
TMA4220: Programming project - part 1 TMA4220 - Numerical solution of partial differential equations using the finite element method The programming project will be split into two parts. This is the first
More informationComputational Methods for Engineering Applications II. Homework Problem Sheet 4
S. Mishra K. O. Lye L. Scarabosio Autumn Term 2015 Computational Methods for Engineering Applications II ETH Zürich D-MATH Homework Problem Sheet 4 This assignment sheet contains some tasks marked as Core
More informationAn introduction to plotting data
An introduction to plotting data Eric D. Black California Institute of Technology v2.0 1 Introduction Plotting data is one of the essential skills every scientist must have. We use it on a near-daily basis
More informationSolution of Non-Homogeneous Dirichlet Problems with FEM
Master Thesis Solution of Non-Homogeneous Dirichlet Problems with FEM Francesco Züger Institut für Mathematik Written under the supervision of Prof. Dr. Stefan Sauter and Dr. Alexander Veit August 27,
More informationAn Introduction to Netgen and NGSolve Session 1 Using Netgen/NGSolve
An Introduction to Netgen and NGSolve Session 1 Using Netgen/NGSolve Christoph Lehrenfeld Center for Computational Engineering Sciences (CCES) RWTH Aachen University University of British Columbia, November
More informationCMU CS 462/662 (INTRO TO COMPUTER GRAPHICS) HOMEWORK 0.0 MATH REVIEW/PREVIEW LINEAR ALGEBRA
CMU CS 462/662 (INTRO TO COMPUTER GRAPHICS) HOMEWORK 0.0 MATH REVIEW/PREVIEW LINEAR ALGEBRA Andrew ID: ljelenak August 25, 2018 This assignment reviews basic mathematical tools you will use throughout
More informationStar Cluster Photometry and the H-R Diagram
Star Cluster Photometry and the H-R Diagram Contents Introduction Star Cluster Photometry... 1 Downloads... 1 Part 1: Measuring Star Magnitudes... 2 Part 2: Plotting the Stars on a Colour-Magnitude (H-R)
More informationExercise 3: GIS data on the World Wide Web
Exercise 3: GIS data on the World Wide Web These web sites are a few examples of sites that are serving free GIS data. Many other sites exist. Search in Google or other search engine to find GIS data for
More informationIntegration using Gaussian Quadrature
Integration using Gaussian Quadrature Tutorials November 25, 2017 Department of Aeronautics, Imperial College London, UK Scientific Computing and Imaging Institute, University of Utah, USA 2 Chapter 1
More informationngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem)
ngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem) Christoph Lehrenfeld (with contributions from F. Heimann, J. Preuß, M. Hochsteger,...) lehrenfeld@math.uni-goettingen.de
More informationAdvanced Higher Order Finite Elements From Theory to Application with Netgen/NGSolve
Advanced Higher Order Finite Elements From Theory to Application with Netgen/NGSolve Christoph Lehrenfeld, Joachim Schöberl Center for Computational Engineering Sciences (CCES) RWTH Aachen University Vancouver
More informationCityGML XFM Application Template Documentation. Bentley Map V8i (SELECTseries 2)
CityGML XFM Application Template Documentation Bentley Map V8i (SELECTseries 2) Table of Contents Introduction to CityGML 1 CityGML XFM Application Template 2 Requirements 2 Finding Documentation 2 To
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationUsing Microsoft Excel
Using Microsoft Excel Objective: Students will gain familiarity with using Excel to record data, display data properly, use built-in formulae to do calculations, and plot and fit data with linear functions.
More informationLife Cycle of Stars. Photometry of star clusters with SalsaJ. Authors: Daniel Duggan & Sarah Roberts
Photometry of star clusters with SalsaJ Authors: Daniel Duggan & Sarah Roberts Photometry of star clusters with SalsaJ Introduction Photometry is the measurement of the intensity or brightness of an astronomical
More informationLeaf Spring (Material, Contact, geometric nonlinearity)
00 Summary Summary Nonlinear Static Analysis - Unit: N, mm - Geometric model: Leaf Spring.x_t Leaf Spring (Material, Contact, geometric nonlinearity) Nonlinear Material configuration - Stress - Strain
More informationM E R C E R W I N WA L K T H R O U G H
H E A L T H W E A L T H C A R E E R WA L K T H R O U G H C L I E N T S O L U T I O N S T E A M T A B L E O F C O N T E N T 1. Login to the Tool 2 2. Published reports... 7 3. Select Results Criteria...
More informationVPython Class 2: Functions, Fields, and the dreaded ˆr
Physics 212E Classical and Modern Physics Spring 2016 1. Introduction VPython Class 2: Functions, Fields, and the dreaded ˆr In our study of electric fields, we start with a single point charge source
More informationLesson 3: Advanced Factoring Strategies for Quadratic Expressions
Advanced Factoring Strategies for Quadratic Expressions Student Outcomes Students develop strategies for factoring quadratic expressions that are not easily factorable, making use of the structure of the
More informationIntroduction. Name: Basic Features of Sunspots. The Solar Rotational Period. Sunspot Numbers
PHYS-1050 Tracking Sunspots Spring 2013 Name: 1 Introduction Sunspots are regions on the solar surface that appear dark because they are cooler than the surrounding photosphere, typically by about 1500
More informationTask 1: Open ArcMap and activate the Spatial Analyst extension.
Exercise 10 Spatial Analyst The following steps describe the general process that you will follow to complete the exercise. Specific steps will be provided later in the step-by-step instructions component
More informationThis format is the result of tinkering with a mixed lecture format for 3 terms. As such, it is still a work in progress and we will discuss
Version 1, August 2016 1 This format is the result of tinkering with a mixed lecture format for 3 terms. As such, it is still a work in progress and we will discuss adaptations both to the general format
More informationLAB Exercise #4 - Answers The Traction Vector and Stress Tensor. Introduction. Format of lab. Preparation reading
LAB Exercise #4 - Answers The Traction Vector and Stress Tensor Due: Thursday, 26 February 2009 (Special Thanks to D.D. Pollard who pioneered this exercise in 1991) Introduction Stress concentrations in
More informationDepartment of Chemical Engineering University of California, Santa Barbara Spring Exercise 2. Due: Thursday, 4/19/09
Department of Chemical Engineering ChE 210D University of California, Santa Barbara Spring 2012 Exercise 2 Due: Thursday, 4/19/09 Objective: To learn how to compile Fortran libraries for Python, and to
More informationComp 11 Lectures. Mike Shah. July 12, Tufts University. Mike Shah (Tufts University) Comp 11 Lectures July 12, / 33
Comp 11 Lectures Mike Shah Tufts University July 12, 2017 Mike Shah (Tufts University) Comp 11 Lectures July 12, 2017 1 / 33 Please do not distribute or host these slides without prior permission. Mike
More informationTask 1: Start ArcMap and add the county boundary data from your downloaded dataset to the data frame.
Exercise 6 Coordinate Systems and Map Projections The following steps describe the general process that you will follow to complete the exercise. Specific steps will be provided later in the step-by-step
More informationDistributed Memory Parallelization in NGSolve
Distributed Memory Parallelization in NGSolve Lukas Kogler June, 2017 Inst. for Analysis and Scientific Computing, TU Wien From Shared to Distributed Memory Shared Memory Parallelization via threads (
More informationCS 7B - Fall Final Exam (take-home portion). Due Wednesday, 12/13 at 2 pm
CS 7B - Fall 2017 - Final Exam (take-home portion). Due Wednesday, 12/13 at 2 pm Write your responses to following questions on separate paper. Use complete sentences where appropriate and write out code
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More informationEXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3)
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3) 121 Textbook Reference: Knight, Chapter 13.1-3, 6. SYNOPSIS In
More informationLab 1: Numerical Solution of Laplace s Equation
Lab 1: Numerical Solution of Laplace s Equation ELEC 3105 last modified August 27, 2012 1 Before You Start This lab and all relevant files can be found at the course website. You will need to obtain an
More informationMultiphysics Modeling
11 Multiphysics Modeling This chapter covers the use of FEMLAB for multiphysics modeling and coupled-field analyses. It first describes the various ways of building multiphysics models. Then a step-by-step
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More informationUnit 15 Section 3 : Symmetry
Unit 15 Section 3 : Symmetry In this section we revise the symmetry of objects and examine the symmetry of regular polygons. There are two types of symmetry: reflective symmetry and rotational symmetry.
More informationExample problem: Adaptive solution of the 2D advection diffusion equation
Chapter 1 Example problem: Adaptive solution of the 2D advection diffusion equation In this example we discuss the adaptive solution of the 2D advection-diffusion problem Solve Two-dimensional advection-diffusion
More informationprogeny. Observe the phenotypes of the F1 progeny flies resulting from this reciprocal cross.
Name Fruit Fly Exercise 8 Goal In this exercise, you will use the StarGenetics, a software tool that simulates mating experiments, to perform your own simulated genetic crosses to analyze the mode of inheritance
More informationFoundations of Computation
The Australian National University Semester 2, 2018 Research School of Computer Science Tutorial 1 Dirk Pattinson Foundations of Computation The tutorial contains a number of exercises designed for the
More informationCalculating NMR Chemical Shifts for beta-ionone O
Calculating NMR Chemical Shifts for beta-ionone O Molecular orbital calculations can be used to get good estimates for chemical shifts. In this exercise we will calculate the chemical shifts for beta-ionone.
More informationApplied C Fri
Applied C++11 2013-01-25 Fri Outline Introduction Auto-Type Inference Lambda Functions Threading Compiling C++11 C++11 (formerly known as C++0x) is the most recent version of the standard of the C++ Approved
More informationStudying Topography, Orographic Rainfall, and Ecosystems (STORE)
Introduction Studying Topography, Orographic Rainfall, and Ecosystems (STORE) Lesson: Using ArcGIS Explorer to Analyze the Connection between Topography, Tectonics, and Rainfall GIS-intensive Lesson This
More informationLESSON 9.1 ROOTS AND RADICALS
LESSON 9.1 ROOTS AND RADICALS LESSON 9.1 ROOTS AND RADICALS 67 OVERVIEW Here s what you ll learn in this lesson: Square Roots and Cube Roots a. Definition of square root and cube root b. Radicand, radical
More informationInstructor Notes for Chapters 3 & 4
Algebra for Calculus Fall 0 Section 3. Complex Numbers Goal for students: Instructor Notes for Chapters 3 & 4 perform computations involving complex numbers You might want to review the quadratic formula
More informationLivePhoto Physics Activity 3. Velocity Change. Motion Detector. Sample
LivePhoto Physics Activity 3 Name: Date: Analyzing Position vs. Time Graphs: The most fundamental measurements of motion involve the determination of an object s location at a series of times. A very effective
More informationExercise in gas chromatography
AM0925 Exercise in gas chromatography Introduction The purpose of this exercise is to investigate the effect of different carrier gases (helium and nitrogen) on the Golay equation and the optimal carrier
More informationChemistry 14CL. Worksheet for the Molecular Modeling Workshop. (Revised FULL Version 2012 J.W. Pang) (Modified A. A. Russell)
Chemistry 14CL Worksheet for the Molecular Modeling Workshop (Revised FULL Version 2012 J.W. Pang) (Modified A. A. Russell) Structure of the Molecular Modeling Assignment The molecular modeling assignment
More information10.1.1Theorem: Module 4 : Local / Global Maximum / Minimum and Curve Sketching
\ Module 4 : Local / Global Maximum / Minimum and Curve Sketching Lecture 10 : Sufficient conditions for increasing / decreasing [Section 10.1] Objectives In this section you will learn the following :
More informationLesson 5b Solving Quadratic Equations
Lesson 5b Solving Quadratic Equations In this lesson, we will continue our work with Quadratics in this lesson and will learn several methods for solving quadratic equations. The first section will introduce
More informationSection 2.1 Differential Equation and Solutions
Section 2.1 Differential Equation and Solutions Key Terms: Ordinary Differential Equation (ODE) Independent Variable Order of a DE Partial Differential Equation (PDE) Normal Form Solution General Solution
More informationCore Mathematics 2 Algebra
Core Mathematics 2 Algebra Edited by: K V Kumaran Email: kvkumaran@gmail.com Core Mathematics 2 Algebra 1 Algebra and functions Simple algebraic division; use of the Factor Theorem and the Remainder Theorem.
More informationLAB 3: WORK AND ENERGY
1 Name Date Lab Day/Time Partner(s) Lab TA (CORRECTED /4/05) OBJECTIVES LAB 3: WORK AND ENERGY To understand the concept of work in physics as an extension of the intuitive understanding of effort. To
More informationName Date Class. Figure 1. The Google Earth Pro drop-down menu.
GIS Student Walk-Through Worksheet Procedure 1. Import historical tornado and hurricane data into Google Earth Pro by following these steps: A. In the Google Earth Pro drop-down menu > click File > Import
More informationConformational Analysis of n-butane
Conformational Analysis of n-butane In this exercise you will calculate the Molecular Mechanics (MM) single point energy of butane in various conformations with respect to internal rotation around the
More informationC++ For Science and Engineering Lecture 17
C++ For Science and Engineering Lecture 17 John Chrispell Tulane University Monday October 4, 2010 Functions and C-Style Strings Three methods for representing the C-style string: An array of char A quoted
More informationMachine Learning: Homework 5
0-60 Machine Learning: Homework 5 Due 5:0 p.m. Thursday, March, 06 TAs: Travis Dick and Han Zhao Instructions Late homework policy: Homework is worth full credit if submitted before the due date, half
More informationIOP2601. Some notes on basic mathematical calculations
IOP601 Some notes on basic mathematical calculations The order of calculations In order to perform the calculations required in this module, there are a few steps that you need to complete. Step 1: Choose
More informationLab 15 Taylor Polynomials
Name Student ID # Instructor Lab Period Date Due Lab 15 Taylor Polynomials Objectives 1. To develop an understanding for error bound, error term, and interval of convergence. 2. To visualize the convergence
More informationCreating Axisymmetric Models in FEMAP
Creating Axisymmetric Models in FEMAP 1. Introduction NE/Nastran does not support 2-d axisymmetric elements. 3-d axisymmetric models are supported, and can be generated with a few additional steps. The
More informationSignal Processing First Lab 11: PeZ - The z, n, and ˆω Domains
Signal Processing First Lab : PeZ - The z, n, and ˆω Domains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations to be made when
More informationLAB 2 - ONE DIMENSIONAL MOTION
Name Date Partners L02-1 LAB 2 - ONE DIMENSIONAL MOTION OBJECTIVES Slow and steady wins the race. Aesop s fable: The Hare and the Tortoise To learn how to use a motion detector and gain more familiarity
More informationNumerical solution of the time-independent 1-D Schrödinger equation. Matthias E. Möbius. September 24, 2010
Numerical solution of the time-independent 1-D Schrödinger equation Matthias E. Möbius September 24, 2010 1 Aim of the computational lab Numerical solution of the one-dimensional stationary Schrödinger
More informationProofs, Strings, and Finite Automata. CS154 Chris Pollett Feb 5, 2007.
Proofs, Strings, and Finite Automata CS154 Chris Pollett Feb 5, 2007. Outline Proofs and Proof Strategies Strings Finding proofs Example: For every graph G, the sum of the degrees of all the nodes in G
More informationDate: Summer Stem Section:
Page 1 of 7 Name: Date: Summer Stem Section: Summer assignment: Build a Molecule Computer Simulation Learning Goals: 1. Students can describe the difference between a molecule name and chemical formula.
More informationSolving Differential Equations on 2-D Geometries with Matlab
Solving Differential Equations on 2-D Geometries with Matlab Joshua Wall Drexel University Philadelphia, PA 19104 (Dated: April 28, 2014) I. INTRODUCTION Here we introduce the reader to solving partial
More informationAssignment 1 Physics/ECE 176
Assignment 1 Physics/ECE 176 Made available: Thursday, January 13, 211 Due: Thursday, January 2, 211, by the beginning of class. Overview Before beginning this assignment, please read carefully the part
More informationModule 10: Free Vibration of an Undampened 1D Cantilever Beam
Module 10: Free Vibration of an Undampened 1D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry 4 Preprocessor 6 Element Type 6 Real Constants and Material Properties 7
More informationModeling a Composite Slot Cross-Section for Torsional Analysis
Modeling a Composite Slot Cross-Section for Torsional Analysis The cross-section in question is shown below (see also p. 432 in the textbook). Due to double symmetry, only one-quarter of the cross-section
More informationthrough any three given points if and only if these points are not collinear.
Discover Parabola Time required 45 minutes Teaching Goals: 1. Students verify that a unique parabola with the equation y = ax + bx+ c, a 0, exists through any three given points if and only if these points
More informationIntroduction to lambda calculus Part 6
Introduction to lambda calculus Part 6 Antti-Juhani Kaijanaho 2017-02-16 1 Untyped lambda calculus 2 Typed lambda calculi 2.1 Dynamically typed lambda calculus with integers 2.2 A model of Lisp 2.3 Simply
More informationLab 1: Handout GULP: an Empirical energy code
Lab 1: Handout GULP: an Empirical energy code We will be using the GULP code as our energy code. GULP is a program for performing a variety of types of simulations on 3D periodic solids, gas phase clusters,
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
More informationSynteny Portal Documentation
Synteny Portal Documentation Synteny Portal is a web application portal for visualizing, browsing, searching and building synteny blocks. Synteny Portal provides four main web applications: SynCircos,
More informationDOC SOLVING SYSTEMS OF EQUATIONS BY SUBSTITUTION CALCULATOR
02 July, 2018 DOC SOLVING SYSTEMS OF EQUATIONS BY SUBSTITUTION CALCULATOR Document Filetype: PDF 471.03 KB 0 DOC SOLVING SYSTEMS OF EQUATIONS BY SUBSTITUTION CALCULATOR In case you have to have assistance
More informationHandling Raster Data for Hydrologic Applications
Handling Raster Data for Hydrologic Applications Prepared by Venkatesh Merwade Lyles School of Civil Engineering, Purdue University vmerwade@purdue.edu January 2018 Objective The objective of this exercise
More informationProject IV Fourier Series
Project IV Fourier Series Robert Jerrard Goal of the project To develop understanding of how many terms of a Fourier series are required in order to well-approximate the original function, and of the differences
More informationDoes Air Have Mass? 1.3 Investigate
Does Air Have Mass? A water bottle with no water in it is sitting on the table. Your friend says the bottle is empty. What would you tell your friend about the bottle? Would you describe the bottle as
More informationA short introduction to FreeFem++
A short introduction to FreeFem++ M. Ersoy Basque Center for Applied Mathematics 12 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 HOW TO solve steady pdes : e.g. Laplace equation
More informationPhotometry of a Globular Cluster
Name(s): Date: Course/Section: Grade: Photometry of a Globular Cluster Learning Objectives: Students will determine the age of a globular cluster using photometry. Checklist: Complete the pre- lab quiz
More informationUniversity of TN Chattanooga Physics 1040L 8/18/2012 PHYSICS 1040L LAB LAB 4: R.C. TIME CONSTANT LAB
PHYSICS 1040L LAB LAB 4: R.C. TIME CONSTANT LAB OBJECT: To study the discharging of a capacitor and determine the time constant for a simple circuit. APPARATUS: Capacitor (about 24 μf), two resistors (about
More informationLab Partner(s) TA Initials (on completion) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE 117 Textbook Reference: Walker, Chapter 10-1,2, Chapter 11-1,3 SYNOPSIS
More informationData Structures & Database Queries in GIS
Data Structures & Database Queries in GIS Objective In this lab we will show you how to use ArcGIS for analysis of digital elevation models (DEM s), in relationship to Rocky Mountain bighorn sheep (Ovis
More informationRelative Photometry with data from the Peter van de Kamp Observatory D. Cohen and E. Jensen (v.1.0 October 19, 2014)
Relative Photometry with data from the Peter van de Kamp Observatory D. Cohen and E. Jensen (v.1.0 October 19, 2014) Context This document assumes familiarity with Image reduction and analysis at the Peter
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationTemperature measurement
Luleå University of Technology Johan Carlson Last revision: July 22, 2009 Measurement Technology and Uncertainty Analysis - E7021E Lab 3 Temperature measurement Introduction In this lab you are given a
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More information1 Solution of Electrostatics Problems with COM- SOL
1 Solution of Electrostatics Problems with COM- SOL This section gives examples demonstrating how Comsol can be used to solve some simple electrostatics problems. 1.1 Laplace s Equation We start with a
More informationConnect the Vernier spectrometer to your lap top computer and power the spectrometer if necessary. Start LoggerPro on your computer.
Connect the Vernier spectrometer to your lap top computer and power the spectrometer if necessary. Start LoggerPro on your computer. The screen shown in Fig. 1 may be displayed. If status line displays
More informationUsing the Stock Hydrology Tools in ArcGIS
Using the Stock Hydrology Tools in ArcGIS This lab exercise contains a homework assignment, detailed at the bottom, which is due Wednesday, October 6th. Several hydrology tools are part of the basic ArcGIS
More informationAlgebra III and Trigonometry Summer Assignment
Algebra III and Trigonometry Summer Assignment Welcome to Algebra III and Trigonometry! This summer assignment is a review of the skills you learned in Algebra II. Please bring this assignment with you
More informationDSP First Lab 11: PeZ - The z, n, and ωdomains
DSP First Lab : PeZ - The, n, and ωdomains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations to be made when using the PeZ GUI.
More informationWMS 10.1 Tutorial GSSHA Applications Precipitation Methods in GSSHA Learn how to use different precipitation sources in GSSHA models
v. 10.1 WMS 10.1 Tutorial GSSHA Applications Precipitation Methods in GSSHA Learn how to use different precipitation sources in GSSHA models Objectives Learn how to use several precipitation sources and
More informationAMS 132: Discussion Section 2
Prof. David Draper Department of Applied Mathematics and Statistics University of California, Santa Cruz AMS 132: Discussion Section 2 All computer operations in this course will be described for the Windows
More informationEOSC 110 Reading Week Activity, February Visible Geology: Building structural geology skills by exploring 3D models online
EOSC 110 Reading Week Activity, February 2015. Visible Geology: Building structural geology skills by exploring 3D models online Geological maps show where rocks of different ages occur on the Earth s
More informationNewton's Laws and Atwood's Machine
Newton's Laws and Atwood's Machine Purpose: In this lab we will verify Newton's Second Law of Motion within estimated uncertainty and propose an explanation if verification is not within estimated uncertainty.
More informationThe Quantizing functions
The Quantizing functions What is quantizing? Quantizing in its fundamental form is a function that automatically moves recorded notes, positioning them on exact note values: For example, if you record
More informationLED Lighting Facts: Product Submission Guide
LED Lighting Facts: Product Submission Guide NOVEMBER 2017 1 P a g e L E D L i g h t i n g F a c t s : M a n u f a c t u r e r P r o d u c t S u b m i s s i o n G u i d e TABLE OF CONTENTS Section 1) Accessing
More informationLAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION
LAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION In this lab you will learn how to use Excel to display the relationship between two quantitative variables, measure the strength and direction of the
More informationScientific Computing Software Concepts for Solving Partial Differential Equations
Scientific Computing Software Concepts for Solving Partial Differential Equations Joachim Schöberl WS04/05 Abstract Many problems in science and engineering are described by partial differential equations.
More information