fsi/adlib-amroc/DetonatingCanister/src/SolidProblem.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@amroc.net

#ifndef ADLIB_SOLIDPROBLEM_H
#define ADLIB_SOLIDPROBLEM_H

#define OWN_ELC_COUPLED_ADLIBSOLVER
#include "external_coupling/AdlibStdELCCoupledProblem.h"

class SolidSolverSpecific : public AdlibELCCoupledSolver<DIM> {
  typedef AdlibELCCoupledSolver<DIM> base;
public:
  SolidSolverSpecific(MPI_Comm solidComm, int numFluidNodes, int firstFluidNode) : 
  base(solidComm, numFluidNodes, firstFluidNode) {}

  virtual void Update(const double& dt) {
    // Boundary conditions 
    // No deformations behind initial detonation location,
    // in z-direction at far away upper end, 
    // and in x-,y-direction at nodes at outer radius
    for (register int n=0; n<adlib::nodes; n++) {
      int offset=3*n;
      double &x = adlib::coordinates[offset];
      double &y = adlib::coordinates[offset+1];
      double &z = adlib::coordinates[offset+2];
      double radius = std::sqrt(x*x+y*y);
      if (radius>=0.017) {
        for (register int m=0; m<2; m++) {
          adlib::boundary[offset+m] = 1;
          adlib::velocities[offset+m] = 0;
          adlib::displacements[offset+m] = 0;
        }
      }
      if (z>=0.095) {
        adlib::boundary[offset+2] = 1;
        adlib::velocities[offset+2] = 0;
        adlib::displacements[offset+2] = 0;
      }
      if (z<=0.015) {
        for (register int m=0; m<3; m++) {
          adlib::boundary[offset+m] = 1;
          adlib::velocities[offset+m] = 0;
          adlib::displacements[offset+m] = 0;
        }
      }
    }
    base::Update(dt);
  }
};


#endif

---