GCC Code Coverage Report

Directory: ../../../builds/dumux-repositories/
File: /builds/dumux-repositories/dumux/test/freeflow/navierstokes/channel/pipe/momentum/problem.hh
Date: 2024-09-21 20:52:54
Exec Total Coverage
Lines: 35 45 77.8%
Functions: 6 12 50.0%
Branches: 63 90 70.0%
Line Branch Exec Source
1 // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2 // vi: set et ts=4 sw=4 sts=4:
3 //
4 // SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
5 // SPDX-License-Identifier: GPL-3.0-or-later
6 //
7
8 #ifndef DUMUX_TEST_FREEFLOW_PIPE_PROBLEM_HH
9 #define DUMUX_TEST_FREEFLOW_PIPE_PROBLEM_HH
10
11 #include <dumux/common/parameters.hh>
12 #include <dumux/common/properties.hh>
13
14 namespace Dumux {
15 /*!
16 * \ingroup NavierStokesTests
17 * \brief Freeflow problem for pipe flow
18 * Simulation of a radially-symmetric pipe flow with circular cross-section
19 */
20 template <class TypeTag, class BaseProblem>
21 class FreeFlowPipeProblem : public BaseProblem
22 {
23 using ParentType = BaseProblem;
24 using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
25 using FVElementGeometry = typename GridGeometry::LocalView;
26 using SubControlVolume = typename FVElementGeometry::SubControlVolume;
27 using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace;
28 using GridView = typename GridGeometry::GridView;
29 using Scalar = GetPropType<TypeTag, Properties::Scalar>;
30 using ModelTraits = GetPropType<TypeTag, Properties::ModelTraits>;
31 using Element = typename GridView::template Codim<0>::Entity;
32 using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
33 using InitialValues = typename ParentType::InitialValues;
34 using Sources = typename ParentType::Sources;
35 using DirichletValues = typename ParentType::DirichletValues;
36 using BoundaryFluxes = typename ParentType::BoundaryFluxes;
37 using BoundaryTypes = typename ParentType::BoundaryTypes;
38
39 public:
40 using Indices = typename ModelTraits::Indices;
41
42 6 FreeFlowPipeProblem(std::shared_ptr<const GridGeometry> gridGeometry)
43
7/20
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✓ Branch 8 taken 6 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 6 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 6 times.
✓ Branch 15 taken 6 times.
✗ Branch 16 not taken.
✓ Branch 18 taken 6 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
 18 : ParentType(gridGeometry)
44 {
45
3/6
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
✗ Branch 7 not taken.
✓ Branch 8 taken 6 times.
 12 name_ = getParamFromGroup<std::string>(this->paramGroup(), "Problem.Name");
46
1/2
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
 6 meanInletVelocity_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.MeanInletVelocity");
47
1/2
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
 6 const Scalar nu = getParam<Scalar>("Component.LiquidKinematicViscosity");
48
1/2
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
 6 rho_ = getParam<Scalar>("Component.LiquidDensity");
49 6 mu_ = nu*rho_;
50
51 36 pipeRadius_ = this->gridGeometry().bBoxMax()[0] - this->gridGeometry().bBoxMin()[0];
52 36 pipeLength_ = this->gridGeometry().bBoxMax()[1] - this->gridGeometry().bBoxMin()[1];
53 6 eps_ = 1e-7*pipeRadius_;
54
55
2/4
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 6 times.
✗ Branch 6 not taken.
 12 std::cout << "-- Reynolds number: " << 2*pipeRadius_*meanInletVelocity_/nu << std::endl;
56 6 }
57
58 const std::string& name() const
59 {
60
1/2
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
 12 return name_;
61 }
62
63 /*!
64 * \brief Specifies which kind of boundary condition should be
65 * used for which equation on a given boundary segment.
66 */
67 9216 BoundaryTypes boundaryTypesAtPos(const GlobalPosition& globalPos) const
68 {
69 9216 BoundaryTypes values;
70
71 if constexpr (ParentType::isMomentumProblem())
72 {
73 36024 if (onLowerBoundary_(globalPos) || onOuterBoundary_(globalPos))
74 values.setAllDirichlet();
75 9204 else if (onInnerBoundary_(globalPos))
76 {
77 4194 values.setDirichlet(Indices::velocityXIdx);
78 4194 values.setNeumann(Indices::momentumYBalanceIdx);
79 }
80 else
81 values.setAllNeumann();
82 }
83 else
84 values.setAllNeumann();
85
86 9216 return values;
87 }
88
89 /*!
90 * \brief Evaluates the boundary conditions for a Neumann control volume.
91 *
92 * \param element The element for which the Neumann boundary condition is set
93 * \param fvGeometry The fvGeometry
94 * \param elemVolVars The element volume variables
95 * \param elemFaceVars The element face variables
96 * \param scvf The boundary sub control volume face
97 */
98 template<class ElementVolumeVariables, class ElementFluxVariablesCache>
99 47608 BoundaryFluxes neumann(const Element& element,
100 const FVElementGeometry& fvGeometry,
101 const ElementVolumeVariables& elemVolVars,
102 const ElementFluxVariablesCache& elemFluxVarsCache,
103 const SubControlVolumeFace& scvf) const
104 {
105 47608 BoundaryFluxes values(0.0);
106
2/2
✓ Branch 0 taken 43334 times.
✓ Branch 1 taken 4274 times.
 47608 const auto& globalPos = scvf.ipGlobal();
107
108 if constexpr (ParentType::isMomentumProblem())
109 {
110 95216 if (onInnerBoundary_(globalPos))
111 values = 0.0; // zero shear stress at symmetry axis
112 8548 else if (onUpperBoundary_(globalPos))
113 4274 values.axpy(analyticalPressure(globalPos), scvf.unitOuterNormal());
114 }
115 else
116 {
117 const auto insideDensity = elemVolVars[scvf.insideScvIdx()].density();
118 values[Indices::conti0EqIdx] = this->faceVelocity(element, fvGeometry, scvf) * insideDensity * scvf.unitOuterNormal();
119 }
120
121 47608 return values;
122 }
123
124 DirichletValues dirichletAtPos(const GlobalPosition& globalPos) const
125 8808 { return analyticalSolution(globalPos); }
126
127 ✗ DirichletValues analyticalSolution(const GlobalPosition& globalPos, Scalar time = 0.0) const
128 {
129 ✗ DirichletValues values(0.0);
130
131 // paraboloid velocity profile
132 if constexpr (ParentType::isMomentumProblem())
133 {
134 ✗ const auto r = globalPos[0] - this->gridGeometry().bBoxMin()[0];
135 ✗ values[Indices::velocityXIdx] = 0.0;
136 ✗ values[Indices::velocityYIdx] = 2.0*meanInletVelocity_*(1.0 - r*r/(pipeRadius_*pipeRadius_));
137 }
138 else
139 values[Indices::pressureIdx] = analyticalPressure(globalPos);
140
141 ✗ return values;
142 }
143
144 ✗ Scalar densityAtPos(const GlobalPosition& globalPos) const
145 ✗ { return rho_; }
146
147 ✗ Scalar effectiveViscosityAtPos(const GlobalPosition& globalPos) const
148 ✗ { return mu_; }
149
150 Scalar pressureAtPos(const GlobalPosition& globalPos) const
151 18060000 { return analyticalPressure(globalPos); }
152
153 Scalar analyticalPressure(const GlobalPosition& globalPos) const
154 {
155 5254274 const auto y = globalPos[1];
156 5254274 return (pipeLength_-y)*meanInletVelocity_*8.0*mu_/(pipeRadius_*pipeRadius_);
157 }
158
159 Scalar analyticalMassFlux() const
160 {
161 return meanInletVelocity_ * M_PI * pipeRadius_*pipeRadius_ * rho_;
162 }
163
164 private:
165 bool onInnerBoundary_(const GlobalPosition &globalPos) const
166
20/20
✓ Branch 0 taken 43334 times.
✓ Branch 1 taken 4274 times.
✓ Branch 2 taken 43334 times.
✓ Branch 3 taken 4274 times.
✓ Branch 4 taken 43334 times.
✓ Branch 5 taken 4274 times.
✓ Branch 6 taken 43334 times.
✓ Branch 7 taken 4274 times.
✓ Branch 8 taken 43334 times.
✓ Branch 9 taken 4274 times.
✓ Branch 10 taken 408 times.
✓ Branch 11 taken 4194 times.
✓ Branch 12 taken 408 times.
✓ Branch 13 taken 4194 times.
✓ Branch 14 taken 408 times.
✓ Branch 15 taken 4194 times.
✓ Branch 16 taken 408 times.
✓ Branch 17 taken 4194 times.
✓ Branch 18 taken 408 times.
✓ Branch 19 taken 4194 times.
 261050 { return globalPos[0] < this->gridGeometry().bBoxMin()[0] + eps_; }
167
168 bool onOuterBoundary_(const GlobalPosition &globalPos) const
169
10/10
✓ Branch 0 taken 4194 times.
✓ Branch 1 taken 4602 times.
✓ Branch 2 taken 4194 times.
✓ Branch 3 taken 4602 times.
✓ Branch 4 taken 4194 times.
✓ Branch 5 taken 4602 times.
✓ Branch 6 taken 4194 times.
✓ Branch 7 taken 4602 times.
✓ Branch 8 taken 4194 times.
✓ Branch 9 taken 4602 times.
 43980 { return globalPos[0] > this->gridGeometry().bBoxMax()[0] - eps_; }
170
171 bool onLowerBoundary_(const GlobalPosition &globalPos) const
172
10/10
✓ Branch 0 taken 8796 times.
✓ Branch 1 taken 420 times.
✓ Branch 2 taken 8796 times.
✓ Branch 3 taken 420 times.
✓ Branch 4 taken 8796 times.
✓ Branch 5 taken 420 times.
✓ Branch 6 taken 8796 times.
✓ Branch 7 taken 420 times.
✓ Branch 8 taken 8796 times.
✓ Branch 9 taken 420 times.
 46080 { return globalPos[1] < this->gridGeometry().bBoxMin()[1] + eps_; }
173
174 bool onUpperBoundary_(const GlobalPosition &globalPos) const
175
5/10
✓ Branch 0 taken 4274 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 4274 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4274 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 4274 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 4274 times.
✗ Branch 9 not taken.
 21370 { return globalPos[1] > this->gridGeometry().bBoxMax()[1] - eps_; }
176
177 std::string name_;
178 Scalar initialPressure_;
179 Scalar meanInletVelocity_;
180 Scalar mu_;
181 Scalar rho_;
182 Scalar pipeRadius_, pipeLength_;
183 Scalar eps_;
184 };
185
186 } // end namespace Dumux
187
188 #endif
189