Acceleration and Deceleration Phase Nonlinear Rayleigh-Taylor Growth at Spherical Interfaces

Similar documents
Equations from The Four Principal Kinetic States of Material Bodies. Copyright 2005 Joseph A. Rybczyk

Invert and multiply. Fractions express a ratio of two quantities. For example, the fraction

10.3 The Quadratic Formula

ME 141. Engineering Mechanics

Physics 201, Lecture 5

How dark matter, axion walls, and graviton production lead to observable Entropy generation in the Early Universe. Dr.

Reinforcement learning

v T Pressure Extra Molecular Stresses Constitutive equations for Stress v t Observation: the stress tensor is symmetric

Some algorthim for solving system of linear volterra integral equation of second kind by using MATLAB 7 ALAN JALAL ABD ALKADER

Ch.4 Motion in 2D. Ch.4 Motion in 2D

() t. () t r () t or v. ( t) () () ( ) = ( ) or ( ) () () () t or dv () () Section 10.4 Motion in Space: Velocity and Acceleration

Science Advertisement Intergovernmental Panel on Climate Change: The Physical Science Basis 2/3/2007 Physics 253

Reinforcement learning

2.25 Advanced Fluid Mechanics

LECTURE 5. is defined by the position vectors r, 1. and. The displacement vector (from P 1 to P 2 ) is defined through r and 1.

Limit Analysis of Stability of Circular Foundation Pit

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Equations and Boundary Value Problems

Numerical Simulation of Natural Gas Flow in Transmission Lines through CFD Method Considering Effective Parameters

Fluids & Bernoulli s Equation. Group Problems 9

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

Control Volume Derivation

Additional Exercises for Chapter What is the slope-intercept form of the equation of the line given by 3x + 5y + 2 = 0?

D zone schemes

Computer Aided Geometric Design

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Photographing a time interval

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

Lecture 2 - Thermodynamics Overview

Department of Chemical Engineering University of Tennessee Prof. David Keffer. Course Lecture Notes SIXTEEN

In the previous section we considered problems where the

F (u) du. or f(t) = t

P a g e 5 1 of R e p o r t P B 4 / 0 9

BINOMIAL THEOREM OBJECTIVE PROBLEMS in the expansion of ( 3 +kx ) are equal. Then k =

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

_ J.. C C A 551NED. - n R ' ' t i :. t ; . b c c : : I I .., I AS IEC. r '2 5? 9

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

Molecular Dynamics Simulations (Leach )

Motion on a Curve and Curvature

Numerical solution of fuzzy differential equations by Milne s predictor-corrector method and the dependency problem

COSC 3361 Numerical Analysis I Ordinary Differential Equations (I) - Introduction

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Compressive modulus of adhesive bonded rubber block

Problem Set 5: Universal Law of Gravitation; Circular Planetary Orbits

0 for t < 0 1 for t > 0

Classification of Equations Characteristics

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997)

The sudden release of a large amount of energy E into a background fluid of density

4/3/2017. PHY 712 Electrodynamics 9-9:50 AM MWF Olin 103

A L A BA M A L A W R E V IE W

III. Direct evolution of the density: The Liouville Operator

PHYS PRACTICE EXAM 2

Comparison between the Discrete and Continuous Time Models

Computer Propagation Analysis Tools

148 CIVIL ENGINEERING

Answers to test yourself questions

4.8 Improper Integrals

Overview. Overview Page 1 of 8

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

Data Structures. Element Uniqueness Problem. Hash Tables. Example. Hash Tables. Dana Shapira. 19 x 1. ) h(x 4. ) h(x 2. ) h(x 3. h(x 1. x 4. x 2.

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

4. Modelling of oscillating-water-column wave energy converters

STUDY OF THE STRESS-STRENGTH RELIABILITY AMONG THE PARAMETERS OF GENERALIZED INVERSE WEIBULL DISTRIBUTION

Final Exam. Tuesday, December hours, 30 minutes

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

Homework 5 for BST 631: Statistical Theory I Solutions, 09/21/2006

1 Using Integration to Find Arc Lengths and Surface Areas

02. MOTION. Questions and Answers

Technical Vibration - text 2 - forced vibration, rotational vibration

Non-Ideal Gas Behavior P.V.T Relationships for Liquid and Solid:

û s L u t 0 s a ; i.e., û s 0

Summary: Binomial Expansion...! r. where

The Global Trade and Environment Model: GTEM

Section 11. Timescales Radiation transport in stars

T h e C S E T I P r o j e c t

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions

Previously. Extensions to backstepping controller designs. Tracking using backstepping Suppose we consider the general system

The Production of Polarization

Experiment 1 Electric field and electric potential

Pythagorean Theorem and Trigonometry

Mathematical Reflections, Issue 5, INEQUALITIES ON RATIOS OF RADII OF TANGENT CIRCLES. Y.N. Aliyev

Water Tunnel Experiment MAE 171A/175A. Objective:

«ACTA»2002. Introduction. Theory. Dissipative Systems. I. D. Chueshov Introduction to the Theory of Infinite-Dimensional Dissipative Systems

Module 4: Moral Hazard - Linear Contracts

Physics Courseware Electromagnetism

Lect. 7 Discrete Fourier Transforms Discrete Fourier Transform: discrete representation of the Fourier Integral

Lecture 10. Solution of Nonlinear Equations - II

Radial geodesics in Schwarzschild spacetime

2D Motion WS. A horizontally launched projectile s initial vertical velocity is zero. Solve the following problems with this information.

On Control Problem Described by Infinite System of First-Order Differential Equations

1. Viscosities: μ = ρν. 2. Newton s viscosity law: 3. Infinitesimal surface force df. 4. Moment about the point o, dm

Illustrating the space-time coordinates of the events associated with the apparent and the actual position of a light source

Lecture 17: Kinetics of Phase Growth in a Two-component System:

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Pasadena Police Department Crime Laboratory Corrective Action

Addition & Subtraction of Polynomials

Transcription:

Aeleion nd Deeleion Pse Nonline Ryleig-Tylo Gow Speil Inefes Pesened o: HEDP Summe Sool Dniel S. Clk nd Mx Tbk Lwene Livemoe Nionl Lbooy Augus 5 Tis wok ws pefomed unde e uspies of e U.S. Depmen of Enegy by e Univesiy of Clifoni Lwene Livemoe Nionl Lbooy unde Con No. W-745-Eng-48. Lwene Livemoe Nionl Lbooy, P.O. Box 88, Livemoe, CA 9455-88 UCRL-POST-434

Abs Te lssi Lyze model fo e nonline evoluion of bubbles in e Ryleig-Tylo insbiliy s eenly been genelized o e se of eleing, speilly imploding inefes [D.S.Clk nd M.Tbk, Pys.Rev.E 7, 553R 5.], e onfiguion mos elevn o ineil onfinemen fusion. Hee, e model is fue exended o e se of bubble gow duing e deeleion sgnion pse of speil sell implosion nd o e gow of spikes duing bo e eleion nd deeleion pses. Diffeenes in e nonline gow es fo bo bubbles nd spikes e found wen omped wi e pln esuls. Te model pediions e veified by ompison wi numeil ydodynmis simulions. UCRL-POST-434

Speil onvegene modifies nonline RT gow fom Lyze model expeions pln speil g Lyze s [955] model is e lssi model fo desibing nonline RT gow bu is sily vlid only fo pln inefes We ve ombined Lyze-like RT model nd Kidde-ype self-simil implosion o develop model fo speil inefes moe elevn o insbiliies on ICF psules Model evels iniil gow like e Lyze model, followed by le ime modifiions Compison wi -D simulions onfims pedied bubble spike gow UCRL-POST-434 3

Moivions Undesnding e nonline pse of RT gow is impon fo IFE wee opimized psule designs my pus e limis of ydodynmi insbiliy. Lyze s model simply quie uely desibes nonline RT gow bu is sily vlid only in pln geomey. Is ee n nlogous model fo speilly onveging sysem moe elevn o ICF ges? Finding n nlogous model mig evel ineesing sling popeies of nonline bubble gow. Su model ould povide igoous elevn nonline es poblem o vlide ydo. odes. Te speil poblem s no been solved, bu ppes solvble UCRL-POST-434 4

Ouline Review of Lyze model Modifiions neessy fo speilly onveging inefes Conneion o self-simil implosions of [Kidde`74] RT bubble model fo onveging geomey: eleion pse gow sgnion deeleion pse gow Lyze-ype models of RT spike gow RT spike model fo onveging geomey: eleion pse gow sgnion deeleion pse gow Appliion o non-self-simil implosions Summy/Conlusions UCRL-POST-434 5

Lyze model is e lssi model fo nonline pse of RT gow unsble inefe g evy fluid line RT gow nonline gow Lyze model z lig fluid spike bubble g u S Appoxime nonline RT pse by bubble ising wiin ylindil boundy UCRL-POST-434 6

Model is nonline bu lol soluion of e ydo. eqns. oninuiy eqn. inompessible/ ioionl Benoulli eqn. momenum inefe moion eqn. φ d d F φ S lnρ Lyze nsz dominn mode: φ φ S ; v φ p gz ρ ; S S Z z nonline in mpliude g u z S k z / φ φ A e J k / / sisfies ydo. eqns. o lowes ode in bu lol UCRL-POST-434 7

Bubbles ise wi onsn veloiies pln inefes inompessible keep only s mode k z / φ φ A e J k / subsiue Lyze model F φ S p φ gz ρ φ S ; S S Z z bubble eig line pse -D simulion Tylo expnd in ububble g / k wee J k ime bu model is sily vlid only fo pln inefes wile ICF ges e speil nd ompessible UCRL-POST-434 8

Subsnil modifiions e equied o oun fo speil onvegene z Reple ylindil boundies wi onil boundies in speil geomey g u S u β Reple unifom gviy g wi speilly oning oodines / g UCRL-POST-434 9

UCRL-POST-434 Hydo. eqns. mus be solved in oning fme fo onveging se Lyze nsz inompessible fo in oning fme, uomilly leds o ompessing flow fixed fme! S S p F S Φ φ φ ρ γ γ φ φ φ φ / exp os 3 ρ ρ P A - φ φ φ oninuiy eqn. in fixed fme oninuiy eqn. in oning fme Lyze nsz fo mode φ effeive g oning fme eqns.

UCRL-POST-434 Inefe eqn. is solved fo mpliude nd uvue subsiue expnd inefe spe in oning fme sepe odes in, ;, b R R S { } 4 3 φ O A b A b A b A S S 4 4 : : b A b A b O A A O mpliude uvue inegion ons. O

UCRL-POST-434 Two nul soluions fo inefe eqn. bubbles spikes 4, 4 O R bubble-ype inefe: inefe spe, 4, >> R,, >> R pidly inesing uvue ses iniil uvue self-simil /ons. uvue bubble, R Flow sould sympoe o one of ese spes depending on iniil vlue of spike-ype inefe: > b < b, R

UCRL-POST-434 3 Fo bubbles, Benoulli eqn. yields oupled ODE s fo mpliude nd sle fo os φ φ P A subsiue inompessible in oning fme R S p F S 4, φ ρ γ γ φ φ Benoulli eqn. on inefe Tylo expnd in 3 : ons. : O O γ onion sle fo keep only s mode mpliude oupled ODE s

In limi >>, O eqn. beomes Kidde s eqn. fo self-simil implosion >> Fo ypil se, eqns. deouple nd s ODE edues o eqn. fo Kidde s self-simil implosion Kidde s eqn. 3γ ons. / / Kidde implosion ime ; γ 5/3 implosion digm self-simil isenopi pofiles fee fll /.7 /.5 /. /.7 /.5 /. UCRL-POST-434 4

Wi speified, nd eqn. n be solved by WKB meod Kidde / Assuming, WKB ppoximion yields bubble eig ~ exp ± d / ~ / / / / / 3/ in good geemen wi numeil soluion of oupled ODEs / ; 3/ 8 numeil soln. WKB soln. u Lyze d u Lyze / Iniilly line gow wi ime s in Lyze model followed by enned bubble gow le in implosion UCRL-POST-434 5

Enned bubble gow onfimed by -D simulions ig pessue lig fluid spike one wlls 8 -D simulion numeil soln. R, ompessing dense fluid u Lyze WKB soln. u Lyze d bubble /. /.5 / onvegene/implosion /.5 Simil esuls found fo nge of mode numbes ; equied subsnil ALE /.75 UCRL-POST-434 6

Kidde s eqn. lso desibes self-simil sgnion By evesing sign of ons. in s ODE flipping dieion of pessue gdien, se desibes self-simil sgnion Kidde s eqn. 3γ >> ons. / oose sgning implosion veloiy / sgnion ; γ 5/3 implosion digm self-simil isenopi pofiles fee fll /.7 /.7 /.5 /.7 /.7 /.5 /. /. /.5 /.5 /. /. UCRL-POST-434 7

Enned bubble gow lso ous duing sgnion Fo sgnion pse, eigenmodes i.e.,, so WKB soln. wi is φ A P sgnion os e ppopie sgnion ~ ~ 5 5 exp ± exp n n d / ; / implosion 8 WKB soln. u Lyze explosion u Lyze d Lyze uve uses pek deeleion / u Lyze Lyze UCRL-POST-434 8

-D simulions onfim enned gow bu sow influene of KH R, ompessing dense fluid ig pessue lig fluid bubble one wlls 8 -D simulion WKB soln. u Lyze KH jeing u Lyze d /. spike /.5 onvegene/implosion KH / Simil esuls found fo nge of mode numbes /.9 /. inefeing voies je ino elee bubble UCRL-POST-434 9

Lyze models n lso desibe spike gow [Zng`98] s sown Lyze-ype models n lso desibe sympoi pse of spike gow: bubble z spike z g u φ φ g u z Z, find spikes fee-fll fom inefe nd sympoe o infinie uvue Z, ~ g Z, ; UCRL-POST-434

Model fo spike gow wi onvegene n lso be onsued Fom genel soln. of inefe eqn. ose spike-ype inefe eleion pse: b R, O 4 deeleion pse: b R, R, pidly inesing uvue, >>, >> 4 R, R, Flow sould sympoe o spike-ype spe if iniil uvue is lge enoug UCRL-POST-434

Benoulli eqn. gin yields oupled ODE s fo mpliude nd onion sle fo spike spe Lyze nsz b φ A P os subsiue ino Benoulli eqn. on inefe nd expnd oupled ODE s O : O : ons. 3γ b 4 b b 4 b / O d d Keeping only leding ode in, eqn. edues o mpliude onion sle fo UCRL-POST-434

Spikes fee-sem in fixed fme O Wi deemined s befoe by, spike mpliude is simply eleion pse: deeleion pse: ; > i.e., in fixed fme, spikes oninue fee-seming moion s in [Zng`98] [ ] ; < < b / ~ ± m Cuvue ip lso s found by [Zng`98] Spike moion is independen of mode numbe, i.e.,., boundy ondiions onil wlls e ignoed s sould be. UCRL-POST-434 3

Model pediions gee wi -D simulions fo eleion pse ig pessue lig fluid spike KH 8 KH model soln. ompessing dense fluid bubble ½ -D simulion u Lyze d /. /.5 onvegene/implosion /.5 /.75 one wlls / Idenil esuls found fo nge of, i.e., spike gow independen of bounding wlls UCRL-POST-434 4

Model pediions gee wi -D simulions fo sgnion pse ompessing dense fluid bubble one wlls 8 model soln. ½ -D simulion KH/wll effe u Lyze d ig pessue lig fluid /. spike /.5 onvegene/implosion /.5 / Gow sues by KH nd ollision wi one wlls /.75 UCRL-POST-434 5

Model pediions gee well wi simulions fo nge of mode numbes bubbles spikes deeleion pse eleion pse UCRL-POST-434 6

Typil ICF implosions e lose o n self-similiy g ons. fee fll : g ons. O Wiou espeing eqn., model esuls fo bubbles o spikes my sill pply fo neby 's Kidde: self-simil bubble spike ~ ~ 3/ exp [ ] d UCRL-POST-434 7

WKB esul my pply o non-self-simil implosions Seing g, WKB esul beomes ~ 3/ ~ exp exp exp ± sin sin g R g R /3/- /3/- in oug geemen wi simulion d u Lyze -D simulion numeil soln. WKB soln. u Lyze d g / R Simil esuls found fo oe mode numbes ; efinemen of mes is on-going Fuue wok: invesige modifiions of Hn suion model due o onvegene. UCRL-POST-434 8

Speil onvegene modifies nonline RT gow fom Lyze model expeions pln speil g Lyze s [955] model is e lssi model fo desibing nonline RT gow bu is sily vlid only fo pln inefes We ve ombined Lyze-like RT model nd Kidde-ype self-simil implosion o develop model fo speil inefes moe elevn o insbiliies on ICF psules Model evels iniil gow like e Lyze model, followed by le ime modifiions Compison wi -D simulions onfims pedied bubble spike gow UCRL-POST-434 9