single_layer_spine_mesh.template.h
Go to the documentation of this file.
1 // LIC// ====================================================================
2 // LIC// This file forms part of oomph-lib, the object-oriented,
3 // LIC// multi-physics finite-element library, available
4 // LIC// at http://www.oomph-lib.org.
5 // LIC//
6 // LIC// Copyright (C) 2006-2023 Matthias Heil and Andrew Hazel
7 // LIC//
8 // LIC// This library is free software; you can redistribute it and/or
9 // LIC// modify it under the terms of the GNU Lesser General Public
10 // LIC// License as published by the Free Software Foundation; either
11 // LIC// version 2.1 of the License, or (at your option) any later version.
12 // LIC//
13 // LIC// This library is distributed in the hope that it will be useful,
14 // LIC// but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // LIC// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // LIC// Lesser General Public License for more details.
17 // LIC//
18 // LIC// You should have received a copy of the GNU Lesser General Public
19 // LIC// License along with this library; if not, write to the Free Software
20 // LIC// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
21 // LIC// 02110-1301 USA.
22 // LIC//
23 // LIC// The authors may be contacted at oomph-lib@maths.man.ac.uk.
24 // LIC//
25 // LIC//====================================================================
26 #ifndef OOMPH_SINGLE_LAYER_SPINE_MESH_HEADER
27 #define OOMPH_SINGLE_LAYER_SPINE_MESH_HEADER
28 
29 // oomph-lib includes
30 #include "../generic/spines.h"
32 
33 namespace oomph
34 {
35  //======================================================================
36  /// Single-layer spine mesh class derived from standard 2D mesh.
37  /// The mesh contains a layer of spinified fluid elements (of type ELEMENT;
38  /// e.g SpineElement<QCrouzeixRaviartElement<2>).
39  /// It constructs the spines and contains the information on how to update
40  /// the nodal positions within the mesh as a function of the spine lengths.
41  /// Equations that determine the spine heights (even if they are pinned)
42  /// must be specified externally or else there will be problems.
43  //======================================================================
44  template<class ELEMENT>
45  class SingleLayerSpineMesh : public RectangularQuadMesh<ELEMENT>,
46  public SpineMesh
47  {
48  public:
49  /// Constructor: Pass number of elements in x-direction, number of
50  /// elements in y-direction, axial length, height of layer, and pointer
51  /// to timestepper (defaults to Steady timestepper)
53  const unsigned& nx,
54  const unsigned& ny,
55  const double& lx,
56  const double& h,
57  TimeStepper* time_stepper_pt = &Mesh::Default_TimeStepper);
58 
59 
60  /// Constructor: Pass number of elements in x-direction, number of
61  /// elements in y-direction, axial length, height of layer, a boolean
62  /// flag to make the mesh periodic in the x-direction, and a pointer
63  /// to timestepper (defaults to Steady timestepper)
65  const unsigned& nx,
66  const unsigned& ny,
67  const double& lx,
68  const double& h,
69  const bool& periodic_in_x,
70  TimeStepper* time_stepper_pt = &Mesh::Default_TimeStepper);
71 
72  /// General node update function implements pure virtual function
73  /// defined in SpineMesh base class and performs specific node update
74  /// actions: along vertical spines
75  virtual void spine_node_update(SpineNode* spine_node_pt)
76  {
77  // Get fraction along the spine
78  double W = spine_node_pt->fraction();
79  // Get spine height
80  double H = spine_node_pt->h();
81  // Set the value of y
82  spine_node_pt->x(1) = this->Ymin + W * H;
83  }
84 
85 
86  protected:
87  /// Helper function to actually build the single-layer spine mesh
88  /// (called from various constructors)
89  virtual void build_single_layer_mesh(TimeStepper* time_stepper_pt);
90  };
91 
92 } // namespace oomph
93 
94 #endif
RectangularQuadMesh is a two-dimensional mesh of Quad elements with Nx elements in the "x" (horizonal...
const unsigned & ny() const
Return number of elements in y direction.
const unsigned & nx() const
Return number of elements in x direction.
double Ymin
Minimum value of y coordinate.
Single-layer spine mesh class derived from standard 2D mesh. The mesh contains a layer of spinified f...
SingleLayerSpineMesh(const unsigned &nx, const unsigned &ny, const double &lx, const double &h, TimeStepper *time_stepper_pt=&Mesh::Default_TimeStepper)
Constructor: Pass number of elements in x-direction, number of elements in y-direction,...
virtual void spine_node_update(SpineNode *spine_node_pt)
General node update function implements pure virtual function defined in SpineMesh base class and per...
virtual void build_single_layer_mesh(TimeStepper *time_stepper_pt)
Helper function to actually build the single-layer spine mesh (called from various constructors)
////////////////////////////////////////////////////////////////////// //////////////////////////////...