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 SettingsSSOLIIIAXI_H__
#define SettingsSSOLIIIAXI_H__
 
#include "constants.hpp"
#include <deal.II/base/types.h>
 
using namespace dealii;
 
class SettingsSSOLIIIAXI : public Constants::Physics
{
public:
  SettingsSSOLIIIAXI(){};
 
  const unsigned int nr_threads_max = 0;
 
  const double mu_0 = permeability_fs;
 
  const double mu_r = 4.0;
 
  const double mu = mu_r * mu_0;
 
  const double K_0 = 1.0;
 
  const double d1 = 0.1;
 
  const double rd1 = sqrt(2) * d1;
 
  const double a1 = 0.3;
 
  const double b1 = 0.6;
 
  const double a2 = 0.9;
 
  const double b2 = 1.2;
 
  const double d2 = 1.5;
 
  const double d3 = 3.0;
 
  const double H_0 = (1.0 / 3.0) * K_0 * (pow(b2, 2) - pow(a2, 2));
 
  const types::material_id mid_1 = 1;
 
  const types::material_id mid_2 = 2;
 
  const types::material_id mid_3 = 3;
 
  const types::boundary_id bid_infty = 4;
 
  const types::boundary_id bid_axi = 1;
 
  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 = SettingsSSOLIIIAXI;
 
#endif