Login / Get an account Logout

sandbox/ghigo/README

This ‘README’ provides on overview of the content of my sandbox.

Blood Flow Modeling

Code

• 1D blood flow model with three 1st-order well-balanced fluxes
• 1D blood flow model with 2nd-order hydrostatic reconstruction
• Adaptive mesh refinement compatible with hydrostatic reconstruction

1D blood flow examples

The following examples are compatible with the three 1D blood flow models.

• Inviscid wave propagation in an artery
• Inviscid tourniquet
• Inviscid analytic solution by O. Delestre
• Viscous wave propagation
• Thacker solution - Oscillations in an aneurysm
• Steady solution in a stenosis
• Inviscid wave propagation in a tapered artery

1D blood flow examples with 2nd-order hydrostatic reconstruction and mesh adaptation

The following examples are specific to the 1D blood flow model with 2nd-order hydrostatic reconstruction.

• Inviscid wave propagation in an artery:
  • 1st-order;
  • 2nd-order;
• Inviscid tourniquet:
  • 1st-order;
  • 2nd-order;
  • 1st-order with adaptivity;
• Inviscid analytic solution by O. Delestre:
  • 1st-order;
  • 2nd-order;
• Thacker solution - Oscillations in an aneurysm:
  • 1st-order;
  • 2nd-order;
• Steady solution in a stenosis:
  • 1st-order;
  • 2nd-order;
  • 1st-order with adaptivity;
  • 2nd-order with adaptivity;
    
• Inviscid wave propagation in a tapered artery:
  • 1st-order;
  • 2nd-order;
  • 1st-order with adaptivity;
  • 2nd-order with adaptivity;

3D blood flow example in a network of rigid arteries using embedded boundaries

• Incompressible flow at Re \approx 10 in a network of rigid arteries:
  • steady;
  • pulsatile.

Particle-Laden Flows Using A Cartesian Grid Embedded Boundary Method

Code

Embedded boundaries

• myembed.h: modified embed.h.
• myembed-tree-moving.h: modified embed-tree.h.

Navier-Stokes equations

• mypoisson.h: modified poisson.h to take into account non-zero face-velocity uf boundary conditions on embedded boundary.
• myembed-moving.h: functions and hooks to the centered.h events to account for moving embedded boundaries.

Advection-diffusion of a tracer field

• mytracer.h: modified tracer.h.
• mydiffusion.h: modified diffusion.h.

Fluid-particle coupling

• myembed-particle.h: functions and hooks to the centered.h events to account for moving particles.

Test cases

Poisson equation

• Poisson equation on a 2D complex domain with Neumann boundary conditions;
• Poisson equation on a 3D complex domain with Neumann boundary conditions;

Heat equation

• Heat equation in a 3D sphere with time-dependent Neumann boundary conditions;

Stokes equations

• Poiseuille flow in a periodic domain
• Poiseuille flow in a periodic channel inclined at 45 degrees
• Couette flow between rotating cylinders
• Wannier flow between rotating excentric cylinders
• Flow past a periodic array of cylinders
• Flow past a periodic array of spheres
• Stokes flow through a complex 2D porous medium
• Stokes flow through a complex 3D porous medium
• Stokes flow through a complex 2D porous medium, randomly refined
• Flow past cylinder moving at the same speed:
  • fixed embedded boundaries, uniform mesh
  • moving embedded boundaries, uniform mesh
• Stability of the embedded face velocity interpolation
• Couette flow past a sphere
• Sphere moving (but fixed) towards a wall in a quiescent fluid
• Spheroid moving (but fixed) parallel to a wall in a quiescent fluid

Navier-Stokes equations

• Hydrostatic balance with refined embedded boundaries
• Hydrostatic balance with refined embedded boundaries in 3D
• Stream flow past cylinder moving at the same speed:
  • fixed embedded boundaries, uniform mesh
  • moving embedded boundaries, uniform mesh
• Flow past a circular cylinder for different Reynolds numbers Re
• Flow past a spheroid for different Reynolds numbers Re
• Vortex shedding behind different cylinders for different Reynolds numbers Re
• Flow past a confined circular cylinder for different Reynolds numbers Re
• Starting flow around a cylinder at the Reynolds number Re=1000:
  • fixed embedded boundaries, adaptive mesh
  • moving embedded boundaries, adaptive mesh
• Rotating cylinder in a steam flow at the Reynolds number Re=200
• Rotating sphere in a steam flow at the Reynolds number Re=200
• Accelerating cylinder in a quiescent fluid at the Reynolds number Re=300
• Accelerating sphere in a quiescent fluid at the Reynolds number Re=300
  
• [Oscillating cylinder in a quiescent fluid at the Reynolds number Re=100 and KC=5] (src/test-navier-stokes/cylinder-oscillating.c)
• Oscillating sphere in a quiescent fluid at the Reynolds number Re=40 and Strouhal number St=3.2
  
• Vertically oscillating cylinder in a stream flow at Re=185
• Starting vortex of a NACA2414 airfoil:
  • 2D: naca2414-starting.c
  • 3D: naca2414-starting3D.c
  • correction: naca-corrected.c
• Pitching NACA0015 airfoil
• Incompressible flow at the Reynolds number Re \approx 10 in a network of rigid arteries:
  • steady;
  • pulsatile.
• Flow past a wind turbine at the Reynolds number Re=1000

Freely moving heavy particles

• Buoyant cylinder advected in a Stokes flow
• Buoyant cylinder advected in an inviscid flow
• Heavy cylinder advected by a pressure-driven flow at a Reynolds number Re=20
• Heavy sphere advected by a pressure-driven flow at a Reynolds number Re=20
• Sphere of near-unity density settling in a closed box at low to moderate Reynolds numbers
• Heavy sphere settling in a large box at high Reynolds or Galileo numbers:
  • sphere-settling-large-domain.c;
  • sphere-settling-large-domain-2.c;
• Cylinder settling in a channel;
• Cylinder drifting in a channel;

Viscoplastic Flows

Code

• myviscosity.h: modified viscosity.h to account for embedded boundaries and variable viscosity.
• myviscosity-viscoplastic.h: regularization method for yield stress fluids.

Test cases

• Stokes flow in a lid-driven cavity at different Bingham numbers Bi;
• Poiseuille flow in a channel;
• Stokes flow past a cylinder at different Bingham numbers Bi;

Potential Bugs?

• Spurious velocity at intersections of non-periodic faces in rotational stokes flow

References

inspired by gitit, powered by darcsit