Walking with coffee: when and why coffee spills

Similar documents
Irrigation &Hydraulics Department lb / ft to kg/lit.

In memory of Dr. Kevin P. Granata, my graduate advisor, who was killed protecting others on the morning of April 16, 2007.

Dynamic Component of Ship s Heeling Moment due to Sloshing vs. IMO IS-Code Recommendations

Sample Solution for Problem 1.a

How Do You Swing? You should be working with new lab partners starting with this lab.

Anatomy of a Homer. Purpose. Required Equipment/Supplies. Optional Equipment/Supplies. Discussion

Wave phenomena in a ripple tank

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS

Experimental study on path instability of rising bubbles

Finite Element Analysis on Trapezoidal Tank to Suppress Sloshing Effect

Quantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges

Acoustical approach to analysis of energy conversions in an oscillating bubble

COURSE NUMBER: ME 321 Fluid Mechanics I Fluid statics. Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET

Development of Technology to Estimate the Flow Field around Ship Hull Considering Wave Making and Propeller Rotating Effects

Putting Report Details: Key and Diagrams: This section provides a visual diagram of the. information is saved in the client s database

Numerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water

Current issues regarding induced acceleration analysis of walking using the integration method to decompose the GRF

b. What is the x-distance from the foot of the cliff to the point of impact in the lake?

ZIN Technologies PHi Engineering Support. PHi-RPT CFD Analysis of Large Bubble Mixing. June 26, 2006

G.C.E (A/L) Examination March In collaboration with

Bioreactor System ERT 314. Sidang /2011

1. A tendency to roll or heel when turning (a known and typically constant disturbance) 2. Motion induced by surface waves of certain frequencies.

PHYS 101 Previous Exam Problems

Body Stabilization of PDW toward Humanoid Walking

Vapour pressure of liquids SURFACE TENSION

. In an elevator accelerating upward (A) both the elevator accelerating upward (B) the first is equations are valid

Unsteady airfoil experiments

Pre AP Physics: Unit 7 Vibrations, Waves, and Sound. Clear Creek High School

The Usage of Propeller Tunnels For Higher Efficiency and Lower Vibration. M. Burak Şamşul


Optimization of an off-road bicycle with four-bar linkage rear suspension

Activity P07: Acceleration of a Cart (Acceleration Sensor, Motion Sensor)

Transverse waves cause particles to vibrate perpendicularly to the direction of the wave's motion (e.g. waves on a string, ripples on a pond).

What is a wave? Even here the wave more or less keeps it s shape and travelled at a constant speed. YouTube. mexicanwave.mov

Stability and Flight Controls

Wind tunnel effects on wingtip vortices

Figure 1 Schematic of opposing air bearing concept

Naval Postgraduate School, Operational Oceanography and Meteorology. Since inputs from UDAS are continuously used in projects at the Naval

Ball Toss. Vernier Motion Detector

Two dimensional kinematics. Projectile Motion

The Starting Point. Prosthetic Alignment in the Transtibial Amputee. Outline. COM Motion in the Coronal Plane

Humanoid Robots and biped locomotion. Contact: Egidio Falotico

The Discussion of this exercise covers the following points: Pumps Basic operation of a liquid pump Types of liquid pumps The centrifugal pump.

Wind Flow Model of Area Surrounding the Case Western Reserve University Wind Turbine

A NOVEL FLOATING OFFSHORE WIND TURBINE CONCEPT: NEW DEVELOPMENTS

Conditions for occurrence of vortex shedding on a large cable stayed bridge. Full scale data from monitoring system

Waves. harmonic wave wave equation one dimensional wave equation principle of wave fronts plane waves law of reflection

Toward a Human-like Biped Robot with Compliant Legs

INSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation

Flight Corridor. The speed-altitude band where flight sustained by aerodynamic forces is technically possible is called the flight corridor.

Proceedings of the International Conference on Emerging Trends in Engineering and Management (ICETEM14) 30-31,December, 2014, Ernakulam, India

Flutter Testing. Wind Tunnel Testing (excerpts from Reference 1)

EFFECT OF CORNER CUTOFFS ON FLOW CHARACTERISTICS AROUND A SQUARE CYLINDER

PHYSICS 12 NAME: Kinematics and Projectiles Review

Garrett McNamara, Portugal, 30 Jan What is a wave?

ITTC Recommended Procedures and Guidelines

Computer Simulation Helps Improve Vertical Column Induced Gas Flotation (IGF) System

Bubble dynamics in drinks

Dynamically stepping over large obstacle utilizing PSO optimization in the B4LC system

Advanced Applications in Naval Architecture Beyond the Prescriptions in Class Society Rules

A numerical Euler-Lagrange method for bubble tower CO2 dissolution modeling

Research on Small Wind Power System Based on H-type Vertical Wind Turbine Rong-Qiang GUAN a, Jing YU b

Core Concept. PowerPoint Lectures Physical Science, 8e. Chapter 5 Wave Motions and Sound. New Symbols for this Chapter 2/20/2011

Research on Goods and the Ship Interaction Based on ADAMS

and its weight (in newtons) when located on a planet with an acceleration of gravity equal to 4.0 ft/s 2.

Biomechanics and Models of Locomotion

Positioning and Control of Boom Crane Luffing with Double-Pendulum Payloads

AP Physics B Summer Homework (Show work)

a. Determine the sprinter's constant acceleration during the first 2 seconds. b. Determine the sprinters velocity after 2 seconds have elapsed.

Lesson 5. Section 2.2: Trigonometric Functions of an Acute Angle 1 = 1

REPORT. A comparative study of the mechanical and biomechanical behaviour of natural turf and hybrid turf for the practise of sports

Proceedings of Meetings on Acoustics

Aerodynamic Analysis of a Symmetric Aerofoil

Influence of rounding corners on unsteady flow and heat transfer around a square cylinder

Analysis of Movement

Gravity wave effects on the calibration uncertainty of hydrometric current meters

Kinematic Differences between Set- and Jump-Shot Motions in Basketball

A NEW GOLF-SWING ROBOT MODEL UTILIZING SHAFT ELASTICITY

Workshop 1: Bubbly Flow in a Rectangular Bubble Column. Multiphase Flow Modeling In ANSYS CFX Release ANSYS, Inc. WS1-1 Release 14.

Experimental Studies on the Instabilities of Viscous Fingering in a Hele-Shaw Cell

Wind Flow Validation Summary

Chapter 15 Wave Motion. Copyright 2009 Pearson Education, Inc.

Sound waves... light waves... water waves...

Effect of Nozzle Twisted Lance on Jet Behavior and Spitting Rate in Top Blown Process

Friction properties of the face of a hand-held tennis racket

THE BRIDGE COLLAPSED IN NOVEMBER 1940 AFTER 4 MONTHS OF ITS OPENING TO TRAFFIC!

Simulation of Free Surface Flows with Surface Tension with ANSYS CFX

Aeroelasticity of Turbomachines

Pre-swirl mechanism in front of a centrifugal compressor: effects on surge line and on unsteady phenomena in surge area

Chapter 4 Surface Tension

Chapter 14. Vibrations and Waves

Application of Simulation Technology to Mitsubishi Air Lubrication System

Dillon Thorse Flow Visualization MCEN 4047 Team Poject 1 March 14th, 2013

Available online at Procedia Engineering 200 (2010) (2009) In situ drag measurements of sports balls

EXPERIMENTAL STUDY ON THE HYDRODYNAMIC BEHAVIORS OF TWO CONCENTRIC CYLINDERS

CFD AND EXPERIMENTAL STUDY OF AERODYNAMIC DEGRADATION OF ICED AIRFOILS

AERODYNAMIC CHARACTERISTICS OF SPIN PHENOMENON FOR DELTA WING

Vertical Wind Energy Engineering Design and Evaluation of a Twisted Savonius Wind Turbine

Conceptual Physics. Chapter 25: Vibrations and Waves Mr. Miller

Lesson 14: Simple harmonic motion, Waves (Sections )

Transcription:

Walking with coffee: when and why coffee spills Hans C. Mayer and Rouslan Krechetnikov Department of Mechanical Engineering University of California at Santa Barbara February 20-24, 2012 Page 1/25

Motivation and introduction It is common everyday knowledge to each of us that any small container filled with liquid must be moved or carried very carefully to avoid spills. - H. Norman Abramson (NASA SP-106, 1966) Spilling coffee when walking a problem for a wide audience Analysis of when and why coffee spills is a confluence of many subject including: (i) biomechanics, (ii) sloshing engineering, and (iii) dynamical systems. The motions of the human body, while seemingly regular, are quite complex and coupled to a container/liquid. As we have all walked with coffee (or liquid filled containers), we have all experienced this problem it is a familiar problem. Page 2/25

Introduction: Sloshing vs. Biomechanics Biomechanics of walking Large body of literature (primarily medical) dedicated to the study of human walking: Motion of body CM, energy expenditure, efficiency of process. Quantifying gait patterns related to gender, age, health, etc. Motion of appendages (arms and legs) in natural positions. Regularity of walking over long periods of time. Sloshing and engineering Large body of handbooks, technical documents, and journal articles (NASA SP-106, NASA SP 8301, Ibrahim, R.A. (2005) ) concerned with Large liquid filled structures (tanks). Structures subjected to considerable accelerations/vibrations. Forces and torques due to sloshing liquid. Suppression of sloshing necessary for vehicle control. Page 3/25

Motivation and introduction: Sloshing vs. Biomechanics Sloshing in cylindrical containers (a) Lateral sloshing: response to translational or pitching motions (antisymmetric modes) (b) Vertical sloshing: response to vertical motions normal to equilibrium surface (symmetric modes) (c) Rotational sloshing: swirling motion of liquid Lateral Back-Forth Vertical Rotational (Abramson et. al, NASA SP-106) Page 4/25

Experimental setup: bits and pieces Coffee mug geometry Filming walking Image analysis Page 5/25

Experimental setup: coffee mug geometry Properties of coffee similar to water (neglect surface tension) (neglect viscosity) Page 6/25

Experimental setup: capturing mug motion top view side view MATLAB code (using image analysis toolbox): along a straight line path motion analyzed in the xy - plane finds markers determines mug coordinates Manual determination of spill (using LED signal) Page 7/25

Experimental setup: movie files Page 8/25

Results from image analysis: coffee mug kinematics (x direction) Typical output of image analysis Note the following: Page 9/25

Results from image analysis: coffee mug kinematics (y direction) Typical output of image analysis Page 10/25

Results from image analysis: coffee mug kinematics (θ direction) Image analysis of coffee mug motion during walking Typical output of image analysis Although wrist is flexible, one can anticipate that the human subject will not allow the mug to pitch substantially (small θ values) Pitching motion captured is really a combination of arm swing and flexible wrist. Demonstrates walking motion produces back-and-forth, vertical, and pitching oscillations. Page 11/25

Results from data analysis: back-and-forth excitation Due to the oscillatory nature of the x direction motion, backand-forth forcing of the coffee mug exists, but it is apparent that the magnitude of the oscillations change with walking speed. Standard deviation of the velocity about the average decreases with increasing walking speed. Page 12/25

Results from data analysis: vertical excitation Magnitude of y - motion, y - velocity, and y acceleration increase with walking speed (although average y, y - velocity, and y - acceleration are essentially zero). Page 13/25

Results from data analysis: comparing all data What can we observe about walking danger zone Maximum accelerations compared to those needed to initiate nearly immediate spilling (based on α) Page 14/25

Results from data analysis: step frequency Page 15/25

Comparison of components of body motion (biomechanics) Back-and-Forth Lateral Vertical M. W. Whittle Human Movement Science 16 pg. 347 (1997) Page 16/25

Results from data analysis: sloshing natural oscillations frequency Cylindrical vessel Antisymmetric mode (natural frequency) Calculated for our coffee mug approximately the same frequency for all coffee mugs! R. A. Ibrahim Liquid Sloshing Dynamics (2005) Conclusion: back-and-forth motion should excite the antisymmetric sloshing mode Page 17/25

Results from data analysis: liquid motion and observation of swirl Closer look suggests that back-and-forth and swirling liquid motions can be excited by walking at typical speeds. (a) back-and-forth sloshing motion (Dodge (2005)) Question: why is rotational motion generated? (b) rotational sloshing motion Page 18/25

Results from data analysis: how long does it take to spill? Time to spill is extracted from all walking movies interpret as nearly instant spill Page 19/25

Results from data analysis: how long does it take to spill? (and thus walking speed) Page 20/25

A more complex picture: dynamical systems point of view y motion of mug smooth oscillatory motion noise Decomposition of motion Effect of noise: Looking more closely at the mug motion during constant acceleration reveals that it is composed of noise and smooth oscillations. Comparing focused and unfocused regimes suggests that noise substantially intensifies the spilling and noise intensity is lower in the focused regime. Unsteady effects: Since walking is usually non-uniform, the frequency of step changes and thus when excitation approaches the natural frequency one can expect resonance to occur much earlier (jump frequency phenomenon) Page 21/25

Modeling of spilling as an equivalent spherical pendulum: steps to spill (a) (b) Figure (a) shows steps to spill N/P versus the normalized initial deflection of the interface for three different coffee levels. Each data point shown in the figure corresponds to the average number of steps N of all model runs that yield a spill divided by the fraction of runs P that result in a spill. Figure (b) reflects the influence of the noise amplitude an on the number of steps to spill for a fixed angle. Page 22/25

Modeling of spilling as an equivalent spherical pendulum: noise effects (c) (d) Figure (c) shows the effect of time-dependent (chirped) forcing frequency on coffee spilling over short time (inset) and longer times. Figure (d) shows example trajectory of spherical pendulum projected onto the x-y plane. A critical spill radius is denoted by the solid red line. Page 23/25

Conclusions: how do we spill coffee? * (a) Either by accelerating too much for a given coffee level (fluid statics) (a) Or, through more complicated dynamical phenomena: Initial acceleration sets an initial sloshing amplitude, which is analogous to the main antisymmetric mode of sloshing. This initial perturbation is amplified by the back-and-forth and pitching excitations since their frequency is close to the natural one because of the choice of normal mug dimensions. Vertical motion also does not lead to resonance as it is a subharmonic excitation (Faraday phenomena). Noise has higher frequency, which makes the antisymmetric mode unstable thus generating a swirl. Time to spill depends on focused / unfocused regime and increases with decreasing maximum acceleration (walking speed). * Scientific American 305, 22 (2011) Page 24/25

Conclusions: how can we prevent spilling? Lessons learned from sloshing dynamics may suggest strategies to control spilling, e.g. via using a flexible container to act as a sloshing absorber in suppressing liquid oscillations. a series of concentric rings (baffles) arranged around the inner wall of a container. (Abramson, NASA SP-8031 (1969)) (Langner SRI Tech Report 7 (1963)) Page 25/25