settings.hpp

/******************************************************************************
 * Copyright (C) Siarhei Uzunbajakau, 2023.
 *
 * This program is free software. You can use, modify, and redistribute it under
 * the terms of the GNU Lesser General Public License as published by the Free
 * Software Foundation, either version 3 or (at your option) any later version.
 * This program is distributed without any warranty.
 *
 * Refer to COPYING.LESSER for more details.
 ******************************************************************************/
 
#ifndef SettingsSLDII_H__
#define SettingsSLDII_H__
 
#include "constants.hpp"
#include <deal.II/base/types.h>
 
using namespace dealii;
 
enum BoundaryConditionType
{
  Dirichlet,
  Exact
};
 
class SettingsSLDII : public Constants::Physics
{
public:
  SettingsSLDII()
    : j_hat({ 1.0, 0.0 })
    , k_hat({ 0.0, 0.0, 1.0 }){};
 
  const unsigned int nr_threads_max = 0;
 
  const double mu_0 = permeability_fs;
 
  const double d1 = (DIMENSION__ == 2) ? 0.07 : 0.05;
 
  const double rd1 = sqrt(DIMENSION__) * d1;
 
  const double a = 0.2;
 
  const double b = 0.4;
 
  const double d_2 = 0.8;
 
  const double d_3 = 2.0;
 
  const Tensor<1, 2> j_hat;
 
  const Tensor<1, 3> k_hat;
 
  const types::boundary_id bid = 1;
 
  const BoundaryConditionType type_of_bc = Exact;
 
  const types::material_id mid_1 = 1;
 
  const types::material_id mid_2 = 2;
 
  const types::material_id mid_3 = 3;
 
  const double mur_1 = 1.0;
 
  const double mur_2 = 4.0;
 
  const double mur_3 = 1.0;
 
  const double mu_1 = mur_1 * mu_0;
 
  const double mu_2 = mur_2 * mu_0;
 
  const double mu_3 = mur_3 * mu_0;
 
  const double H_0 = 1.0;
 
  const double eps = 1e-12;
 
  const bool print_time_tables = false;
 
  const bool project_exact_solution = false;
 
  const bool log_cg_convergence = false;
};
 
using Settings = SettingsSLDII;
 
#endif