[wayne14]

Hi everyone,

I am working on a 1D steady state darcy flow through different porous materials. The solver reads like this:

    while (simple.loop())
    {
        Info<< "Time = " << runTime.timeName() << nl << endl;

        while (simple.correctNonOrthogonal())
        {
            solve
            (
                fvm::laplacian(k/mu, p)
            );
        }


		U = -k/mu*fvc::grad(p);


        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

The permeability k (four different porous materials) along the distance is shown in the first figure. The piecewise linear distribution of pressure p in the second figure seems reasonable to me. But

U = -k/mu*fvc::grad(p);

which causes velocity pulses at the interface of adjacent porous materials in the third figure. Is this result physical?

Please give me some advice.

Thank you!
Yan

[wayne14]

Originally Posted by wayne14

Hi everyone,

I am working on a 1D steady state darcy flow through different porous materials. The solver reads like this:

    while (simple.loop())
    {
        Info<< "Time = " << runTime.timeName() << nl << endl;

        while (simple.correctNonOrthogonal())
        {
            solve
            (
                fvm::laplacian(k/mu, p)
            );
        }


        U = -k/mu*fvc::grad(p);


        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

The permeability k (four different porous materials) along the distance is shown in the first figure. The piecewise linear distribution of pressure p in the second figure seems reasonable to me. But

U = -k/mu*fvc::grad(p);

which causes velocity pulses at the interface of adjacent porous materials in the third figure. Is this result physical?

Please give me some advice.

Thank you!
Yan

Problem solved.

The velocity should be reconstructed from surface flux instead of -k/mu*fvc::grad(p), and fvSchemes need to be set carefully.