aslFDElChemKinetics.h

Go to the documentation of this file.

/*
 * Advanced Simulation Library <http://asl.org.il>
 * 
 * Copyright 2015 Avtech Scientific <http://avtechscientific.com>
 *
 *
 * This file is part of Advanced Simulation Library (ASL).
 *
 * ASL is free software: you can redistribute it and/or modify it
 * under the terms of the GNU Affero General Public License as
 * published by the Free Software Foundation, version 3 of the License.
 *
 * ASL is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with ASL. If not, see <http://www.gnu.org/licenses/>.
 *
 */
 
 
#ifndef ASLFDELCHEMKINETICS_H
#define ASLFDELCHEMKINETICS_H
 
#include "aslSingleKernelNM.h"
#include "acl/aclMath/aclVectorOfElementsDef.h"
 
 
namespace acl
{
  class VectorOfElementsData;
  class VectorOfElements;
}
 
namespace asl
{
  class VectorTemplate;
  template <typename V> class DataWithGhostNodes;
  typedef DataWithGhostNodes<acl::VectorOfElementsData> DataWithGhostNodesACLData;
  typedef std::shared_ptr<DataWithGhostNodesACLData> SPDataWithGhostNodesACLData;
  

  class FDBVKinetics: public SingleKernelMapNM
  {
    public:
      typedef SPDataWithGhostNodesACLData Data;
      typedef SPAbstractDataWithGhostNodes Field;
      typedef acl::VectorOfElements Param;
      
    private:
      acl::SPKernel kernelJ;
      std::vector<Data> aI;     
      Data efSPhi;
      Field phi;
 
      vector<double> nI;
      Param j0;
      Param beta;
      double n;
 
      virtual void init0();
    public:     
      FDBVKinetics(Data a0,
                   double n0,
                   Data aI,
                   double nI,
                   Data phiS,
                   Field phi,
                   const Param & j0,
                   const Param & beta,
                   double n);
      void setElectricFieldSource(Field phi);
      Field getElectricFieldSource() const;
      inline Data & getAI(unsigned int i=0);
      void addAI(Data ai, double ni);     
      void executeJ();
  };
 
  typedef std::shared_ptr<FDBVKinetics> SPFDBVKinetics;
 
 
  SPFDBVKinetics generateFDBVKinetics(SPDataWithGhostNodesACLData a0,
                                      double n0,
                                      SPDataWithGhostNodesACLData aI,
                                      double nI,
                                      SPDataWithGhostNodesACLData phiS,
                                      SPAbstractDataWithGhostNodes phi,
                                      double j0,
                                      double beta,
                                      double n);
  
 
  
// ------------------------- Implementation ------------------------
  
  inline FDBVKinetics::Data & FDBVKinetics::getAI(unsigned int i)
  {
    return aI[i];
  }
    
} // asl
#endif // ASLFDADVECTIONDIFFUSION_H