DOCSLIB.ORG

The Structure of Vorticity in Cumulonimbus Convection

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
The Structure of Vorticity in Cumulonimbus Convection 9_9:3 THE STRUCTURE OF VORTICITY IN CUMULONIMBUS CONVECTION: A NUMERICAL STUDY by Martin Johnson Atmospheric Physics Group Department of Physics Imperial College of Science and Technology A thesis submitted for the Degree of Doctor of Philosophy in the University of London August 1978 • THE STRUCTURE OF VORTICITY IN CUMULONIMBUS CONVECTION : A NUMERICAL STUDY Martin Johnson ABSTRACT This thesis investigates the vorticity and potential vorticity fields generated by a three-dimensional primitive equation numerical model of deep convection. Using the momentum and continuity of mass equations solved by Miller's deep. convection model, a vorticity equation is derived in both continuous and finite-difference forms. The latter is derived in such a way as to be consistent with both its continuous analogue and the primitive equations.. It is shown that the primitive equation model time integration scheme introduces second-order time errors into the horizontal components of the vorticity equation, and a change in the scheme is suggested which would remove these errors. Using the finite-difference form of the vorticity equation, the fields of the three components of (i) the vorticity vector and (ii) the diagnostic terms in the vorticity equation are computed for three different simulations: a squall line, a stationary storm, and a splitting storm system. In all cases the vorticity field is surprisingly well organised. It is found that the horizontal components of vorticity reach values of order of magnitude 10-2 s-1 in all the storms around the edges of the updraught cores and downdraught 'cold pools'. The vertical vorticity attains this magnitude over a relatively large area in the first and third cases, and over a very restricted area in the second. This corresponds to a significantly larger amount of (horizontal) vorticity in the mean flow of the first and third cases. In all the storms, the vertical compohent appears as a 'vortex doublet', the positive (cyclonic) value being about twice the negative; the area of positive vertical vorticity lies just inside the updraught on the side nearest to the downdraught, while the negative area lies on the opposite side of the updraught. Both are contained mainly in the 500-800 mb layer. A simplified model of the vorticity field is described in terms of vortex lines, and using the results of the diagnostic calculations, a theory is suggested for their development. Comparisons are made between the results of the above, and dual- Doppler radar observations of severe storms of the United States and computer models of tornadoes. Equations are obtained for potential vorticity, liquid water potential vorticity, and equivalent potential vorticity, and it is shown that due to baroclinity, diffusion and diabatic effects, these quantities are not conserved. The results are presented of the computation of the first of these quantities. CONTENTS ABSTRACT 1 NOTATION 5 CHAPTER ONE - INTRODUCTION 6 1.1 Observations of cumulonimbus convection 7 1.2 Analytical models of cumulonimbus convection 10 1.3 Numerical models of convection 13 1.4 Tornadoes 17 CHAPTER TWO - DEFINITIONS OF TERMS AND DERIVATION OF VORTICITY AND CIRCULATION EQUATIONS IN CONTINUOUS FORM 20 2.1 Definitions 20 2.2 Derivation of the vorticity equation 21 2.2.1 Primitive equations 22 2.2.2 Vorticity equation 23 2.3 Circulation equation 27 2.3.1 Derivation of circulation equation 27 2.3.2 Consequences of circulation equation 27 CHAPTER THREE - DERIVATION OF VORTICITY EQUATION - FINITE DIFFERENCE FORM 29 3.1 Primitive equation model of cumulonimbus convection 29 3.1.1 Definition of symbols 29 3.1.2 Configuration of model grid 30 3.1.3 Primitive equations in finite difference form; time integration scheme 32 3.2 Derivation of vorticity equation 34 3.2.1 Choice of finite-difference scheme and points in grid for calculations 34 3.2.2 The vorticity vector in finite difference form; further simplification of the vertical momentum equation 37 3.2.3 Derivation of the vorticity equation: x-component 38 3.2.4 Derivation of the vorticity equation: y-component 44 3. 3.2.5 Derivation of the vorticity equation: vertical component 45 CHAPTER FOUR - THE VORTICAL STRUCTURE OF SIMULATED CUMULONIMBUS 48 4.1 Method of presentation of results 48 4.2 Case A : Tropical cumulonimbus 49 4.3 Case B : The Hampstead storm 51 -1 4.3.1 • Results for 32-44 minute means 52 4.3.2 Results for 64-76 minute means 53 4.4 Case C : Splitting storm system 53 Key to diagrams 55 Diagrams for case A 56 Diagrams for case B (i) 68 Diagrams for case B (ii) 80 Diagrams for case C 92 CHAPTER FIVE - INTERPRETATION OF RESULTS 104 5.1 Model of vorticity field in cumulonimbus 104 5.2 Physical interpretation of results 108 5.3 Comparison of model and simulation results 116 5.4 Generation of vertical component; comparison with dual-Doppler observations 121 5.5 Tornadoes 123 CHAPTER SIX - POTENTIAL VORTICITY 125 6.1 Equation of general potential vorticity of deep convection 125 6.2 Special cases 126 6.2.1 'Ordinary' potential vorticity 126 6.2.2 Equivalent potential vorticity 128 6.2.3 Liquid water potential vorticity 129 6.2.4 Potential vorticity of deep convection 129 6.3 Potential vorticity equations_insimplified models of cumulonimbus 131 6.3.1 Moist inviscid cumulonimbus model 131 6.3.2 Dry inviscid cumulonimbus model 132 6.4 Potential vorticity fields in simulated cumulonimbus 133 4. CHAPTER SEVEN - SUMMARY AND CONCLUDING REMARKS 143 7.1 Summary 143 7.2 Suggestions for future research 145 APPENDIX - PRIMITIVE EQUATION MODEL TIME INTEGRATION SCHEME 147 REFERENCES 152 ACKNOWLEDGEMENTS 156 5. NOTATION Cp Specific heat of air at constant pressure F Initiating heat source F F F Acceleration components due to subgrid scale oe 3 P momentum transfer Vector (Fe Fy,- ygps) g Acceleration due to gravity g(1 + Z ) g* h Height of isobaric surfaces K K8' q Diffusion coefficients Total mass of liquid water per unit mass of air Z Mass of cloudwater per unit mass of air c Mass of rainwater per unit mass of air Latent heat of condensation of water p Pressure q Specific humidity Q Sink of specific humidity due to condensation T Temperature u E (uw,w) Velocity vector ✓ Terminal velocity of raindrops Circulation Vorticity vector Potential temperature 0 Virtual potential temperature eE Equivalent potential temperature W Liquid water potential temperature Potential vorticity of cumulonimbus convection TICB Equivalent potential vorticity HE LW Liquid water potential vorticity P Air density In 6 $* Zn 64 Zn e (I)E E cLW in OLW 0 Flux of vorticity along vortex tube; strength of vortex tube X (e*ye s ) - gh' + 1/2 (u2 v2) Dp/Dt SZ Instantaneous angular velocity of fluid element CORRIGENDA AND ADDENDA Page 21 line 6 Delete 'It can be shown,' .... 'that'. Page 37 line 14 For 'differerence' read 'difference' Page 38 line 15 For 'AE ' read 'AE 1/2,1/2,1/2 0,0,1/2 Page 38 line 16 For '6Y.F s ' read I SYPX ' P - Page 49 line 25 It should be made clear that the horizontal dimensions of the vorticity model were restricted to 20km x 20 km due to limits on computer storage area, Although slight problems do exist near the primitive equation model boundaries, the fact that these are avoided in the vorticity model is merely a co-incidental advantage _ of the restriction. Page 50 line 15 For 'Fig. 4.7(a)' read 'Fig. 4.9(a)' Page 104 line 15 Delete 'i.e. the dominant wavenumber of the vorticity field is the same as that of the storm system'. Pages 106-7, Figs. 5.1-5.3 Vortex lines are assumed to be moving from left to right in these figures. Pages 108, 137, Figs. 5.4, 6.4 The trajectory shown in these figures are computed according to the method presented in Miller and Betts (1977). Pages 108-110 In addition to the information presented on Fig.5.4 and in the text, the data in the following Table may be found useful. The Table illustrates in more detail how fl and the terms in the equation for ri vary-along the trajectory in the immediate vicinity of the updraught. Position on n Dn/dt Stret- Tilt- Comp- Brclnty Diff- trajectory ching ing ress. usion relative to 29kmx29km grid 10-7s-2 I J P 10-3s-1 10-7s-2 10-7s-2 10-7s- 10-7s-2 10-7s-2 11.6 8.6 939 -6.57 -99 -16 2 3 -510 390 12.4 8.6 904 -2.21 221 24 -6 5 159 39 13.0 9.2 816 3.30 53 34 -14 -35 76 -8 13.2 9.3 606 7..57 322 77 7 -12 763 -513 14.6 9.2 460 5.75 262 108 17 -76 399 -186 15.6 9.5 447 6.31 450 100 -9 -36 270 125 Corrigenda and Addenda (continued) Page 110 Line 2 For '-88' read Page 155 The following reference should be added: Miller, M.J. and Betts, A.K. 1977 Traveling Convective Storms over Venezuela, Mon.Wea.Rev., 105, pp.833-848. 6. CHAPTER ONE INTRODUCTION 'Convection' is the name given to any motion arising directly from the conversion of available potential energy and which is involved in the exchange of heat, moisture and momentum. This type of motion occurs on all scales, and there is a close interaction between them; thus precise knowledge of the features of one particular scale of motion is required in order to understand fully the other scales of motion with which it interacts.
Load more

Recommended publications
  • Energy, Vorticity and Enstrophy Conserving Mimetic Spectral Method for the Euler Equation
    Master of Science Thesis Energy, vorticity and enstrophy conserving mimetic spectral method for the Euler equation D.J.D. de Ruijter, BSc 5 September 2013 Faculty of Aerospace Engineering · Delft University of Technology Energy, vorticity and enstrophy conserving mimetic spectral method for the Euler equation Master of Science Thesis For obtaining the degree of Master of Science in Aerospace Engineering at Delft University of Technology D.J.D. de Ruijter, BSc 5 September 2013 Faculty of Aerospace Engineering · Delft University of Technology Copyright ⃝c D.J.D. de Ruijter, BSc All rights reserved. Delft University Of Technology Department Of Aerodynamics, Wind Energy, Flight Performance & Propulsion The undersigned hereby certify that they have read and recommend to the Faculty of Aerospace Engineering for acceptance a thesis entitled \Energy, vorticity and enstro- phy conserving mimetic spectral method for the Euler equation" by D.J.D. de Ruijter, BSc in partial fulfillment of the requirements for the degree of Master of Science. Dated: 5 September 2013 Head of department: prof. dr. F. Scarano Supervisor: dr. ir. M.I. Gerritsma Reader: dr. ir. A.H. van Zuijlen Reader: P.J. Pinto Rebelo, MSc Summary The behaviour of an inviscid, constant density fluid on which no body forces act, may be modelled by the two-dimensional incompressible Euler equations, a non-linear system of partial differential equations. If a fluid whose behaviour is described by these equations, is confined to a space where no fluid flows in or out, the kinetic energy, vorticity integral and enstrophy integral within that space remain constant in time. Solving the Euler equations accompanied by appropriate boundary and initial conditions may be done analytically, but more often than not, no analytical solution is available.
    [Show full text]
  • Soaring Weather
    Chapter 16 SOARING WEATHER While horse racing may be the "Sport of Kings," of the craft depends on the weather and the skill soaring may be considered the "King of Sports." of the pilot. Forward thrust comes from gliding Soaring bears the relationship to flying that sailing downward relative to the air the same as thrust bears to power boating. Soaring has made notable is developed in a power-off glide by a conven­ contributions to meteorology. For example, soar­ tional aircraft. Therefore, to gain or maintain ing pilots have probed thunderstorms and moun­ altitude, the soaring pilot must rely on upward tain waves with findings that have made flying motion of the air. safer for all pilots. However, soaring is primarily To a sailplane pilot, "lift" means the rate of recreational. climb he can achieve in an up-current, while "sink" A sailplane must have auxiliary power to be­ denotes his rate of descent in a downdraft or in come airborne such as a winch, a ground tow, or neutral air. "Zero sink" means that upward cur­ a tow by a powered aircraft. Once the sailcraft is rents are just strong enough to enable him to hold airborne and the tow cable released, performance altitude but not to climb. Sailplanes are highly 171 r efficient machines; a sink rate of a mere 2 feet per second. There is no point in trying to soar until second provides an airspeed of about 40 knots, and weather conditions favor vertical speeds greater a sink rate of 6 feet per second gives an airspeed than the minimum sink rate of the aircraft.
    [Show full text]
  • Effect of Wind on Near-Shore Breaking Waves
    EFFECT OF WIND ON NEAR-SHORE BREAKING WAVES by Faydra Schaffer A Thesis Submitted to the Faculty of The College of Engineering and Computer Science in Partial Fulfillment of the Requirements for the Degree of Master of Science Florida Atlantic University Boca Raton, Florida December 2010 Copyright by Faydra Schaffer 2010 ii ACKNOWLEDGEMENTS The author wishes to thank her mother and family for their love and encouragement to go to college and be able to have the opportunity to work on this project. The author is grateful to her advisor for sponsoring her work on this project and helping her to earn a master’s degree. iv ABSTRACT Author: Faydra Schaffer Title: Effect of wind on near-shore breaking waves Institution: Florida Atlantic University Thesis Advisor: Dr. Manhar Dhanak Degree: Master of Science Year: 2010 The aim of this project is to identify the effect of wind on near-shore breaking waves. A breaking wave was created using a simulated beach slope configuration. Testing was done on two different beach slope configurations. The effect of offshore winds of varying speeds was considered. Waves of various frequencies and heights were considered. A parametric study was carried out. The experiments took place in the Hydrodynamics lab at FAU Boca Raton campus. The experimental data validates the knowledge we currently know about breaking waves. Offshore winds effect is known to increase the breaking height of a plunging wave, while also decreasing the breaking water depth, causing the wave to break further inland. Offshore winds cause spilling waves to react more like plunging waves, therefore increasing the height of the spilling wave while consequently decreasing the breaking water depth.
    [Show full text]
  • NWS Unified Surface Analysis Manual
    Unified Surface Analysis Manual Weather Prediction Center Ocean Prediction Center National Hurricane Center Honolulu Forecast Office November 21, 2013 Table of Contents Chapter 1: Surface Analysis – Its History at the Analysis Centers…………….3 Chapter 2: Datasets available for creation of the Unified Analysis………...…..5 Chapter 3: The Unified Surface Analysis and related features.……….……….19 Chapter 4: Creation/Merging of the Unified Surface Analysis………….……..24 Chapter 5: Bibliography………………………………………………….…….30 Appendix A: Unified Graphics Legend showing Ocean Center symbols.….…33 2 Chapter 1: Surface Analysis – Its History at the Analysis Centers 1. INTRODUCTION Since 1942, surface analyses produced by several different offices within the U.S. Weather Bureau (USWB) and the National Oceanic and Atmospheric Administration’s (NOAA’s) National Weather Service (NWS) were generally based on the Norwegian Cyclone Model (Bjerknes 1919) over land, and in recent decades, the Shapiro-Keyser Model over the mid-latitudes of the ocean. The graphic below shows a typical evolution according to both models of cyclone development. Conceptual models of cyclone evolution showing lower-tropospheric (e.g., 850-hPa) geopotential height and fronts (top), and lower-tropospheric potential temperature (bottom). (a) Norwegian cyclone model: (I) incipient frontal cyclone, (II) and (III) narrowing warm sector, (IV) occlusion; (b) Shapiro–Keyser cyclone model: (I) incipient frontal cyclone, (II) frontal fracture, (III) frontal T-bone and bent-back front, (IV) frontal T-bone and warm seclusion. Panel (b) is adapted from Shapiro and Keyser (1990) , their FIG. 10.27 ) to enhance the zonal elongation of the cyclone and fronts and to reflect the continued existence of the frontal T-bone in stage IV.
    [Show full text]
  • An Overview of Impellers, Velocity Profile and Reactor Design
    An Overview of Impellers, Velocity Profile and Reactor Design Praveen Patel1, Pranay Vaidya1, Gurmeet Singh2 1Indian Institute of Technology Bombay, India 1Indian Oil Corporation Limited, R&D Centre Faridabad Abstract: This paper presents a simulation operation and hence in estimating the approach to develop a model for understanding efficiency of the operating system. The the mixing phenomenon in a stirred vessel. The involved processes can be analysed to optimize mixing in the vessel is important for effective the products of a technology, through defining chemical reaction, heat transfer, mass transfer the model with adequate parameters. and phase homogeneity. In some cases, it is Reactors are always important parts of any very difficult to obtain experimental working industry, and for a detailed study using information and it takes a long time to collect several parameters many readings are to be the necessary data. Such problems can be taken and the system has to be observed for all solved using computational fluid dynamics adversities of environment to correct for any (CFD) model, which is less time consuming, inexpensive and has the capability to visualize skewness of values. As such huge amount of required data, it takes a long time to collect the real system in three dimensions. enough to build a model. In some cases, also it As reactor constructions and impeller is difficult to obtain experimental information configurations were identified as the potent also. Pilot plant experiments can be considered variables that could affect the macromixing as an option, but this conventional method has phenomenon and hydrodynamics, these been left far behind by the advent of variables were modelled.
    [Show full text]
  • Vorticity and Strain Analysis Using Mohr Diagrams FLOW
    Journal of Structural Geology, Vol, 10, No. 7, pp. 755 to 763, 1988 0191-8141/88 $03.00 + 0.00 Printed in Great Britain © 1988 Pergamon Press plc Vorticity and strain analysis using Mohr diagrams CEES W. PASSCHIER and JANOS L. URAI Instituut voor Aardwetenschappen, P.O. Box 80021,3508TA, Utrecht, The Netherlands (Received 16 December 1987; accepted in revisedform 9 May 1988) Abstract--Fabric elements in naturally deformed rocks are usually of a highly variable nature, and measurements contain a high degree of uncertainty. Calculation of general deformation parameters such as finite strain, volume change or the vorticity number of the flow can be difficult with such data. We present an application of the Mohr diagram for stretch which can be used with poorly constrained data on stretch and rotation of lines to construct the best fit to the position gradient tensor; this tensor describes all deformation parameters. The method has been tested on a slate specimen, yielding a kinematic vorticity number of 0.8 _+ 0.1. INTRODUCTION Two recent papers have suggested ways to determine this number in naturally deformed rocks (Ghosh 1987, ONE of the aims in structural geology is the reconstruc- Passchier 1988). These methods rely on the recognition tion of finite deformation parameters and the deforma- that certain fabric elements in naturally deformed rocks tion history for small volumes of rock from the geometry have a memory for sense of shear and flow vorticity and orientation of fabric elements. Such data can then number. Some examples are rotated porphyroblast- be used for reconstruction of large-scale deformation foliation systems (Ghosh 1987); sets of folded and patterns and eventually in tracing regional tectonics.
    [Show full text]
  • QUICK REFERENCE GUIDE XT and XE Ovens
    QUICK REFERENCE GUIDE XT AND XE OVENS SYMBOL OVEN FUNCTION FUNCTION SELECTION SETTING BEST FOR Select on the desired Activates the upper and the temperature through the This cooking mode is suitable for any Bake either oven control knob lower heating elements or the touch control smart kind of dishes and it is great for baking botton. 100°F - 500°F and roasting. Food probe allowed. Select on the desired Reduced cooking time (up to 10%). Activates the upper and the temperature through the Ideal for multi-level baking and roasting Convection bake either oven control knob lower heating elements or the touch control smart of meat and poultry. botton. 100°F - 500°F Food probe allowed. Select on the desired Reduced cooking time (up to 10%). Activates the circular heating temperature through the Ideal for multi-level baking and roasting, Convection elements and the convection either oven control knob especially for cakes and pastry. Food fan or the touch control smart probe allowed. botton. 100°F - 500°F Activates the upper heating Select on the desired Reduced cooking time (up to 10%). temperature through the Best for multi-level baking and roasting, Turbo element, the circular heating either oven control knob element and the convection or the touch control smart especially for pizza, focaccia and fan botton. 100°F - 500°F bread. Food probe allowed. Ideal for searing and roasting small cuts of beef, pork, poultry, and for 4 power settings – LOW Broil Activates the broil element (1) to HIGH (4) grilling vegetables. Food probe allowed. Activates the broil element, Ideal for browning fish and other items 4 power settings – LOW too delicate to turn and thicker cuts of Convection broil the upper heating element (1) to HIGH (4) and the convection fan steaks.
    [Show full text]
  • Numerical Study of a Mixing Layer in Spatial and Temporal Development of a Stratified Flow M. Biage, M. S. Fernandes & P. Lo
    Transactions on Engineering Sciences vol 18, © 1998 WIT Press, www.witpress.com, ISSN 1743-3533 Numerical study of a mixing layer in spatial and temporal development of a stratified flow M. Biage, M. S. Fernandes & P. Lopes Jr Department of Mechanical Engineering Center of Exact Science and Tecnology Federal University ofUberldndia, UFU- MG-Brazil [email protected] Abstract The two-dimensional plane shear layers for high Reynolds numbers are simulated using the transient explicit Spectral Collocation method. The fast and slow air streams are submitted to different densities. A temporally growing mixing layer, developing from a hyperbolic tangent velocity profile to which is superposed an infinitesimal white-noise perturbation over the initial condition, is studied. In the same way, a spatially growing mixing layer is also study. In this case, the flow in the inlet of the domain is perturbed by a time dependent white noise of small amplitude. The structure of the vorticity and the density are visualized for high Reynolds number for both cases studied. Broad-band kinetic energy and temperature spectra develop after the first pairing are observed. The results showed that similar conclusions are obtained for the temporal and spatial mixing layers. The spreading rate of the layer in particular is in good agreement with the works presented in the literature. 1 Introduction The intensely turbulent region formed at the boundary between two parallel fluid streams of different velocity has been studied for many years. This kind of flow is called of mixing layer. The mixing layer of a newtonian fluid, non-reactive and compressible flow, practically, is one of the most simple turbulent flow.
    [Show full text]
  • Scaling the Vorticity Equation
    Vorticity Equation in Isobaric Coordinates To obtain a version of the vorticity equation in pressure coordinates, we follow the same procedure as we used to obtain the z-coordinate version: ∂ ∂ [y-component momentum equation] − [x-component momentum equation] ∂x ∂y Using the p-coordinate form of the momentum equations, this is: ∂ ⎡∂v ∂v ∂v ∂v ∂Φ ⎤ ∂ ⎡∂u ∂u ∂u ∂u ∂Φ ⎤ ⎢ + u + v +ω + fu = − ⎥ − ⎢ + u + v +ω − fv = − ⎥ ∂x ⎣ ∂t ∂x ∂y ∂p ∂y ⎦ ∂y ⎣ ∂t ∂x ∂y ∂p ∂x ⎦ which yields: d ⎛ ∂u ∂v ⎞ ⎛ ∂ω ∂u ∂ω ∂v ⎞ vorticity equation in ()()ζ p + f = − ζ p + f ⎜ + ⎟ − ⎜ − ⎟ dt ⎝ ∂x ∂y ⎠ p ⎝ ∂y ∂p ∂x ∂p ⎠ isobaric coordinates Scaling The Vorticity Equation Starting with the z-coordinate form of the vorticity equation, we begin by expanding the total derivative and retaining only nonzero terms: ∂ζ ∂ζ ∂ζ ∂ζ ∂f ⎛ ∂u ∂v ⎞ ⎛ ∂w ∂v ∂w ∂u ⎞ 1 ⎛ ∂p ∂ρ ∂p ∂ρ ⎞ + u + v + w + v = − ζ + f ⎜ + ⎟ − ⎜ − ⎟ + ⎜ − ⎟ ()⎜ ⎟ ⎜ ⎟ 2 ⎜ ⎟ ∂t ∂x ∂y ∂z ∂y ⎝ ∂x ∂y ⎠ ⎝ ∂x ∂z ∂y ∂z ⎠ ρ ⎝ ∂y ∂x ∂x ∂y ⎠ Next, we will evaluate the order of magnitude of each of these terms. Our goal is to simplify the equation by retaining only those terms that are important for large-scale midlatitude weather systems. 1 Scaling Quantities U = horizontal velocity scale W = vertical velocity scale L = length scale H = depth scale δP = horizontal pressure fluctuation ρ = mean density δρ/ρ = fractional density fluctuation T = time scale (advective) = L/U f0 = Coriolis parameter β = “beta” parameter Values of Scaling Quantities (midlatitude large-scale motions) U 10 m s-1 W 10-2 m s-1 L 106 m H 104 m δP (horizontal) 103 Pa ρ 1 kg m-3 δρ/ρ 10-2 T 105 s -4 -1 f0 10 s β 10-11 m-1 s-1 2 ∂ζ ∂ζ ∂ζ ∂ζ ∂f ⎛ ∂u ∂v ⎞ ⎛ ∂w ∂v ∂w ∂u ⎞ 1 ⎛ ∂p ∂ρ ∂p ∂ρ ⎞ + u + v + w + v = − ζ + f ⎜ + ⎟ − ⎜ − ⎟ + ⎜ − ⎟ ()⎜ ⎟ ⎜ ⎟ 2 ⎜ ⎟ ∂t ∂x ∂y ∂z ∂y ⎝ ∂x ∂y ⎠ ⎝ ∂x ∂z ∂y ∂z ⎠ ρ ⎝ ∂y ∂x ∂x ∂y ⎠ ∂ζ ∂ζ ∂ζ U 2 , u , v ~ ~ 10−10 s−2 ∂t ∂x ∂y L2 ∂ζ WU These inequalities appear because w ~ ~ 10−11 s−2 the two terms may partially offset ∂z HL one another.
    [Show full text]
  • A STUDY of ELEVATED CONVECTION and ITS IMPACTS on SURFACE WEATHER CONDITIONS a Dissertat
    A STUDY OF ELEVATED CONVECTION AND ITS IMPACTS ON SURFACE WEATHER CONDITIONS _______________________________________ A Dissertation presented to the Faculty of the Graduate School at the University of Missouri-Columbia _______________________________________________________ In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy _____________________________________________________ by Joshua Kastman Dr. Patrick Market, Dissertation Supervisor December 2017 © copyright by Joshua S. Kastman 2017 All Rights Reserved The undersigned, appointed by the dean of the Graduate School, have examined the dissertation entitled A Study of Elevated Convection and its Impacts on Surface Weather Conditions presented by Joshua Kastman, a candidate for the degree of doctor of philosophy, Soil, Environmental, and Atmospheric Sciences and hereby certify that, in their opinion, it is worthy of acceptance. ____________________________________________ Professor Patrick Market ____________________________________________ Associate Professor Neil Fox ____________________________________________ Professor Anthony Lupo ____________________________________________ Associate Professor Sonja Wilhelm Stannis ACKNOWLEDGMENTS I would like to begin by thanking Dr. Patrick Market for all of his guidance and encouragement throughout my time at the University of Missouri. His advice has been invaluable during my graduate studies. His mentorship has meant so much to me and I look forward his advice and friendship in the years to come. I would also like to thank Anthony Lupo, Neil Fox and Sonja Wilhelm-Stannis for serving as committee members and for their advice and guidance. I would like to thank the National Science Foundation for funding the project. I would l also like to thank my wife Anna for her, support, patience and understanding while I pursued this degree. Her love and partnerships means everything to me.
    [Show full text]
  • Ductile Deformation - Concepts of Finite Strain
    327 Ductile deformation - Concepts of finite strain Deformation includes any process that results in a change in shape, size or location of a body. A solid body subjected to external forces tends to move or change its displacement. These displacements can involve four distinct component patterns: - 1) A body is forced to change its position; it undergoes translation. - 2) A body is forced to change its orientation; it undergoes rotation. - 3) A body is forced to change size; it undergoes dilation. - 4) A body is forced to change shape; it undergoes distortion. These movement components are often described in terms of slip or flow. The distinction is scale- dependent, slip describing movement on a discrete plane, whereas flow is a penetrative movement that involves the whole of the rock. The four basic movements may be combined. - During rigid body deformation, rocks are translated and/or rotated but the original size and shape are preserved. - If instead of moving, the body absorbs some or all the forces, it becomes stressed. The forces then cause particle displacement within the body so that the body changes its shape and/or size; it becomes deformed. Deformation describes the complete transformation from the initial to the final geometry and location of a body. Deformation produces discontinuities in brittle rocks. In ductile rocks, deformation is macroscopically continuous, distributed within the mass of the rock. Instead, brittle deformation essentially involves relative movements between undeformed (but displaced) blocks. Finite strain jpb, 2019 328 Strain describes the non-rigid body deformation, i.e. the amount of movement caused by stresses between parts of a body.
    [Show full text]
  • Solution to Two-Dimensional Incompressible Navier-Stokes Equations with SIMPLE, SIMPLER and Vorticity-Stream Function Approaches
    Solution to two-dimensional Incompressible Navier-Stokes Equations with SIMPLE, SIMPLER and Vorticity-Stream Function Approaches. Driven-Lid Cavity Problem: Solution and Visualization. by Maciej Matyka Computational Physics Section of Theoretical Physics University of Wrocław in Poland Department of Physics and Astronomy Exchange Student at University of Link¨oping in Sweden [email protected] http://panoramix.ift.uni.wroc.pl/∼maq 30 czerwca 2004 roku Streszczenie In that report solution to incompressible Navier - Stokes equations in non - dimensional form will be presented. Standard fundamental methods: SIMPLE, SIMPLER (SIMPLE Revised) and Vorticity-Stream function approach are compared and results of them are analyzed for standard CFD test case - Drived Cavity flow. Different aspect ratios of cavity and different Reynolds numbers are studied. 1 Introduction will be solved on rectangular, staggered grid. Then, solu- tion on non-staggered grid with vorticity-stream function The main problem is to solve two-dimensional Navier- form of NS equations will be shown. Stokes equations. I will consider two different mathemati- cal formulations of that problem: ¯ u,v,p primitive variables formulation 2 Math background ¯ ζ, ψ vorticity-stream function approach We will consider two-dimensional Navier-Stokes equations I will provide full solution with both of these methods. in non-dimensional form1: First we will consider three standard, primitive component formulations, where fundamental Navier-Stokes equation 1We consider flow without external forces i.e. without gravity. −→ Guess: Solve (3),(4) for: ∂ u −→ −→ 1 2−→ = −( u ∇) u − ∇ϕ + ∇ u (1) n n n n+1 n+1 ∂t Re (P*) ,(U*) ,(V*) (U*) ,(V*) D = ∇−→u = 0 (2) Where equation (2) is a continuity equation which has Solve (6) for: to be true for the final result.
    [Show full text]