41 velocityComponentLaplacianFvMotionSolver,
49 Foam::velocityComponentLaplacianFvMotionSolver::
50 velocityComponentLaplacianFvMotionSolver
63 "pointMotionU" + cmptName_,
64 fvMesh_.time().timeName(),
69 pointMesh::New(fvMesh_)
75 "cellMotionU" + cmptName_,
78 IOobject::READ_IF_PRESENT,
85 pointMotionU_.dimensions(),
88 cellMotionBoundaryTypes<scalar>(pointMotionU_.boundaryField())
92 motionDiffusivity::New(*
this, lookup(
"diffusivity"))
99 else if (cmptName_ ==
"y")
103 else if (cmptName_ ==
"z")
111 "velocityComponentLaplacianFvMotionSolver::"
112 "velocityComponentLaplacianFvMotionSolver"
113 "(const polyMesh& mesh, Istream& msData)"
114 ) <<
"Given component name " << cmptName_ <<
" should be x, y or z"
144 + fvMesh_.time().deltaT().value()*pointMotionU_.internalField()
147 twoDCorrectPoints(tcurPoints());
157 movePoints(fvMesh_.points());
159 diffusivityPtr_->correct();
160 pointMotionU_.boundaryField().updateCoeffs();
166 diffusivityPtr_->operator()(),
168 "laplacian(diffusivity,cellMotionU)"
183 diffusivityPtr_.reset(NULL);