View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0001795OpenFOAMBugpublic2015-07-23 08:022015-07-25 01:03
Reporteruser1044Assigned Tohenry  
PrioritylowSeverityminorReproducibilityalways
Status closedResolutionno change required 
PlatformLinuxOSCentOSOS Version6.5
Summary0001795: Explode channel395 with vectorSemiImplicitSource
DescriptionI did Allrun of the tutorial of channel395 with vectorSemiImplicitSource in stead of pressureGradientExplicitSource. The calculation of the cectorSemiImplict has been bursted at Time=7.0 with timeStep=0.2.
Steps To Reproduce1.Copy the channel395 to user-dir.

2.Change the fvOption as follows,
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version 2.0;
    format ascii;
    class dictionary;
    location "system";
    object fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

momentumSource
{
    type vectorSemiImplicitSource;
    active on; //on/off switch
    selectionMode all; //cellSet // points //cellZone

    vectorSemiImplicitSourceCoeffs
    {
        volumeMode specific;
        injectionRateSuSp
        {
            U (( 0.1335 0 0 ) 0);
        }
    }
}

3.Just Allrun.

4.The flow keeps developing as the following figure of channel395_semiImplicitSource_Time6.8.png

5.Terminate the calculation at the next time step of Time=7.0 as the uploaded log file.
Additional InformationI might misunderstand the specification of vectorSemiImplicitSource, in which the term of Su of injectionRateSuSp is the same as Ubar in pressureGradientExplicitSource, that is, I might have to change the Su to lower value. If then, I'd like you to tell me what value of Su.
On the otherhand, I found the member function of fv::SemiImplicitSource<Type>::correct not implemented in SemiImplicitSource.C. In that case, the "do nothing" function of correct to be implemented in fvOption.C must be called, so the fvOptions.correct(U) does nothing in pimpleFaom solver. Is this correct? the "SemiImplict" means the correction of U is not needed?
Tagscyclic, fvOptions, vectorSemiImplicitSource

Activities

user1044

2015-07-23 08:02

 

channel395_semiImplicitSource_Time6.8.png (180,099 bytes)   
channel395_semiImplicitSource_Time6.8.png (180,099 bytes)   

user1044

2015-07-23 08:05

 

log.pimpleFoam (43,086 bytes)

user1044

2015-07-23 08:12

  ~0005123

Last edited: 2015-07-23 11:34

 I uploaded the same figure of pressureGradientExplicitSource with the same legend-scale at Time=308, in which the flow is fully developed. As you can see, Ux exceeds the term of Su, 0.1335, quite much. In order to confirm the volumeMode of specific is correct, I changed the mode to absolute, but the courant number exceeds over up to 2.6, Ux also exceeds 12.0, so the calculation has been bursted.

user1044

2015-07-23 08:13

 

channel395_pressureGradientExplicitSource.png (275,799 bytes)

user1044

2015-07-23 12:16

  ~0005125

Last edited: 2015-07-23 12:20

 I tried the calculation of the same model with simpleFoam of kOmegaSST... Its calculation is supposed to be correct up to 5000 iteration, as the uploaded png file, channel395_specific_simpleFoam_with_SST.png. So, I'm sure the term of Su of injectionRateSuSp in vectorSemiImplicitSource is exactly the same as Ubar of pressureGradientExplicitSource. It's funny the result of steady analysis of simpleFoam is supposed to be unsteady rather than that of pimpleFoam.

user1044

2015-07-23 12:16

 

channel395_specific_simpleFoam_with_SST.png (60,039 bytes)   
channel395_specific_simpleFoam_with_SST.png (60,039 bytes)   

henry

2015-07-23 17:22

manager   ~0005129

This is a request for user support rather than a bug report.

Issue History

Date Modified Username Field Change
2015-07-23 08:02 user1044 New Issue
2015-07-23 08:02 user1044 File Added: channel395_semiImplicitSource_Time6.8.png
2015-07-23 08:05 user1044 Tag Attached: cyclic
2015-07-23 08:05 user1044 Tag Attached: fvOptions
2015-07-23 08:05 user1044 Tag Attached: vectorSemiImplicitSource
2015-07-23 08:05 user1044 File Added: log.pimpleFoam
2015-07-23 08:12 user1044 Note Added: 0005123
2015-07-23 08:13 user1044 File Added: channel395_pressureGradientExplicitSource.png
2015-07-23 11:34 user1044 Note Edited: 0005123
2015-07-23 12:16 user1044 Note Added: 0005125
2015-07-23 12:16 user1044 File Added: channel395_specific_simpleFoam_with_SST.png
2015-07-23 12:20 user1044 Note Edited: 0005125
2015-07-23 17:22 henry Note Added: 0005129
2015-07-23 17:22 henry Status new => closed
2015-07-23 17:22 henry Assigned To => henry
2015-07-23 17:22 henry Resolution open => no change required