www.mooseframework.org
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
CNSFVNoSlopeLimiting Class Reference

A user object that does no slope limiting in multi-dimensions. More...

#include <CNSFVNoSlopeLimiting.h>

Inheritance diagram for CNSFVNoSlopeLimiting:
[legend]

Public Member Functions

 CNSFVNoSlopeLimiting (const InputParameters &parameters)
 
virtual std::vector< RealGradient > limitElementSlope () const
 compute the limited slope of the cell More...
 
virtual void initialize ()
 
virtual void finalize ()
 
virtual void computeElement ()
 
virtual const std::vector< RealGradient > & getElementSlope (dof_id_type elementid) const
 accessor function call More...
 
virtual void execute ()
 
virtual void pre ()
 
virtual void onElement (const Elem *elem)
 
virtual void onBoundary (const Elem *elem, unsigned int side, BoundaryID bnd_id)
 
virtual void onInternalSide (const Elem *elem, unsigned int side)
 
virtual void onInterface (const Elem *elem, unsigned int side, BoundaryID bnd_id)
 
virtual void post ()
 
virtual void subdomainChanged ()
 
virtual bool keepGoing ()
 
virtual void meshChanged ()
 
void join (const ElementLoopUserObject &)
 

Protected Member Functions

virtual void serialize (std::string &serialized_buffer)
 
virtual void deserialize (std::vector< std::string > &serialized_buffers)
 
virtual void caughtMooseException (MooseException &e)
 
virtual void computeBoundary ()
 
virtual void computeInternalSide ()
 
virtual void computeInterface ()
 

Protected Attributes

std::map< dof_id_type, std::vector< RealGradient > > _lslope
 store the updated slopes into this map indexed by element ID More...
 
bool _include_bc
 option whether to include BCs More...
 
const SlopeReconstructionBase_rslope
 slope reconstruction user object More...
 
const MooseArray< Point > & _q_point_face
 required data for face assembly More...
 
QBase *& _qrule_face
 
const MooseArray< Real > & _JxW_face
 
const MooseArray< Point > & _normals_face
 
unsigned int & _side
 current side of the current element More...
 
const Elem *& _side_elem
 
const Real & _side_volume
 
const Elem *& _neighbor_elem
 the neighboring element More...
 
MooseMesh & _mesh
 
const Elem * _current_elem
 
const Real & _current_elem_volume
 
unsigned int _current_side
 
const Elem * _current_neighbor
 
const MooseArray< Point > & _q_point
 
QBase *& _qrule
 
const MooseArray< Real > & _JxW
 
const MooseArray< Real > & _coord
 
bool _have_interface_elems
 true if we have cached interface elements, false if they need to be cached. We want to (re)cache only when mesh changed More...
 
std::set< dof_id_type > _interface_elem_ids
 List of element IDs that are on the processor boundary and need to be send to other processors. More...
 
SubdomainID _subdomain
 The subdomain for the current element. More...
 
SubdomainID _old_subdomain
 The subdomain for the last element. More...
 

Detailed Description

A user object that does no slope limiting in multi-dimensions.

Definition at line 22 of file CNSFVNoSlopeLimiting.h.

Constructor & Destructor Documentation

CNSFVNoSlopeLimiting::CNSFVNoSlopeLimiting ( const InputParameters &  parameters)

Definition at line 19 of file CNSFVNoSlopeLimiting.C.

20  : SlopeLimitingBase(parameters)
21 {
22 }
SlopeLimitingBase(const InputParameters &parameters)

Member Function Documentation

void ElementLoopUserObject::caughtMooseException ( MooseException &  e)
protectedvirtualinherited

Definition at line 253 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

254 {
255  std::string what(e.what());
256  _fe_problem.setException(what);
257 }
void ElementLoopUserObject::computeBoundary ( )
protectedvirtualinherited

Definition at line 231 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::onBoundary().

232 {
233 }
void SlopeLimitingBase::computeElement ( )
virtualinherited

Reimplemented from ElementLoopUserObject.

Definition at line 68 of file SlopeLimitingBase.C.

69 {
70  dof_id_type _elementID = _current_elem->id();
71 
72  _lslope[_elementID] = limitElementSlope();
73 }
virtual std::vector< RealGradient > limitElementSlope() const =0
compute the slope of the cell
std::map< dof_id_type, std::vector< RealGradient > > _lslope
store the updated slopes into this map indexed by element ID
void ElementLoopUserObject::computeInterface ( )
protectedvirtualinherited

Definition at line 241 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::onInterface().

242 {
243 }
void ElementLoopUserObject::computeInternalSide ( )
protectedvirtualinherited

Definition at line 236 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::onInternalSide().

237 {
238 }
void SlopeLimitingBase::deserialize ( std::vector< std::string > &  serialized_buffers)
protectedvirtualinherited

Definition at line 95 of file SlopeLimitingBase.C.

Referenced by SlopeLimitingBase::finalize().

96 {
97  // The input string stream used for deserialization
98  std::istringstream iss;
99 
100  mooseAssert(serialized_buffers.size() == _app.n_processors(),
101  "Unexpected size of serialized_buffers: " << serialized_buffers.size());
102 
103  for (auto rank = decltype(_app.n_processors())(0); rank < serialized_buffers.size(); ++rank)
104  {
105  if (rank == processor_id())
106  continue;
107 
108  iss.str(serialized_buffers[rank]); // populate the stream with a new buffer
109  iss.clear(); // reset the string stream state
110 
111  // Load the communicated data into all of the other processors' slots
112 
113  unsigned int size = 0;
114  iss.read((char *)&size, sizeof(size));
115 
116  for (unsigned int i = 0; i < size; i++)
117  {
118  dof_id_type key;
119  loadHelper(iss, key, this);
120 
121  std::vector<RealGradient> value;
122  loadHelper(iss, value, this);
123 
124  // merge the data we received from other procs
125  _lslope.insert(std::pair<dof_id_type, std::vector<RealGradient>>(key, value));
126  }
127  }
128 }
std::map< dof_id_type, std::vector< RealGradient > > _lslope
store the updated slopes into this map indexed by element ID
void ElementLoopUserObject::execute ( )
virtualinherited

Definition at line 69 of file ElementLoopUserObject.C.

70 {
71  ConstElemRange & elem_range = *_mesh.getActiveLocalElementRange();
72 
73  try
74  {
75  pre();
76 
77  _subdomain = std::numeric_limits<SubdomainID>::max();
78  ConstElemRange::const_iterator el = elem_range.begin();
79  for (el = elem_range.begin(); el != elem_range.end(); ++el)
80  {
81  if (!keepGoing())
82  break;
83 
84  const Elem * elem = *el;
85  unsigned int cur_subdomain = elem->subdomain_id();
86 
88  _subdomain = cur_subdomain;
89 
90  if (this->hasBlocks(_subdomain))
91  {
94 
95  onElement(elem);
96 
97  for (unsigned int side = 0; side < elem->n_sides(); side++)
98  {
99  std::vector<BoundaryID> boundary_ids = _mesh.getBoundaryIDs(elem, side);
100 
101  if (boundary_ids.size() > 0)
102  for (std::vector<BoundaryID>::iterator it = boundary_ids.begin();
103  it != boundary_ids.end();
104  ++it)
105  onBoundary(elem, side, *it);
106 
107  if (elem->neighbor(side) != NULL)
108  {
109  if (this->hasBlocks(elem->neighbor(side)->subdomain_id()))
110  onInternalSide(elem, side);
111  if (boundary_ids.size() > 0)
112  for (std::vector<BoundaryID>::iterator it = boundary_ids.begin();
113  it != boundary_ids.end();
114  ++it)
115  onInterface(elem, side, *it);
116  }
117  } // sides
118  }
119  } // range
120 
121  post();
122  }
123  catch (MooseException & e)
124  {
126  }
127 }
virtual void onInterface(const Elem *elem, unsigned int side, BoundaryID bnd_id)
virtual void onInternalSide(const Elem *elem, unsigned int side)
virtual void caughtMooseException(MooseException &e)
virtual void onBoundary(const Elem *elem, unsigned int side, BoundaryID bnd_id)
SubdomainID _subdomain
The subdomain for the current element.
virtual void onElement(const Elem *elem)
SubdomainID _old_subdomain
The subdomain for the last element.
void SlopeLimitingBase::finalize ( )
virtualinherited

Reimplemented from ElementLoopUserObject.

Definition at line 131 of file SlopeLimitingBase.C.

132 {
134 
135  if (_app.n_processors() > 1)
136  {
137  std::vector<std::string> send_buffers(1);
138  std::vector<std::string> recv_buffers;
139 
140  recv_buffers.reserve(_app.n_processors());
141  serialize(send_buffers[0]);
142  comm().allgather_packed_range((void *)(nullptr),
143  send_buffers.begin(),
144  send_buffers.end(),
145  std::back_inserter(recv_buffers));
146  deserialize(recv_buffers);
147  }
148 }
virtual void serialize(std::string &serialized_buffer)
virtual void deserialize(std::vector< std::string > &serialized_buffers)
const std::vector< RealGradient > & SlopeLimitingBase::getElementSlope ( dof_id_type  elementid) const
virtualinherited

accessor function call

Definition at line 56 of file SlopeLimitingBase.C.

Referenced by CNSFVMaterial::computeQpProperties(), and AEFVMaterial::computeQpProperties().

57 {
58  Threads::spin_mutex::scoped_lock lock(_mutex);
59  std::map<dof_id_type, std::vector<RealGradient>>::const_iterator pos = _lslope.find(elementid);
60 
61  if (pos == _lslope.end())
62  mooseError("Limited slope is not cached for element id '", elementid, "' in ", __FUNCTION__);
63 
64  return pos->second;
65 }
static Threads::spin_mutex _mutex
std::map< dof_id_type, std::vector< RealGradient > > _lslope
store the updated slopes into this map indexed by element ID
void SlopeLimitingBase::initialize ( )
virtualinherited

Reimplemented from ElementLoopUserObject.

Definition at line 48 of file SlopeLimitingBase.C.

49 {
51 
52  _lslope.clear();
53 }
std::map< dof_id_type, std::vector< RealGradient > > _lslope
store the updated slopes into this map indexed by element ID
void ElementLoopUserObject::join ( const ElementLoopUserObject )
inherited

Definition at line 221 of file ElementLoopUserObject.C.

222 {
223 }
virtual bool ElementLoopUserObject::keepGoing ( )
inlinevirtualinherited

Definition at line 81 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::execute().

81 { return true; }
std::vector< RealGradient > CNSFVNoSlopeLimiting::limitElementSlope ( ) const
virtual

compute the limited slope of the cell

Implements SlopeLimitingBase.

Definition at line 25 of file CNSFVNoSlopeLimiting.C.

26 {
28 }
const SlopeReconstructionBase & _rslope
slope reconstruction user object
virtual const std::vector< RealGradient > & getElementSlope(dof_id_type elementid) const
accessor function call to get element slope values
void ElementLoopUserObject::meshChanged ( )
virtualinherited

Reimplemented in SlopeReconstructionBase.

Definition at line 246 of file ElementLoopUserObject.C.

Referenced by SlopeReconstructionBase::meshChanged().

247 {
248  _interface_elem_ids.clear();
249  _have_interface_elems = false;
250 }
std::set< dof_id_type > _interface_elem_ids
List of element IDs that are on the processor boundary and need to be send to other processors...
bool _have_interface_elems
true if we have cached interface elements, false if they need to be cached. We want to (re)cache only...
void ElementLoopUserObject::onBoundary ( const Elem *  elem,
unsigned int  side,
BoundaryID  bnd_id 
)
virtualinherited

Definition at line 153 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

154 {
155  _current_side = side;
156  computeBoundary();
157 }
void ElementLoopUserObject::onElement ( const Elem *  elem)
virtualinherited

Definition at line 146 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

147 {
148  _current_elem = elem;
149  computeElement();
150 }
void ElementLoopUserObject::onInterface ( const Elem *  elem,
unsigned int  side,
BoundaryID  bnd_id 
)
virtualinherited

Definition at line 188 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

189 {
190  _current_elem = elem;
191  // Pointer to the neighbor we are currently working on.
192  _current_neighbor = elem->neighbor(side);
193 
194  // Get the global id of the element and the neighbor
195  const dof_id_type elem_id = elem->id();
196  const dof_id_type neighbor_id = _current_neighbor->id();
197 
198  // TODO: add if-statement to check if this needs to be executed
199  if ((_current_neighbor->active() && (_current_neighbor->level() == elem->level()) &&
200  (elem_id < neighbor_id)) ||
201  (_current_neighbor->level() < elem->level()))
202  {
204  }
205 
206  if (!_have_interface_elems &&
207  (_current_elem->processor_id() != _current_neighbor->processor_id()))
208  {
209  // if my current neighbor is on another processor store the current element
210  // ID for later communication
211  _interface_elem_ids.insert(_current_elem->id());
212  }
213 }
std::set< dof_id_type > _interface_elem_ids
List of element IDs that are on the processor boundary and need to be send to other processors...
bool _have_interface_elems
true if we have cached interface elements, false if they need to be cached. We want to (re)cache only...
void ElementLoopUserObject::onInternalSide ( const Elem *  elem,
unsigned int  side 
)
virtualinherited

Definition at line 160 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

161 {
162  _current_elem = elem;
163  // Pointer to the neighbor we are currently working on.
164  _current_neighbor = elem->neighbor(side);
165 
166  // Get the global id of the element and the neighbor
167  const dof_id_type elem_id = elem->id();
168  const dof_id_type neighbor_id = _current_neighbor->id();
169 
170  // TODO: add if-statement to check if this needs to be executed
171  if ((_current_neighbor->active() && (_current_neighbor->level() == elem->level()) &&
172  (elem_id < neighbor_id)) ||
173  (_current_neighbor->level() < elem->level()))
174  {
176  }
177 
178  if (!_have_interface_elems &&
179  (_current_elem->processor_id() != _current_neighbor->processor_id()))
180  {
181  // if my current neighbor is on another processor store the current element ID for later
182  // communication
183  _interface_elem_ids.insert(_current_elem->id());
184  }
185 }
virtual void computeInternalSide()
std::set< dof_id_type > _interface_elem_ids
List of element IDs that are on the processor boundary and need to be send to other processors...
bool _have_interface_elems
true if we have cached interface elements, false if they need to be cached. We want to (re)cache only...
void ElementLoopUserObject::post ( )
virtualinherited

Definition at line 216 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

217 {
218 }
void ElementLoopUserObject::pre ( )
virtualinherited

Definition at line 136 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

137 {
138 }
void SlopeLimitingBase::serialize ( std::string &  serialized_buffer)
protectedvirtualinherited

Definition at line 76 of file SlopeLimitingBase.C.

Referenced by SlopeLimitingBase::finalize().

77 {
78  std::ostringstream oss;
79 
80  // First store the number of elements to send
81  unsigned int size = _interface_elem_ids.size();
82  oss.write((char *)&size, sizeof(size));
83 
84  for (auto it = _interface_elem_ids.begin(); it != _interface_elem_ids.end(); ++it)
85  {
86  storeHelper(oss, *it, this);
87  storeHelper(oss, _lslope[*it], this);
88  }
89 
90  // Populate the passed in string pointer with the string stream's buffer contents
91  serialized_buffer.assign(oss.str());
92 }
std::set< dof_id_type > _interface_elem_ids
List of element IDs that are on the processor boundary and need to be send to other processors...
std::map< dof_id_type, std::vector< RealGradient > > _lslope
store the updated slopes into this map indexed by element ID
void ElementLoopUserObject::subdomainChanged ( )
virtualinherited

Definition at line 141 of file ElementLoopUserObject.C.

Referenced by ElementLoopUserObject::execute().

142 {
143 }

Member Data Documentation

const MooseArray<Real>& ElementLoopUserObject::_coord
protectedinherited

Definition at line 100 of file ElementLoopUserObject.h.

const Elem* ElementLoopUserObject::_current_elem
protectedinherited
const Real& ElementLoopUserObject::_current_elem_volume
protectedinherited
const Elem* ElementLoopUserObject::_current_neighbor
protectedinherited
unsigned int ElementLoopUserObject::_current_side
protectedinherited

Definition at line 94 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::onBoundary().

bool ElementLoopUserObject::_have_interface_elems
protectedinherited

true if we have cached interface elements, false if they need to be cached. We want to (re)cache only when mesh changed

Definition at line 103 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::finalize(), ElementLoopUserObject::meshChanged(), ElementLoopUserObject::onInterface(), and ElementLoopUserObject::onInternalSide().

bool SlopeLimitingBase::_include_bc
protectedinherited

option whether to include BCs

Definition at line 48 of file SlopeLimitingBase.h.

Referenced by CNSFVMinmaxSlopeLimiting::limitElementSlope().

std::set<dof_id_type> ElementLoopUserObject::_interface_elem_ids
protectedinherited

List of element IDs that are on the processor boundary and need to be send to other processors.

Definition at line 105 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::meshChanged(), ElementLoopUserObject::onInterface(), ElementLoopUserObject::onInternalSide(), SlopeLimitingBase::serialize(), and SlopeReconstructionBase::serialize().

const MooseArray<Real>& ElementLoopUserObject::_JxW
protectedinherited

Definition at line 99 of file ElementLoopUserObject.h.

const MooseArray<Real>& SlopeLimitingBase::_JxW_face
protectedinherited

Definition at line 56 of file SlopeLimitingBase.h.

std::map<dof_id_type, std::vector<RealGradient> > SlopeLimitingBase::_lslope
protectedinherited
MooseMesh& ElementLoopUserObject::_mesh
protectedinherited
const Elem*& SlopeLimitingBase::_neighbor_elem
protectedinherited

the neighboring element

Definition at line 66 of file SlopeLimitingBase.h.

const MooseArray<Point>& SlopeLimitingBase::_normals_face
protectedinherited

Definition at line 57 of file SlopeLimitingBase.h.

SubdomainID ElementLoopUserObject::_old_subdomain
protectedinherited

The subdomain for the last element.

Definition at line 111 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::execute().

const MooseArray<Point>& ElementLoopUserObject::_q_point
protectedinherited

Definition at line 97 of file ElementLoopUserObject.h.

const MooseArray<Point>& SlopeLimitingBase::_q_point_face
protectedinherited

required data for face assembly

Definition at line 54 of file SlopeLimitingBase.h.

QBase*& ElementLoopUserObject::_qrule
protectedinherited

Definition at line 98 of file ElementLoopUserObject.h.

QBase*& SlopeLimitingBase::_qrule_face
protectedinherited

Definition at line 55 of file SlopeLimitingBase.h.

const SlopeReconstructionBase& SlopeLimitingBase::_rslope
protectedinherited

slope reconstruction user object

Definition at line 51 of file SlopeLimitingBase.h.

Referenced by limitElementSlope(), CNSFVMinmaxSlopeLimiting::limitElementSlope(), and CNSFVWENOSlopeLimiting::limitElementSlope().

unsigned int& SlopeLimitingBase::_side
protectedinherited

current side of the current element

Definition at line 60 of file SlopeLimitingBase.h.

const Elem*& SlopeLimitingBase::_side_elem
protectedinherited

Definition at line 62 of file SlopeLimitingBase.h.

const Real& SlopeLimitingBase::_side_volume
protectedinherited

Definition at line 63 of file SlopeLimitingBase.h.

SubdomainID ElementLoopUserObject::_subdomain
protectedinherited

The subdomain for the current element.

Definition at line 108 of file ElementLoopUserObject.h.

Referenced by ElementLoopUserObject::execute().


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