The base class for boundaries between one-dimensional spatial domains. More...

#include <Boundary1D.h>

Inheritance diagram for Boundary1D:

[legend]

Detailed Description

The base class for boundaries between one-dimensional spatial domains.

The boundary may have its own internal variables, such as surface species coverages.

The boundary types are an inlet, an outlet, a symmetry plane, and a surface.

The public methods are all virtual, and the base class implementations throw exceptions.

Definition at line 42 of file Boundary1D.h.

Public Member Functions
	Boundary1D ()
	Default constructor.

void	init () override
	Initialize.

string	domainType () const override
	Domain type flag.

bool	isConnector () override
	True if the domain is a connector domain.

virtual void	setTemperature (double t)
	Set the temperature.

virtual double	temperature ()
	Temperature [K].

virtual size_t	nSpecies ()
	Get the number of species.

virtual void	setMoleFractions (const string &xin)
	Set the mole fractions by specifying a string.

virtual void	setMoleFractions (const double *xin)
	Set the mole fractions by specifying an array.

virtual double	massFraction (size_t k)
	Mass fraction of species k.

virtual void	setMdot (double mdot)
	Set the total mass flow rate [kg/m²/s].

virtual void	setSpreadRate (double V0)
	Set tangential velocity gradient [1/s] at this boundary.

virtual double	spreadRate ()
	Tangential velocity gradient [1/s] at this boundary.

virtual double	mdot ()
	The total mass flow rate [kg/m2/s].

void	setupGrid (size_t n, const double *z) override
	called to set up initial grid, and after grid refinement

void	fromArray (SolutionArray &arr, double *soln) override
	Restore the solution for this domain from a SolutionArray.

Public Member Functions inherited from Domain1D
	Domain1D (size_t nv=1, size_t points=1, double time=0.0)
	Constructor.

	Domain1D (const Domain1D &)=delete

Domain1D &	operator= (const Domain1D &)=delete

virtual string	domainType () const
	Domain type flag.

string	type () const
	String indicating the domain implemented.

size_t	domainIndex ()
	The left-to-right location of this domain.

virtual bool	isConnector ()
	True if the domain is a connector domain.

void	setSolution (shared_ptr< Solution > sol)
	Set the solution manager.

virtual void	setKinetics (shared_ptr< Kinetics > kin)
	Set the kinetics manager.

virtual void	setTransport (shared_ptr< Transport > trans)
	Set transport model to existing instance.

const OneDim &	container () const
	The container holding this domain.

void	setContainer (OneDim *c, size_t index)
	Specify the container object for this domain, and the position of this domain in the list.

void	setBandwidth (int bw=-1)
	Set the Jacobian bandwidth. See the discussion of method bandwidth().

size_t	bandwidth ()
	Set the Jacobian bandwidth for this domain.

virtual void	init ()
	Initialize.

virtual void	setInitialState (double *xlocal=0)

virtual void	setState (size_t point, const double state, double x)

virtual void	resetBadValues (double *xg)
	When called, this function should reset "bad" values in the state vector such as negative species concentrations.

virtual void	resize (size_t nv, size_t np)
	Resize the domain to have nv components and np grid points.

Refiner &	refiner ()
	Return a reference to the grid refiner.

size_t	nComponents () const
	Number of components at each grid point.

void	checkComponentIndex (size_t n) const
	Check that the specified component index is in range.

void	checkComponentArraySize (size_t nn) const
	Check that an array size is at least nComponents().

size_t	nPoints () const
	Number of grid points in this domain.

void	checkPointIndex (size_t n) const
	Check that the specified point index is in range.

void	checkPointArraySize (size_t nn) const
	Check that an array size is at least nPoints().

virtual string	componentName (size_t n) const
	Name of component `n`. May be overloaded.

void	setComponentName (size_t n, const string &name)
	Set the name of the component `n` to `name`.

virtual size_t	componentIndex (const string &name) const
	index of component with name `name`.

void	setBounds (size_t n, double lower, double upper)
	Set the upper and lower bounds for a solution component, n.

void	setTransientTolerances (double rtol, double atol, size_t n=npos)
	Set tolerances for time-stepping mode.

void	setSteadyTolerances (double rtol, double atol, size_t n=npos)
	Set tolerances for steady-state mode.

double	rtol (size_t n)
	Relative tolerance of the nth component.

double	atol (size_t n)
	Absolute tolerance of the nth component.

double	steady_rtol (size_t n)
	Steady relative tolerance of the nth component.

double	steady_atol (size_t n)
	Steady absolute tolerance of the nth component.

double	transient_rtol (size_t n)
	Transient relative tolerance of the nth component.

double	transient_atol (size_t n)
	Transient absolute tolerance of the nth component.

double	upperBound (size_t n) const
	Upper bound on the nth component.

double	lowerBound (size_t n) const
	Lower bound on the nth component.

void	initTimeInteg (double dt, const double *x0)
	Performs the setup required before starting a time-stepping solution.

void	setSteadyMode ()
	Set the internally-stored reciprocal of the time step to 0.0, which is used to indicate that the problem is in steady-state mode.

bool	steady ()
	True if in steady-state mode.

bool	transient ()
	True if not in steady-state mode.

void	needJacUpdate ()
	Set this if something has changed in the governing equations (for example, the value of a constant has been changed, so that the last-computed Jacobian is no longer valid.

virtual void	eval (size_t j, double x, double r, integer *mask, double rdt=0.0)
	Evaluate the residual function at point j.

size_t	index (size_t n, size_t j) const
	Returns the index of the solution vector, which corresponds to component n at grid point j.

double	value (const double *x, size_t n, size_t j) const
	Returns the value of solution component n at grid point j of the solution vector x.

virtual void	setJac (MultiJac *jac)

virtual shared_ptr< SolutionArray >	asArray (const double *soln) const
	Save the state of this domain as a SolutionArray.

shared_ptr< SolutionArray >	toArray (bool normalize=false) const
	Save the state of this domain to a SolutionArray.

virtual void	fromArray (SolutionArray &arr, double *soln)
	Restore the solution for this domain from a SolutionArray.

void	fromArray (const shared_ptr< SolutionArray > &arr)
	Restore the solution for this domain from a SolutionArray.

shared_ptr< Solution >	solution () const
	Return thermo/kinetics/transport manager used in the domain.

size_t	size () const
	Return the size of the solution vector (the product of m_nv and m_points).

void	locate ()
	Find the index of the first grid point in this domain, and the start of its variables in the global solution vector.

virtual size_t	loc (size_t j=0) const
	Location of the start of the local solution vector in the global solution vector.

size_t	firstPoint () const
	The index of the first (that is, left-most) grid point belonging to this domain.

size_t	lastPoint () const
	The index of the last (that is, right-most) grid point belonging to this domain.

void	linkLeft (Domain1D *left)
	Set the left neighbor to domain 'left.

void	linkRight (Domain1D *right)
	Set the right neighbor to domain 'right.'.

void	append (Domain1D *right)
	Append domain 'right' to this one, and update all links.

Domain1D *	left () const
	Return a pointer to the left neighbor.

Domain1D *	right () const
	Return a pointer to the right neighbor.

double	prevSoln (size_t n, size_t j) const
	Value of component n at point j in the previous solution.

void	setID (const string &s)
	Specify an identifying tag for this domain.

string	id () const
	Returns the identifying tag for this domain.

virtual void	show (std::ostream &s, const double *x)
	Print the solution.

virtual void	show (const double *x)
	Print the solution.

double	z (size_t jlocal) const
	Get the coordinate [m] of the point with local index `jlocal`

double	zmin () const
	Get the coordinate [m] of the first (leftmost) grid point in this domain.

double	zmax () const
	Get the coordinate [m] of the last (rightmost) grid point in this domain.

void	setProfile (const string &name, double values, double soln)
	Set initial values for a component at each grid point.

vector< double > &	grid ()
	Access the array of grid coordinates [m].

const vector< double > &	grid () const
	Access the array of grid coordinates [m].

double	grid (size_t point) const

virtual void	setupGrid (size_t n, const double *z)
	called to set up initial grid, and after grid refinement

virtual void	_getInitialSoln (double *x)
	Writes some or all initial solution values into the global solution array, beginning at the location pointed to by x.

virtual double	initialValue (size_t n, size_t j)
	Initial value of solution component n at grid point j.

virtual void	_finalize (const double *x)
	In some cases, a domain may need to set parameters that depend on the initial solution estimate.

void	forceFullUpdate (bool update)
	In some cases, for computational efficiency some properties (such as transport coefficients) may not be updated during Jacobian evaluations.

void	setData (shared_ptr< vector< double > > &data)
	Set shared data pointer.

Protected Member Functions
void	_init (size_t n)
	Initialize member variables based on the adjacent domains.

Protected Member Functions inherited from Domain1D
virtual AnyMap	getMeta () const
	Retrieve meta data.

virtual void	setMeta (const AnyMap &meta)
	Retrieve meta data.

Protected Attributes
Flow1D *	m_flow_left = nullptr
	Flow domain to the left of this boundary.

Flow1D *	m_flow_right = nullptr

size_t	m_left_nv = 0
	Flow domain to the right of this boundary.

size_t	m_right_nv = 0
	Number of state vector components in right flow domain.

size_t	m_left_nsp = 0
	Number of species in left flow domain.

size_t	m_right_nsp = 0
	Number of species in right flow domain.

ThermoPhase *	m_phase_left = nullptr
	Thermo object used by left flow domain.

ThermoPhase *	m_phase_right = nullptr
	Thermo object used by right flow domain.

double	m_temp = 0.0
	Temperature of the boundary.

double	m_mdot = 0.0
	Mass flow rate at the boundary.

Protected Attributes inherited from Domain1D
shared_ptr< vector< double > >	m_state
	data pointer shared from OneDim

double	m_rdt = 0.0
	Reciprocal of the time step.

size_t	m_nv = 0
	Number of solution components.

size_t	m_points
	Number of grid points.

vector< double >	m_slast
	Solution vector at the last time step.

vector< double >	m_max
	Upper bounds on solution components.

vector< double >	m_min
	Lower bounds on solution components.

vector< double >	m_rtol_ss
	Relative tolerances for steady mode.

vector< double >	m_rtol_ts
	Relative tolerances for transient mode.

vector< double >	m_atol_ss
	Absolute tolerances for steady mode.

vector< double >	m_atol_ts
	Absolute tolerances for transient mode.

vector< double >	m_z
	1D spatial grid coordinates

OneDim *	m_container = nullptr
	Parent OneDim simulation containing this and adjacent domains.

size_t	m_index
	Left-to-right location of this domain.

size_t	m_iloc = 0
	Starting location within the solution vector for unknowns that correspond to this domain.

size_t	m_jstart = 0
	Index of the first point in this domain in the global point list.

Domain1D *	m_left = nullptr
	Pointer to the domain to the left.

Domain1D *	m_right = nullptr
	Pointer to the domain to the right.

string	m_id
	Identity tag for the domain.

unique_ptr< Refiner >	m_refiner
	Refiner object used for placing grid points.

vector< string >	m_name
	Names of solution components.

int	m_bw = -1
	See bandwidth()

bool	m_force_full_update = false
	see forceFullUpdate()

shared_ptr< Solution >	m_solution
	Composite thermo/kinetics/transport handler.

Constructor & Destructor Documentation

◆ Boundary1D()

Boundary1D ( )

Default constructor.

Definition at line 16 of file Boundary1D.cpp.

Member Function Documentation

◆ init()

void init ( )

inlineoverridevirtual

Initialize.

This method is called by OneDim::init() for each domain once at the beginning of a simulation. Base class method does nothing, but may be overloaded.

Reimplemented from Domain1D.

Reimplemented in Inlet1D, Empty1D, Symm1D, Outlet1D, OutletRes1D, Surf1D, and ReactingSurf1D.

Definition at line 48 of file Boundary1D.h.

◆ domainType()

string domainType ( ) const

inlineoverridevirtual

Domain type flag.

Since: Starting in Cantera 3.1, the return type is a string.

Reimplemented from Domain1D.

Reimplemented in Inlet1D, Empty1D, Symm1D, Outlet1D, OutletRes1D, Surf1D, and ReactingSurf1D.

Definition at line 52 of file Boundary1D.h.

◆ isConnector()

bool isConnector ( )

inlineoverridevirtual

True if the domain is a connector domain.

Reimplemented from Domain1D.

Definition at line 56 of file Boundary1D.h.

◆ setTemperature()

virtual void setTemperature ( double t )

inlinevirtual

Set the temperature.

Reimplemented in Inlet1D.

Definition at line 61 of file Boundary1D.h.

◆ temperature()

virtual double temperature ( )

inlinevirtual

Temperature [K].

Definition at line 66 of file Boundary1D.h.

◆ nSpecies()

virtual size_t nSpecies ( )

inlinevirtual

Get the number of species.

Reimplemented in Inlet1D, and OutletRes1D.

Definition at line 71 of file Boundary1D.h.

◆ setMoleFractions() [1/2]

virtual void setMoleFractions ( const string & xin )

inlinevirtual

Set the mole fractions by specifying a string.

Reimplemented in Inlet1D, and OutletRes1D.

Definition at line 76 of file Boundary1D.h.

◆ setMoleFractions() [2/2]

virtual void setMoleFractions ( const double * xin )

inlinevirtual

Set the mole fractions by specifying an array.

Reimplemented in Inlet1D, and OutletRes1D.

Definition at line 81 of file Boundary1D.h.

◆ massFraction()

virtual double massFraction ( size_t k )

inlinevirtual

Mass fraction of species k.

Reimplemented in Inlet1D, and OutletRes1D.

Definition at line 86 of file Boundary1D.h.

◆ setMdot()

virtual void setMdot ( double mdot )

inlinevirtual

Set the total mass flow rate [kg/m²/s].

Definition at line 91 of file Boundary1D.h.

◆ setSpreadRate()

virtual void setSpreadRate ( double V0 )

inlinevirtual

Set tangential velocity gradient [1/s] at this boundary.

Reimplemented in Inlet1D.

Definition at line 96 of file Boundary1D.h.

◆ spreadRate()

virtual double spreadRate ( )

inlinevirtual

Tangential velocity gradient [1/s] at this boundary.

Reimplemented in Inlet1D.

Definition at line 101 of file Boundary1D.h.

◆ mdot()

virtual double mdot ( )

inlinevirtual

The total mass flow rate [kg/m2/s].

Definition at line 106 of file Boundary1D.h.

◆ setupGrid()

void setupGrid	(	size_t	n,
		const double *	z
	)

inlineoverridevirtual

called to set up initial grid, and after grid refinement

Reimplemented from Domain1D.

Definition at line 110 of file Boundary1D.h.

◆ fromArray()

void fromArray	(	SolutionArray &	arr,
		double *	soln
	)

overridevirtual

Restore the solution for this domain from a SolutionArray.

Parameters

[in]	arr	SolutionArray defining the state of this domain
[out]	soln	Value of the solution vector, local to this domain

Since: New in Cantera 3.0.

Reimplemented from Domain1D.

Reimplemented in Inlet1D, OutletRes1D, Surf1D, and ReactingSurf1D.

Definition at line 76 of file Boundary1D.cpp.

◆ _init()

void _init ( size_t n )

protected

Initialize member variables based on the adjacent domains.

Parameters

n	Number of state variables associated with the boundary object

Definition at line 20 of file Boundary1D.cpp.

Member Data Documentation

◆ m_flow_left

Flow1D* m_flow_left = nullptr

protected

Flow domain to the left of this boundary.

Definition at line 119 of file Boundary1D.h.

◆ m_flow_right

Flow1D* m_flow_right = nullptr

protected

Definition at line 120 of file Boundary1D.h.

◆ m_left_nv

size_t m_left_nv = 0

protected

Flow domain to the right of this boundary.

Number of state vector components in left flow domain

Definition at line 121 of file Boundary1D.h.

◆ m_right_nv

size_t m_right_nv = 0

protected

Number of state vector components in right flow domain.

Definition at line 122 of file Boundary1D.h.

◆ m_left_nsp

size_t m_left_nsp = 0

protected

Number of species in left flow domain.

Definition at line 123 of file Boundary1D.h.

◆ m_right_nsp

size_t m_right_nsp = 0

protected

Number of species in right flow domain.

Definition at line 124 of file Boundary1D.h.

◆ m_phase_left

ThermoPhase* m_phase_left = nullptr

protected

Thermo object used by left flow domain.

Definition at line 125 of file Boundary1D.h.

◆ m_phase_right

ThermoPhase* m_phase_right = nullptr

protected

Thermo object used by right flow domain.

Definition at line 126 of file Boundary1D.h.

◆ m_temp

double m_temp = 0.0

protected

Temperature of the boundary.

Definition at line 129 of file Boundary1D.h.

◆ m_mdot

double m_mdot = 0.0

protected

Mass flow rate at the boundary.

Definition at line 131 of file Boundary1D.h.

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

Detailed Description

Public Member Functions

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ Boundary1D()

Member Function Documentation

◆ init()

◆ domainType()

◆ isConnector()

◆ setTemperature()

◆ temperature()

◆ nSpecies()

◆ setMoleFractions() [1/2]

◆ setMoleFractions() [2/2]

◆ massFraction()

◆ setMdot()

◆ setSpreadRate()

◆ spreadRate()

◆ mdot()

◆ setupGrid()

◆ fromArray()

◆ _init()

Member Data Documentation

◆ m_flow_left

◆ m_flow_right

◆ m_left_nv

◆ m_right_nv

◆ m_left_nsp

◆ m_right_nsp

◆ m_phase_left

◆ m_phase_right

◆ m_temp

◆ m_mdot