CFLOW H. Massol, C. Jaupart GOALS AND INTENT

CFLOW H. Massol, C. Jaupart GOALS AND INTENT

CFLOW H. Massol, C. Jaupart GOALS AND INTENT OF CFLOW 2-D PRESSURE STATE IN THE CONDUIT EXPLOSIVITY OF LAVA DOMES ESTIMATE OF GAS OVERPRESSURE MT UNZEN, JAPAN HETEROGENEITY OF GAS CONTENT IN FLOWS AND DOMES GAS LOSS THROUGH CONDUIT WALLS MT ST HELENS, USA

SCHEMATIC VIEW OF A VOLCANIC CONDUIT BUBBLY SECTION INTEGRATION DOMAIN Exsolution level VISCOUS AND COMPRESSIBLE FLOW

INCOMPRESSIBLE FLOW ORIGINALITY: 2-D METHOD: Finite Element MAGMA CHAMBER RHEOLOGY (1) = -2 e + 2 ( 3 .v) + Pg - K ( Pg = P + K ( .v)

.v) Pg: Gas pressure Shear viscosity K: Bulk viscosity VISCOSITY+COMPRESSIBILITY GAS OVERPRESSURE DOME EXPLOSIVITY RHEOLOGY (2) o , p f P

g b Pm l K, R [ r r ] r=R

2 RR 2 = pg - 4 l 3 b b ( 4 1- K = l 3 b3 1- 3

R ) [ rr ]r=R = pm - 3K 2 pm = pb b R R z

zz = p atm u=0 orrz = 0 DOMAIN AND BC H u=0 w=0 u=0

rz = 0 h 0 zz = ps r u=0 a Criteres Criteres darret

BASIC EQUATIONS darret 10-2 10-210-8 10-8 Conservation of momentum [ D ] U = SU [ D ] W = SW Conservation of mass [ M ] = S Artificial time Mass lumping

Petrov Galerkin weighting CAPABILITIES VARIABLE MELT VISCOSITY VARIABLE CONDUIT GEOMETRY HORIZONTAL AND VERTICAL VELOCITY COMPONENTS VARIABLE COMPRESSIBILITY ASSUMPTIONS AND LIMITS EQUILIBRIUM DEGASSING ONLY VALID BEFORE FRAG. LEVEL ANALYTICAL SOLUTION HYPOTHESES

- No horizontal velocity - Constant compressibility - Constant viscosity BUT: Gas pressure varies in both directions Numerical model benchmark RESULTS PARABOLIC PRESSURE PROFILE ACROSS THE CONDUIT DIMENSIONLESS NUMBER, D PH-Pa K+4/3 a2 = P

=D 2 P0-pa H EXAMPLE RESULT (1) U=0 = 106 Pa.s x0 = 0.5 Wt% P(0,H) = 0.56 MPa OVERPRESSURE AT THE CENTER OF THE CONDUIT EXIT

Compatible with the analytical solution EXAMPLE RESULT (2) rz =0 x0 = 0.5 Wt%

P(0,H) = 0.56 MPa P(a,H) = 1.5 MPa Variable viscosity (Hess and Dingwell, 96) EXAMPLE RESULT (3) x0 = 0.5 Wt% P(0,H) = 0.34 MPa P(a,H) = 1.1 MPa CONCLUSIONS GAS PHASE IS OVERPRESSURED / DOME EXPLOSIVITY

HORIZONTAL PRESSURE GRADIENT / VITRIFIED MARGINS, HETEROGENEITY IN GAS CONTENT IN FLOWS AND DOMES IMPORTANCE OF THE EXIT BOUNDARY CONDITIONS / CREASE STRUCTURE FUTURE WORK BOUNDARY CONDITIONS - CONDUIT WALLS AND - COUPLING WITH FLOW CRYSTALS NUCLASCENT (H. Massol, T. Koyaguchi)

Turbulent Flow 1-D Finite difference Cylindrical Geometry Steady state Variable viscosity Non-equilibrium degassing Fragmentation Level Laminar Flow Nucleation of Bubbles EVOLUTION OF DISSOLVED WATER IN THE MELT

H=5000 m a = 50 m 0 = 106 Pa.s x0 = 4wt% = 0.02 N m-1 D = 10-11 m2 s-1 EVOLUTION OF PRESSURE AND NUMBER OF BUBBLES H=5000 m a = 50 m 0 = 106 Pa.s x0 = 4wt% = 0.02 N m-1 D = 10-11 m2 s-1

OUTPUT OF THE MODEL BUBBLE SIZES BUBBLE DENSITY PRESSURE INSIDE BUBBLES NEXT STEP: CONTINUOUS BUBBLE SIZE DISTRIBUTION

Recently Viewed Presentations

  • Circle1- Limbo Dantes First Circle of Hell is

    Circle1- Limbo Dantes First Circle of Hell is

    The Eight Circle of Hell is resided by the fraudulent. Dante and Virgil reach it on the back of Geryon, a flying monster with different natures, just like the fraudulent. This circle of Hell is divided into 10 Bolgias or...
  • Airedale and Partners Enhanced Health in Care Homes

    Airedale and Partners Enhanced Health in Care Homes

    Teleconsultation by secure video link between nursing and residential homes and the Airedale Digital Hub. Providing clinical consultation not a logarithm based approach like 111.
  • 1. This destructive force can be either chemical or physical.

    1. This destructive force can be either chemical or physical.

    This destructive force can be either chemical or physical. a. Erosionb. Depositionc. ... How do farmers use contour plowing and terracing to help their crops? a. Contour plowing and terracing ... This destructive force can be either chemical or physical.
  • C++ Programming: Program Design Including Data Structures ...

    C++ Programming: Program Design Including Data Structures ...

    Some compilers initialize global variables to default values Scope resolution operator in C++ is :: By using the scope resolution operator A global variable declared before the definition of a function (or block) can be accessed by the function (or...
  • CH 103: SOLUBILITY AND SOLUTIONS  In todays experiment

    CH 103: SOLUBILITY AND SOLUTIONS In todays experiment

    What is mass percent? Which, if any, of these concentrations will change with temperature? Only molarity (M) will change with temperature. The volume of solution changes with temperature. In contrast, the moles of solute, mass of solute, moles of solvent,...
  • Communication Skills - Loudoun County Public Schools

    Communication Skills - Loudoun County Public Schools

    Formal Communication. A formal communication network is the system of official channels that carry organizationally approved messages.. Communication channels generally follow the reporting relationships in the firm. Communication flows upward, downward, and across the organization in a prescribed manner.
  • Excess Reactant - teachnlearnchem.com

    Excess Reactant - teachnlearnchem.com

    Notice: Ratio is the same in both… 2 mol H2(g) + 1 mol O2(g) 2 mol H2O(g) limiting reactant = H2 excess reactant = O2 2 mol H2(g) 3 mol O2(g) 2 mol H2O(g) + 2 mol O2 left over...
  • The Market System and the Circular Flow

    The Market System and the Circular Flow

    Households are the owners of the resources (selling their resources in the resource market) and income flows into the households. Households are the ones who buy the goods and services in the products market. The prices that are paid in...