parametres.hpp

Go to the documentation of this file.

//Copyright 2017 Laurent Monasse

/*
  This file is part of CELIA3D.
  
  CELIA3D is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
  
  CELIA3D 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 General Public License for more details.
  
  You should have received a copy of the GNU General Public License
  along with CELIA3D.  If not, see <http://www.gnu.org/licenses/>.
*/

#include <cstdlib>
#include <iostream>
#include <fstream>
#include <stdio.h>
#include <sstream>
#ifndef PARAMETRES_H
#define PARAMETRES_H

using std::string;

int rep = false;
int numrep = 18;  



bool explicite = true;
bool flag_2d = false;
bool couplage1d=false;
const bool exact_swap=false;

//Fluid parameters
const double gam = 1.4;                   
const double eps =  0.00000000000001;     
const double epsa = 0.5;                  
const int ordremax = 11;                  
const int marge = 6;                      
const double eps_vide =  0.0000000001;     
const int N_dim=3;

const double X0 = 0;              
const double Y0 = 0;
const double Z0 = 0;

const int Nx =18;                 
const int Ny =18;                 
const int Nz =18;                 

const double domainex = 0.3;            
const double domainey = 0.3;          
const double domainez = 0.3;          

const double deltax = domainex/Nx;      
const double deltay = domainey/Ny;       
const double deltaz = domainez/Nz;       

//Solid parameters
const double rhos = 2698.9; 
const double nu = 0.; 
const double E = 5.; 
const double k_max = 0.01; 

//Parametres temprels
const double T = 0.1;             
const double cfl = 0.5;            
const double cfls = 0.5;           
const int nimp = 10;                
const double dtimp = T/nimp;        
const int Nmax = 1000000;           


const int BC_x_in =  2;                 
const int BC_x_out = 2;                 
const int BC_y_in =  2;                 
const int BC_y_out = 2;                 
const int BC_z_in =  2;                 
const int BC_z_out = 2;                 

double Rho(double x = 0.,double y = 0., double z = 0.);

double U(double x = 0.,double y = 0., double z = 0.);

double V(double x = 0.,double y = 0., double z = 0.); 

double W(double x = 0.,double y = 0., double z = 0.);

double P(double x = 0.,double y = 0., double z = 0., double dx = deltax, double dy = deltay, double dz = deltaz); 

#endif