[Samarth_2412]

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::ddt(U)
          + 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
(
    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