Introduction to Parallelism in CASTEP
|
|
- Lionel Carson
- 6 years ago
- Views:
Transcription
1 to ism in CASTEP Phil Hasnip August 2016
2 Recap: what does CASTEP do? CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei: Ĥ k [ρ]ψ bk = E bk ψ bk where particle b has the bth solution ( band ) at the Brilliouin zone sampling point k, and Ĥ k [ρ] = 2 2m 2 + ˆV HXC [ρ] + ˆV ext.
3 Bloch s theorem and plane-waves Recall that Bloch s theorem let us write: ψ bk (r) = e ik r u bk (r), where u k (r) is periodic and e ik.r is an arbitrary phase factor. We express u k (r) as a Fourier series: u bk (r) = G c Gbk e ig r ψ bk (r) = e ik r G c Gbk e ig.r = G c Gbk e i(g+k) r where c Gbk are complex Fourier coefficients, and the sum is over all the reciprocal lattice vectors, or G-vectors.
4 The wavefunction The wavefunction is one of the main data objects in Castep: ψ bk (r) = G c Gbk e i(g+k) r The complex coefficients c Gbk are what Castep is trying to compute, and take up a lot of the computer s memory. G : a reciprocal lattice vector ( G-vector ) b : band index k : a Brillouin zone sampling point ( k-point )
5 k-point sampling The bands for different k-points are independent of each other, so we get a different set of Kohn-Sham equations at each: Ĥ k [ρ]ψ bk = E bk ψ bk where ρ(r) = bk ψ bk (r) 2
6 Where does CASTEP spend its time? Applying Ĥk to ψ bk The kinetic energy is applied in reciprocal-space The local potential is applied in real-space We need to Fourier transform between the two spaces. Orthogonalisation of ψ bk We need to ensure our trial bands are orthogonal to each other. We compute the overlap matrix between all pairs of bands, and invert it.
7 Fourier transforms A 3D Fourier transform can be performed as 3 separate 1D transformations one in each direction (x, y and z). Time to transform ψ bk (G) ψ bk (r) scales as N G log(n G ). Every band at every k-point has to be transformed, so total time is N G N b N k log(n G ).
8 Orthogonalisation We construct the band-overlap matrix for each k-point S nmk = ψ nk ψ mk. Total time scales as N G N 2 b N k. Invert S k to find an orthogonalising transformation at each k-point. Total time scales as N 3 b N k. Apply transformation to get orthogonal bands. Total time scales as N G N 2 b N k.
9 Large calculations As we simulate larger and larger systems, N G and N b increase and N k decreases. Time for Fourier transforms scales as N G log N G N b N k. Time for orthogonalisation scales as N G N 2 b N k. Orthogonalisation dominates in large calculations.
10 As the simulation system gets bigger and bigger, the orthogonalisation time dominates. We want to be able to use more computer cores in our calculation to speed it up how can we do this?
11 k-point parallelism Bands at different k-points are almost entirely independent of each other give each core a subset of the k-points. each core solves a subset of Kohn-Sham equations Cores only communicate when constructing the density ρ(r) = bk ψ bk (r) 2
12 TiN Benchmark The TiN simulation is a small standard benchmark 33 atoms 8 k-points 164 bands 10,972 G-vectors
13 k-point parallelism in action
14 k-points and large systems k-point parallelism is almost perfect. As simulations get bigger, N k gets smaller. the bigger the simulation, the fewer the cores we can use!
15 G-vector parallelism S nmk = ψ nk ψ mk = G c Gnk c Gmk Orthogonalisation is a sum over G-vectors give each core a subset of G-vectors. Contributions to S are summed over cores. N G is large so can use lots of cores. As simulation size increases, N G also increases.
16 G-vector parallelism in action
17 G-vector parallelism Each core only has some of the G-vectors. Each core only has some of the real-space r-vectors. Fourier transform: all G-vectors contribute at all points in real-space.
18 G-vector parallelism 3D transform can be performed as 3 1D transforms Give each core all G-vectors in a column in z Each core does transform in z All cores swap data so they have y-column data Each core does transform in y All cores swap data so they have x-column data Each core does transform in x Each core ends up with real-space data in x
19 G-vector parallelism Start: G-vectors inside cut-off sphere put on grid.
20 G-vector parallelism Now perform FFT in z-direction...
21 G-vector parallelism Transpose (swap) data into y-columns.
22 G-vector parallelism Now perform FFT in y-direction...
23 G-vector parallelism Transpose data into x-columns.
24 G-vector parallelism Now perform FFT in x-direction...
25 G-vector parallelism Now have real-space data in x-columns.
26 G-vector parallelism Actual transforms distribute well Transpositions are a problem Every core has to communicate with every other core! time scales as N 2 core. as N core increases, Fourier transform will dominate. when the communication time is comparable to the compute time, there s no point using more cores (it might even make CASTEP slower). There are ways to optimise the FFT, but the basic problem remains the same.
27 G-vector parallelism in action
28 k and G parallelism k-point and G-vector parallelism is independent Can combine both to improve scaling E.g. if N k = 2, N G = 9, 000 and N core = 6: Data k-point 1 k-point 2 G-vectors 1-3,000 core 1 core 4 G-vectors 3,001-6,001 core 2 core 5 G-vectors 6,001-9,000 core 3 core 6 For any k-point, the G-vector data is split across 3 cores this is 3-way G-vector parallelism. For any subset of G-vectors the k-point data is split between 2 cores this is 2-way k-point parallelism
29 k+g parallelism in action
30 Al 2 O 3 Benchmark The Al 2 O 3 surface simulation is a larger standard benchmark 270 atoms 2 k-points 778 bands 88,184 G-vectors
31 k+g parallelism in action
32 Multi-core nodes Common to have multi-core processors May have several processors per node Each core on the same node can access shared memory (RAM) Communications can use this shared RAM instead of the network Access via the parameter: num_proc_in_smp : <integer>
33 k+g+smp parallelism in action
34 Optimal performance It is always worth exploiting k-point parallelism when you can, but not all computers let you run on any number of cores. If you can t use N core = N k then try to have a high common factor between them. E.g. if N k = 35, N core = 35 will give an excellent speed-up, but N core = 5 or 7 will also be very efficient. N core = 20 or 21 would use some G-vector parallelism, but also give good efficiency. (Note that 2-way G-vector parallelism is not very quick, so N core = 10 or 14 might not be the best choices.)
35 Optimal performance Remember that as you increase the G-vector parallelism, the communication time increases. Eventually your calculation will scale poorly, and if you keep increasing N core it will even start to run slower. CASTEP always defaults to using as much k-point parallelism as it can, and then uses G-vector parallelism across any other cores.
36 Very large calculations and the For isolated or very large simulation systems you only need 1 k-point. For well-isolated or extremely large simulations this k-point can be the special point k=(0,0,0), called the. Why do we care?
37 calculations Bands at Γ are real in real-space, not complex. the Fourier coefficients for -G are the complex conjugate of those at G don t bother with -G; only need half the G-vectors Inner products are real, not complex we don t need to bother computing the imaginary parts. Bands take up only half as much memory. FFT about 2x faster, orthogonalisation 8x faster. CASTEP detects if you re only using the and uses these optimisations automatically.
38 calculations No k-point parallelism (N k = 1). Orthogonalisation speed-up better than FFT one Can show poorer scaling. Worth using if this k-point sampling is sufficiently accurate. Occasionally worth using a bigger simulation cell if it allows accurate sampling.
39 Wavefunction revision Recall that the wavefunction is: ψ bk (r) = G c Gbk e i(g+k) r G : a reciprocal lattice vector ( G-vector ) distributed in G-vector parallelism k : a Brillouin zone sampling point ( k-point ) distributed in k-point parallelism b : band index can we distribute bands? Yes!
40 Hamiltonian is the same for all bands at the same k-point Fourier transforms of different bands are independent perfect scaling with band-parallelism when applying the Hamiltonian.
41 Orthogonalisation Need to construct overlap matrix S at each k-point S nm = ψ n ψ m Inner product is between all pairs of bands Need all-to-all communications as band-parallelism increases, communication dominates
42 Band distribution We distribute the bands in a round-robin fashion, e.g. if N b = 11 and N core = 3: Core Bands 1 1,4,7,10 2 2,5,8,11 3 3,6,9 For the FFTs each core just transforms its own bands.
43 Band parallelism performance FFTs scale perfectly Orthogonalisation requires all-to-all amongst the cores communication time scales as N 2 core communications dominate as N core increases
44 k, G and B parallelism k-point, G-vector and band-parallelism are independent Can combine all three to improve scaling Define: kp-group: group of cores with same G-vectors and bands, but different k-points gv-group: group of cores with same k-points and bands, but different G-vectors bnd-group: group of cores with same G-vectors and k-points, but different bands
45 k+g+smp+band parallelism in action
46 CASTEP performance Everything scales well with k-point parallelism As the number of cores in the gv-group increases, the communication time in the FFT dominates As the number of cores in the band-group increases, the communication time in the orthogonalisation dominates need to find the right balance between gv- and bnd-parallelism.
47 Using band-parallelism Relatively new functionality Will compute ground-state energies, forces and stresses Accessed via a devel_code string, e.g.: devel_code : bandpar=2 in your.param file would use 2-way band-parallelism.
48 Why a devel_code? ism and G-vector parallelism use the network differently Low-latency networks good G-vector parallelism High-bandwidth networks good band-parallelism Difficult for CASTEP to decide the best parallelisation strategy Some limitations to current implementation
49 Limitations Not all tasks supported; currently supports tasks: energy geometryoptimisation moleculardynamics DFT+U not supported These limitations are temporary!
50 k-point parallelism is very efficient, but you eventually run out of k-points G-vector parallelism is good but becomes worse as you use more and more cores. ism available for many calculations; fairly good but becomes worse with more and more cores Combining these parallelisms allows Castep to scale well to many cores, but careful choice of the no. cores can improve performance considerably: Know how many k-points you re using (-dryrun) For phonons, use the phonon_kpoints tool NB Path-integral MD has an extra level of parallelism task-farming which is very efficient.
Introducing k-point Parallelism Into VASP. Asimina Maniopoulou Numerical Algorithms Group Ltd, HECToR CSE
Introducing k-point Parallelism Into VASP Asimina Maniopoulou Numerical Algorithms Group Ltd, HECToR CSE Plan for the talk VASP Implementation of k-point parallelization Benchmarks -performance enhancement
More informationIntroducing k-point Parallelism into VASP
Introducing k-point Parallelism into VASP Asimina Maniopoulou HECToR CSE Team, NAG Ltd., UK support@nag.co.uk August 24, 2011 Contents 1 Introduction 1 1.1 VASP..............................................
More informationPUV Wave Directional Spectra How PUV Wave Analysis Works
PUV Wave Directional Spectra How PUV Wave Analysis Works Introduction The PUV method works by comparing velocity and pressure time series. Figure 1 shows that pressure and velocity (in the direction of
More informationFast Floating Point Compression on the Cell BE Processor
Fast Floating Point Compression on the Cell BE Processor Ajith Padyana, T.V. Siva Kumar, P.K.Baruah Sri Satya Sai University Prasanthi Nilayam - 515134 Andhra Pradhesh, India ajith.padyana@gmail.com, tvsivakumar@gmail.com,
More informationWave Motion. interference destructive interferecne constructive interference in phase. out of phase standing wave antinodes resonant frequencies
Wave Motion Vocabulary mechanical waves pulse continuous periodic wave amplitude period wavelength period wave velocity phase transverse wave longitudinal wave intensity displacement amplitude phase velocity
More informationSeismic Survey Designs for Converted Waves
Seismic Survey Designs for Converted Waves James A. Musser* GMG/AXIS Inc., Denver, CO 1720 Red Cloud Road, Longmont, CO, 80501, USA jmusser@gmg.com ABSTRACT Designing converted wave 3D surveys is considerably
More informationSeismic waves. Seismic waves, like all waves, transfer energy from one place to another without moving material. Seismic Waves 1 Author Paul Denton
Seismic waves When an earthquake happens deep underground a crack will start to open on a pre-existing line of weakness in the Earth s brittle crust. This crack will then grow larger and larger, relieving
More informationROSE-HULMAN INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering. Mini-project 3 Tennis ball launcher
Mini-project 3 Tennis ball launcher Mini-Project 3 requires you to use MATLAB to model the trajectory of a tennis ball being shot from a tennis ball launcher to a player. The tennis ball trajectory model
More informationFirst Steps Towards the AEI 10m Prototype Single Arm Test Auto Alignment
First Steps Towards the AEI 10m Prototype Single Arm Test Auto Alignment Sean Leavey and the AEI 10m Prototype Team Albert Einstein Institute Hanover Germany 26th March 2014 The AEI 10m Prototype Test
More informationThe Evolution of Transport Planning
The Evolution of Transport Planning On Proportionality and Uniqueness in Equilibrium Assignment Michael Florian Calin D. Morosan Background Several bush-based algorithms for computing equilibrium assignments
More informationSection 1 Types of Waves. Distinguish between mechanical waves and electromagnetic waves.
Section 1 Types of Waves Objectives Recognize that waves transfer energy. Distinguish between mechanical waves and electromagnetic waves. Explain the relationship between particle vibration and wave motion.
More informationConvection Current Exploration:
Heat on Earth 8.10A RECOGNIZE THAT THE SUN PROVIDES THE ENERGY THAT DRIVES CONVECTION WITHIN THE ATMOSPHERE AND OCEANS, PRODUCING WINDS AND OCEAN CURRENTS [INCORPORATE 6.6B INTO CONVECTION] A few reminders
More information6/16/2010 DAG Execu>on Model, Work and Depth 1 DAG EXECUTION MODEL, WORK AND DEPTH
6/16/2010 DAG Execu>on Model, Work and Depth 1 DAG EXECUTION MODEL, WORK AND DEPTH 6/16/2010 DAG Execu>on Model, Work and Depth 2 Computa>onal Complexity of (Sequen>al) Algorithms Model: Each step takes
More informationAdvanced pre and post-processing in Windsim
University of Perugia Department of Industrial Engineering Francesco Castellani Advanced pre and post-processing in Windsim CONTENTS Pre-processing 1) Domain control: *.gws construction 2) Advanced grid
More informationChapter 15 Wave Motion. Copyright 2009 Pearson Education, Inc.
Chapter 15 Wave Motion 15-1 Characteristics of Wave Motion All types of traveling waves transport energy. Study of a single wave pulse shows that it is begun with a vibration and is transmitted through
More informationAlgorithm for Line Follower Robots to Follow Critical Paths with Minimum Number of Sensors
International Journal of Computer (IJC) ISSN 2307-4523 (Print & Online) Global Society of Scientific Research and Researchers http://ijcjournal.org/ Algorithm for Line Follower Robots to Follow Critical
More informationMassey Method. Introduction. The Process
Massey Method Introduction Massey s Method, also referred to as the Point Spread Method, is a rating method created by mathematics professor Kenneth Massey. It is currently used to determine which teams
More informationWaves. harmonic wave wave equation one dimensional wave equation principle of wave fronts plane waves law of reflection
Waves Vocabulary mechanical wave pulse continuous periodic wave amplitude wavelength period frequency wave velocity phase transverse wave longitudinal wave intensity displacement wave number phase velocity
More informationWorkshop 1: Bubbly Flow in a Rectangular Bubble Column. Multiphase Flow Modeling In ANSYS CFX Release ANSYS, Inc. WS1-1 Release 14.
Workshop 1: Bubbly Flow in a Rectangular Bubble Column 14. 5 Release Multiphase Flow Modeling In ANSYS CFX 2013 ANSYS, Inc. WS1-1 Release 14.5 Introduction This workshop models the dispersion of air bubbles
More informationRanking teams in partially-disjoint tournaments
Ranking teams in partially-disjoint tournaments Alex Choy Mentor: Chris Jones September 16, 2013 1 Introduction Throughout sports, whether it is professional or collegiate sports, teams are ranked. In
More informationCS 4649/7649 Robot Intelligence: Planning
CS 4649/7649 Robot Intelligence: Planning Differential Kinematics, Probabilistic Roadmaps Sungmoon Joo School of Interactive Computing College of Computing Georgia Institute of Technology S. Joo (sungmoon.joo@cc.gatech.edu)
More information1. Predict what will happen in the following situation. Sketch below your prediction of the interference pattern when the waves overlap:
Name P Physics1 Waves Unit WVE EHVIOR Date INTRODUCTION y the end of these activities you should have: an understanding of interference as a behavior which occurs as two waves interfere while moving through
More informationDevelopment of Fluid-Structure Interaction Program for the Mercury Target
Chapter 4 Epoch Making Simulation Development of Fluid-Structure Interaction Program for the Mercury Target Project Representative Chuichi Arakawa Authors Chuichi Arakawa Takuma Kano Ryuta Imai Japan Atomic
More informationGOLOMB Compression Technique For FPGA Configuration
GOLOMB Compression Technique For FPGA Configuration P.Hema Assistant Professor,EEE Jay Shriram Group Of Institutions ABSTRACT Bit stream compression is important in reconfigurable system design since it
More informationCAM Final Report John Scheele Advisor: Paul Ohmann I. Introduction
CAM Final Report John Scheele Advisor: Paul Ohmann I. Introduction Herds are a classic complex system found in nature. From interactions amongst individual animals, group behavior emerges. Historically
More informationChapter # 08 Waves. [WAVES] Chapter # 08
Chapter # 08 Waves Q2) Write short answers of the following questions. i) What is the difference between progressive and stationary waves? Answer: Progressive Waves 1 Progressive waves are the result of
More informationCFD ANALYSIS AND COMPARISON USING ANSYS AND STAR-CCM+ OF MODEL AEROFOIL SELIG 1223
International Journal of Mechanical Engineering and Technology (IJMET) Volume 8, Issue 11, November 2017, pp. 312 318, Article ID: IJMET_08_11_034 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=8&itype=11
More informationThe Usage of Propeller Tunnels For Higher Efficiency and Lower Vibration. M. Burak Şamşul
The Usage of Propeller Tunnels For Higher Efficiency and Lower Vibration M. Burak Şamşul ITU AYOC 2014 - Milper Pervane Teknolojileri Company Profile MILPER is established in 2011 as a Research and Development
More informationPIG MOTION AND DYNAMICS IN COMPLEX GAS NETWORKS. Dr Aidan O Donoghue, Pipeline Research Limited, Glasgow
PIG MOTION AND DYNAMICS IN COMPLEX GAS NETWORKS Dr Aidan O Donoghue, Pipeline Research Limited, Glasgow A model to examine pigging and inspection of gas networks with multiple pipelines, connections and
More informationInvestigating the Problems of Ship Propulsion on a Supercomputer
Ship Elbrus powered by What I See is What I Understand. Так мы его задумали Investigating the Problems of Ship Propulsion on a Supercomputer A.A. Aksenov, S.V. Zhluktov, D.P. Silaev, S.A. Kharchenko, E.A.
More informationMotion Control of a Bipedal Walking Robot
Motion Control of a Bipedal Walking Robot Lai Wei Ying, Tang Howe Hing, Mohamed bin Hussein Faculty of Mechanical Engineering Universiti Teknologi Malaysia, 81310 UTM Skudai, Johor, Malaysia. Wylai2@live.my
More informationEffect of Diameter on the Aerodynamics of Sepaktakraw Balls, A Computational Study
ISSN 1750-9823 (print) International Journal of Sports Science and Engineering Vol. 03 (2009) No. 01, pp. 017-021 Effect of Diameter on the Aerodynamics of Sepaktakraw Balls, A Computational Study Zahari
More informationVibrations are the sources of waves. A vibration creates a disturbance in a given medium, that disturbance travels away from the source, carrying
Vibrations are the sources of waves. A vibration creates a disturbance in a given medium, that disturbance travels away from the source, carrying energy with it, we call this traveling disturbance a wave.
More informationEE 364B: Wind Farm Layout Optimization via Sequential Convex Programming
EE 364B: Wind Farm Layout Optimization via Sequential Convex Programming Jinkyoo Park 1 Introduction In a wind farm, the wakes formed by upstream wind turbines decrease the power outputs of downstream
More informationRunning WIEN2K on Ranger
Running WIEN2K on Ranger with both coarse and fine parallelism Hang Liu Texas Advanced Computing Center May 8, 2012 Outline Introduction Setting WIEN2K in User s Account 2 Introduction WIEN2K WIEN2K Software
More informationAuthor s Name Name of the Paper Session. Positioning Committee. Marine Technology Society. DYNAMIC POSITIONING CONFERENCE September 18-19, 2001
Author s Name Name of the Paper Session PDynamic Positioning Committee Marine Technology Society DYNAMIC POSITIONING CONFERENCE September 18-19, 2001 POWER PLANT SESSION A New Concept for Fuel Tight DP
More informationOut-of-Core Cholesky Factorization Algorithm on GPU and the Intel MIC Co-processors
Out-of-Core Cholesky Factorization Algorithm on GPU and the Intel MIC Co-processors Ben Chan (Chinese University of Hong Kong) Nina Qian (Chinese University of Hong Kong) Mentors: Ed D Azevedo (ORNL) Shiquan
More informationAnalysis and Research of Mooring System. Jiahui Fan*
nd International Conference on Computer Engineering, Information Science & Application Technology (ICCIA 07) Analysis and Research of Mooring System Jiahui Fan* School of environment, North China Electric
More informationg L Agenda Chapter 13 Problem 28 Equations of Motion for SHM: What if we have friction or drag? Driven Oscillations; Resonance 4/30/14 k m f = 1 2π
Agenda Today: HW quiz, More simple harmonic motion and waves Thursday: More waves Midterm scores will be posted by Thursday. Chapter 13 Problem 28 Calculate the buoyant force due to the surrounding air
More informationTransverse waves cause particles to vibrate perpendicularly to the direction of the wave's motion (e.g. waves on a string, ripples on a pond).
Waves Introduction A vibration must be the source of a wave. Waves in turn also cause vibrations. They are intrinsically connected. Waves transmit energy. There are different ways in which waves can be
More informationAGA Swiss McMahon Pairing Protocol Standards
AGA Swiss McMahon Pairing Protocol Standards Final Version 1: 2009-04-30 This document describes the Swiss McMahon pairing system used by the American Go Association (AGA). For questions related to user
More informationThis portion of the piping tutorial covers control valve sizing, control valves, and the use of nodes.
Piping Tutorial A piping network represents the flow of fluids through several pieces of equipment. If sufficient variables (flow rate and pressure) are specified on the piping network, CHEMCAD calculates
More informationAnnual Unit Extended Team Championships. Sanctioning and Reporting Guide. January, 2017
Annual Unit Extended Team Championships Sanctioning and Reporting Guide January, 2017-1 - Table of Contents 1) General Principles page 3 Masterpoint Awards page 3 Sanction Application and Fees page 3 Conditions
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 9, 2010 http://acousticalsociety.org/ 159th Meeting Acoustical Society of America/NOISE-CON 2010 Baltimore, Maryland 19-23 April 2010 Session 1pBB: Biomedical
More informationChapter 11 Waves. Waves transport energy without transporting matter. The intensity is the average power per unit area. It is measured in W/m 2.
Energy can be transported by particles or waves: Chapter 11 Waves A wave is characterized as some sort of disturbance that travels away from a source. The key difference between particles and waves is
More information*MAT GAS MIXTURE, a new gas mixture model for airbag applications
4 th European LS-DYNA Users Conference Occupant III / Airbag *MAT GAS MIXTURE, a new gas mixture model for airbag applications Author: Lars Olovsson LSTC Correspondence: tel/fax: +46-852 030 881 e-mail:
More informationLECTURE 5 TRAVELING WAVES. Instructor: Kazumi Tolich
LECTURE 5 TRAVELING WAVES Instructor: Kazumi Tolich Lecture 5 2 15.1 The wave model Mechanical waves Electromagnetic and matter waves The transverse and longitudinal waves Traveling waves Waves on a string
More informationHorse Farm Management s Report Writer. User Guide Version 1.1.xx
Horse Farm Management s Report Writer User Guide Version 1.1.xx August 30, 2001 Before you start 3 Using the Report Writer 4 General Concepts 4 Running the report writer 6 Creating a new Report 7 Opening
More informationWaves, Light, and Sound
CHAPTER 14 Waves, Light, and Sound LESSON 1 Waves What do you think? Read the two statements below and decide whether you agree or disagree with them. Place an A in the Before column if you agree with
More informationAn Analysis of Reducing Pedestrian-Walking-Speed Impacts on Intersection Traffic MOEs
An Analysis of Reducing Pedestrian-Walking-Speed Impacts on Intersection Traffic MOEs A Thesis Proposal By XIAOHAN LI Submitted to the Office of Graduate Studies of Texas A&M University In partial fulfillment
More informationOcean Fishing Fleet Scheduling Path Optimization Model Research. Based On Improved Ant Colony Algorithm
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 205) Ocean Fishing Fleet Scheduling Path Optimization Model Research Based On Improved Ant Colony Algorithm Li Jia-lin2,
More informationGolf Ball Impact: Material Characterization and Transient Simulation
12 th International LS-DYNA Users Conference Blast/Impact(1) Golf Ball Impact: Material Characterization and Transient Simulation Xiaohu Liu, David Quinn and Jorgen Bergström Veryst Engineering, LLC 47A
More informationElectronic Structure Workshop, Spring 2017
Electronic Structure Workshop, Spring 2017 Weine Olovsson @Linköping University, Wednesday 29th March 2017 National Supercomputer Centre in Linköping Sweden Introduction Welcome! Presentations from SNIC
More informationPhysics 11. Unit 7 (Part 1) Wave Motion
Physics 11 Unit 7 (Part 1) Wave Motion 1. Introduction to wave Wave motion is a popular phenomenon that we observe often in our daily lives. For example, light waves, sound waves, radio waves, water waves,
More informationPaul Burkhardt. May 19, 2016
GraphEx Symposium 2016 U.S. National Security Agency Research Directorate May 19, 2016 Graphs are everywhere! Why graph? Graph-theoretic approaches are appealing to many fields simple data abstraction
More informationSolids, Liquids, and Gases
chapter 14 Solids, Liquids, and Gases section 3 Behavior of Gases What You ll Learn how a gas exerts pressure on its container how changing pressure, temperature, or volume affect a gas Before You Read
More informationA STUDY OF THE LOSSES AND INTERACTIONS BETWEEN ONE OR MORE BOW THRUSTERS AND A CATAMARAN HULL
A STUDY OF THE LOSSES AND INTERACTIONS BETWEEN ONE OR MORE BOW THRUSTERS AND A CATAMARAN HULL L Boddy and T Clarke, Austal Ships, Australia SUMMARY CFD analysis has been conducted on a 100m catamaran hull
More informationCOMPLEX NUMBERS. Powers of j. Consider the quadratic equation; It has no solutions in the real number system since. j 1 ie. j 2 1
COMPLEX NUMBERS Consider the quadratic equation; x 0 It has no solutions in the real number system since x or x j j ie. j Similarly x 6 0 gives x 6 or x 6 j Powers of j j j j j. j j j j.j jj j j 5 j j
More informationTitle: 4-Way-Stop Wait-Time Prediction Group members (1): David Held
Title: 4-Way-Stop Wait-Time Prediction Group members (1): David Held As part of my research in Sebastian Thrun's autonomous driving team, my goal is to predict the wait-time for a car at a 4-way intersection.
More informationA Novel Decode-Aware Compression Technique for Improved Compression and Decompression
A Novel Decode-Aware Compression Technique for Improved Compression and Decompression J. Suresh Babu, K. Tirumala Rao & P. Srinivas Department Of Electronics & Comm. Engineering, Nimra College of Engineering
More informationThe Estimation Of Compressor Performance Using A Theoretical Analysis Of The Gas Flow Through the Muffler Combined With Valve Motion
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering The Estimation Of Compressor Performance Using A Theoretical Analysis Of The Gas Flow Through
More informationWhat kind of gamer are you?
What kind of gamer are you? Guillaume DECOOPMAN - Louis PIOTROWSKI - Jean-Marie SEGUY ENAC May 20, 2015 G. DECOOPMAN - L. PIOTROWSKI - J.M. SEGUY What kind of gamer are you? 1 / 26 1 Introduction 2 3 4
More informationTraffic circles. February 9, 2009
Traffic circles February 9, 2009 Abstract The use of a traffic circle is a relatively common means of controlling traffic in an intersection. Smaller Traffic circles can be especially effective in routing
More informationMECHANICAL WAVES AND SOUND
MECHANICAL WAVES AND SOUND Waves Substances have a stable equilibrium state Uniform pressure everywhere throughout the substance Atomic springs are at their equilibrium length Can make a wave by disturbing
More informationBesides the reported poor performance of the candidates there were a number of mistakes observed on the assessment tool itself outlined as follows:
MATHEMATICS (309/1) REPORT The 2013 Mathematics (309/1) paper was of average standard. The paper covered a wide range of the syllabus. It was neither gender bias nor culture bias. It did not have language
More informationA HYBRID METHOD FOR CALIBRATION OF UNKNOWN PARTIALLY/FULLY CLOSED VALVES IN WATER DISTRIBUTION SYSTEMS ABSTRACT
A HYBRID METHOD FOR CALIBRATION OF UNKNOWN PARTIALLY/FULLY CLOSED VALVES IN WATER DISTRIBUTION SYSTEMS Nhu Cuong Do 1,2, Angus Simpson 3, Jochen Deuerlein 4, Olivier Piller 5 1 University of Saskatchewan,
More informationOpen Research Online The Open University s repository of research publications and other research outputs
Open Research Online The Open University s repository of research publications and other research outputs Developing an intelligent table tennis umpiring system Conference or Workshop Item How to cite:
More informationLUMAZOTE T ECHNICAL BRO CHU RE
by Chemiluminescence Standards AST D4629 / ASTM D6069 / ASTM D 5176 / ASTM D 5762 / NF EN 12260 / NF M 07-058 Licence TOTAL LUMAZOTE T ECHNICAL BRO CHU RE ERALY & Associés Phone: +33 (0) 1 77 04 80 97
More informationOrganize information about waves. Differentiate two main types of waves.
Lesson 1 Waves Scan Lesson 1. Read the lesson titles and bold words. Look at the pictures. Identify three facts you discovered about waves. Record your facts in your Science Journal. What are waves? Organize
More informationLOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12
LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/ This experiment will introduce you to the kinetic properties of low-pressure gases. You will make observations on the rates with which selected
More informationCover Sheet-Block 6 Wave Properties
Cover Sheet-Block 6 Wave Properties Name Standards-Physics 4 a b c d 4a. Students know waves carry energy from one place to another. 4. b. Students know how to identify transverse and longitudinal waves
More informationBayesian Optimized Random Forest for Movement Classification with Smartphones
Bayesian Optimized Random Forest for Movement Classification with Smartphones 1 2 3 4 Anonymous Author(s) Affiliation Address email 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
More informationA MODEL FOR ANALYSIS OF THE IMPACT BETWEEN A TENNIS RACKET AND A BALL
A MODEL FOR ANALYSIS OF THE IMPACT BETWEEN A TENNIS RACKET AND A BALL Hiroshi Maeda, Masaaki Okauchi and Yoshio Shimada Faculty of Engineering, Oita University, Oita City, Oita, Japan The purpose of this
More informationU S F O S B u o y a n c y And Hydrodynamic M a s s
1 U S F O S B u o y a n c y And Hydrodynamic M a s s 2 CONTENTS: 1 INTRODUCTION... 3 2 ACCURACY LEVELS... 3 2.1 LEVEL-0... 3 2.2 LEVEL-1... 3 2.3 PANEL MODEL... 3 3 EX 1. SINGLE PIPE. NON FLOODED... 4
More informationInternational Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.
DESIGN AND ANALYSIS OF FEED CHECK VALVE AS CONTROL VALVE USING CFD SOFTWARE R.Nikhil M.Tech Student Industrial & Production Engineering National Institute of Engineering Mysuru, Karnataka, India -570008
More informationOPTIMIZATION OF A WAVE CANCELLATION MULTIHULL SHIP USING CFD TOOLS
OPTIMIZATION OF A AVE CANCELLATION MULTIHULL SHIP USING CFD TOOLS C. Yang, R. Löhner and O. Soto School of Computational Sciences, George Mason University Fairfax VA 030-4444, USA ABSTRACT A simple CFD
More informationUnits of Chapter 14. Types of Waves Waves on a String Harmonic Wave Functions Sound Waves Standing Waves Sound Intensity The Doppler Effect
Units of Chapter 14 Types of Waves Waves on a String Harmonic Wave Functions Sound Waves Standing Waves Sound Intensity The Doppler Effect Units of Chapter 14 Optional Superposition and Interference Beats
More informationWind turbine Varying blade length with wind speed
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, PP 01-05 www.iosrjournals.org Wind turbine Varying blade length with wind speed Mohammed Ashique
More informationISOLATION OF NON-HYDROSTATIC REGIONS WITHIN A BASIN
ISOLATION OF NON-HYDROSTATIC REGIONS WITHIN A BASIN Bridget M. Wadzuk 1 (Member, ASCE) and Ben R. Hodges 2 (Member, ASCE) ABSTRACT Modeling of dynamic pressure appears necessary to achieve a more robust
More informationSoundCast Design Intro
SoundCast Design Intro Basic Design SoundCast and Daysim 3 Land use attributes Households & Individuals SoundCast DaySim Travel demand simulator Trips and Households, Excel Summary Sheets, EMME network
More informationFluids, Pressure and buoyancy
Fluids, Pressure and buoyancy Announcements: CAPA due Friday at 10pm. Comment on the hint in Problem 5. CAPA solutions from previous sets can be found by logging onto CAPA and selecting View Previous Set
More informationChapter 10. Physical Characteristics of Gases
Chapter 10 Physical Characteristics of Gases Kinetic Molecular Theory An understanding of the behavior of atoms that make up matter Ideal gas: an imaginary gas that perfectly fits all assumptions of the
More informationAP Physics 1 Summer Packet Review of Trigonometry used in Physics
AP Physics 1 Summer Packet Review of Trigonometry used in Physics For some of you this material will seem pretty familiar and you will complete it quickly. For others, you may not have had much or any
More informationOlympus Production Contour Plot
Olympus Production Contour Plot Overview An alternative to the Scale Scenario calculation is the Contour Diagram. The conditions defined in the Scale Scenario are specific to production locations (wellhead,
More informationANSWERS TO QUESTIONS IN THE NOTES AUTUMN 2018
ANSWERS TO QUESTIONS IN THE NOTES AUTUMN 2018 Section 1.2 Example. The discharge in a channel with bottom width 3 m is 12 m 3 s 1. If Manning s n is 0.013 m -1/3 s and the streamwise slope is 1 in 200,
More informationApplied Econometrics with. Time, Date, and Time Series Classes. Motivation. Extension 2. Motivation. Motivation
Applied Econometrics with Extension 2 Time, Date, and Time Series Classes Time, Date, and Time Series Classes Christian Kleiber, Achim Zeileis 2008 2017 Applied Econometrics with R Ext. 2 Time, Date, and
More informationPhysics Workbook WALCH PUBLISHING
Physics Workbook WALCH PUBLISHING Table of Contents To the Student.............................. vii Unit 1: Forces Activity 1 Distance and Displacement................ 1 Activity 2 Vector and Scalar Quantities...............
More informationSimulation of Arterial Traffic Using Cell Transmission Model
1 Simulation of Arterial Traffic Using Cell Transmission Model Dongyan Su University of California, Berkeley 2167c Etcheverry Hall, University of California, Berkeley, CA 94720-1740 Tel: 1-510-928-8771,
More informationStatic Fluids. **All simulations and videos required for this package can be found on my website, here:
DP Physics HL Static Fluids **All simulations and videos required for this package can be found on my website, here: http://ismackinsey.weebly.com/fluids-hl.html Fluids are substances that can flow, so
More informationChapter 16. Waves-I Types of Waves
Chapter 16 Waves-I 16.2 Types of Waves 1. Mechanical waves. These waves have two central features: They are governed by Newton s laws, and they can exist only within a material medium, such as water, air,
More informationPolynomial DC decompositions
Polynomial DC decompositions Georgina Hall Princeton, ORFE Joint work with Amir Ali Ahmadi Princeton, ORFE 7/31/16 DIMACS Distance geometry workshop 1 Difference of convex (dc) programming Problems of
More informationDynamic Programming: The Matrix Chain Algorithm
Dynamic Programming: The Matrix Chain Algorithm Andreas Klappenecker [partially based on slides by Prof. Welch] Matrix Chain Problem Suppose that we want to multiply a sequence of rectangular matrices.
More informationNote! In this lab when you measure, round all measurements to the nearest meter!
Distance and Displacement Lab Note! In this lab when you measure, round all measurements to the nearest meter! 1. Place a piece of tape where you will begin your walk outside. This tape marks the origin.
More informationVibration of floors and footfall analysis
Webinar Autodesk Robot Structural Analysis Professional 20/04/2016 Vibration of floors and footfall analysis Artur Kosakowski Rafał Gawęda Webinar summary In this webinar we will focus on the theoretical
More information2 When Some or All Labels are Missing: The EM Algorithm
CS769 Spring Advanced Natural Language Processing The EM Algorithm Lecturer: Xiaojin Zhu jerryzhu@cs.wisc.edu Given labeled examples (x, y ),..., (x l, y l ), one can build a classifier. If in addition
More information6. An oscillator makes four vibrations in one second. What is its period and frequency?
Period and Frequency 19.1 The period of a pendulum is the time it takes to move through one cycle. As the ball on the string is pulled to one side and then let go, the ball moves to the side opposite the
More informationITW RANSBURG REA AND VECTOR
ITW RANSBURG REA AND VECTOR APPLICATORS TECHNICAL MANUAL (Replaces IL-246-A) January - 2008 ELECTROSTATIC TIC SPRAY TECHNIQUES IMPORTANT ANT: : Before using this equipment, carefully read SAFETY PRECAUTIONS,
More informationSection 2 Multiphase Flow, Flowing Well Performance
Section 2 Multiphase Flow, Flowing Well Performance Multiphase Vertical Flow When producing an oil or gas well, the flow of the fluids up the tubing will be in most cases be 2 phase, liquid and gas. The
More informationDesigning a Traffic Circle By David Bosworth For MATH 714
Designing a Traffic Circle By David Bosworth For MATH 714 Abstract We all have had the experience of sitting in a traffic jam, or we ve seen cars bunched up on a road for some no apparent reason. The use
More informationMarch Madness Basketball Tournament
March Madness Basketball Tournament Math Project COMMON Core Aligned Decimals, Fractions, Percents, Probability, Rates, Algebra, Word Problems, and more! To Use: -Print out all the worksheets. -Introduce
More information