CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 3, 2021, 10:58
Default Generation of velocity field by solving LNSE around base flow generated by icoFoam
  #1
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)
Code:
fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phiB, U)
          + fvc::div(phi, UB)
          - fvm::laplacian(nu, U)
         );
In createFields.H, added the following

Code:
Info<< "Reading field UB\n" << endl;
volVectorField UB
(
    IOobject
    (
        "UB",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

surfaceScalarField phiB
(
    IOobject
    (
        "phiB",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    fvc::flux(UB)
);
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

Reply

Tags
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 17:36
chtMultiRegionSimpleFoam: maximum number of iterations excedeed. Nkl OpenFOAM Running, Solving & CFD 19 October 10, 2019 02:42
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 13:40
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 18:57
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 06:37


All times are GMT -4. The time now is 05:36.