Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 418
    • Issues 418
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 7
    • Merge requests 7
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #2329
Closed
Open
Issue created Jan 13, 2022 by Kutalmış Berçin@kutiMaintainer

BUG: turbulentDigitalFilterInlet: Forward-stepwise option is inoperative

Summary

In turbulentDigitalFilterInletFvPatchVectorField.C, the following assignment operation is inoperative:

for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
{
    U.component(dir) =
        U0_.component(dir)*coeffs1FSM_[dir]
      + U.component(dir)*coeffs2FSM_[dir];
}

The LHS U.component(dir) = is unassignable; therefore, the content of U has never been altered by the intended RHS operations.

For this reason, DFM was acting as a single-cell box rather than the FSM simplifications are in effect.

This situation should have adverse effects on integral scale syntheses (which is luckily not considerably important), and should not have any effect on Reynolds stresses or mean quantities.

Environment information

base0    = base
base1    = develop
api      = 2112
patch    = 0
HEAD     = befbcfce24
version  = com
compiler = Clang (system)
         = clang version 9.0.1 
mpi      = SYSTEMOPENMPI
         = mpirun (Open MPI) 1.10.7.0.5e373bf1fd
OS       = Description:     openSUSE Leap 15.3
opts     = linux64ClangDPInt32Opt

Possible fixes

    for (direction dir = 0; dir < pTraits<vector>::nComponents; ++dir)
    {
        U.replace
        (
            dir,
            U0.component(dir)*coeffs1FSM_[dir] + U.component(dir)*coeffs2FSM_[dir]
        );
    }
Edited Jan 13, 2022 by Kutalmış Berçin
Assignee
Assign to
Time tracking