Dear users
I need to run a thermochemical 3D convection in a box.
However, I receive this error message:
An error occurred in line <2360> of file </home/francesco-radica/Documenti/aspect/source/simulator/helper_functions.cc> in function
void aspect::Simulator<dim>::check_consistency_of_boundary_conditions() const [with int dim = 3]
The violated condition was:
is_element( pb.first.first, boundary_temperature_manager.get_fixed_temperature_boundary_indicators() ) == false && is_element( pb.first.second, boundary_temperature_manager.get_fixed_temperature_boundary_indicators() ) == false && is_element( pb.first.first, boundary_composition_manager.get_fixed_composition_boundary_indicators() ) == false && is_element( pb.first.second, boundary_composition_manager.get_fixed_composition_boundary_indicators() ) == false && is_element( pb.first.first, boundary_indicator_lists[0] ) == false && is_element( pb.first.second, boundary_indicator_lists[0] ) == false && is_element( pb.first.first, boundary_indicator_lists[1] ) == false && is_element( pb.first.second, boundary_indicator_lists[1] ) == false && is_element( pb.first.first, boundary_indicator_lists[3] ) == false && is_element( pb.first.second, boundary_indicator_lists[3] ) == false
Additional information:
Periodic boundaries must not have boundary conditions set.
Stacktrace:
-----------
#0 aspect: aspect::Simulator<3>::check_consistency_of_boundary_conditions() const
#1 aspect: aspect::Simulator<3>::Simulator(ompi_communicator_t*, dealii::ParameterHandler&)
#2 aspect: void run_simulator<3>(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, bool, bool, bool)
#3 aspect: main
--------------------------------------------------------
Aborting!
After some tests, the error is given by these lines:
subsection Geometry model
set Model name = box
subsection Box
set X extent = 6000.e3
set Y extent = 2890.e3
set Z extent = 6000.e3
set X periodic = true
set X repetitions = 1
set Z periodic = true
set Z repetitions = 1
end
end
subsection Boundary temperature model
set List of model names = box
set Fixed temperature boundary indicators = right, left
subsection Box
set Right temperature = 3700
set Top temperature = 300
end
end
The error disappears if I delete Z periodic and also if I delete the boundary temperature model.
The point is that I checked the manual, cookbooks, and tests, and there isn’t an example related, so I have no idea how to solve the problem. I think that, even if not explicit, the right and left boundaries are automatically set to 0, and so I get the error for Z.
How can I solve the problem?
Francesco.