Main Page: Difference between revisions

From openpipeflow.org
Jump to navigation Jump to search
 
(91 intermediate revisions by the same user not shown)
Line 1: Line 1:
{{#ev:youtube|http://youtu.be/hSj0VnmSQro|280|right|Slow streaks (blue) and vortex structures (yellow)}}
<!--
[[File:6551noshiftloop.mp4|280|right|Slow streaks (blue) and vortex structures (yellow)]]
[[File:OpenPFlogo1.png|thumb|right|Slow streaks (blue) and vortex structures (yellow)]]
-->
==== Quicklinks ====
[[Download]] | [[Manual]] | [[Tutorial]] | [[Database]]
==Overview==
==Overview==
* Primitive variable pipe flow code for incompressible flow.
'''openpipeflow.org is a free resource for researchers, engineers, educators and the interested public.'''
* Double-Fourier (theta,z) + finite difference (r).
 
Pipe flow is a simple and familiar set up, yet the flow patterns exhibit [[Fun_stuff|rich chaotic dynamics]].  This provides a setting for investigating the principles of simulation at one level, and at another, for developing new methods designed to probe fundamental properties of dynamical systems. 
 
The majority of mathematical techniques described on these pages are applicable to a huge range of problems, and [[Manual|subroutines for well-known methods]] are designed to be callable from any code.  The core [[Manual|pipe flow code]] is designed to be flexible yet very fast.
 
'''Aims'''
* To make accessible a range of modelling techniques.
* To facilitate rapid entry into the world of numerical simulation and fluid dynamics.
* To provide flexible modules for more the use and development of advanced techniques in research.
 
'''Code features'''
* Primitive-variable pipe-flow code for incompressible flow.
* Simple scripts for visualisation with Matlab/Octave/Visit.
* Readable Fortran 90, uses modules and derived types, no esoteric extensions.
* Core program <3000 lines.
* Spatial discretisation: double-Fourier (theta,z) + finite difference (r).
* PPE formulation; influence matrix corrects boundary conditions to machine precision.
* PPE formulation; influence matrix corrects boundary conditions to machine precision.
* Second-order predictor-corrector method, automatic timestep control.
* Second-order predictor-corrector method, automatic timestep control.
* Serial or parallel (with MPI).
* May be run on a single core or in parallel (with MPI).  Essentially linear scaling with number of cores.
* Fortran 90, uses modules and derived types, no esoteric extensions.
* '2-dimensional' parallelisation, radial+axial split.
* Core program <3000 lines.
* Jacobian-Free Newton-Krylov (JFNK) solver.  Now at https://github.com/apwillis1/JFNK-Hookstep .
Details for the above can be found in the [[Manual]].
This article, [[File:TheOpenpipeflowSolver.pdf]], provides an overview of the code and its context.
 
==Database==
The [[Database]] provides sample parameters and initial conditions from which to launch new simulations.  In general, simulations start most reliably from an initial state computed for similar parameters.  A range of starting points are provided.
 
==Features to appear/wishlist==
* Module for the immersed boundary method (IBM).
* More FAQ + documentation.
<!--
* tutorial, inc. matlab example.
* Utilities for Krylov methods -- Newton-Krylov, Arnoldi.
* phys-statefiles for direct upload to matlab/visit etc.
 
==Related codes==
Closely related codes are currently obtained by request.
* Taylor--Couette flow.
* Flow in a periodic box.
* Pipeflow, potential formulation.


==Version 1.02b==
==Problem-independent codes==
* Documentation [[Media:Openpipeflow-1.02b-doc.pdf]]
To be uploaded
* Tarball [[Media:Openpipeflow-1.02b.tgz]]
* GMRES
* Newton-Hookstep
* Arnoldi
-->


==Citation==
==Citation==
*        In talks:
Please cite this article: [[File:TheOpenpipeflowSolver.pdf]]
                openpipeflow.org
*        In articles:
                1) URL above if possible and
                2) "based on code described in \cite{WK09}."
                  @article{WK09,
                        Author = {A. P. Willis and R. R. Kerswell},
                        Year = {2009},
                        Title = {Turbulent dynamics of pipe flow captured in a
                                reduced model: puff relaminarisation and
                                localised `edge' states},
                        Journal = {J.\ Fluid Mech.},
                        Volume = {619},
                        Pages = {213-233}  }


==Author==
==Author==
:      Ashley P. Willis,
:      [http://maths.dept.shef.ac.uk/maths/staff_info_402.html Ashley P. Willis],
:      School of Mathematics and Statistics,
:      [http://www.sheffield.ac.uk/maths School of Mathematics and Statistics (SoMaS)],
:      University of Sheffield,
:      [http://www.sheffield.ac.uk/maths University of Sheffield, U.K.]
:      South Yorkshire,
:      S3 7RH,
:      U.K.
:      a.p.willis/at/sheffield.ac.uk
:      a.p.willis/at/sheffield.ac.uk
==Thanks==
* John Gibson ([http://channelflow.org/gibson channelflow.org])
* Predrag Cvitanović ([http://www.cns.gatech.edu/~predrag/ GaTech] [http://chaosbook.org chaosbook.org])
* Rich Kerswell ([http://www.damtp.cam.ac.uk/user/rrk26/ Cambridge])
* many other people!
* EPSRC GR/S76144/01, EP/K03636X/1
* The [http://www.sheffield.ac.uk/maths University of Sheffield].

Latest revision as of 13:24, 27 March 2024

Slow streaks (blue) and vortex structures (yellow)

Quicklinks

Download | Manual | Tutorial | Database

Overview

openpipeflow.org is a free resource for researchers, engineers, educators and the interested public.

Pipe flow is a simple and familiar set up, yet the flow patterns exhibit rich chaotic dynamics. This provides a setting for investigating the principles of simulation at one level, and at another, for developing new methods designed to probe fundamental properties of dynamical systems.

The majority of mathematical techniques described on these pages are applicable to a huge range of problems, and subroutines for well-known methods are designed to be callable from any code. The core pipe flow code is designed to be flexible yet very fast.

Aims

  • To make accessible a range of modelling techniques.
  • To facilitate rapid entry into the world of numerical simulation and fluid dynamics.
  • To provide flexible modules for more the use and development of advanced techniques in research.

Code features

  • Primitive-variable pipe-flow code for incompressible flow.
  • Simple scripts for visualisation with Matlab/Octave/Visit.
  • Readable Fortran 90, uses modules and derived types, no esoteric extensions.
  • Core program <3000 lines.
  • Spatial discretisation: double-Fourier (theta,z) + finite difference (r).
  • PPE formulation; influence matrix corrects boundary conditions to machine precision.
  • Second-order predictor-corrector method, automatic timestep control.
  • May be run on a single core or in parallel (with MPI). Essentially linear scaling with number of cores.
  • '2-dimensional' parallelisation, radial+axial split.
  • Jacobian-Free Newton-Krylov (JFNK) solver. Now at https://github.com/apwillis1/JFNK-Hookstep .

Details for the above can be found in the Manual. This article, File:TheOpenpipeflowSolver.pdf, provides an overview of the code and its context.

Database

The Database provides sample parameters and initial conditions from which to launch new simulations. In general, simulations start most reliably from an initial state computed for similar parameters. A range of starting points are provided.

Features to appear/wishlist

  • Module for the immersed boundary method (IBM).
  • More FAQ + documentation.

Citation

Please cite this article: File:TheOpenpipeflowSolver.pdf

Author

Ashley P. Willis,
School of Mathematics and Statistics (SoMaS),
University of Sheffield, U.K.
a.p.willis/at/sheffield.ac.uk

Thanks