From ab6aac74e20531556e8d7a2a0f3a637ecfd5a960 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 9 Sep 2013 17:34:37 +0100
Subject: [PATCH] ENH: curvatureSeparation film model - check that gravity is
set/valid
---
.../curvatureSeparation/curvatureSeparation.C | 20 +++++++++++++++++--
1 file changed, 18 insertions(+), 2 deletions(-)
diff --git a/src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/curvatureSeparation/curvatureSeparation.C b/src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/curvatureSeparation/curvatureSeparation.C
index 27376b14d5e..0d8011404cf 100644
--- a/src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/curvatureSeparation/curvatureSeparation.C
+++ b/src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/curvatureSeparation/curvatureSeparation.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -229,8 +229,24 @@ curvatureSeparation::curvatureSeparation
deltaByR1Min_(coeffs().lookupOrDefault<scalar>("deltaByR1Min", 0.0)),
definedPatchRadii_(),
magG_(mag(owner.g().value())),
- gHat_(owner.g().value()/magG_)
+ gHat_(vector::zero)
{
+ if (magG_ < ROOTVSMALL)
+ {
+ FatalErrorIn
+ (
+ "curvatureSeparation::curvatureSeparation"
+ "("
+ "const surfaceFilmModel&, "
+ "const dictionary&"
+ ")"
+ )
+ << "Acceleration due to gravity must be non-zero"
+ << exit(FatalError);
+ }
+
+ gHat_ = owner.g().value()/magG_;
+
List<Tuple2<word, scalar> > prIn(coeffs().lookup("definedPatchRadii"));
const wordList& allPatchNames = owner.regionMesh().boundaryMesh().names();
--
GitLab