PHYS Tutorial 7: Random Walks & Monte Carlo Integration
|
|
- Abel Anderson
- 5 years ago
- Views:
Transcription
1 PHYS Tutorial 7: Random Walks & Monte Carlo Integration The goal of this tutorial is to model a random walk in two dimensions and observe a phase transition as parameters are varied. Additionally, some Monte Carlo integration exercises are included at the end for practice. 1. Zombie Apocalypse Rumour has it that mad scientists from UBC have been working on a dangerous virus that turns anyone it infects into a zombie. Alarmed, the City of Vancouver has contacted you because of your expertise in producing fast and reliable simulations. Their request is simple: they want you to model the potential consequences of a zombie invasion and to find ways to mitigate the damages it would cause to the population. The setup The simplest model of a zombie invasion is that of a two-dimensional random walk. For the sake of the argument, we will assume that both zombies and terrified humans walk randomly on a Cartesian lattice, and that the infection is transmitted when walkers find themselves on the same site as a zombie. We will also assume that the virus is so potent that zombies eventually die of organ failure, and thus stop infecting humans. An effective simulation therefore needs to keep track of the health status of all walkers on the lattice at each time step, updating it as necessary. In order to assess the consequences of a pandemic outbreak, many parameters have to be taken into account. Among those are the density of the population, the survival time of a zombie, the probability of transmitting the disease, and the fraction of people immune to the virus, to name only a few. Definitions Luckily for you, the officials at City Hall already have an incomplete version of a zombie invasion code. In it, they define the following parameters: a grid size N, such that the lattice possesses N N sites; a density variable that specifies the fraction of sites occupied by people; a time T which measures the number of iterations (hours) it takes for a zombie to die; a number of walker M, expressed in terms of N and density. 1
2 The advantage of vectorization Intuitively, it would make sense to code the zombie invasion with a series of embedded for and if loops that cover the realm of all possibilities for the walkers: do we take a step in x? in y? does the infection get transmitted? and so forth. However, coding that way is extremely inefficient. The City of Vancouver is in a state of emergency and time is a limited resource, so you need to find a way to get results quickly. Since MATLAB excels at manipulating matrices and vectors, the best solution is to minimize the quantity of loops, which act on individual elements, in favour of operations that act on arrays all at once. This technique is called vectorization. 1 Note: In a vectorized random walk, all walkers take random steps simultaneously. A subtle problem caused by this is that neighbours can never find themselves on the same site after a random step, which implies that zombies can never infect adjacent walkers. An easy fix is to infect the walkers who are on the same site as a zombie as well as those who step on a site where a zombie used to be at the previous iteration. Question 1: Initialization To initialize the random walk, the first step is to randomly assign positions to all M walkers. To keep track of the location of every walker, we can define two vectors x and y with the statement randi(n, M, 1), which creates a column vector whose M elements take random integer values between 1 and N. Thus, walker j has position {x(j), y(j)} on the lattice. The second step is to turn a random walker into a zombie. After this step, we need to be able to distinguish between healthy, infected, and dead walkers. To do so, let us initialize a vector, which we will call infect and whose elements take either the value 0 (healthy), 1 (zombie), or 2 (dead). Therefore, we would like infect to be 0 everywhere except at the position of the infected walker, where it takes the value 1. Remember that zombies can only wander around for T hours before dying, which means that we also have to keep track of the time that passes after every infection. An easy way to do so is through a vector hours initialized to the value T for everyone, but which decreases every iteration for all zombies. When an element hours(j) reaches 0, we will need to update infect(j) from 1 to 2 (where 2 means death; the walker stops moving altogether). Additionally, we will initialize three counters to keep track of casualties and time: ndeaths = 0, nzombies = 1, and iter = 1. 1 On my 2012 MacBook Pro, a simulation at high population density constructed of loops takes more than 200 seconds to run, compared to 1 second on the vectorized version of the same code. 2
3 Question 2: Random walk The random walk will take place in a while loop that keeps on running as long as the conditions nzombies > 0 and iter < Nitermax are met. To walk randomly, we need random directions. Convince yourself that the statement >> d e l t a x = 2 floor (2 rand) 1 ; where rand is a number randomly chosen in the open interval (0, 1), defines a variable deltax which exclusively takes the values -1 or 1. As mentioned earlier, the walkers in a vectorized random walk all take steps simultaneously. Look at the Random Walk part of the incomplete code, and add a section for the random walk in the y direction. Be sure to understand how to use array-wide statements like deltax(infect == 2) = 0, or what type of variable walkx and walky are, as knowing this is the key to eschew the use of loops. Notice how some problems might arise in xnew if walkers take steps outside of the N N lattice. Implement periodic boundary conditions in both the x and y directions to prevent walkers from leaving the simulation. Question 3: Infection and death counts There are a few tasks your code needs to do before each random step: If a walker is infected, subtract 1 from its hours counter; Update the status of all zombies who die on this iteration; Update the counters nzombies and ndeaths. Do not use if loops! Use logical operations instead, as in Question 2. For instance, the number of dead walkers is simply ndeaths = sum(infect == 2) in vectorized notation. Question 4: Contamination The contamination step occurs on two separate occasions: a healthy walker and a zombie land on the same site; a walker steps on a site that a zombie just left. We need to compare the x and y coordinates of the healthy walkers to those of the infected ones and, if there are matches, proceed with contamination. The statement C = ismember(a, B) compares the arrays A and B and outputs a logical array C taking the value 1 for all elements of A also found in B, and 0 for the unique ones. 3
4 Therefore, all that you need to do is to define a variable called I coord which contains the x and y coordinates (old and new) of the infected zombies. Also, make sure to understand why H coord has been initialized the way it has. Question 5: Additional modifications Make sure that your code works properly. When it does, add these two features to make it more realistic: There is a possibility that some people are naturally immune to the virus. Assume that a fraction pimmune of the population (chosen at random) cannot be contaminated, and modify the initialization part of your code to take this into account. Hint: in the case of immunization, inf ect should take an additional value. Some viruses are more contagious than others. Add a probability of transmission that determines whether contamination occurs or not. The variable infection status should be modified accordingly. Hint: The health status of walkers corresponds to the binary true/infected (1) and false/healthy (0). Therefore, a simple way to implement contagion is to make infection status a logical array obtained from a logical comparison to a transmission probability. If k walkers have been in contact with zombies after some iteration, then infection status should be a logical vector of length k where each element determines whether one of these walker randomly becomes a zombie or not. Question 6: Critical phenomena Now that you have a working model of a zombie invasion, you can answer some questions that the mayor Gregor Robertson has for you: Given a density of 0.5 and a transmission probability of 0.8, what is the minimum value of pimmune needed so that most people stay alive? What if the density is 0.25? Given a density of 0.5 and assuming that 30% of the population is immune against the virus, how contagious does the virus need to be in order to become pandemic? How would these results change if the zombies died after 12 hours (instead of 24)? You should observe critical phenomena as you play with various parameters. Indeed, there seems to be a phase transition between the two population states mostly alive and mostly dead. (1) Given the random nature of your simulation, you should run your code many times with the same parameters and average your results. How would you estimate the error associated with your predictions? 4
5 2. Monte Carlo Integration Question 1: Estimating π using Monte Carlo Explain why the following lines calculate π: >> k = 1 e7 ; >> p i e s t = 4 sum(sum(rand(k,2).ˆ2,2) <=1)/ k ; Question 2: Hyperspheres Monte Carlo integration shines the most in higher dimensional applications. Use the logic of Question 1 to calculate the volume inside S d, the d-sphere, defined by for d = 2, i.e. for a ball in R 3 ; for d = 5. d+1 x 2 i = 1, (2) i=1 Question 3: Importance sampling Compute the integral I = by using Monte Carlo integration with Uniform sampling; 1 0 ( x 1/3 + x ) dx (3) 10 Importance sampling with the sampling function w(x) = 2 3 x 1/3. The actual answer is Plot the error for k going from 10 2 to 10 6 on a loglog plot; what do you observe? Question 4: Importance sampling (again) Repeat the above exercice with the integral I = 1 0 e x dx (4) with the sampling function w(x) = 1 + x, since e x 1 + x on the interval (0, 1). 5
ROSE-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 informationCalculation of Trail Usage from Counter Data
1. Introduction 1 Calculation of Trail Usage from Counter Data 1/17/17 Stephen Martin, Ph.D. Automatic counters are used on trails to measure how many people are using the trail. A fundamental question
More informationORF 201 Computer Methods in Problem Solving. Final Project: Dynamic Programming Optimal Sailing Strategies
Princeton University Department of Operations Research and Financial Engineering ORF 201 Computer Methods in Problem Solving Final Project: Dynamic Programming Optimal Sailing Strategies Due 11:59 pm,
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 informationTaking Your Class for a Walk, Randomly
Taking Your Class for a Walk, Randomly Daniel Kaplan Macalester College Oct. 27, 2009 Overview of the Activity You are going to turn your students into an ensemble of random walkers. They will start at
More informationEpidemics and zombies
Epidemics and zombies (Sethna, "Entropy, Order Parameters, and Complexity", ex. 6.21) 2018, James Sethna, all rights reserved. This exercise is based on Alemi and Bierbaum's class project, published in
More informationLab 4: Root Locus Based Control Design
Lab 4: Root Locus Based Control Design References: Franklin, Powell and Emami-Naeini. Feedback Control of Dynamic Systems, 3 rd ed. Addison-Wesley, Massachusetts: 1994. Ogata, Katsuhiko. Modern Control
More informationUNIVERSITY OF WATERLOO
UNIVERSITY OF WATERLOO Department of Chemical Engineering ChE 524 Process Control Laboratory Instruction Manual January, 2001 Revised: May, 2009 1 Experiment # 2 - Double Pipe Heat Exchanger Experimental
More informationSoftware Reliability 1
Software Reliability 1 Software Reliability What is software reliability? the probability of failure-free software operation for a specified period of time in a specified environment input sw output We
More informationOpleiding Informatica
Opleiding Informatica Determining Good Tactics for a Football Game using Raw Positional Data Davey Verhoef Supervisors: Arno Knobbe Rens Meerhoff BACHELOR THESIS Leiden Institute of Advanced Computer Science
More informationA Hare-Lynx Simulation Model
1 A Hare- Simulation Model What happens to the numbers of hares and lynx when the core of the system is like this? Hares O Balance? S H_Births Hares H_Fertility Area KillsPerHead Fertility Births Figure
More informationAssignment #3 Breakout! Due: 12pm on Wednesday, February 6th This assignment should be done individually (not in pairs)
Chris Piech Assn #3 CS 106A January 28, 2019 Assignment #3 Breakout! Due: 12pm on Wednesday, February 6th This assignment should be done individually (not in pairs) Based on a handout by Eric Roberts with
More informationAssignment #3 Breakout!
Eric Roberts Handout #25 CS 106A January 22, 2010 Assignment #3 Breakout! Due: Wednesday, February 3, 5:00P.M. Your job in this assignment is to write the classic arcade game of Breakout, which was invented
More information1.2 Example 1: A simple hydraulic system
Note: It is possible to use more than one fluid in the Hydraulic library. This is important because you can model combined cooling and lubrication systems of a library. The hydraulic library assumes a
More informationIteration: while, for, do while, Reading Input with Sentinels and User-defined Functions
Iteration: while, for, do while, Reading Input with Sentinels and User-defined Functions This programming assignment uses many of the ideas presented in sections 6 and 7 of the course notes. You are advised
More informationESTIMATION OF THE DESIGN WIND SPEED BASED ON
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 2009, Taipei, Taiwan ESTIMATION OF THE DESIGN WIND SPEED BASED ON UNCERTAIN PARAMETERS OF THE WIND CLIMATE Michael Kasperski 1 1
More informationAPPENDIX A COMPUTATIONALLY GENERATED RANDOM DIGITS 748 APPENDIX C CHI-SQUARE RIGHT-HAND TAIL PROBABILITIES 754
IV Appendices APPENDIX A COMPUTATIONALLY GENERATED RANDOM DIGITS 748 APPENDIX B RANDOM NUMBER TABLES 750 APPENDIX C CHI-SQUARE RIGHT-HAND TAIL PROBABILITIES 754 APPENDIX D LINEAR INTERPOLATION 755 APPENDIX
More informationAssignment 3: Breakout!
CS106A Winter 2011-2012 Handout #16 February 1, 2011 Assignment 3: Breakout! Based on a handout by Eric Roberts and Mehran Sahami Your job in this assignment is to write the classic arcade game of Breakout,
More informationDeer Population Student Guide
Deer Population Student Guide In many places, deer have become nuisance animals because they are so numerous. In some areas, a hunting season has been introduced or lengthened to reduce the number of deer.
More informationBased on a handout by Eric Roberts
Mehran Sahami Handout #19 CS 106A October 15, 2018 Assignment #3 Breakout! Due: 1:30pm on Wednesday, October 24th This assignment should be done individually (not in pairs) Your Early Assignment Help (YEAH)
More informationDesign of Experiments Example: A Two-Way Split-Plot Experiment
Design of Experiments Example: A Two-Way Split-Plot Experiment A two-way split-plot (also known as strip-plot or split-block) design consists of two split-plot components. In industry, these designs arise
More informationWhen zombies attack! The maths of infection rates
TEACHERS RESOURCE MEOW Epidemiology When zombies attack! The maths of infection rates Game 1 1. Give each student a sheet of paper with a letter at the top of a column with plenty of room below to keep
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 informationExcel Solver Case: Beach Town Lifeguard Scheduling
130 Gebauer/Matthews: MIS 213 Hands-on Tutorials and Cases, Spring 2015 Excel Solver Case: Beach Town Lifeguard Scheduling Purpose: Optimization under constraints. A. GETTING STARTED All Excel projects
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 informationLaser Tag Pro: Battle Rifle
Laser Tag Pro: Battle Rifle Congratulations and thank you for your purchase of the Battle Rifle Pro Bundle! We hope you are excited to take advantage of everything this unit has to offer. Whether for business
More informationUnderstanding safety life cycles
Understanding safety life cycles IEC/EN 61508 is the basis for the specification, design, and operation of safety instrumented systems (SIS) Fast Forward: IEC/EN 61508 standards need to be implemented
More informationThree-Dimensional Ray-Cast Pong
Three-Dimensional Ray-Cast Pong Richard Hughes Elizabeth Power! Overview What is Pong? Traditional Pong is a two-dimensional game that simulates table tennis. The player controls a paddle by moving it
More information5.1 Introduction. Learning Objectives
Learning Objectives 5.1 Introduction Statistical Process Control (SPC): SPC is a powerful collection of problem-solving tools useful in achieving process stability and improving capability through the
More informationCHEMICAL ENGINEERING LABORATORY CHEG 239W. Control of a Steam-Heated Mixing Tank with a Pneumatic Process Controller
CHEMICAL ENGINEERING LABORATORY CHEG 239W Control of a Steam-Heated Mixing Tank with a Pneumatic Process Controller Objective The experiment involves tuning a commercial process controller for temperature
More informationVerification and Validation Pathfinder Release 0730 x64
403 Poyntz Avenue, Suite B Manhattan, KS 66502 USA +1.785.770.8511 www.thunderheadeng.com Verification and Validation Pathfinder 2014.2 Release 0730 x64 Disclaimer Thunderhead Engineering makes no warranty,
More information2600T Series Pressure Transmitters Plugged Impulse Line Detection Diagnostic. Pressure Measurement Engineered solutions for all applications
Application Description AG/266PILD-EN Rev. C 2600T Series Pressure Transmitters Plugged Impulse Line Detection Diagnostic Pressure Measurement Engineered solutions for all applications Increase plant productivity
More informationChapter 14 Coal Boiler Flue Gas Scrubber II
Chapter 14 Coal Boiler Flue Gas Scrubber II The Application The following case study is a continuation of the base case scenario created previous chapter. The base case scenario contained four process
More informationApplication of the probabilistic-fuzzy method of assessing the risk of a ship manoeuvre in a restricted area
Application of the probabilistic-fuzzy method of assessing the risk of a ship manoeuvre in a restricted area L. ~ucrna', Z. pietrzykowski2 Maritime University of Szczecin ul. Waly Chrobrego 1/2 70-500
More informationECO 199 GAMES OF STRATEGY Spring Term 2004 Precept Materials for Week 3 February 16, 17
ECO 199 GAMES OF STRATEGY Spring Term 2004 Precept Materials for Week 3 February 16, 17 Illustration of Rollback in a Decision Problem, and Dynamic Games of Competition Here we discuss an example whose
More informationQueue analysis for the toll station of the Öresund fixed link. Pontus Matstoms *
Queue analysis for the toll station of the Öresund fixed link Pontus Matstoms * Abstract A new simulation model for queue and capacity analysis of a toll station is presented. The model and its software
More informationThe risk assessment of ships manoeuvring on the waterways based on generalised simulation data
Safety and Security Engineering II 411 The risk assessment of ships manoeuvring on the waterways based on generalised simulation data L. Gucma Maritime University of Szczecin, Poland Abstract This paper
More informationMIL-STD-883G METHOD
STEADY-STATE LIFE 1. PURPOSE. The steady-state life test is performed for the purpose of demonstrating the quality or reliability of devices subjected to the specified conditions over an extended time
More information4. Please Do Break the Crystal
4. Please Do Break the Crystal Tell the broken plate you are sorry. Mary Robertson. Programming constructs and algorithmic paradigms covered in this puzzle: Break statements, radix representations. You
More informationSection 1: Multiple Choice Explained EXAMPLE
CFSP Process Applications Section 1: Multiple Choice Explained EXAMPLE Candidate Exam Number (No Name): Please write down your name in the above provided space. Only one answer is correct. Please circle
More informationintended velocity ( u k arm movements
Fig. A Complete Brain-Machine Interface B Human Subjects Closed-Loop Simulator ensemble action potentials (n k ) ensemble action potentials (n k ) primary motor cortex simulated primary motor cortex neuroprosthetic
More informationLab # 03: Visualization of Shock Waves by using Schlieren Technique
AerE545 Lab # 03: Visualization of Shock Waves by using Schlieren Technique Objectives: 1. To get hands-on experiences about Schlieren technique for flow visualization. 2. To learn how to do the optics
More informationCFD AND EXPERIMENTAL STUDY OF AERODYNAMIC DEGRADATION OF ICED AIRFOILS
Colloquium FLUID DYNAMICS 2008 Institute of Thermomechanics AS CR, v.v.i., Prague, October 22-24, 2008 p.1 CFD AND EXPERIMENTAL STUDY OF AERODYNAMIC DEGRADATION OF ICED AIRFOILS Vladimír Horák 1, Dalibor
More informationBackground Summary Kaibab Plateau: Source: Kormondy, E. J. (1996). Concepts of Ecology. Englewood Cliffs, NJ: Prentice-Hall. p.96.
Assignment #1: Policy Analysis for the Kaibab Plateau Background Summary Kaibab Plateau: Source: Kormondy, E. J. (1996). Concepts of Ecology. Englewood Cliffs, NJ: Prentice-Hall. p.96. Prior to 1907, the
More informationVerification and Validation Pathfinder
403 Poyntz Avenue, Suite B Manhattan, KS 66502 USA +1.785.770.8511 www.thunderheadeng.com Verification and Validation Pathfinder 2015.1 Release 0504 x64 Disclaimer Thunderhead Engineering makes no warranty,
More informationLecture 1: Knot Theory
Math 7H Professor: Padraic Bartlett Lecture 1: Knot Theory Week 1 UCSB 015 1 Introduction Outside of mathematics, knots are ways to loop a single piece of string around itself: In mathematics, we mean
More informationBlack Box testing Exercises. Lecturer: Giuseppe Santucci
Black Box testing Exercises Lecturer: Giuseppe Santucci Exercise 1 A software program for the simulation of scuba diving with special gas mixture gives indication about maximum time of stay on the bottom,
More informationChapter 12 Practice Test
Chapter 12 Practice Test 1. Which of the following is not one of the conditions that must be satisfied in order to perform inference about the slope of a least-squares regression line? (a) For each value
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 informationSafety Critical Systems
Safety Critical Systems Mostly from: Douglass, Doing Hard Time, developing Real-Time Systems with UML, Objects, Frameworks And Patterns, Addison-Wesley. ISBN 0-201-49837-5 1 Definitions channel a set of
More informationZombies and Differential Equations
May 11, 2015 Topics Covered Basic Predator-Prey Equation Basic Equation Adjusted for Zombies Basic Equation with Latency The Quarantine Model The Treatment Equation Intelligent Human Resistance Introduction
More informationLegendre et al Appendices and Supplements, p. 1
Legendre et al. 2010 Appendices and Supplements, p. 1 Appendices and Supplement to: Legendre, P., M. De Cáceres, and D. Borcard. 2010. Community surveys through space and time: testing the space-time interaction
More informationSteam-Boiler Control Specification Problem
Steam-Boiler Control Specification Problem Jean-Raymond Abrial 1 Problem Statement 1.1 ntroduction This text constitutes an informal specification of a program which serves to control the level of water
More informationAppendix 12 Consequences
Appendix 12 Consequences Overview All of the consequence data used in the risk analysis was provided by the IPET Consequence team whose work is reported in Volume VII of the IPET report. The following
More information1 SE/P-02. Experimental and Analytical Studies on Thermal-Hydraulic Performance of a Vacuum Vessel Pressure Suppression System in ITER
1 SE/P-2 Experimental and Analytical Studies on Thermal-Hydraulic Performance of a Vacuum Vessel Pressure Suppression System in ITER K. Takase 1), H. Akimoto 1) 1) Japan Atomic Energy Research Institute,
More informationE. Agu, M. Kasperski Ruhr-University Bochum Department of Civil and Environmental Engineering Sciences
EACWE 5 Florence, Italy 19 th 23 rd July 29 Flying Sphere image Museo Ideale L. Da Vinci Chasing gust fronts - wind measurements at the airport Munich, Germany E. Agu, M. Kasperski Ruhr-University Bochum
More informationNew generation of Electronic Card Systems: The 4-D Card
New generation of Electronic Card Systems: The 4-D Card Capt. Alain Richard Project director Maritime Innovation arichard@imar.ca Context In the present context of globalization and competitiveness, it
More informationRewriting the Zombie Apocalypse The Day the Zombies Came
Rewriting the Zombie Apocalypse The Day the Zombies Came 826 National, Inc. All rights reserved. It Was Just a Normal Day What is a normal week like for you? Who are the people you come in contact with?
More informationSaturated-Unsaturated Consolidation
Introduction Saturated-Unsaturated Consolidation This example looks at the volume change behavior of a saturated-unsaturated column under loading, wetting and drying conditions. Feature Highlights GeoStudio
More informationPrior Knowledge: Students should have an understanding that plants and animals compete for resources such as food, space, water, air and shelter.
Science Lesson Plan Form Teacher: 4 th Grade Lesson: Predator/Prey SPI: 2.1 Science Goal: Recognize the impact of predation and competition on an ecosystem. What is the big idea of this standard? All life
More informationDetermining Occurrence in FMEA Using Hazard Function
Determining Occurrence in FMEA Using Hazard Function Hazem J. Smadi Abstract FMEA has been used for several years and proved its efficiency for system s risk analysis due to failures. Risk priority number
More informationRESEARCH OPPURTINITIES ON AIRCRAFT EMERGENCY EVACUATION. Presented by: Dr. Minesh POUDEL
RESEARCH OPPURTINITIES ON AIRCRAFT EMERGENCY EVACUATION Presented by: Dr. Minesh POUDEL Table of Content 1. Introduction and Objective 2. Aircraft Emergency Evacuation: Analysis of main Issues and regulations
More informationAmmonia Synthesis with Aspen Plus V8.0
Ammonia Synthesis with Aspen Plus V8.0 Part 2 Closed Loop Simulation of Ammonia Synthesis 1. Lesson Objectives Review Aspen Plus convergence methods Build upon the open loop Ammonia Synthesis process simulation
More informationA Research on the Airflow Efficiency Analysis according to the Variation of the Geometry Tolerance of the Sirocco Fan Cut-off for Air Purifier
A Research on the Airflow Efficiency Analysis according to the Variation of the Geometry Tolerance of the Sirocco Fan Cut-off for Air Purifier Jeon-gi Lee*, Choul-jun Choi*, Nam-su Kwak*, Su-sang Park*
More informationA SEMI-PRESSURE-DRIVEN APPROACH TO RELIABILITY ASSESSMENT OF WATER DISTRIBUTION NETWORKS
A SEMI-PRESSURE-DRIVEN APPROACH TO RELIABILITY ASSESSMENT OF WATER DISTRIBUTION NETWORKS S. S. OZGER PhD Student, Dept. of Civil and Envir. Engrg., Arizona State Univ., 85287, Tempe, AZ, US Phone: +1-480-965-3589
More informationARTIFICIAL NEURAL NETWORK BASED DESIGN FOR DUAL LATERAL WELL APPLICATIONS
The Pennsylvania State University the Graduate School Department of Energy and Mineral Engineering ARTIFICIAL NEURAL NETWORK BASED DESIGN FOR DUAL LATERAL WELL APPLICATIONS Thesis in Energy and Mineral
More informationAdvanced LOPA Topics
11 Advanced LOPA Topics 11.1. Purpose The purpose of this chapter is to discuss more complex methods for using the LOPA technique. It is intended for analysts who are competent with applying the basic
More informationThe benefits of the extended diagnostics feature. Compact, well-proven, and flexible
ABB MEASUREMENT & ANALYTICS TECHNICAL INFORMATION PositionMaster EDP300 Extended Diagnostics Compact, well-proven, and flexible The benefits of the extended diagnostics feature The PositionMaster EDP300
More informationNumerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water
Numerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water Hong Xu, Chokri Guetari ANSYS INC. Abstract Transient numerical simulations of the rise of a train of gas bubbles in a liquid
More informationScience Skills Station
Science Skills Station Objective 1. Interpret and analyze data so to determine the relationship between resource availability and carrying capacity of a population. 2. Identify biotic and abiotic factors
More informationAgilent Dimension Software for ELSD User Manual
Agilent Dimension Software for ELSD User Manual Agilent Dimension Software for ELSD User Manual Agilent Technologies Notices Agilent Technologies, Inc. 2011 No part of this manual may be reproduced in
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 informationCVEN Computer Applications in Engineering and Construction. Programming Assignment #4 Analysis of Wave Data Using Root-Finding Methods
CVEN 30-501 Computer Applications in Engineering and Construction Programming Assignment #4 Analysis of Wave Data Using Root-Finding Methods Date distributed: 9/30/016 Date due: 10/14/016 at 3:00 PM (electronic
More informationSection 1: Multiple Choice
CFSP Process Applications Section 1: Multiple Choice EXAMPLE Candidate Exam Number (No Name): Please write down your name in the above provided space. Only one answer is correct. Please circle only the
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 informationFUBA RULEBOOK VERSION
FUBA RULEBOOK 15.07.2018 2.1 VERSION Introduction FUBA is a board game, which simulates football matches from a tactical view. The game includes its most important details. The players take roles of the
More informationRULES AND REGULATIONS OF FIXED ODDS BETTING GAMES
RULES AND REGULATIONS OF FIXED ODDS BETTING GAMES Royalhighgate Public Company Ltd. 04.04.2014 Table of contents SECTION I: GENERAL RULES... 6 ARTICLE 1 - GENERAL REGULATIONS... 6 ARTICLE 2 - THE HOLDING
More informationEVOLVING HEXAPOD GAITS USING A CYCLIC GENETIC ALGORITHM
Evolving Hexapod Gaits Using a Cyclic Genetic Algorithm Page 1 of 7 EVOLVING HEXAPOD GAITS USING A CYCLIC GENETIC ALGORITHM GARY B. PARKER, DAVID W. BRAUN, AND INGO CYLIAX Department of Computer Science
More informationPolicy Gradient RL to learn fast walk
Policy Gradient RL to learn fast walk Goal: Enable an Aibo to walk as fast as possible Policy Gradient RL to learn fast walk Goal: Enable an Aibo to walk as fast as possible Start with a parameterized
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 informationThe system design must obey these constraints. The system is to have the minimum cost (capital plus operating) while meeting the constraints.
Computer Algorithms in Systems Engineering Spring 2010 Problem Set 6: Building ventilation design (dynamic programming) Due: 12 noon, Wednesday, April 21, 2010 Problem statement Buildings require exhaust
More informationINTRODUCTION Page 2 SETTINGS
Page 1 of 11 SCRIMMAGE RULES AND INSTRUCTIONS INTRODUCTION ------------------------------------------------- Page 2 SETTINGS OBJECTIVE OF THE GAME --------------------------------------- Page 3 THE FIELD
More informationGlobal Journal of Engineering Science and Research Management
SIMULATION AND OPTIMZING TRAFFIC FLOW AT SIGNALIZED INTERSECTION USING MATLAB Dr Mohammed B. Abduljabbar*, Dr Amal Ali, Ruaa Hameed * Assist Prof., Civil Engineering Department, Al-Mustansiriayah University,
More informationCS Problem Solving and Object-Oriented Programming Lab 2 - Methods, Variables and Functions in Alice Due: September 23/24
CS 101 - Problem Solving and Object-Oriented Programming Lab 2 - Methods, Variables and Functions in Alice Due: September 23/24 Pre-lab Preparation Before coming to lab, you are expected to have: Read
More informationIntroduction to Pattern Recognition
Introduction to Pattern Recognition Jason Corso SUNY at Buffalo 12 January 2009 J. Corso (SUNY at Buffalo) Introduction to Pattern Recognition 12 January 2009 1 / 28 Pattern Recognition By Example Example:
More informationSafety Manual. Process pressure transmitter IPT-1* 4 20 ma/hart. Process pressure transmitter IPT-1*
Safety Manual Process pressure transmitter IPT-1* 4 20 ma/hart Process pressure transmitter IPT-1* Contents Contents 1 Functional safety 1.1 General information... 3 1.2 Planning... 4 1.3 Instrument parameter
More informationInfrared Thermal Imaging, Inc.
www.itimaging.com Infrared Inspection Data and Addressing Discrepancies in A Fired Heater Program By Ty Keeth Infrared inspection of fired heaters has quickly become the standard for fired heater inspection.
More informationINSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation
Lab 1 Standing Waves on a String Learning Goals: To distinguish between traveling and standing waves To recognize how the wavelength of a standing wave is measured To recognize the necessary conditions
More informationBlack Sea Bass Encounter
Black Sea Bass Encounter Below is an adaptation of the Shark Encounter (Lawrence Hall of Science: MARE 2002) lesson plan to be about Black Sea Bass and to incorporate information learned from Dr. Jensen
More informationChapter 20. Planning Accelerated Life Tests. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 20 Planning Accelerated Life Tests William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationEnergy capture performance
Energy capture performance Cost of energy is a critical factor to the success of marine renewables, in order for marine renewables to compete with other forms of renewable and fossil-fuelled power generation.
More informationCarrying Capacity Activity. 5 th Grade PSI. Teacher s Notes Procedure: Simulation 1 Regular herds
Carrying Capacity Activity 5 th Grade PSI Teacher s Notes Procedure: Simulation 1 Regular herds 1. Before the start of the activity, count out 4 beans per student. Place the beans in the center of the
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 informationAmendments to the International Convention on Maritime Search and Rescue of 27 April 1979
Downloaded on July 27, 2018 Amendments to the International Convention on Maritime Search and Rescue of 27 April 1979 Region United Nations (UN) Subject Maritime Sub Subject Type Amendments Reference Number
More informationGuided Study Program in System Dynamics System Dynamics in Education Project System Dynamics Group MIT Sloan School of Management 1
Guided Study Program in System Dynamics System Dynamics in Education Project System Dynamics Group MIT Sloan School of Management 1 Solutions to Assignment #17 Friday, March 12, 1999 Reading Assignment:
More informationUsing MATLAB with CANoe
Version 2.0 2017-03-09 Application Note AN-IND-1-007 Author Restrictions Abstract Vector Informatik GmbH Public Document This application note describes the usage of MATLAB /Simulink combined with CANoe.
More informationRanger Walking Initiation Stephanie Schneider 5/15/2012 Final Report for Cornell Ranger Research
1 Ranger Walking Initiation Stephanie Schneider sns74@cornell.edu 5/15/2012 Final Report for Cornell Ranger Research Abstract I joined the Biorobotics Lab this semester to gain experience with an application
More informationDigital Level Control One and Two Loops Proportional and Integral Control Single-Loop and Cascade Control
Digital Level Control One and Two Loops Proportional and Integral Control Single-Loop and Cascade Control Introduction This experiment offers a look into the broad field of process control. This area of
More informationPURE SUBSTANCE. Nitrogen and gaseous air are pure substances.
CLASS Third Units PURE SUBSTANCE Pure substance: A substance that has a fixed chemical composition throughout. Air is a mixture of several gases, but it is considered to be a pure substance. Nitrogen and
More informationCS 351 Design of Large Programs Zombie House
CS 351 Design of Large Programs Zombie House Instructor: Joel Castellanos e-mail: joel@unm.edu Web: http://cs.unm.edu/~joel/ Office: Electrical and Computer Engineering building (ECE). Room 233 2/23/2017
More information