Bridging the Gap: Harnessing advanced Computational Libraries for coupled multi-physics problems
|
|
- Leo McCormick
- 5 years ago
- Views:
Transcription
1 Bridging the Gap: Harnessing advanced Computational Libraries for coupled multi-physics problems Marc Spiegelman 1,2 and Cian R Wilson 1 Plus enormous contributions from PETSc, FEniCS and AMCG groups 1 Lamont-Doherty Earth Observatory, Columbia University, New York 2 DEES/DAPAM, Columbia University, New York BTG, Princeton, 1 2 October 2012
2 Introduction Motivation Computation is essential for exploring non-linear multi-physics problems in Geosciences Major advances in hardware, software and algorithms make complex problems more accessible However, still a considerable gap between Geosciences and Computational Sciences/Math Some Existing Barriers: Overall Complexity of both software libraries and models Lack of transparency/reproducibility/reusability of current models Language/Training issues between computation and Earth Sciences Spiegelman and Wilson TerraFERMA BTG
3 Introduction Motivation Computation is essential for exploring non-linear multi-physics problems in Geosciences Major advances in hardware, software and algorithms make complex problems more accessible However, still a considerable gap between Geosciences and Computational Sciences/Math Some Existing Barriers: Overall Complexity of both software libraries and models Lack of transparency/reproducibility/reusability of current models Language/Training issues between computation and Earth Sciences We need better tools for exploring complex (and simple) models and making advanced computation accessible to more general users Spiegelman and Wilson TerraFERMA BTG
4 Target Application : multi-physics models of magmatic plate boundaries image courtesy of the Neptune project (wwwneptunewashingtonedu) Central to Solid Earth Geosciences: Essential for understanding Seismic, Volcanic and Tsunamagenic Natural Hazards global tectonics and geochemical cycling Spiegelman and Wilson TerraFERMA BTG
5 Target Application : multi-physics models of magmatic plate boundaries image courtesy of the Neptune project (wwwneptunewashingtonedu) Central to Solid Earth Geosciences: Essential for understanding Seismic, Volcanic and Tsunamagenic Natural Hazards global tectonics and geochemical cycling Multi-physics: Tightly coupled non-linear systems Coupled fluid-solid mechanics Complex solid rheologies Coupled thermodynamics/geodynamics Spiegelman and Wilson TerraFERMA BTG
6 Target Application : multi-physics models of magmatic plate boundaries image courtesy of the Neptune project (wwwneptunewashingtonedu) Central to Solid Earth Geosciences: Essential for understanding Seismic, Volcanic and Tsunamagenic Natural Hazards global tectonics and geochemical cycling Multi-physics: Tightly coupled non-linear systems Coupled fluid-solid mechanics Complex solid rheologies Coupled thermodynamics/geodynamics Considerable uncertainty in equations, constitutive relations and coupling Spiegelman and Wilson TerraFERMA BTG
7 Model Requirements User Flexibility: Choice of equations, geometry, elements, constitutive relations, coupling Wide choice of solvers/preconditioners for coupled non-linear problems Ability to rapidly compose a range of models from simple process models to regional geodynamic models All choices available at or near run time Spiegelman and Wilson TerraFERMA BTG
8 Model Requirements Infrastructure: Residual monitoring for the full coupled problem Model reproducibility: Transparent options system Regression tested Generic model services: Standardized I/O Checkpointing Monitoring and detectors Parallel, scaleable, open source, version controlled, regression tested free Spiegelman and Wilson TerraFERMA BTG
9 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation Au = f F(u) = 0 u n+1 = f(u n ) Spiegelman and Wilson TerraFERMA BTG
10 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation Au = f F(u) = 0 u n+1 = f(u n ) Spiegelman and Wilson TerraFERMA BTG
11 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation Au = f F(u) = 0 u n+1 = f(u n ) Spiegelman and Wilson TerraFERMA BTG
12 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation Au = f F(u) = 0 u n+1 = f(u n ) Spiegelman and Wilson TerraFERMA BTG
13 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation Au = f F(u) = 0 u n+1 = f(u n ) P E T S c Spiegelman and Wilson TerraFERMA BTG
14 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation F E ni C S Au = f F(u) = 0 u n+1 = f(u n ) P E T S c Spiegelman and Wilson TerraFERMA BTG
15 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation F E ni C S G m sh Au = f F(u) = 0 u n+1 = f(u n ) P E T S c Spiegelman and Wilson TerraFERMA BTG
16 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation S P u D F E ni C S G m sh Au = f F(u) = 0 u n+1 = f(u n ) P E T S c Spiegelman and Wilson TerraFERMA BTG
17 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation S P u D F E ni C S G m sh Au = f F(u) = 0 u n+1 = f(u n ) P E T S c TerraFERMA Transparent Finite Element Rapid Model Assembler Spiegelman and Wilson TerraFERMA BTG
18 Application: Fluid Migration in Subduction Zones ( ( )) vs + v T s 2η = p +ϕ 0 ϕk 2 v s = ϕ 0 P ζ K µ P + P ( ) K ζ = µ k p + Γ ρ ρ f ϕ 0 Dϕ Dt = (1 ϕ 0ϕ) P ζ + Γ ϕ 0 Spiegelman and Wilson TerraFERMA BTG
19 Example: Infinite Prandtl Thermal Convection 2 µ v + v T 2 + p T k = 0, v = 0, T v T + T =0 t Ra Spiegelman and Wilson TerraFERMA BTG
20 Example: Infinite Prandtl Thermal Convection 2 µ v + v T 2 + p T k = 0, v = 0, T v T + T =0 t Ra Spiegelman and Wilson TerraFERMA BTG
21 Example: Infinite Prandtl Thermal Convection 2 µ v + v T 2 + p T k = 0, v = 0, T v T + T =0 t Ra Spiegelman and Wilson TerraFERMA BTG
22 Example: Infinite Prandtl Thermal Convection 2 µ v + v T 2 + p T k = 0, v = 0, T v T + T =0 t Ra Let u = (v, p, T ) u i+1 = u i αj(u i ) 1 r (u i ) Spiegelman and Wilson TerraFERMA BTG
23 Example: FEniCS Equation Description Given function u = (v, p, T ), test function u t = (v t, p t, T t ) and trial function u a = (v a, p a, T a ): Weak form of residual, r(u): rv = r p = r T = Ω Ω Ω [( v t + v T t 2 ) ( ) ] v + v T : 2µ v t p (v t ) z T, 2 p t v, [ T t ((T T n ) + tv θ T θ ) + t ] Ra T t T θ r =rv + r p + r T Weak form of Jacobian, J (u): J = r (u) Spiegelman and Wilson TerraFERMA BTG
24 Example: FEniCS Equation Description Given function u = (v, p, T ), test function u t = (v t, p t, T t ) and trial function u a = (v a, p a, T a ): Weak form of residual, r(u): r v = (inner(sym(grad(v t)), 2*mu*sym(grad(v))) - div(v t)*p - T*v t[1] )*dx r p = p t*div(v)*dx r T = (T t*((t - T n) + dt*inner(v theta, grad(t theta))) + (dt/ra)*inner(grad(t t), grad(t theta)) )*dx r = r v + r p + r T Weak form of Jacobian, J (u): J = derivative(r, u, u a) Spiegelman and Wilson TerraFERMA BTG
25 Example: PETSc Block Preconditioning & Solvers J(u i )δu = r(u i ) J = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
26 Example: PETSc Block Preconditioning & Solvers J(u i )δu = r(u i ) J = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Stokes Advection-Diffusion Spiegelman and Wilson TerraFERMA BTG
27 Example: PETSc Block Preconditioning & Solvers J PC (u i ) 1 J(u i )δu = J PC (u i ) 1 r(u i ) J = J PC = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T 2 M 0 B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
28 Example: PETSc Block Preconditioning & Solvers J PC (u i ) 1 J(u i )δu = J PC (u i ) 1 r(u i ) J = J PC = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T 2 M 0 B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
29 Example: PETSc Block Preconditioning & Solvers J PC (u i ) 1 J(u i )δu = J PC (u i ) 1 r(u i ) J = J PC = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T 2 M 0 B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
30 Example: PETSc Block Preconditioning & Solvers J PC (u i ) 1 J(u i )δu = J PC (u i ) 1 r(u i ) J = J PC = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T 2 M 0 B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
31 Example: PETSc Block Preconditioning & Solvers J PC (u i ) 1 J(u i )δu = J PC (u i ) 1 r(u i ) J = J PC = K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T B 1 B 2 0 A K 11 K 12 G 1 C 1 K 21 K 22 G 2 C 2 G T 1 G T 2 M 0 B 1 B 2 0 A System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Velocity 1 (preonly,mg) Velocity 2 (preonly,mg) Spiegelman and Wilson TerraFERMA BTG
32 Generic Structure of PDE Based Numerical Models Options System Constitutive PDEs Eqs Meshes/ Geometry Discrete Eqs Solvers Visualization Error Est Validation S P u D F E ni C S G m sh Au = f F(u) = 0 u n+1 = f(u n ) P E T S c TerraFERMA Transparent Finite Element Rapid Model Assembler Spiegelman and Wilson TerraFERMA BTG
33 Example: TerraFERMA Options Spiegelman and Wilson TerraFERMA BTG
34 Example: TerraFERMA Options Spiegelman and Wilson TerraFERMA BTG
35 Example: Setting Equations in TerraFERMA Spiegelman and Wilson TerraFERMA BTG
36 Example: Setting Equations in TerraFERMA Spiegelman and Wilson TerraFERMA BTG
37 Example: Solver Options in TerraFERMA Spiegelman and Wilson TerraFERMA BTG
38 Example: Solver Options in TerraFERMA System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Spiegelman and Wilson TerraFERMA BTG
39 Example: Solver Options in TerraFERMA System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Spiegelman and Wilson TerraFERMA BTG
40 Example: Solver Options in TerraFERMA System (fgmres, fieldsplit) Temperature (gmres, ilu) Stokes (fgmres, fieldsplit) Pressure (cg,sor) Velocity (preonly,fieldsplit) Spiegelman and Wilson TerraFERMA BTG
41 Example: Benchmark (Blankenbach, 1989) Nu v rms q 1 q 2 T e z e Ra = 10 4 isoviscous Ra = 10 5 isoviscous Ra = 10 6 isoviscous Ra = 10 4 T -dep viscosity Benchmark Benchmark Benchmark Benchmark Spiegelman and Wilson TerraFERMA BTG
42 Current Benchmarks Incompressible 2D Convection (Blankenbach et al, 1989) Incompressible 3D Convection (Busse et al, 1994) Cylindrical Laminar Plumes (Vatteville et al, 2009) Compressible 2D Convection (King et al, 2009) Kinematic Subduction Zones (van Keken et al, 2008) Linearized Free Surface Evolution (Kramer et al, 2012) Non-linear magma waves (Simpson and Spiegelman, 2011) Cylindrical Convection (Rhodri Davies: Spiegelman and Wilson TerraFERMA BTG
43 Fluid Migration Buoyancy Buoyancy and Pore Pressure Spiegelman and Wilson TerraFERMA BTG
44 Fluid Migration Buoyancy Buoyancy and Pore Pressure Spiegelman and Wilson TerraFERMA BTG
45 Fluid Migration Buoyancy Buoyancy and Pore Pressure Spiegelman and Wilson TerraFERMA BTG
46 Fluid Migration Buoyancy Buoyancy and Pore Pressure Spiegelman and Wilson TerraFERMA BTG
47 Fluid Migration Buoyancy Buoyancy and Pore Pressure Spiegelman and Wilson TerraFERMA BTG
48 Conclusions Many important problems in Earth Sciences can be described by non-linear coupled systems of partial differential equations Much of the complexity modelling these systems stems from the nature of multi-physics where small changes in the coupling between fields or constitutive relations can lead to radical changes in behavior and the resulting demands on discretizations and solvers Many established single-physics models have locked in solution strategies that are unsuited or difficult to extend to the more challenging non-linear behaviour of a multi-physics system We require a significant increase in flexibility for users for both composing problems (from simple model problems to more realistic coupled simulations) and changing solver strategies Several advanced, open-source computational libraries mean that it is now possible to provide much greater flexibility to the user TerraFERMA is a model building framework that gives access to these libraries while also providing a transparent description of the model Spiegelman and Wilson TerraFERMA BTG
49 Spiegelman and Wilson TerraFERMA BTG
50 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p + ϕ 0 ϕ k 2 v s = ϕ 0 P ζ Fluid flow for pore pressure, P, and porosity, ϕ: ( K µ P + P ζ = K µ k p Dϕ Dt = ( ) P 1 ϕ 0 ϕ ζ + Γ ϕ 0 ) + Γ ρ ρ f ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
51 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p + ϕ 0 ϕ k 2 v s = ϕ 0 P ζ Fluid flow for pore pressure, P, and porosity, ϕ: ( K µ P + P ζ = K µ k p Dϕ Dt = ( ) P 1 ϕ 0 ϕ ζ + Γ ϕ 0 ) + Γ ρ ρ f ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
52 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p + ϕ 0 ϕ k 2 v s = ϕ 0 P ζ Fluid flow for pore pressure, P, and porosity, ϕ: ( K µ P + P ζ = K µ k p Dϕ Dt = ( ) P 1 ϕ 0 ϕ ζ + Γ ϕ 0 ) + Γ ρ ρ f ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
53 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p + ϕ 0 ϕ k 2 v s = ϕ 0 P ζ Fluid flow for pore pressure, P, and porosity, ϕ: ( K µ P + P ζ = K µ k p Dϕ Dt = ( ) P 1 ϕ 0 ϕ ζ + Γ ϕ 0 ) + Γ ρ ρ f ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
54 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p + ϕ 0 ϕ k 2 v s = ϕ 0 P ζ Fluid flow for pore pressure, P, and porosity, ϕ: ( K µ P + P ζ = K µ k p Dϕ Dt = ( ) P 1 ϕ 0 ϕ ζ + Γ ϕ 0 ) + Γ ρ ρ f ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
55 Equations (Katz et al, 2007) Solid flow for solid velocity, v s, and dynamic pressure, p : ( ( )) vs + v T s 2 η = p 2 v s = 0 Fluid flow for pore pressure, P, and porosity, ϕ: K µ P + P ζ = K µ k + Γ ρ ρ f ϕ 0 Dϕ Dt = P ζ + Γ ϕ 0 Constitutive relations: η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
56 Solid Flow Velocity streamlines 50km T = T s T = zt 0 /50 v i = 0 600km v = v 0 T = T 0 T = T s + (T 0 T s)erf ( z/ κt age) Spiegelman and Wilson TerraFERMA BTG
57 Equations Buoyancy and Pore Pressure: Constitutive relations: K µ P + P ζ = K µ k + Γ ρ ρ f ϕ 0 Dϕ Dt = P ζ + Γ ϕ 0 η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
58 Equations Buoyancy (zero compaction length approximation): Dϕ Dt = K µ k + Γ ( ) ρ + 1 ϕ 0 ρ f Buoyancy and Pore Pressure: Constitutive relations: K µ P + P ζ = K µ k + Γ ρ ρ f ϕ 0 Dϕ Dt = P ζ + Γ ϕ 0 η = η (T, ɛ, ϕ), ζ = ηϕ m, K = ϕ n Spiegelman and Wilson TerraFERMA BTG
An Introduction and Tutorial to the McKenzie Equations for magma migration
An Introduction and Tutorial to the McKenzie Equations for magma migration Marc Spiegelman 1,2, Richard Katz 1,3 and Gideon Simpson 2 1 Lamont-Doherty Earth Obs. of Columbia University 2 Dept. of Applied
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationNonlinear Preconditioning in PETSc
Nonlinear Preconditioning in PETSc Matthew Knepley PETSc Team Computation Institute University of Chicago Department of Molecular Biology and Physiology Rush University Medical Center Algorithmic Adaptivity
More informationNonlinear Preconditioning in PETSc
Nonlinear Preconditioning in PETSc Matthew Knepley PETSc Team Computation Institute University of Chicago Challenges in 21st Century Experimental Mathematical Computation ICERM, Providence, RI July 22,
More informationBenchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003
Introduction Benchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003 In early October 2002 a group of researchers met at the University of Michigan at Ann Arbor,
More informationA benchmark for the modeling of subduction zone thermal structure
A benchmark for the modeling of subduction zone thermal structure Peter van Keken (UM) Claire Currie (Dalhousie) Scott King (PU) Richard Katz (LDEO) Marc Spiegelman (LDEO) History Gre out of 2002 Margins
More informationNumerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation
Physics of the Earth and Planetary Interiors 163 (2007) 52 68 Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation R.F. Katz a,, M.G. Knepley b,
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationThe effect of different geometries on the thermal mantle convection
The effect of different geometries on the thermal mantle convection Mátyás Herein Eötvös Loránd University Faculty of Science Department of Geophysics and Space Sciences 2011 I. Introduction Outline II.
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationSubduction II Fundamentals of Mantle Dynamics
Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011 Rheology Elasticity vs. viscous deformation η
More informationBeall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/ /g
GSA Data Repository 2018160 Beall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/10.1130/g39943.1. Numerical Modeling Methodology Stokes
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu Univ. of Texas at Austin January 12, 2007 1 Introduction 2 Weighted Residuals 3 Poisson Equation 4 Other Examples 5 Essential
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationWhat is the Rate-Limiting Step In furthering our Understanding of Subduction Dynamics?
What is the Rate-Limiting Step In furthering our Understanding of Subduction Dynamics? Magali Billen - Katrina Arredondo, - Juan Rodriguez, - Erin Burkett, --Margerete Jadamec Processes in Series are Rate-Limited
More information1 Scaling analysis and non-dimensional numbers
1 Scaling analysis and non-dimensional numbers While this is a textbook on numerical analysis, it is crucial to keep the nature of the physical processes which we would like to model in mind. This will
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationGeodynamics of MARGINS
Geodynamics of MARGINS Scott King (Virginia Tech), Peter van Keken (Univ. Michigan), Mike Gurnis (Caltech), Marc Spiegelman (Lamont), Shijie Zhong (Univ. Colorado), Magali Billen (U.C. Davis) contact:
More informationFinding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems
Finding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems AMSC 663-664 Final Report Minghao Wu AMSC Program mwu@math.umd.edu Dr. Howard Elman Department of Computer
More informationES265 Order of Magnitude Phys & Chem Convection
ES265 Order of Magnitude Phys & Chem Convection Convection deals with moving fluids in which there are spatial variations in temperature or chemical concentration. In forced convection, these variations
More informationFinite element models for coupled magma/mantle dynamics with anisotropic viscosity
Finite element models for coupled magma/mantle dynamics with anisotropic viscosity Sjoerd Hack October 3, 2014 1 Introduction Partially molten rocks consist of a matrix of solid grains with magma flowing
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationNumerical Simulation of the Thermal Convection and Subduction Process in the Mantle
Chapter 1 Earth Science Numerical Simulation of the Thermal Convection and Subduction Process in the Mantle Project Representative Yoshio Fukao Institute for Research on Earth Evolution, Japan Agency for
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationComposing Nonlinear Solvers
Composing Nonlinear Solvers Matthew Knepley Computational and Applied Mathematics Rice University MIT Aeronautics and Astronautics Boston, MA May 10, 2016 Matt (Rice) PETSc MIT 1 / 69 What is PETSc? PETSc
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationLaminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid
Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationIterative Methods for Incompressible Flow
Iterative Methods for Incompressible Flow Melanie McKay Thesis submitted to the Faculty of Graduate and Postdoctoral Studies in partial fulfillment of the requirements for the degree of Master s of Science
More informationarxiv: v1 [physics.geo-ph] 23 Oct 2018
submitted to Geophys. J. Int. A new formulation for coupled magma/mantle dynamics arxiv:1810.10105v1 [physics.geo-ph] 23 Oct 2018 Juliane Dannberg 1, Rene Gassmöller 1, Ryan Grove 2, Timo Heister 3 1 University
More informationAn Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD
An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface
More informationCHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION
CHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION OBJECTIVE The objective of the experiment is to compare the heat transfer characteristics of free and forced convection.
More informationarxiv: v2 [physics.geo-ph] 4 Sep 2009
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, 1 2 A Multiscale Model of Partial Melts 2: Numerical Results G. Simpson arxiv:0903.0437v2 [physics.geo-ph] 4 Sep 2009 3 4 5 6 7 8 Department
More informationTectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop
.. Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Politecnico di Milano - MOX, Dipartimento di
More informationB.1 NAVIER STOKES EQUATION AND REYNOLDS NUMBER. = UL ν. Re = U ρ f L μ
APPENDIX B FLUID DYNAMICS This section is a brief introduction to fluid dynamics. Historically, a simplified concept of the boundary layer, the unstirred water layer, has been operationally used in the
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationINDIAN INSTITUTE OF TECHNOOGY, KHARAGPUR Date: -- AN No. of Students: 5 Sub. No.: ME64/ME64 Time: Hours Full Marks: 6 Mid Autumn Semester Examination Sub. Name: Convective Heat and Mass Transfer Instructions:
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationNumerical Simulation of the Mantle Convection and Subduction Process
Numerical Simulation of the Mantle Convection and Subduction Process Project Representative Yoshio Fukao Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology
More informationEntropic Evaluation of Dean Flow Micromixers
COMSOL Conference, Boston, 2013 Brian Vyhnalek, Petru S. Fodor and Miron Kaufman Physics Department Cleveland State University Entropic Evaluation of Dean Flow Micromixers ABSTRACT We study the use of
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationThe Use of CFD Simulations in Learning Fluid Mechanics at the Undergraduate Level
Presented at the COMSOL Conference 2009 Boston The Use of CFD Simulations in Learning Fluid Mechanics at the Undergraduate Level Marc K. Smith Woodruff School of Mechanical Engineering Georgia Institute
More informationHEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1
HEAT TRANSFER BY CONVECTION Dr. Şaziye Balku 1 CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in the
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationStudies of Arc Volcanism and Mantle Behavior in Subduction Zones
1 Studies of Arc Volcanism and Mantle Behavior in Subduction Zones Using COMSOL Lee, Changyeol Faculty of fearth and denvironmental lsciences Chonnam National University, Gwangju, Republic of Korea Sponsor:
More informationIterative calculation of the heat transfer coefficient
Iterative calculation of the heat transfer coefficient D.Roncati Progettazione Ottica Roncati, Ferrara - Italy Aim The plate temperature of a cooling heat sink is an important parameter that has to be
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More information6.2 Modeling of Systems and Components
Chapter 6 Modelling of Equipment, Processes, and Systems 61 Introduction Modeling is one of the most important elements of thermal system design Most systems are analyzed by considering equations which
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationNumerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Numerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux Ekkehard Holzbecher Weierstrass Institute for Applied
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationDynamic processes in the mantle wedge
Mantle Processes and Geodynamics Dynamic processes in the mantle wedge Ikuko Wada Woods Hole Oceanographic Institution Something Amazing: Despite the apparent complexity of the mantle wedge dynamics, the
More informationRAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8
HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More information4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium
Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4-1 4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4.1 Problem Statement The stress and pore pressure changes due to the expansion
More informationModeling Myths. Modeling Facts
Spiegelman: Columbia University, February 8, 22 Modeling Myths Modeling is Difficult Numerical models can do everything Numerical models will solve your problems Modeling Facts Modeling is simple book-keeping
More informationToward less synchronous composable multilevel methods for implicit multiphysics simulation
Toward less synchronous composable multilevel methods for implicit multiphysics simulation Jed Brown 1, Mark Adams 2, Peter Brune 1, Matt Knepley 3, Barry Smith 1 1 Mathematics and Computer Science Division,
More informationNumerical Methods in Geophysics. Introduction
: Why numerical methods? simple geometries analytical solutions complex geometries numerical solutions Applications in geophysics seismology geodynamics electromagnetism... in all domains History of computers
More informationExpressive Environments and Code Generation for High Performance Computing
Expressive Environments and Code Generation for High Performance Computing Garth N. Wells University of Cambridge SIAM Parallel Processing, 20 February 2014 Collaborators Martin Alnæs, Johan Hake, Richard
More informationApplication of He s homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate
Physics Letters A 37 007) 33 38 www.elsevier.com/locate/pla Application of He s homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate M. Esmaeilpour, D.D. Ganji
More informationTransport equation cavitation models in an unstructured flow solver. Kilian Claramunt, Charles Hirsch
Transport equation cavitation models in an unstructured flow solver Kilian Claramunt, Charles Hirsch SHF Conference on hydraulic machines and cavitation / air in water pipes June 5-6, 2013, Grenoble, France
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationAnswer ALL questions in Section A, and TWO questions from Section B.
UNIVERSITY OF EAST ANGLIA School of Environmental Sciences Main Series Undergraduate Examination 2012-2013 Candidate s no.: GEODYNAMICS: EARTH S ENGINE ENV-2A43 Time allowed: 2 hours. Answer ALL questions
More informationAn Overview of Impellers, Velocity Profile and Reactor Design
An Overview of s, Velocity Profile and Reactor Design Praveen Patel 1, Pranay Vaidya 1, Gurmeet Singh 2 1 Indian Institute of Technology Bombay, India 1 Indian Oil Corporation Limited, R&D Centre Faridabad
More informationA Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationFrictional rheologies have a wide range of applications in engineering
A liquid-crystal model for friction C. H. A. Cheng, L. H. Kellogg, S. Shkoller, and D. L. Turcotte Departments of Mathematics and Geology, University of California, Davis, CA 95616 ; Contributed by D.
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationUnderstanding Hot-Wire Anemometry
Thermal Minutes Understanding Hot-Wire Anemometry Introduction Hot-wire anemometry is a technique for measuring the velocity of fluids, and can be used in many different fields. A hot-wire anemometer consists
More informationNumerical Modelling in Fortran: day 10. Paul Tackley, 2016
Numerical Modelling in Fortran: day 10 Paul Tackley, 2016 Today s Goals 1. Useful libraries and other software 2. Implicit time stepping 3. Projects: Agree on topic (by final lecture) (No lecture next
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationUsing PETSc Solvers in PyLith
Using PETSc Solvers in PyLith Matthew Knepley, Brad Aagaard, and Charles Williams Computational and Applied Mathematics Rice University CIG All-Hands PyLith Tutorial 2016 UC Davis June 19, 2016 M. Knepley
More informationRecycling the lid: Effects of subduction and stirring on boundary layer dynamics in bottom-heated planetary mantle convection
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L20318, doi:10.1029/2006gl027668, 2006 Recycling the lid: Effects of subduction and stirring on boundary layer dynamics in bottom-heated
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationSolution of Partial Differential Equations
Solution of Partial Differential Equations Introduction and Problem Statement Combination of Variables R. Shankar Subramanian We encounter partial differential equations routinely in transport phenomena.
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationCoupled calculations in OpenFOAM -
Coupled calculations in OpenFOAM - Multiphysics handling, structures and solvers, Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg Chalmers University of Technology November 14, 2012 Outline
More informationHeat and Mass Transfer
1 Comments on six papers published by S.P. Anjali Devi and R. Kandasamy in Heat and Mass Transfer, ZAMM, Mechanics Research Communications, International Communications in Heat and Mass Transfer, Communications
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationTransient Heat Transfer Experiment. ME 331 Introduction to Heat Transfer. June 1 st, 2017
Transient Heat Transfer Experiment ME 331 Introduction to Heat Transfer June 1 st, 2017 Abstract The lumped capacitance assumption for transient conduction was tested for three heated spheres; a gold plated
More informationAnalysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications
Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications P. Steinert* 1, I. Schaarschmidt 1, R. Paul 1, M. Zinecker 1, M. Hackert-Oschätzchen 1, Th. Muschalek
More information