|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Anders Simonsen
Join Date: Sep 2013
Posts: 18
Rep Power: 12 ![]() |
Hi
I'm facing a problem with implementing a linearized source term for the momentum equations. The U-equation in my solver states (Turbulence+Steady state): Code:
tmp<fvVectorMatrix> tUEqn ( fvm::div(phi, U) + turbulence->divDevRhoReff(U) ); // [N/m^3] Code:
p[0] = p_Su[0] + p_Sp[0] * (U[0] - U_Ref[0]) p[1] = p_Su[1] + p_Sp[1] * (U[1] - U_Ref[1]) p[2] = p_Su[2] + p_Sp[2] * (U[2] - U_Ref[2]) Code:
tmp<fvVectorMatrix> tUEqn ( fvm::div(phi, U) + turbulence->divDevRhoReff(U) == p_Su - fvc::Sp(p_Sp,U) - fvm::SuSp(-p_Sp,U) ); // [N/m^3] The format shown above works for scalar equations, where I've successfully implemented it for scalar species as shown below: Code:
fvScalarMatrix YiEqn ( mvConvection->fvmDiv(phi, Yi) - fvm::laplacian(turbulence->muEff(), Yi) == Yi_Su - fvc::Sp(Yi_Sp,Yi) - fvm::SuSp(-Yi_Sp,Yi) ); // [kg/s-m^3] Can someone help me with my problem? ![]() I'd like to avoid using fvOptions, as I can't specify the source terms for each cell individually. Best regards Anders |
|
![]() |
![]() |
![]() |
Tags |
momentum equation, source term |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
polynomial BC | srv537 | OpenFOAM Pre-Processing | 4 | December 3, 2016 09:07 |
[Other] Adding solvers from DensityBasedTurbo to foam-extend 3.0 | Seroga | OpenFOAM Community Contributions | 9 | June 12, 2015 17:18 |
[swak4Foam] Swak4FOAM 0.2.3 / OF2.2.x installation error | FerdiFuchs | OpenFOAM Community Contributions | 27 | April 16, 2014 15:14 |
friction forces icoFoam | ofslcm | OpenFOAM | 3 | April 7, 2012 10:57 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 17:51 |