diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/YEqn.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/YEqn.H index 687acd4e70..8f77f23bc2 100644 --- a/applications/solvers/lagrangian/trackedReactingParcelFoam/YEqn.H +++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/YEqn.H @@ -28,6 +28,7 @@ tmp > mvConvection == reactingParcels.Srho(i) + kappa*chemistry.RR(i)().dimensionedInternalField() + + pointMassSources.Su(i) ); Yi.max(0.0); diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/createMulticomponentPointSources.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/createMulticomponentPointSources.H new file mode 100644 index 0000000000..358d7f0a14 --- /dev/null +++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/createMulticomponentPointSources.H @@ -0,0 +1,8 @@ +Info<< "\nConstructing multi-component mass flow rate point sources" << endl; +timeActivatedExplicitMulticomponentPointSource pointMassSources +( + "pointMassSources", + mesh, + Y, + dimMass/dimTime +); diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/trackedReactingParcelFoam.C b/applications/solvers/lagrangian/trackedReactingParcelFoam/trackedReactingParcelFoam.C index 77a67502da..d730ca44b3 100644 --- a/applications/solvers/lagrangian/trackedReactingParcelFoam/trackedReactingParcelFoam.C +++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/trackedReactingParcelFoam.C @@ -37,6 +37,7 @@ Description #include "ReactingCloudThermoTypes.H" #include "radiationModel.H" #include "porousZones.H" +#include "timeActivatedExplicitMulticomponentPointSource.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -51,6 +52,7 @@ int main(int argc, char *argv[]) #include "createFields.H" #include "createRadiationModel.H" #include "createClouds.H" + #include "createMulticomponentPointSources.H" #include "readPISOControls.H" #include "initContinuityErrs.H" #include "readTimeControls.H"