Point source (or sink) that adds (removes) fluid at a constant mass flux rate for times between the specified start and end times. More...

#include <PorousFlowSquarePulsePointSource.h>

Inheritance diagram for PorousFlowSquarePulsePointSource:

[legend]

Public Member Functions
	PorousFlowSquarePulsePointSource (const InputParameters &parameters)

virtual void	addPoints () override

virtual Real	computeQpResidual () override

Protected Attributes
const Real	_mass_flux
	The constant mass flux (kg/s) More...

const Point	_p
	The location of the point source (sink) More...

const Real	_start_time
	The time at which the point source (sink) starts operating. More...

const Real	_end_time
	The time at which the point source (sink) stops operating. More...

Detailed Description

Point source (or sink) that adds (removes) fluid at a constant mass flux rate for times between the specified start and end times.

If no start and end times are specified, the source (sink) starts at the start of the simulation and continues to act indefinitely

Definition at line 24 of file PorousFlowSquarePulsePointSource.h.

Constructor & Destructor Documentation

◆ PorousFlowSquarePulsePointSource()

PorousFlowSquarePulsePointSource::PorousFlowSquarePulsePointSource ( const InputParameters & parameters )

Definition at line 32 of file PorousFlowSquarePulsePointSource.C.

   : DiracKernel(parameters),
     _mass_flux(getParam<Real>("mass_flux")),
     _p(getParam<Point>("point")),
     _start_time(getParam<Real>("start_time")),
     _end_time(getParam<Real>("end_time"))
 {
   // Sanity check to ensure that the end_time is greater than the start_time
   if (_end_time <= _start_time)
     mooseError(name(),
                ": start time for PorousFlowSquarePulsePointSource is ",
                _start_time,
                " but it must be less than end time ",
                _end_time);
 }

Member Function Documentation

◆ addPoints()

void PorousFlowSquarePulsePointSource::addPoints ( )

overridevirtual

Definition at line 50 of file PorousFlowSquarePulsePointSource.C.

 {
   addPoint(_p, 0);
 }

◆ computeQpResidual()

Real PorousFlowSquarePulsePointSource::computeQpResidual ( )

overridevirtual

There are six cases for the start and end time in relation to t-dt and t. If the interval (t-dt,t) is only partly but not fully within the (start,end) interval, then the mass_flux is scaled so that the total mass added (or removed) is correct

Definition at line 56 of file PorousFlowSquarePulsePointSource.C.

 {
   Real factor = 0.0;
  
   if (_t < _start_time || _t - _dt >= _end_time)
     factor = 0.0;
   else if (_t - _dt < _start_time)
   {
     if (_t <= _end_time)
       factor = (_t - _start_time) / _dt;
     else
       factor = (_end_time - _start_time) / _dt;
   }
   else
   {
     if (_t <= _end_time)
       factor = 1.0;
     else
       factor = (_end_time - (_t - _dt)) / _dt;
   }
  
   // Negative sign to make a positive mass_flux in the input file a source
   return -_test[_i][_qp] * factor * _mass_flux;
 }