CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Generation of velocity field by solving LNSE around base flow generated by icoFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   May 3, 2021, 11:58
Default Generation of velocity field by solving LNSE around base flow generated by icoFoam
New Member
Samarth Acharya
Join Date: May 2021
Posts: 2
Rep Power: 0
Samarth_2412 is on a distinguished road
Hi all,
I am trying to study the instability of boundary layer flows by modifying icoFoam solver to create the LNSE solver, I have made the following changes to the standard icoFoam solver, where U is the perturbation term and UB the base flow velocity

dirIcoFoam.c: (modified part in icoFoam)
fvVectorMatrix UEqn
          + fvm::div(phiB, U)
          + fvc::div(phi, UB)
          - fvm::laplacian(nu, U)
In createFields.H, added the following

Info<< "Reading field UB\n" << endl;
volVectorField UB

surfaceScalarField phiB
I have successfully compiled the modified solver and also generated results of the base flow using icoFoam at the required reynolds number.
Current methodology to approach the problem:
  1. generate base flow using icoFoam
  2. use the solutions obtained at the final time step as initial file for "UB" in the "0" directory for the LNSE solver run
  3. analysis of the perturbed flow "U"
My confusions: I am confused about step 2 in my methodology, as I'm unsure whether it will be acceptable to simply take the final step in the base flow sim and neglect the previous time steps, or should i be looking to incorporate the previous time steps as well, and if so, can anyone guide me how to execute this?
Also, i would appreciate if anyone found an error in my code and made me aware of it!

Thank you
Samarth_2412 is offline   Reply With Quote


boundary layer flow, icofoam, instability, linearized navier stokes, perturbation

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with chtMultiregionFoam radiation boundary condition baran_foam OpenFOAM Running, Solving & CFD 10 December 17, 2019 18:36
chtMultiRegionSimpleFoam: maximum number of iterations excedeed. Nkl OpenFOAM Running, Solving & CFD 19 October 10, 2019 03:42
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 14:40
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 19:57
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37

All times are GMT -4. The time now is 03:58.