biharmonic.h

Go to the documentation of this file.

#ifndef BIHARMONIC_H
#define BIHARMONIC_H
 
#include "libmesh/equation_systems.h"
#include "libmesh/replicated_mesh.h"
#include "libmesh/exodusII_io.h"
#include "libmesh/mesh_refinement.h"
 
// Bring in bits from the libMesh namespace.
// Just the bits we're using, since this is a header.
using libMesh::EquationSystems;
using libMesh::ExodusII_IO;
using libMesh::Point;
using libMesh::Real;
using libMesh::ReplicatedMesh;
 
#ifdef LIBMESH_ENABLE_AMR
using libMesh::MeshRefinement;
#endif
 
class Biharmonic : public EquationSystems
{
public:
  // Initial state enumeration
  enum InitialStateEnum {STRIP = 0,
                         ROD   = 1,
                         BALL  = 2};
 
  // Free energy enumeration
  enum FreeEnergyEnum {DOUBLE_WELL         = 1,
                       DOUBLE_OBSTACLE     = 2,
                       LOG_DOUBLE_WELL     = 3,
                       LOG_DOUBLE_OBSTACLE = 4};
 
  Biharmonic(ReplicatedMesh & mesh);
 
 
  ~Biharmonic()
  {
    // delete _meshRefinement;
  }
 
 
  // Misc. getters
  bool verbose()         { return _verbose; }
  Real dt0()             { return _dt0; }
  Real dt()              { return _dt; }
 
 
  // Public interface functions
  void viewParameters();
  void init();
  void step(const Real & dt = -1.0);
  void output(int timestep, const Real & t, Real & o_t, bool force = false);
  void run();
 
private:
  unsigned int  _dim, _N;
  Real _kappa, _theta, _theta_c;
  Real _tol;
  bool _growth, _degenerate, _cahn_hillard, _netforce;
  FreeEnergyEnum  _energy;
  int _log_truncation;
  bool _verbose;
  InitialStateEnum  _initialState;
  Point _initialCenter;
  Real _initialWidth;
  Real _dt0, _dt, _t0, _T, _t1;
  Real _cnWeight;
  //
  std::string  _ofile_base, _ofile;
  std::unique_ptr<ExodusII_IO> _exio;
  Real    _o_dt;
  int     _o_count;
  //
  friend class JR;
  class JR;       // forward
  ReplicatedMesh & _mesh;
  JR * _jr;
};
 
#endif // BIHARMONIC_H