Problem ( GetPot const& data_file ) :
M_Tstart (0.),
M_saveEvery (1),
M_returnValue (EXIT_FAILURE)
{
using namespace LifeV;
FSIOperator::solid_Type::material_Type::StructureMaterialFactory::instance().registerProduct ( "linearVenantKirchhoff", &FSIOperator::createVenantKirchhoffLinear );
FSIOperator::solid_Type::material_Type::StructureMaterialFactory::instance().registerProduct ( "exponential", &FSIOperator::createExponentialMaterialNonLinear );
FSIOperator::solid_Type::material_Type::StructureMaterialFactory::instance().registerProduct ( "neoHookean", &FSIOperator::createNeoHookeanMaterialNonLinear );
FSIOperator::solid_Type::material_Type::StructureMaterialFactory::instance().registerProduct ( "nonLinearVenantKirchhoff", &FSIOperator::createVenantKirchhoffNonLinear );
std::cout << "register MonolithicGE : " << FSIMonolithicGE::S_register << std::endl;
std::cout << "register MonolithicGI : " << FSIMonolithicGI::S_register << std::endl;
M_data = dataPtr_Type ( new data_Type() );
M_data->setup ( data_file );
//M_data->dataSolid()->setTimeData( M_data->dataFluid()->dataTime() ); //Same TimeData for fluid & solid
//M_data->showMe();
M_fsi = fsi_solver_ptr ( new FSISolver( ) );
MPI_Barrier ( MPI_COMM_WORLD );
M_fsi->setData ( M_data );
M_fsi->FSIOper()->setDataFile ( data_file ); //TO BE REMOVED!
// Setting FESpace and DOF
std::string fluidMeshPartitioned = data_file ( "problem/fluidMeshPartitioned", "none" );
std::string solidMeshPartitioned = data_file ( "problem/solidMeshPartitioned", "none" );
#ifdef HAVE_HDF5
if ( fluidMeshPartitioned.compare ( "none" ) )
{
FSIOperator::meshFilter_Type fluidMeshFilter ( data_file, fluidMeshPartitioned );
fluidMeshFilter.setComm ( M_fsi->FSIOper()->worldComm() );
FSIOperator::meshFilter_Type solidMeshFilter ( data_file, solidMeshPartitioned );
solidMeshFilter.setComm ( M_fsi->FSIOper( )->worldComm( ) );
M_fsi->FSIOper( )->partitionMeshes ( fluidMeshFilter, solidMeshFilter );
M_fsi->FSIOper( )->setupFEspace( );
M_fsi->FSIOper( )->setupDOF ( fluidMeshFilter );
fluidMeshFilter.closeFile( );
solidMeshFilter.closeFile( );
}
else
#endif
{
M_fsi->FSIOper( )->partitionMeshes( );
M_fsi->FSIOper( )->setupFEspace( );
M_fsi->FSIOper( )->setupDOF( );
}
#ifdef DEBUG
debugStream ( 10000 ) << "Setting up the FESpace and DOF \n";
#endif
MPI_Barrier ( MPI_COMM_WORLD );
#ifdef DEBUG
debugStream ( 10000 ) << "Setting up the BC \n";
#endif
M_fsi->setFluidBC ( BCh_monolithicFlux ( true ) );
M_fsi->setSolidBC ( BCh_monolithicSolid ( *M_fsi->FSIOper( ) ) );
M_fsi->setup();
M_fsi->setFluidBC ( BCh_monolithicFluid ( *M_fsi->FSIOper( ), true ) );
M_fsi->setHarmonicExtensionBC ( BCh_harmonicExtension ( *M_fsi->FSIOper( ) ) );
dynamic_cast<LifeV::FSIMonolithic*> (M_fsi->FSIOper().get() )->mergeBCHandlers();
#ifdef DEBUG
debugStream ( 10000 ) << "BC set\n";
#endif
std::string const exporterType = data_file ( "exporter/type", "ensight" );
std::string const fluidName = data_file ( "exporter/fluid/filename", "fluid" );
std::string const solidName = data_file ( "exporter/solid/filename", "solid" );
#ifdef HAVE_HDF5
if (exporterType.compare ("hdf5") == 0)
{
M_exporterFluid.reset ( new hdf5Filter_Type ( data_file, fluidName) );
M_exporterSolid.reset ( new hdf5Filter_Type ( data_file, solidName) );
}
else
#endif
{
if (exporterType.compare ("none") == 0)
{
M_exporterFluid.reset ( new ExporterEmpty<RegionMesh<LinearTetra> > ( data_file, M_fsi->FSIOper()->uFESpace().mesh(), fluidName, M_fsi->FSIOper()->uFESpace().map().comm().MyPID() ) );
M_exporterSolid.reset ( new ExporterEmpty<RegionMesh<LinearTetra> > ( data_file, M_fsi->FSIOper()->dFESpace().mesh(), solidName, M_fsi->FSIOper()->uFESpace().map().comm().MyPID() ) );
}
else
{
M_exporterFluid.reset ( new ensightFilter_Type ( data_file, fluidName) );
M_exporterSolid.reset ( new ensightFilter_Type ( data_file, solidName) );
}
}
// load using ensight/hdf5
M_saveEvery = data_file ("exporter/saveEvery", 1);
// load using ensight/hdf5
std::string restartType (data_file ("importer/restartFSI", "false" ) );
std::cout << "The load state is: " << restartType << std::endl;
if ( !restartType.compare ("true") )
{
restartFSI (data_file);
}
else
{
M_fsi->initialize();
M_velAndPressure.reset ( new vector_Type ( M_fsi->FSIOper()->fluid().getMap(), M_exporterFluid->mapType() ) );
M_fluidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->mmFESpace().map(), M_exporterFluid->mapType() ) );
M_solidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->dFESpace().map(), M_exporterSolid->mapType() ) );
}
M_exporterFluid->setMeshProcId (M_fsi->FSIOper()->uFESpace().mesh(), M_fsi->FSIOper()->uFESpace().map().comm().MyPID() );
M_exporterSolid->setMeshProcId (M_fsi->FSIOper()->dFESpace().mesh(), M_fsi->FSIOper()->dFESpace().map().comm().MyPID() );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "f-velocity",
M_fsi->FSIOper()->uFESpacePtr(), M_velAndPressure, UInt (0) );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::ScalarField, "f-pressure",
M_fsi->FSIOper()->pFESpacePtr(), M_velAndPressure,
UInt (3 * M_fsi->FSIOper()->uFESpace().dof().numTotalDof() ) );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "f-displacement",
M_fsi->FSIOper()->mmFESpacePtr(), M_fluidDisp, UInt (0) );
M_exporterSolid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "s-displacement",
M_fsi->FSIOper()->dFESpacePtr(), M_solidDisp, UInt (0) );
//M_fsi->FSIOper()->fluid().setupPostProc(); //this has to be called if we want to initialize the postProcess
FC0.initParameters ( *M_fsi->FSIOper(), 3);
M_data->dataFluid()->dataTime()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataFluid()->dataTime()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataSolid()->dataTime()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataSolid()->dataTime()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->timeDataALE()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->timeDataALE()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
}
Problem ( GetPot const& data_file, boost::shared_ptr<Epetra_Comm> comm) :
M_Tstart (0.),
M_tolSave (1),
M_saveEvery (1),
M_comm ( comm )
{
using namespace LifeV;
// Building the communicator for output purposes
//M_comm.reset( new Epetra_Comm( MPI_COMM_WORLD ) );
M_data = dataPtr_Type ( new data_Type() );
M_data->setup ( data_file );
M_data->showMe();
#ifdef DEBUG
debugStream ( 10000 ) << "creating FSISolver with operator : " << method << "\n";
#endif
M_fsi = fsi_solver_ptr ( new FSISolver( ) );
MPI_Barrier ( MPI_COMM_WORLD );
#ifdef DEBUG
debugStream ( 10000 ) << "Setting up data from GetPot \n";
#endif
M_fsi->setData ( M_data );
M_fsi->FSIOper()->setDataFile ( data_file ); //TO BE REMOVED!
MPI_Barrier ( MPI_COMM_WORLD );
// Setting FESpace and DOF
std::string fluidMeshPartitioned = data_file ( "problem/fluidMeshPartitioned", "none" );
std::string solidMeshPartitioned = data_file ( "problem/solidMeshPartitioned", "none" );
#ifdef HAVE_HDF5
if ( fluidMeshPartitioned.compare ( "none" ) )
{
FSIOperator::meshFilter_Type fluidMeshFilter ( data_file, fluidMeshPartitioned );
fluidMeshFilter.setComm ( M_fsi->FSIOper()->worldComm() );
FSIOperator::meshFilter_Type solidMeshFilter ( data_file, solidMeshPartitioned );
solidMeshFilter.setComm ( M_fsi->FSIOper( )->worldComm( ) );
M_fsi->FSIOper( )->partitionMeshes ( fluidMeshFilter, solidMeshFilter );
M_fsi->FSIOper( )->setupFEspace( );
M_fsi->FSIOper( )->setupDOF ( fluidMeshFilter );
fluidMeshFilter.closeFile( );
solidMeshFilter.closeFile( );
}
else
#endif
{
M_fsi->FSIOper( )->partitionMeshes( );
M_fsi->FSIOper( )->setupFEspace( );
M_fsi->FSIOper( )->setupDOF( );
}
std::cout << "register MonolithicGI : " << FSIMonolithicGI::S_register << std::endl;
debugStream ( 10000 ) << "Setting up the FESpace and DOF \n";
MPI_Barrier ( MPI_COMM_WORLD );
#ifdef DEBUG
debugStream ( 10000 ) << "Setting up the BC \n";
#endif
M_fsi->setFluidBC ( BCh_monolithicFlux ( true ) );
M_fsi->setSolidBC ( BCh_monolithicSolid ( *M_fsi->FSIOper( ) ) );
M_fsi->setup (/*data_file*/);
M_fsi->setFluidBC ( BCh_monolithicFluid ( *M_fsi->FSIOper( ), true ) );
M_fsi->setHarmonicExtensionBC ( BCh_harmonicExtension ( *M_fsi->FSIOper( ) ) );
dynamic_cast<LifeV::FSIMonolithic*> (M_fsi->FSIOper().get() )->mergeBCHandlers();
#ifdef DEBUG
debugStream ( 10000 ) << "BC set\n";
#endif
std::string const exporterType = data_file ( "exporter/type", "ensight" );
std::string const fluidName = data_file ( "exporter/fluid/filename", "fluid" );
std::string const solidName = data_file ( "exporter/solid/filename", "solid" );
#ifdef HAVE_HDF5
if (exporterType.compare ("hdf5") == 0)
{
M_exporterFluid.reset ( new hdf5Filter_Type ( data_file, fluidName) );
M_exporterSolid.reset ( new hdf5Filter_Type ( data_file, solidName) );
}
else
#endif
{
if (exporterType.compare ("none") == 0)
{
M_exporterFluid.reset ( new ExporterEmpty<RegionMesh<LinearTetra> > ( data_file, M_fsi->FSIOper()->uFESpace().mesh(), fluidName, M_fsi->FSIOper()->uFESpace().map().comm().MyPID() ) );
M_exporterSolid.reset ( new ExporterEmpty<RegionMesh<LinearTetra> > ( data_file, M_fsi->FSIOper()->dFESpace().mesh(), solidName, M_fsi->FSIOper()->uFESpace().map().comm().MyPID() ) );
}
else
{
M_exporterFluid.reset ( new ensightFilter_Type ( data_file, fluidName) );
M_exporterSolid.reset ( new ensightFilter_Type ( data_file, solidName) );
}
}
//reading the saveEvery
M_saveEvery = data_file ("exporter/saveEvery", 1);
M_tolSave = data_file ("exporter/tolSave", 1);
// load using ensight/hdf5
std::string restartType (data_file ("importer/restartFSI", "false") );
if (!restartType.compare ("true") )
{
restartFSI (data_file);
}
else if ( !restartType.compare ("vectors") )
{
initializeWithVectors( );
M_velAndPressure.reset ( new vector_Type ( M_fsi->FSIOper()->fluid().getMap(), LifeV::Unique ) );
M_fluidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->mmFESpace().map(), LifeV::Unique ) );
M_solidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->dFESpace().map(), LifeV::Unique ) );
// M_solidVel.reset ( new vector_Type( M_fsi->FSIOper()->dFESpace().map(), M_exporterSolid->mapType() ));
M_WS.reset ( new vector_Type ( M_fsi->FSIOper()->dFESpace().map(), LifeV::Unique ) );
}
else
{
M_fsi->initializeMonolithicOperator();
M_velAndPressure.reset ( new vector_Type ( M_fsi->FSIOper()->fluid().getMap(), LifeV::Unique ) );
M_fluidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->mmFESpace().map(), LifeV::Unique ) );
M_solidDisp.reset ( new vector_Type ( M_fsi->FSIOper()->dFESpace().map(), LifeV::Unique ) );
// M_solidVel.reset ( new vector_Type( M_fsi->FSIOper()->dFESpace().map(), M_exporterSolid->mapType() ));
M_WS.reset ( new vector_Type ( M_fsi->FSIOper()->dFESpace().map(), LifeV::Unique ) );
}
M_exporterFluid->setMeshProcId (M_fsi->FSIOper()->uFESpace().mesh(), M_fsi->FSIOper()->uFESpace().map().comm().MyPID() );
M_exporterSolid->setMeshProcId (M_fsi->FSIOper()->dFESpace().mesh(), M_fsi->FSIOper()->dFESpace().map().comm().MyPID() );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "f-velocity",
M_fsi->FSIOper()->uFESpacePtr(), M_velAndPressure, UInt (0) );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::ScalarField, "f-pressure",
M_fsi->FSIOper()->pFESpacePtr(), M_velAndPressure,
UInt (3 * M_fsi->FSIOper()->uFESpace().dof().numTotalDof() ) );
M_exporterFluid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "f-displacement",
M_fsi->FSIOper()->mmFESpacePtr(), M_fluidDisp, UInt (0) );
M_exporterSolid->addVariable ( ExporterData<FSIOperator::mesh_Type>::VectorField, "s-displacement",
M_fsi->FSIOper()->dFESpacePtr(), M_solidDisp, UInt (0) );
// M_exporterSolid->addVariable( ExporterData<FSIOperator::mesh_Type>::VectorField, "s-velocity",
// M_fsi->FSIOper()->dFESpacePtr(), M_solidVel, UInt(0) );
// M_exporterSolid->addVariable( ExporterData<FSIOperator::mesh_Type>::VectorField, "s-ws",
// M_fsi->FSIOper()->dFESpacePtr(), M_WS, UInt(0) );
M_fsi->FSIOper()->fluid().setupPostProc(); //this has to be called if we want to initialize the postProcess
// Initializing either resistance of absorbing BCs
// At the moment, the resistance BCs are applied explicitly in order to limit the ill-conditioning
// of the linear system
Real resistance = data_file ("fluid/physics/resistance", 0.0);
Real hydrostatic = data_file ("fluid/physics/hydrostatic", 0.0);
R1.initParameters ( OUTLET, resistance, hydrostatic, "outlet-3" );
//FC2.initParameters ( *M_fsi->FSIOper(), OUTLET);
M_data->dataFluid()->dataTime()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataFluid()->dataTime()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataSolid()->dataTime()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->dataSolid()->dataTime()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->timeDataALE()->setInitialTime ( M_data->dataFluid()->dataTime()->initialTime() );
M_data->timeDataALE()->setTime ( M_data->dataFluid()->dataTime()->initialTime() );
}
