mirror of
https://github.com/OPM/ResInsight.git
synced 2025-01-08 23:23:01 -06:00
428 lines
19 KiB
C++
428 lines
19 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) 2023- Equinor ASA
|
|
//
|
|
// ResInsight 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.
|
|
//
|
|
// ResInsight 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 at <http://www.gnu.org/licenses/gpl.html>
|
|
// for more details.
|
|
//
|
|
/////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#include "RifFaultReactivationModelExporter.h"
|
|
|
|
#include "RigFaultReactivationModel.h"
|
|
#include "RigGriddedPart3d.h"
|
|
|
|
#include "RiaApplication.h"
|
|
#include "RiaVersionInfo.h"
|
|
|
|
#include "RifInpExportTools.h"
|
|
|
|
#include <fstream>
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::exportToStream( std::ostream& stream, const RigFaultReactivationModel& model )
|
|
{
|
|
std::string applicationNameAndVersion = STRPRODUCTVER;
|
|
|
|
using PartBorderSurface = RigGriddedPart3d::BorderSurface;
|
|
std::vector<std::pair<PartBorderSurface, std::string>> borders = { { PartBorderSurface::UpperSurface, "top" },
|
|
{ PartBorderSurface::FaultSurface, "fault" },
|
|
{ PartBorderSurface::LowerSurface, "base" } };
|
|
|
|
using FaultGridPart = RigFaultReactivationModel::GridPart;
|
|
std::map<std::pair<FaultGridPart, PartBorderSurface>, int> faces = { { { FaultGridPart::PART1, PartBorderSurface::FaultSurface }, 4 },
|
|
{ { FaultGridPart::PART1, PartBorderSurface::UpperSurface }, 4 },
|
|
{ { FaultGridPart::PART1, PartBorderSurface::LowerSurface }, 4 },
|
|
{ { FaultGridPart::PART2, PartBorderSurface::FaultSurface }, 6 },
|
|
{ { FaultGridPart::PART2, PartBorderSurface::UpperSurface }, 6 },
|
|
{ { FaultGridPart::PART2, PartBorderSurface::LowerSurface }, 6 } };
|
|
|
|
std::map<RigGriddedPart3d::Boundary, std::string> boundaries = {
|
|
{ RigGriddedPart3d::Boundary::Bottom, "bottom" },
|
|
{ RigGriddedPart3d::Boundary::Back, "back" },
|
|
{ RigGriddedPart3d::Boundary::Front, "front" },
|
|
{ RigGriddedPart3d::Boundary::FarSide, "farside" },
|
|
};
|
|
|
|
double faultFriction = 0.0;
|
|
|
|
printHeading( stream, applicationNameAndVersion );
|
|
printParts( stream, model, borders, faces );
|
|
printAssembly( stream, model, boundaries );
|
|
printMaterials( stream );
|
|
printInteractionProperties( stream, faultFriction );
|
|
printBoundaryConditions( stream, boundaries );
|
|
printPredefinedFields( stream );
|
|
printInteractions( stream, borders );
|
|
printSteps( stream );
|
|
|
|
// TODO: improve error handling
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::exportToFile( const std::string& filePath,
|
|
const RigFaultReactivationModel& model )
|
|
{
|
|
std::ofstream stream( filePath );
|
|
return exportToStream( stream, model );
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printHeading( std::ostream& stream, const std::string& applicationNameAndVersion )
|
|
{
|
|
if ( RifInpExportTools::printHeading( stream, "Heading" ) &&
|
|
RifInpExportTools::printComment( stream, std::string( "Generated by: " ).append( applicationNameAndVersion ) ) &&
|
|
RifInpExportTools::printHeading( stream, "Preprint, echo=NO, model=NO, history=NO, contact=NO" ) )
|
|
{
|
|
return { true, "" };
|
|
}
|
|
|
|
return { false, "Failed to write header to fault reactivation INP." };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printParts(
|
|
std::ostream& stream,
|
|
const RigFaultReactivationModel& model,
|
|
const std::vector<std::pair<RigGriddedPart3d::BorderSurface, std::string>>& borders,
|
|
const std::map<std::pair<RigFaultReactivationModel::GridPart, RigGriddedPart3d::BorderSurface>, int>& faces )
|
|
{
|
|
RifInpExportTools::printSectionComment( stream, "PARTS" );
|
|
|
|
auto parts = model.allGridParts();
|
|
int partIndex = 1;
|
|
for ( const auto& part : parts )
|
|
{
|
|
std::string partNameHeading = "Part-" + std::to_string( partIndex );
|
|
RifInpExportTools::printHeading( stream, "Part, name=" + partNameHeading );
|
|
|
|
auto grid = model.grid( part );
|
|
|
|
const std::vector<cvf::Vec3d>& nodes = grid->nodes();
|
|
RifInpExportTools::printNodes( stream, nodes );
|
|
|
|
const std::vector<std::vector<unsigned int>>& elements = grid->elementIndices();
|
|
RifInpExportTools::printElements( stream, elements );
|
|
|
|
std::string partName = "part" + std::to_string( partIndex );
|
|
RifInpExportTools::printNodeSet( stream, partName, 1, nodes.size(), false );
|
|
RifInpExportTools::printElementSet( stream, partName, 1, elements.size() );
|
|
|
|
const std::map<RigGriddedPart3d::BorderSurface, std::vector<unsigned int>>& borderSurfaceElements = grid->borderSurfaceElements();
|
|
|
|
for ( auto [border, borderName] : borders )
|
|
{
|
|
auto elementIt = faces.find( { part, border } );
|
|
CAF_ASSERT( elementIt != faces.end() );
|
|
int elementSide = elementIt->second;
|
|
|
|
std::string sideName = "S" + std::to_string( elementSide );
|
|
auto surfaceElements = borderSurfaceElements.find( border );
|
|
if ( surfaceElements != borderSurfaceElements.end() )
|
|
{
|
|
std::string borderElementName = "_" + borderName + "_" + sideName;
|
|
RifInpExportTools::printElementSet( stream, borderElementName, surfaceElements->second );
|
|
RifInpExportTools::printSurface( stream, borderName, borderElementName, sideName );
|
|
}
|
|
}
|
|
|
|
RifInpExportTools::printComment( stream, "Section: sand" );
|
|
RifInpExportTools::printHeading( stream, "Solid Section, elset=" + partName + ", material=sand" );
|
|
|
|
RifInpExportTools::printLine( stream, "," );
|
|
RifInpExportTools::printHeading( stream, "End Part" );
|
|
|
|
partIndex++;
|
|
}
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string>
|
|
RifFaultReactivationModelExporter::printAssembly( std::ostream& stream,
|
|
const RigFaultReactivationModel& model,
|
|
const std::map<RigGriddedPart3d::Boundary, std::string>& boundaries )
|
|
|
|
{
|
|
// ASSEMBLY part
|
|
RifInpExportTools::printSectionComment( stream, "ASSEMBLY" );
|
|
RifInpExportTools::printHeading( stream, "Assembly, name=Assembly" );
|
|
|
|
auto parts = model.allGridParts();
|
|
|
|
int partIndex = 1;
|
|
for ( const auto& part : parts )
|
|
{
|
|
std::string partName = "Part-" + std::to_string( partIndex );
|
|
std::string instanceName = partName + "-1";
|
|
RifInpExportTools::printHeading( stream, "Instance, name=" + instanceName + ", part=" + partName );
|
|
|
|
std::string nodeSetName = "part_" + std::to_string( partIndex ) + "_PP_";
|
|
auto grid = model.grid( part );
|
|
|
|
const std::vector<cvf::Vec3d>& nodes = grid->nodes();
|
|
RifInpExportTools::printNodeSet( stream, nodeSetName, 1, nodes.size(), true );
|
|
|
|
RifInpExportTools::printHeading( stream, "End Instance" );
|
|
|
|
partIndex++;
|
|
}
|
|
|
|
partIndex = 1;
|
|
for ( const auto& part : parts )
|
|
{
|
|
std::string partName = "Part-" + std::to_string( partIndex );
|
|
std::string instanceName = partName + "-1";
|
|
|
|
for ( auto [boundary, boundaryName] : boundaries )
|
|
{
|
|
// Create boundary condition sets for each side of the parts (except top).
|
|
auto grid = model.grid( part );
|
|
|
|
auto boundaryNodes = grid->boundaryNodes();
|
|
auto boundaryElements = grid->boundaryElements();
|
|
|
|
std::string setName = "Set-" + boundaryName;
|
|
const std::vector<unsigned int>& nodes = boundaryNodes[boundary];
|
|
RifInpExportTools::printNodeSet( stream, setName, instanceName, nodes );
|
|
|
|
const std::vector<unsigned int>& elements = boundaryElements[boundary];
|
|
RifInpExportTools::printElementSet( stream, setName, instanceName, elements );
|
|
}
|
|
partIndex++;
|
|
}
|
|
|
|
RifInpExportTools::printHeading( stream, "End Assembly" );
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printMaterials( std::ostream& stream )
|
|
{
|
|
// MATERIALS PART
|
|
struct Material
|
|
{
|
|
std::string name;
|
|
double density;
|
|
double elastic1;
|
|
double elastic2;
|
|
double permeability1;
|
|
double permeability2;
|
|
};
|
|
|
|
RifInpExportTools::printSectionComment( stream, "MATERIALS" );
|
|
|
|
std::vector<Material> materials = {
|
|
Material{ .name = "sand", .density = 2000.0, .elastic1 = 5e+09, .elastic2 = 0.2, .permeability1 = 1e-09, .permeability2 = 0.3 } };
|
|
for ( Material mat : materials )
|
|
{
|
|
RifInpExportTools::printHeading( stream, "Material, name=" + mat.name );
|
|
RifInpExportTools::printHeading( stream, "Density" );
|
|
RifInpExportTools::printNumber( stream, mat.density );
|
|
RifInpExportTools::printHeading( stream, "Elastic" );
|
|
RifInpExportTools::printNumbers( stream, { mat.elastic1, mat.elastic2 } );
|
|
|
|
RifInpExportTools::printHeading( stream, "Permeability, specific=1." );
|
|
RifInpExportTools::printNumbers( stream, { mat.permeability1, mat.permeability2 } );
|
|
}
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printInteractionProperties( std::ostream& stream, double faultFriction )
|
|
{
|
|
RifInpExportTools::printSectionComment( stream, "INTERACTION PROPERTIES" );
|
|
|
|
// Fault interaction
|
|
RifInpExportTools::printHeading( stream, "Surface Interaction, name=fault" );
|
|
RifInpExportTools::printNumber( stream, 1.0 );
|
|
RifInpExportTools::printHeading( stream, "Friction" );
|
|
RifInpExportTools::printNumber( stream, faultFriction );
|
|
RifInpExportTools::printHeading( stream, "Surface Behavior, no separation, pressure-overclosure=HARD" );
|
|
// Non-fault interaction
|
|
RifInpExportTools::printHeading( stream, "Surface Interaction, name=non-fault" );
|
|
RifInpExportTools::printNumber( stream, 1.0 );
|
|
RifInpExportTools::printHeading( stream, "Cohesive Behavior" );
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string>
|
|
RifFaultReactivationModelExporter::printBoundaryConditions( std::ostream& stream,
|
|
const std::map<RigGriddedPart3d::Boundary, std::string>& boundaries )
|
|
{
|
|
auto printBoundaryCondition =
|
|
[]( std::ostream& stream, const std::string& boundaryName, const std::string& boundarySetName, const std::string& symmetryType )
|
|
{
|
|
RifInpExportTools::printComment( stream, "Name: BC-" + boundaryName + " Type: Symmetry/Antisymmetry/Encastre" );
|
|
RifInpExportTools::printHeading( stream, "Boundary" );
|
|
RifInpExportTools::printLine( stream, boundarySetName + ", " + symmetryType );
|
|
};
|
|
|
|
std::map<RigGriddedPart3d::Boundary, std::string> symmetryTypes = {
|
|
{ RigGriddedPart3d::Boundary::Bottom, "ZSYMM" },
|
|
{ RigGriddedPart3d::Boundary::Back, "YSYMM" },
|
|
{ RigGriddedPart3d::Boundary::Front, "YSYMM" },
|
|
{ RigGriddedPart3d::Boundary::FarSide, "XSYMM" },
|
|
|
|
};
|
|
|
|
RifInpExportTools::printSectionComment( stream, "BOUNDARY CONDITIONS" );
|
|
for ( auto [boundary, boundaryName] : boundaries )
|
|
{
|
|
std::string boundarySetName = "Set-" + boundaryName;
|
|
std::string symmetryType = symmetryTypes[boundary];
|
|
printBoundaryCondition( stream, boundaryName, boundarySetName, symmetryType );
|
|
}
|
|
|
|
std::string partSymmetry = "XSYMM";
|
|
printBoundaryCondition( stream, "z1", "Part-1-1.part1", partSymmetry );
|
|
printBoundaryCondition( stream, "z2", "Part-2-1.part2", partSymmetry );
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printPredefinedFields( std::ostream& stream )
|
|
{
|
|
// PREDEFINED FIELDS
|
|
struct PredefinedField
|
|
{
|
|
std::string type;
|
|
std::string initialConditionType;
|
|
std::string shortName;
|
|
std::string partName;
|
|
double value;
|
|
};
|
|
|
|
std::vector<PredefinedField> fields = {
|
|
PredefinedField{ .type = "Void ratio", .initialConditionType = "RATIO", .shortName = "VR1", .partName = "Part-1-1.part1", .value = 0.3 },
|
|
PredefinedField{ .type = "Pore pressure",
|
|
.initialConditionType = "PORE PRESSURE",
|
|
.shortName = "PP1",
|
|
.partName = "Part-1-1.part1",
|
|
.value = 0.0 },
|
|
PredefinedField{ .type = "Pore pressure",
|
|
.initialConditionType = "PORE PRESSURE",
|
|
.shortName = "PP2",
|
|
.partName = "Part-2-1.part2",
|
|
.value = 0.0 },
|
|
PredefinedField{ .type = "Void ratio", .initialConditionType = "RATIO", .shortName = "VR2", .partName = "Part-2-1.part2", .value = 0.3 },
|
|
|
|
};
|
|
|
|
RifInpExportTools::printSectionComment( stream, "PREDEFINED FIELDS" );
|
|
|
|
for ( auto field : fields )
|
|
{
|
|
RifInpExportTools::printComment( stream, "Name: Predefined Field " + field.shortName + " Type: " + field.type );
|
|
RifInpExportTools::printHeading( stream, "Initial Conditions, TYPE=" + field.initialConditionType );
|
|
RifInpExportTools::printLine( stream, field.partName + ", " + std::to_string( field.value ) );
|
|
}
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string> RifFaultReactivationModelExporter::printSteps( std::ostream& stream )
|
|
{
|
|
int numSteps = 2;
|
|
|
|
for ( int i = 0; i < numSteps; i++ )
|
|
{
|
|
std::string stepName = "Step-" + std::to_string( i + 1 );
|
|
RifInpExportTools::printSectionComment( stream, "STEP: " + stepName );
|
|
|
|
RifInpExportTools::printHeading( stream, "Step, name=" + stepName + ", nlgeom=NO" );
|
|
|
|
std::string stepType = i == 0 ? "Geostatic, utol" : "Soils, utol=1.0";
|
|
RifInpExportTools::printHeading( stream, stepType );
|
|
RifInpExportTools::printNumbers( stream, { 1.0, 1.0, 1e-05, 1.0 } );
|
|
|
|
RifInpExportTools::printSectionComment( stream, "BOUNDARY CONDITIONS" );
|
|
RifInpExportTools::printHeading( stream, "Boundary" );
|
|
RifInpExportTools::printLine( stream, "Part-1-1.part1_PP_, 8, 8" );
|
|
RifInpExportTools::printHeading( stream, "Boundary" );
|
|
std::string extra = i != 0 ? ", 1e+07" : "";
|
|
RifInpExportTools::printLine( stream, "Part-2-1.part2_PP_, 8, 8" + extra );
|
|
|
|
RifInpExportTools::printSectionComment( stream, "OUTPUT REQUESTS" );
|
|
RifInpExportTools::printHeading( stream, "Restart, write, frequency=0" );
|
|
|
|
RifInpExportTools::printSectionComment( stream, "FIELD OUTPUT: F-Output-1" );
|
|
RifInpExportTools::printHeading( stream, "Output, field, variable=PRESELECT" );
|
|
|
|
RifInpExportTools::printSectionComment( stream, "HISTORY OUTPUT: H-Output-1" );
|
|
RifInpExportTools::printHeading( stream, "Output, history, variable=PRESELECT" );
|
|
|
|
RifInpExportTools::printHeading( stream, "End Step" );
|
|
}
|
|
|
|
return { true, "" };
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
std::pair<bool, std::string>
|
|
RifFaultReactivationModelExporter::printInteractions( std::ostream& stream,
|
|
const std::vector<std::pair<RigGriddedPart3d::BorderSurface, std::string>>& borders )
|
|
{
|
|
RifInpExportTools::printSectionComment( stream, "INTERACTIONS" );
|
|
for ( const auto& [border, borderName] : borders )
|
|
{
|
|
RifInpExportTools::printComment( stream, "Interaction: " + borderName );
|
|
|
|
std::string interactionName = "non-fault";
|
|
std::string extra;
|
|
if ( border == RigGriddedPart3d::BorderSurface::FaultSurface )
|
|
{
|
|
interactionName = "fault";
|
|
extra = ", adjust=0.0";
|
|
}
|
|
|
|
RifInpExportTools::printHeading( stream,
|
|
"Contact Pair, interaction=" + interactionName + ", small sliding, type=SURFACE TO SURFACE" + extra );
|
|
|
|
std::string part1Name = "Part-1-1";
|
|
std::string part2Name = "Part-2-1";
|
|
RifInpExportTools::printLine( stream, part1Name + "." + borderName + ", " + part2Name + "." + borderName );
|
|
}
|
|
|
|
return { true, "" };
|
|
}
|