From c7792518b102f35ac060744ba1d1eadee29988d9 Mon Sep 17 00:00:00 2001 From: Henry <Henry> Date: Wed, 4 Feb 2015 23:03:45 +0000 Subject: [PATCH] MRFSource: Update the velocity boundary conditions for changes in rotation speed --- src/fvOptions/sources/derived/MRFSource/MRFSource.C | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/src/fvOptions/sources/derived/MRFSource/MRFSource.C b/src/fvOptions/sources/derived/MRFSource/MRFSource.C index 59fd4ff021f..71343060f93 100644 --- a/src/fvOptions/sources/derived/MRFSource/MRFSource.C +++ b/src/fvOptions/sources/derived/MRFSource/MRFSource.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2012-2015 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -103,6 +103,9 @@ void Foam::fv::MRFSource::addSup const label fieldI ) { + // Update the velocity boundary conditions for changes in rotation speed + mrfPtr_->correctBoundaryVelocity(const_cast<volVectorField&>(eqn.psi())); + // Add to rhs of equation mrfPtr_->addCoriolis(eqn, true); } @@ -115,6 +118,9 @@ void Foam::fv::MRFSource::addSup const label fieldI ) { + // Update the velocity boundary conditions for changes in rotation speed + mrfPtr_->correctBoundaryVelocity(const_cast<volVectorField&>(eqn.psi())); + // Add to rhs of equation mrfPtr_->addCoriolis(rho, eqn, true); } -- GitLab