fsi/sfc-amroc/WaterBlastPlastic/src/ShellEnforceBCFunctor.h

Go to the documentation of this file.

// -*- C++ -*-
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// 
//                           Ralf Deiterding, Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2005 All Rights Reserved
//
// <LicenseText>
//
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#ifndef SHELLENFORCEBCFUNCTOR_H
#define SHELLENFORCEBCFUNCTOR_H

#include "shells/driverCC/SVertexFunctors.h"
#include "shells/driverCC/PrescribeVarFunctor.h"

#include <functional>


namespace shells {
    class ShellEnforceBCFunctor;
}


class shells::ShellEnforceBCFunctor :
    public std::unary_function<shells::SVertexS*, void> {
    
public:
    ShellEnforceBCFunctor() {}
    
    ~ShellEnforceBCFunctor() {}

    void operator()(shells::SVertex *vtx) const {
        SVertexCoordinate extractC;
        SVertexCoordinate::DataType *coor = extractC(vtx);      
        const double xcoor = coor[0];
        const double ycoor = coor[1];
        const double zcoor = coor[2];
            
        // set the velocities and accelerations to zero at the boundaries
        if (std::sqrt(ycoor*ycoor+zcoor*zcoor)>0.041) {
            SVertexVelocity extractV;
            SVertexAcceleration extractA;

            SVertexVelocity::DataType *vel = extractV(vtx);
            SVertexAcceleration::DataType *acc = extractA(vtx);
            for (unsigned i=0; i<SVertexVelocity::numVar; ++i) {
                vel[i] = 0.0;
                acc[i] = 0.0;
            }
        }       

        if (std::sqrt(ycoor*ycoor+zcoor*zcoor)>0.032) {
            SVertexVelocity extractV;
            SVertexAcceleration extractA;

            SVertexVelocity::DataType *vel = extractV(vtx);
            SVertexAcceleration::DataType *acc = extractA(vtx);
            for (unsigned i=0; i<SVertexVelocity::numVar; ++i) {
                vel[0] = 0.0;
                acc[0] = 0.0;
            }
        }       
        return;
    }    
};
#endif

---