CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Darcy flow through different porous materials

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 6, 2016, 10:29
Default Darcy flow through different porous materials
  #1
Member
 
Yan Wang
Join Date: May 2015
Location: Beijing
Posts: 41
Rep Power: 10
wayne14 is on a distinguished road
Hi everyone,

I am working on a 1D steady state darcy flow through different porous materials. The solver reads like this:
Code:
    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
Code:
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
Attached Images
File Type: png k.png (30.5 KB, 13 views)
File Type: png p.png (34.2 KB, 11 views)
File Type: png U.png (41.0 KB, 12 views)
wayne14 is offline   Reply With Quote

Old   January 7, 2016, 00:11
Default
  #2
Member
 
Yan Wang
Join Date: May 2015
Location: Beijing
Posts: 41
Rep Power: 10
wayne14 is on a distinguished road
Quote:
Originally Posted by wayne14 View Post
Hi everyone,

I am working on a 1D steady state darcy flow through different porous materials. The solver reads like this:
Code:
    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
Code:
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.
wayne14 is offline   Reply With Quote

Reply

Tags
darcy flow, darcy's law, porous

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
Turbulent flow in porous medium Chander CFX 4 January 17, 2017 10:50
multiphase cylindrical flow through porous medium; BC adrianabm OpenFOAM Running, Solving & CFD 1 September 7, 2015 11:37
Flow in porous media + phase change styleworker FLOW-3D 1 June 27, 2012 15:51
rotating porous media in a general flow a_dores FLUENT 0 October 31, 2010 05:50
porous media/jump mass flow problem Marijo FLUENT 0 February 17, 2006 12:49


All times are GMT -4. The time now is 01:23.