Direct numerical simulations of convective flows Oceanography - Gas - PowerPoint PPT Presentation

Direct numerical simulations of convective flows Oceanography - Gas Transfer Sam Fredriksson, University of Gothenburg

Outline of this presentation • Aim with this presentation • Background • Route forward • Work so far Carbon Cycle – IPCC (2007)

How are CO 2 fluxes determined? Typically: 1.CO 2 concentration in surface water (pCO 2 ) 2.Flux model: F = k ( C aq – C eq ) 3.Piston velocity/gas transfer koefficient ( k ) calculated from average wind speed at 10 m height.

D. Bade 2009, Encycl. of Inland Waters

Surfactant influence of flow field and gas transfer Clean Oley alcohol II k as a funct. of time Surface velocity (surfactant concentration cholesterol) Surface vorticity Surface divergence McKenna 2000

Gas exchange research strategy Forming a group of researchers Gas exchange and concentration measurements Wind measurements Turbulence measurements Infrared camera recording Numerical modeling

Route forward numerical simulations General • DNS (Direct Numerial Simulations) with OpenFoam • DNS with convection • DNS with surface shear • DNS with surfactant • LES (Large Eddy Simulations) • Waves …

Route forward numerical simulations Lars Davidsson et al. ”DNS of plane vertical channel with and without buoyancy”, Turbulence, Heat and General Mass Transfer, 2003 • DNS with OpenFoam 1.4 1.2 • DNS with convection 1 0.8 • DNS with surface shear Umean/Ubar Umean 0.6 _Lada/Ubar • DNS with surfactant 0.4 0.2 • LES/DNS moving surface 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 • Waves … 0.2 0.18 0.16 0.14 0.12 urms/Ubar vrms/Ubar 0.1 wrms/Ubar urms_Lada/Ubar 0.08 urms_Lada/Ubar wrms_Lada/Ubar 0.06 0.04 0.02 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

Route forward numerical simulations General Robert Handler: ”Direct numerical simulations of free convection • DNS with OpenFoam beneath an air-water interface at low • DNS with convection Rayleigh numbers”, Physics of Fluids, 2003 • DNS with surface shear • DNS with surfactant • LES/DNS moving surface • Waves …

Computational domain Evaporation

Boussinesq formulation Momentum and Continuity �� U ∙ � U � � 1 �� ν � � � � �� ∙ U � 0 Linear approximation of density � �� 1 � � � � � � , � � � �� Temperature field is governed by a convection-diffusion equation given by �� U ∙ � T � � � � T � S

Temperature field, Ra q ≈ 5x10 6 0.076 m Played 10 times real time

IR-measurements at Bornö L z = 0.122 (Side = 0.17 m) Ra ~ QL z 4 ≈ 1 m Magnus Gålfalk

Surfactant - Finite Area Method - fam

Areas for investigation Analysis • Mesh convergence – it seems that there is a need for finer meshes with present formulation compared to Handler et al. • Size of computational box and mesh resolution • More thorough analysis and comparison against experiments etc. • Surfactant formulation, cyclic boundaries for fam • Find correlation between IR and gas transfer • Find a good measure/representation of the temperature field (rms, skewness, flatness, divergence etc.) • Waves etc.

Thank you

Direct numerical simulations of convective flows Oceanography - Gas - PowerPoint PPT Presentation

Direct numerical simulations of convective flows Oceanography - Gas Transfer Sam Fredriksson, University of Gothenburg Outline of this presentation Aim with this presentation Background Route forward Work so far Carbon Cycle

DNS of Multiphase Flows Simple Front Tracking DNS of Multiphase Flows Direct Numerical

Governing Equations 4. The governing equations are mathematical statements of the physical

Volume of Fluid j- 1/2 j+ 1/2 H ( x ) C j I ( x ) ( x ) j -1 j j +1 Interface cell DNS of

Corrected Velocities 3. The incompressible velocity at the new time, u n+1, is the predicted

Surface a small segment of the front. For two-dimensional flows we can use Tension the

NUMERICAL SIMULATIONS OF VISCOELASTIC FLOWS IN FIBROUS POROUS MEDIA H. L. Liu and W. R. Hwang*

Surface Tension 4-1. The key thing about surface tension is that it is a tension. The force is

Governing not change, we have incompressible flow Equations Conservation of momentum. For

A Falling Drop scale. Symbols for the various physical quantities follow the usual convention.

Constructing density. the marker function 4. There are several ways to construct a marker

The boundary Flows-3 conditions A Simple Solver for Variable Density Flow (3 of 3) Gretar

Direct Numerical Simulation of Turbulent Wall Flows at Constant Power Input Y. Hasegawa1, B.

A Domain Decomposition Method for Large Scale Simulations of Two-phase Flows with Moving Contact

Direct Numerical Simulation Large Eddy Simulation TURBULENT FLOWS AND INHERENT STRUCTURES Tony

Direct Numerical Simulation of bubbles with Adaptive Mesh Refinement (AMR) on parallel

1. One of the major challenges in simulating multiphase flows is the advection of the marker

Assimilating visible satellite images for convective scale numerical weather prediction Leonhard

Direct Assimilation of Radar Reflectivity Data using a Convective scale EnKF System Jingyao Luo 1

ON NUMERICAL UPSCALING FOR STOKES AND ON NUMERICAL UPSCALING FOR STOKES AND STOKES- -BRINKMAN

1. When simulating multiphase flows on fixed grids we must update the density and viscosity

Numerical Simulations of CO 2 Geo-Sequestration using PETSc Henrik B using Institute for

1. Thank you for your interest in these lectures. Before getting fully immersed in the

Evolution of convective of convective cloud cloud top top height height: : Evolution

Energy flows and modelling approaches Energy flows in buildings external long-wave local