Base class for a tracer particle solver. This solver can advance particles in a pre-defined velocity field. More...

#include <CD_TracerParticleSolver.H>

Collaboration diagram for TracerParticleSolver:

[legend]

Public Member Functions
	TracerParticleSolver ()
	Default constructor.

	TracerParticleSolver (const TracerParticleSolver &a_other)=delete
	Disallowed constructor.

	TracerParticleSolver (const TracerParticleSolver &&a_other)=delete
	Disallowed move constructor.

	TracerParticleSolver (const RefCountedPtr< AmrMesh > &a_amr, const RefCountedPtr< ComputationalGeometry > &a_compGeom)
	Full constructor.

TracerParticleSolver &	operator= (const TracerParticleSolver &a_other)=delete
	Disallowed assignment operator.

TracerParticleSolver &	operator= (const TracerParticleSolver &&a_other)=delete
	Disallowed move assignment operator.

virtual	~TracerParticleSolver ()
	Destructor.

virtual void	allocate ()
	Allocate storage for this solver.

virtual void	registerOperators () const
	Register operators needed for AMR core functionality.

virtual void	setAmr (const RefCountedPtr< AmrMesh > &a_amrMesh)
	Set AmrMesh.

virtual void	setComputationalGeometry (const RefCountedPtr< ComputationalGeometry > &a_compGeom)
	Set the computational geometry.

virtual void	setName (const std::string &a_name) noexcept
	Set the solver name.

virtual void	setVolumeScale (const bool a_scale) noexcept
	Turn on/off volume scaling.

virtual void	setRealm (const std::string &a_realm)
	Set the solver realm.

virtual void	setPhase (const phase::which_phase &a_phase)
	Set the solver phase.

virtual void	setTime (const int a_step, const Real a_time, const Real a_dt)
	Set the time for this solver.

virtual void	setVelocity (const EBAMRCellData &a_velocityField)
	Set the tracer particle velocity field.

virtual void	preRegrid (const int a_lbase, const int a_oldFinestLevel)
	Perform pre-regrid operations.

virtual void	regrid (const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel)
	Regrid this solver.

virtual void	remap ()
	Remap particles.

virtual void	setDeposition (const DepositionType a_deposition) noexcept
	Set deposition method.

virtual void	deposit (EBAMRCellData &a_phi) const noexcept
	Deposit particle weight on mesh.

virtual void	interpolateWeight (const EBAMRCellData &a_scalar) noexcept
	Interpolate a scalar field onto the particle weight.

virtual void	interpolateVelocities ()
	Interpolate particles velocities.

virtual void	parseOptions ()
	Parse solver options.

virtual void	parseRuntimeOptions ()
	Parse solver run-time options.

virtual void	writePlotFile ()
	Write plot file.

virtual int	getNumberOfPlotVariables () const
	Get the number of plot variables.

virtual Vector< std::string >	getPlotVariableNames () const
	Get plot variable names.

virtual void	writePlotData (LevelData< EBCellFAB > &a_output, int &a_comp, const std::string &a_outputRealm, const int a_level) const noexcept
	Write plot data.

virtual Real	computeDt () const
	Compute dt = dx/max(v_x, v_y, v_z) minimized over all particles.

virtual ParticleContainer< P > &	getParticles ()
	Get all particles.

virtual const ParticleContainer< P > &	getParticles () const
	Get all particles. Const version.

const EBAMRCellData &	getVelocityField () const
	Return the velocity field.

virtual DepositionType	getDepositionType () const
	Get the deposition type.

virtual CoarseFineDeposition	getCoarseFineDepositionType () const
	Get the coarse-fine deposition type.

virtual DepositionType	getInterpolationType () const
	Get the interpolation type.

Protected Member Functions
void	parseDeposition ()
	Parse deposition method.

void	parsePlotVariables ()
	Parse plot variables.

void	parseVerbosity ()
	Parse solver verbosity.

virtual void	writeData (LevelData< EBCellFAB > &a_output, int &a_comp, const EBAMRCellData &a_data, const std::string a_outputRealm, const int a_level, const bool a_interpToCentroids, const bool a_interpGhost) const noexcept
	Write data to output. Convenience function.

Protected Attributes
RefCountedPtr< AmrMesh >	m_amr
	Handle to AMR mesh.

RefCountedPtr< ComputationalGeometry >	m_computationalGeometry
	Handle to computational geometry.

DepositionType	m_deposition
	Deposition type when depositing particles to grid.

DepositionType	m_interpolation
	Deposition type when interpolating to particles.

CoarseFineDeposition	m_coarseFineDeposition
	Coarse-fine deposition strategy.

std::string	m_realm
	Realm where this solver lives.

std::string	m_name
	Solver name.

std::string	m_className
	Class name.

phase::which_phase	m_phase
	Phase where this solver lives.

Real	m_dt
	Time step.

Real	m_time
	Time.

int	m_timeStep
	Time step.

int	m_verbosity
	Verbosity level.

bool	m_plotWeight
	Hook for adding weight to plot file.

bool	m_plotVelocity
	Hook for adding velocity field to plot file.

bool	m_volumeScale
	Do volumetric scaling when depositing (or not).

EBAMRCellData	m_velocityField
	Mesh-based velocity field. Must be set through setVelocity.

ParticleContainer< P >	m_particles
	Computational particles.

Detailed Description

template<typename P>
class TracerParticleSolver

Base class for a tracer particle solver. This solver can advance particles in a pre-defined velocity field.

The user can set the velocity field through public member functions. This class is templated so the user can switch tracer particle implementations.

This is a single-phase solver – i.e. the particles only live on one of the phases. Extensions to multiphase is certainly possible.

The particles are stored Struct-of-Arrays (ParticleContainer). P is the user PAYLOAD type; the position and weight are container-owned columns (not payload members). The template requirements on the payload type P are: 1) It must expose the velocity as SpaceDim scalar payload columns named D_DECL(v_x, v_y, v_z), so the solver can interpolate the mesh velocity field into them and compute the advective time step. 2) It must have a ParticleTraits

specialization listing its columns (including v_x/v_y/v_z).

Constructor & Destructor Documentation

◆ TracerParticleSolver() [1/3]

template<typename P >

TracerParticleSolver::TracerParticleSolver ( const TracerParticleSolver & a_other )

delete

Disallowed constructor.

Parameters

[in] a_other Other solver.

◆ TracerParticleSolver() [2/3]

template<typename P >

TracerParticleSolver::TracerParticleSolver ( const TracerParticleSolver && a_other )

delete

Disallowed move constructor.

Parameters

[in] a_other Other solver.

◆ TracerParticleSolver() [3/3]

template<typename P >

TracerParticleSolver< P >::TracerParticleSolver	(	const RefCountedPtr< AmrMesh > &	a_amr,
		const RefCountedPtr< ComputationalGeometry > &	a_compGeom
	)

inline

Full constructor.

Parameters

[in]	a_amr	Handle to AmrMesh.
[in]	a_compGeom	Computational geometry.

Member Function Documentation

◆ computeDt()

template<typename P >

Real TracerParticleSolver::computeDt ( ) const

virtual

Compute dt = dx/max(v_x, v_y, v_z) minimized over all particles.

Returns: Stable advection time step.

◆ deposit()

template<typename P >

void TracerParticleSolver::deposit ( EBAMRCellData & a_phi ) const

inlinevirtualnoexcept

Deposit particle weight on mesh.

Parameters

[out] a_phi Deposited weight.

◆ getCoarseFineDepositionType()

template<typename P >

CoarseFineDeposition TracerParticleSolver::getCoarseFineDepositionType ( ) const

inlinevirtual

Get the coarse-fine deposition type.

Returns: Active coarse-fine deposition strategy.

◆ getDepositionType()

template<typename P >

DepositionType TracerParticleSolver::getDepositionType ( ) const

inlinevirtual

Get the deposition type.

Returns: Active deposition type.

◆ getInterpolationType()

template<typename P >

DepositionType TracerParticleSolver::getInterpolationType ( ) const

inlinevirtual

Get the interpolation type.

Returns: Active interpolation type.

◆ getNumberOfPlotVariables()

template<typename P >

int TracerParticleSolver::getNumberOfPlotVariables ( ) const

inlinevirtual

Get the number of plot variables.

Returns: Returns number of plot variables.

◆ getParticles() [1/2]

template<typename P >

ParticleContainer & TracerParticleSolver::getParticles ( )

inlinevirtual

Get all particles.

Returns: m_particles

◆ getParticles() [2/2]

template<typename P >

const ParticleContainer & TracerParticleSolver::getParticles ( ) const

inlinevirtual

Get all particles. Const version.

Returns: m_particles

◆ getPlotVariableNames()

template<typename P >

Vector< std::string > TracerParticleSolver::getPlotVariableNames ( ) const

inlinevirtual

Get plot variable names.

Returns: List of plot variable names.

◆ getVelocityField()

template<typename P >

const EBAMRCellData & TracerParticleSolver::getVelocityField ( ) const

inline

Return the velocity field.

Returns: m_velocityField

◆ interpolateWeight()

template<typename P >

void TracerParticleSolver::interpolateWeight ( const EBAMRCellData & a_scalar )

virtualnoexcept

Interpolate a scalar field onto the particle weight.

Parameters

[in] a_scalar Scalar field to interpolate from.

◆ operator=() [1/2]

template<typename P >

TracerParticleSolver & TracerParticleSolver::operator= ( const TracerParticleSolver && a_other )

delete

Disallowed move assignment operator.

Parameters

[in] a_other Other solver.

◆ operator=() [2/2]

template<typename P >

TracerParticleSolver & TracerParticleSolver::operator= ( const TracerParticleSolver & a_other )

delete

Disallowed assignment operator.

Parameters

[in] a_other Other solver.

◆ preRegrid()

template<typename P >

void TracerParticleSolver< P >::preRegrid	(	const int	a_lbase,
		const int	a_oldFinestLevel
	)

inlinevirtual

Perform pre-regrid operations.

Parameters

[in]	a_lbase	Coarsest level that changed during regrid.
[in]	a_oldFinestLevel	Finest grid level before the regrid operation.

Note: This stores m_phi and m_source.

◆ regrid()

template<typename P >

void TracerParticleSolver< P >::regrid	(	const int	a_lmin,
		const int	a_oldFinestLevel,
		const int	a_newFinestLevel
	)

inlinevirtual

Regrid this solver.

Parameters

[in]	a_lmin	Coarsest level where grids did not change.
[in]	a_oldFinestLevel	Finest AMR level before the regrid.
[in]	a_newFinestLevel	Finest AMR level after the regrid.

Note: The velocity is NOT remeshed – the user is responsible for setting this after the regrid.

◆ setAmr()

template<typename P >

void TracerParticleSolver::setAmr ( const RefCountedPtr< AmrMesh > & a_amrMesh )

inlinevirtual

Set AmrMesh.

Parameters

[in] a_amrMesh AmrMesh

◆ setComputationalGeometry()

template<typename P >

void TracerParticleSolver::setComputationalGeometry ( const RefCountedPtr< ComputationalGeometry > & a_compGeom )

inlinevirtual

Set the computational geometry.

Parameters

[in] a_compGeom Computational geometry

◆ setDeposition()

template<typename P >

void TracerParticleSolver::setDeposition ( const DepositionType a_deposition )

inlinevirtualnoexcept

Set deposition method.

Parameters

[in] a_deposition Deposition method

◆ setName()

template<typename P >

void TracerParticleSolver::setName ( const std::string & a_name )

inlinevirtualnoexcept

Set the solver name.

Parameters

[in] a_name Solver name

◆ setPhase()

template<typename P >

void TracerParticleSolver::setPhase ( const phase::which_phase & a_phase )

inlinevirtual

Set the solver phase.

Parameters

[in] a_phase Phase

◆ setRealm()

template<typename P >

void TracerParticleSolver::setRealm ( const std::string & a_realm )

inlinevirtual

Set the solver realm.

Parameters

[in] a_realm Realm name.

◆ setTime()

template<typename P >

void TracerParticleSolver< P >::setTime	(	const int	a_step,
		const Real	a_time,
		const Real	a_dt
	)

virtual

Set the time for this solver.

Parameters

[in]	a_step	Time step number
[in]	a_time	Time (in seconds)
[in]	a_dt	Time step increment

Note: This sets m_step=a_step, m_time=a_time, m_dt=a_dt

◆ setVelocity()

template<typename P >

void TracerParticleSolver::setVelocity ( const EBAMRCellData & a_velocityField )

inlinevirtual

Set the tracer particle velocity field.

Parameters

[in] a_velocityField Velocity field.

◆ setVolumeScale()

template<typename P >

void TracerParticleSolver::setVolumeScale ( const bool a_scale )

inlinevirtualnoexcept

Turn on/off volume scaling.

Parameters

[in] a_scale Volume scale or not.

◆ writeData()

template<typename P >

void TracerParticleSolver< P >::writeData	(	LevelData< EBCellFAB > &	a_output,
		int &	a_comp,
		const EBAMRCellData &	a_data,
		const std::string	a_outputRealm,
		const int	a_level,
		const bool	a_interpToCentroids,
		const bool	a_interpGhost
	)		const

protectedvirtualnoexcept

Write data to output. Convenience function.

Parameters

[in,out]	a_output	Output data holder.
[in,out]	a_comp	Starting component where this solver begins writing the output.
[in]	a_data	Data to write.
[in]	a_outputRealm	Realm where a_output belongs.
[in]	a_level	Grid level.
[in]	a_interpToCentroids	If true, a_data will be interpolated to cell centroids before writing to a_output.
[in]	a_interpGhost	If true, interpolate ghost cells.

◆ writePlotData()

template<typename P >

void TracerParticleSolver< P >::writePlotData	(	LevelData< EBCellFAB > &	a_output,
		int &	a_comp,
		const std::string &	a_outputRealm,
		const int	a_level
	)		const

inlinevirtualnoexcept

Write plot data.

Parameters

[in,out]	a_output	Output data holder.
[in,out]	a_comp	Starting component where this solver begins writing the output.
[in]	a_outputRealm	Realm where a_output belongs.
[in]	a_level	Grid level.

Note: This will write the plot data in this solver to a_output, starting on a_comp.

◆ writePlotFile()

template<typename P >

void TracerParticleSolver::writePlotFile ( )

inlinevirtual

Write plot file.

The name of the plot file is m_name.stepXXXXX.DIM.hdf5

Note: This calls writePlotData(...)

The documentation for this class was generated from the following files:

Source/TracerParticles/CD_TracerParticleSolver.H
Source/TracerParticles/CD_TracerParticleSolverImplem.H

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TracerParticleSolver() [1/3]

◆ TracerParticleSolver() [2/3]

◆ TracerParticleSolver() [3/3]

Member Function Documentation

◆ computeDt()

◆ deposit()

◆ getCoarseFineDepositionType()

◆ getDepositionType()

◆ getInterpolationType()

◆ getNumberOfPlotVariables()

◆ getParticles() [1/2]

◆ getParticles() [2/2]

◆ getPlotVariableNames()

◆ getVelocityField()

◆ interpolateWeight()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ preRegrid()

◆ regrid()

◆ setAmr()

◆ setComputationalGeometry()

◆ setDeposition()

◆ setName()

◆ setPhase()

◆ setRealm()

◆ setTime()

◆ setVelocity()

◆ setVolumeScale()

◆ writeData()

◆ writePlotData()

◆ writePlotFile()