Split set velocity and set acceleration

c-p-schmidt · c-p-schmidt · commit 25c3b366543f · 2025-05-02T11:12:27.000+02:00
into two separate methods and use const reference
diff --git a/src/elch/4C_elch_moving_boundary_algorithm.cpp b/src/elch/4C_elch_moving_boundary_algorithm.cpp
@@ -112,7 +112,8 @@ void ElCh::MovingBoundaryAlgorithm::time_loop()
   {
     // transfer convective velocity = fluid velocity - grid velocity
     scatra_field()->set_convective_velocity(*fluid_field()->convective_vel());
-    scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->convective_vel());
+    scatra_field()->set_velocity_field(*fluid_field()->convective_vel());
+    scatra_field()->set_acceleration_field(*fluid_field()->hist());
   }
 
   // transfer moving mesh data
@@ -266,7 +267,8 @@ void ElCh::MovingBoundaryAlgorithm::solve_scatra()
       {
         // transfer convective velocity = fluid velocity - grid velocity
         scatra_field()->set_convective_velocity(*fluid_field()->convective_vel());
-        scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->convective_vel());
+        scatra_field()->set_velocity_field(*fluid_field()->convective_vel());
+        scatra_field()->set_acceleration_field(*fluid_field()->hist());
       }
     }
     break;
diff --git a/src/fs3i/4C_fs3i_fps3i_partitioned.cpp b/src/fs3i/4C_fs3i_fps3i_partitioned.cpp
@@ -652,7 +652,7 @@ void FS3I::PartFPS3I::set_velocity_fields()
       for (unsigned i = 0; i < scatravec_.size(); ++i)
       {
         scatravec_[i]->scatra_field()->set_convective_velocity(*convel[i]);
-        scatravec_[i]->scatra_field()->set_velocity_field(nullptr, vel[i]);
+        scatravec_[i]->scatra_field()->set_velocity_field(*vel[i]);
       }
       break;
     }
diff --git a/src/fs3i/4C_fs3i_partitioned.cpp b/src/fs3i/4C_fs3i_partitioned.cpp
@@ -704,8 +704,7 @@ void FS3I::PartFS3I::set_velocity_fields() const
   {
     scatravec_[i]->scatra_field()->set_convective_velocity(
         *vol_mortar_master_to_slavei(i, convel[i]));
-    scatravec_[i]->scatra_field()->set_velocity_field(
-        nullptr, vol_mortar_master_to_slavei(i, vel[i]));
+    scatravec_[i]->scatra_field()->set_velocity_field(*vol_mortar_master_to_slavei(i, vel[i]));
   }
 }
 
diff --git a/src/loma/4C_loma_algorithm.cpp b/src/loma/4C_loma_algorithm.cpp
@@ -313,7 +313,7 @@ void LowMach::Algorithm::initial_calculations()
 {
   // set initial velocity field for evaluation of initial scalar time derivative in SCATRA
   scatra_field()->set_convective_velocity(*fluid_field()->velnp());
-  scatra_field()->set_velocity_field(nullptr, fluid_field()->velnp());
+  scatra_field()->set_velocity_field(*fluid_field()->velnp());
   if (scatra_field()->fine_scale_velocity_field_required())
   {
     scatra_field()->set_fine_scale_velocity(*fluid_field()->fs_vel());
@@ -543,8 +543,9 @@ void LowMach::Algorithm::set_fluid_values_in_scatra()
   {
     case Inpar::FLUID::timeint_afgenalpha:
     {
+      scatra_field()->set_acceleration_field(*fluid_field()->accam());
       scatra_field()->set_convective_velocity(*fluid_field()->velaf());
-      scatra_field()->set_velocity_field(fluid_field()->accam(), fluid_field()->velaf());
+      scatra_field()->set_velocity_field(*fluid_field()->velaf());
       if (scatra_field()->fine_scale_velocity_field_required() and
           fluid_field()->fs_vel() != nullptr)
       {
@@ -555,8 +556,9 @@ void LowMach::Algorithm::set_fluid_values_in_scatra()
     case Inpar::FLUID::timeint_one_step_theta:
     case Inpar::FLUID::timeint_bdf2:
     {
+      scatra_field()->set_acceleration_field(*fluid_field()->hist());
       scatra_field()->set_convective_velocity(*fluid_field()->velnp());
-      scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->velnp());
+      scatra_field()->set_velocity_field(*fluid_field()->velnp());
       if (scatra_field()->fine_scale_velocity_field_required() and
           fluid_field()->fs_vel() != nullptr)
       {
diff --git a/src/poroelast_scatra/4C_poroelast_scatra_base.cpp b/src/poroelast_scatra/4C_poroelast_scatra_base.cpp
@@ -180,7 +180,7 @@ void PoroElastScaTra::PoroScatraBase::set_velocity_fields()
   }
 
   scatra_->scatra_field()->set_convective_velocity(*convel);
-  scatra_->scatra_field()->set_velocity_field(nullptr, velnp);
+  scatra_->scatra_field()->set_velocity_field(*velnp);
 }
 
 /*----------------------------------------------------------------------*
diff --git a/src/scatra/4C_scatra_algorithm.cpp b/src/scatra/4C_scatra_algorithm.cpp
@@ -234,8 +234,9 @@ void ScaTra::ScaTraAlgorithm::prepare_time_step_convection()
   //(fluid initial field was set inside the constructor of fluid base class)
   if (step() == 1)
   {
+    scatra_field()->set_acceleration_field(*fluid_field()->hist());
     scatra_field()->set_convective_velocity(*fluid_field()->velnp());
-    scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->velnp());
+    scatra_field()->set_velocity_field(*fluid_field()->velnp());
   }
 
   // prepare time step (+ initialize one-step-theta scheme correctly with
@@ -298,9 +299,9 @@ void ScaTra::ScaTraAlgorithm::set_velocity_field()
     case Inpar::FLUID::timeint_npgenalpha:
     case Inpar::FLUID::timeint_afgenalpha:
     {
+      scatra_field()->set_acceleration_field(*fluid_to_scatra(fluid_field()->accam()));
       scatra_field()->set_convective_velocity(*fluid_to_scatra(fluid_field()->velaf()));
-      scatra_field()->set_velocity_field(
-          fluid_to_scatra(fluid_field()->accam()), fluid_to_scatra(fluid_field()->velaf()));
+      scatra_field()->set_velocity_field(*fluid_to_scatra(fluid_field()->velaf()));
       if (scatra_field()->fine_scale_velocity_field_required() and
           fluid_field()->fs_vel() != nullptr)
       {
@@ -312,9 +313,9 @@ void ScaTra::ScaTraAlgorithm::set_velocity_field()
     case Inpar::FLUID::timeint_bdf2:
     case Inpar::FLUID::timeint_stationary:
     {
+      scatra_field()->set_acceleration_field(*fluid_to_scatra(fluid_field()->hist()));
       scatra_field()->set_convective_velocity(*fluid_to_scatra(fluid_field()->velnp()));
-      scatra_field()->set_velocity_field(
-          fluid_to_scatra(fluid_field()->hist()), fluid_to_scatra(fluid_field()->velnp()));
+      scatra_field()->set_velocity_field(*fluid_to_scatra(fluid_field()->velnp()));
       if (scatra_field()->fine_scale_velocity_field_required() and
           fluid_field()->fs_vel() != nullptr)
       {
diff --git a/src/scatra/4C_scatra_timint_implicit.cpp b/src/scatra/4C_scatra_timint_implicit.cpp
@@ -1461,6 +1461,20 @@ void ScaTra::ScaTraTimIntImpl::set_mean_concentration(
   mean_conc_ = MeanConc;
 }
 
+/*----------------------------------------------------------------------*
+ *----------------------------------------------------------------------*/
+void ScaTra::ScaTraTimIntImpl::set_acceleration_field(
+    const Core::LinAlg::Vector<double>& acceleration) const
+{
+  // time measurement
+  TEUCHOS_FUNC_TIME_MONITOR("SCATRA: set acceleration");
+
+  FOUR_C_ASSERT(nds_vel() < discret_->num_dof_sets(), "Too few dof sets on scatra discretization!");
+
+  // provide scatra discretization with acceleration field if required
+  discret_->set_state(nds_vel(), "acceleration field", acceleration);
+}
+
 /*----------------------------------------------------------------------*
  *----------------------------------------------------------------------*/
 void ScaTra::ScaTraTimIntImpl::set_convective_velocity(
@@ -1515,23 +1529,18 @@ bool ScaTra::ScaTraTimIntImpl::fine_scale_velocity_field_required() const
 
 /*----------------------------------------------------------------------*
  *----------------------------------------------------------------------*/
-void ScaTra::ScaTraTimIntImpl::set_velocity_field(
-    std::shared_ptr<const Core::LinAlg::Vector<double>> acc,
-    std::shared_ptr<const Core::LinAlg::Vector<double>> vel)
+void ScaTra::ScaTraTimIntImpl::set_velocity_field(const Core::LinAlg::Vector<double>& velocity)
 {
   // time measurement
-  TEUCHOS_FUNC_TIME_MONITOR("SCATRA: set velocity fields");
+  TEUCHOS_FUNC_TIME_MONITOR("SCATRA: set velocity field");
 
   // checks
   FOUR_C_ASSERT(velocity_field_type_ == Inpar::ScaTra::velocity_Navier_Stokes,
       "Wrong set_velocity_field() called for velocity field type {}!", velocity_field_type_);
   FOUR_C_ASSERT(nds_vel() < discret_->num_dof_sets(), "Too few dof sets on scatra discretization!");
 
   // provide scatra discretization with velocity
-  if (vel != nullptr) discret_->set_state(nds_vel(), "velocity field", *vel);
-
-  // provide scatra discretization with acceleration field if required
-  if (acc != nullptr) discret_->set_state(nds_vel(), "acceleration field", *acc);
+  discret_->set_state(nds_vel(), "velocity field", velocity);
 }
 
 /*----------------------------------------------------------------------*
diff --git a/src/scatra/4C_scatra_timint_implicit.hpp b/src/scatra/4C_scatra_timint_implicit.hpp
@@ -244,6 +244,9 @@ namespace ScaTra
      */
     void set_external_force() const;
 
+    //! set the @acceleration vector to the scalar transport discretization
+    void set_acceleration_field(const Core::LinAlg::Vector<double>& acceleration) const;
+
     //! set the @convective_velocity vector to the scalar transport discretization
     void set_convective_velocity(const Core::LinAlg::Vector<double>& convective_velocity) const;
 
@@ -253,15 +256,8 @@ namespace ScaTra
     //! return whether setting of the fine scale velocity is required
     [[nodiscard]] bool fine_scale_velocity_field_required() const;
 
-    /*!
-     * @brief set velocity field (+ pressure and acceleration field as well as fine-scale velocity
-     * field, if required)
-     *
-     * @param acc     acceleration vector
-     * @param vel     velocity vector
-     */
-    void set_velocity_field(std::shared_ptr<const Core::LinAlg::Vector<double>> acc,
-        std::shared_ptr<const Core::LinAlg::Vector<double>> vel);
+    //! set the @velocity vector to the scalar transport discretization
+    void set_velocity_field(const Core::LinAlg::Vector<double>& velocity);
 
     void set_wall_shear_stresses(std::shared_ptr<const Core::LinAlg::Vector<double>> wss);
 
diff --git a/src/ssi/4C_ssi_coupling.cpp b/src/ssi/4C_ssi_coupling.cpp
@@ -142,7 +142,7 @@ void SSI::SSICouplingMatchingVolume::set_velocity_fields(
     std::shared_ptr<const Core::LinAlg::Vector<double>> vel)
 {
   scatra->scatra_field()->set_convective_velocity(*convvel);
-  scatra->scatra_field()->set_velocity_field(nullptr, vel);
+  scatra->scatra_field()->set_velocity_field(*vel);
 }
 
 /*----------------------------------------------------------------------*/
@@ -290,7 +290,7 @@ void SSI::SSICouplingNonMatchingBoundary::set_velocity_fields(
   scatra->scatra_field()->set_convective_velocity(
       *adaptermeshtying_->master_to_slave(*extractor_->extract_cond_vector(*convvel)));
   scatra->scatra_field()->set_velocity_field(
-      nullptr, adaptermeshtying_->master_to_slave(*extractor_->extract_cond_vector(*vel)));
+      *adaptermeshtying_->master_to_slave(*extractor_->extract_cond_vector(*vel)));
 }
 
 /*----------------------------------------------------------------------*/
@@ -410,7 +410,7 @@ void SSI::SSICouplingNonMatchingVolume::set_velocity_fields(
   scatra->scatra_field()->set_convective_velocity(
       *volcoupl_structurescatra_->apply_vector_mapping21(*convvel));
   scatra->scatra_field()->set_velocity_field(
-      nullptr, volcoupl_structurescatra_->apply_vector_mapping21(*vel));
+      *volcoupl_structurescatra_->apply_vector_mapping21(*vel));
 }
 
 /*----------------------------------------------------------------------*/
@@ -606,7 +606,7 @@ void SSI::SSICouplingMatchingVolumeAndBoundary::set_velocity_fields(
     std::shared_ptr<const Core::LinAlg::Vector<double>> vel)
 {
   scatra->scatra_field()->set_convective_velocity(*convvel);
-  scatra->scatra_field()->set_velocity_field(nullptr, vel);
+  scatra->scatra_field()->set_velocity_field(*vel);
 }
 
 /*----------------------------------------------------------------------*/
diff --git a/src/ssti/4C_ssti_algorithm.cpp b/src/ssti/4C_ssti_algorithm.cpp
@@ -290,9 +290,9 @@ void SSTI::SSTIAlgorithm::distribute_structure_solution() const
   const auto convective_velocity = Core::LinAlg::create_vector(*structure_->dof_row_map());
 
   scatra_field()->set_convective_velocity(*convective_velocity);
-  scatra_field()->set_velocity_field(nullptr, structure_->velnp());
+  scatra_field()->set_velocity_field(*structure_->velnp());
   thermo_field()->set_convective_velocity(*convective_velocity);
-  thermo_field()->set_velocity_field(nullptr, structure_->velnp());
+  thermo_field()->set_velocity_field(*structure_->velnp());
 }
 
 /*----------------------------------------------------------------------*/

Original file line number	Diff line number	Diff line change
`@@ -112,7 +112,8 @@ void ElCh::MovingBoundaryAlgorithm::time_loop()`
`112`	`112`	`{`
`113`	`113`	`// transfer convective velocity = fluid velocity - grid velocity`
`114`	`114`	`scatra_field()->set_convective_velocity(*fluid_field()->convective_vel());`
`115`		`- scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->convective_vel());`
	`115`	`+ scatra_field()->set_velocity_field(*fluid_field()->convective_vel());`
	`116`	`+ scatra_field()->set_acceleration_field(*fluid_field()->hist());`
`116`	`117`	`}`
`117`	`118`
`118`	`119`	`// transfer moving mesh data`
`@@ -266,7 +267,8 @@ void ElCh::MovingBoundaryAlgorithm::solve_scatra()`
`266`	`267`	`{`
`267`	`268`	`// transfer convective velocity = fluid velocity - grid velocity`
`268`	`269`	`scatra_field()->set_convective_velocity(*fluid_field()->convective_vel());`
`269`		`- scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->convective_vel());`
	`270`	`+ scatra_field()->set_velocity_field(*fluid_field()->convective_vel());`
	`271`	`+ scatra_field()->set_acceleration_field(*fluid_field()->hist());`
`270`	`272`	`}`
`271`	`273`	`}`
`272`	`274`	`break;`
Original file line number	Diff line number	Diff line change
`@@ -652,7 +652,7 @@ void FS3I::PartFPS3I::set_velocity_fields()`
`652`	`652`	`for (unsigned i = 0; i < scatravec_.size(); ++i)`
`653`	`653`	`{`
`654`	`654`	`scatravec_[i]->scatra_field()->set_convective_velocity(*convel[i]);`
`655`		`- scatravec_[i]->scatra_field()->set_velocity_field(nullptr, vel[i]);`
	`655`	`+ scatravec_[i]->scatra_field()->set_velocity_field(*vel[i]);`
`656`	`656`	`}`
`657`	`657`	`break;`
`658`	`658`	`}`
Original file line number	Diff line number	Diff line change
`@@ -704,8 +704,7 @@ void FS3I::PartFS3I::set_velocity_fields() const`
`704`	`704`	`{`
`705`	`705`	`scatravec_[i]->scatra_field()->set_convective_velocity(`
`706`	`706`	`*vol_mortar_master_to_slavei(i, convel[i]));`
`707`		`- scatravec_[i]->scatra_field()->set_velocity_field(`
`708`		`- nullptr, vol_mortar_master_to_slavei(i, vel[i]));`
	`707`	`+ scatravec_[i]->scatra_field()->set_velocity_field(*vol_mortar_master_to_slavei(i, vel[i]));`
`709`	`708`	`}`
`710`	`709`	`}`
`711`	`710`
Original file line number	Diff line number	Diff line change
`@@ -313,7 +313,7 @@ void LowMach::Algorithm::initial_calculations()`
`313`	`313`	`{`
`314`	`314`	`// set initial velocity field for evaluation of initial scalar time derivative in SCATRA`
`315`	`315`	`scatra_field()->set_convective_velocity(*fluid_field()->velnp());`
`316`		`- scatra_field()->set_velocity_field(nullptr, fluid_field()->velnp());`
	`316`	`+ scatra_field()->set_velocity_field(*fluid_field()->velnp());`
`317`	`317`	`if (scatra_field()->fine_scale_velocity_field_required())`
`318`	`318`	`{`
`319`	`319`	`scatra_field()->set_fine_scale_velocity(*fluid_field()->fs_vel());`
`@@ -543,8 +543,9 @@ void LowMach::Algorithm::set_fluid_values_in_scatra()`
`543`	`543`	`{`
`544`	`544`	`case Inpar::FLUID::timeint_afgenalpha:`
`545`	`545`	`{`
	`546`	`+ scatra_field()->set_acceleration_field(*fluid_field()->accam());`
`546`	`547`	`scatra_field()->set_convective_velocity(*fluid_field()->velaf());`
`547`		`- scatra_field()->set_velocity_field(fluid_field()->accam(), fluid_field()->velaf());`
	`548`	`+ scatra_field()->set_velocity_field(*fluid_field()->velaf());`
`548`	`549`	`if (scatra_field()->fine_scale_velocity_field_required() and`
`549`	`550`	`fluid_field()->fs_vel() != nullptr)`
`550`	`551`	`{`
`@@ -555,8 +556,9 @@ void LowMach::Algorithm::set_fluid_values_in_scatra()`
`555`	`556`	`case Inpar::FLUID::timeint_one_step_theta:`
`556`	`557`	`case Inpar::FLUID::timeint_bdf2:`
`557`	`558`	`{`
	`559`	`+ scatra_field()->set_acceleration_field(*fluid_field()->hist());`
`558`	`560`	`scatra_field()->set_convective_velocity(*fluid_field()->velnp());`
`559`		`- scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->velnp());`
	`561`	`+ scatra_field()->set_velocity_field(*fluid_field()->velnp());`
`560`	`562`	`if (scatra_field()->fine_scale_velocity_field_required() and`
`561`	`563`	`fluid_field()->fs_vel() != nullptr)`
`562`	`564`	`{`
Original file line number	Diff line number	Diff line change
`@@ -180,7 +180,7 @@ void PoroElastScaTra::PoroScatraBase::set_velocity_fields()`
`180`	`180`	`}`
`181`	`181`
`182`	`182`	`scatra_->scatra_field()->set_convective_velocity(*convel);`
`183`		`- scatra_->scatra_field()->set_velocity_field(nullptr, velnp);`
	`183`	`+ scatra_->scatra_field()->set_velocity_field(*velnp);`
`184`	`184`	`}`
`185`	`185`
`186`	`186`	`/----------------------------------------------------------------------`
Original file line number	Diff line number	Diff line change
`@@ -234,8 +234,9 @@ void ScaTra::ScaTraAlgorithm::prepare_time_step_convection()`
`234`	`234`	`//(fluid initial field was set inside the constructor of fluid base class)`
`235`	`235`	`if (step() == 1)`
`236`	`236`	`{`
	`237`	`+ scatra_field()->set_acceleration_field(*fluid_field()->hist());`
`237`	`238`	`scatra_field()->set_convective_velocity(*fluid_field()->velnp());`
`238`		`- scatra_field()->set_velocity_field(fluid_field()->hist(), fluid_field()->velnp());`
	`239`	`+ scatra_field()->set_velocity_field(*fluid_field()->velnp());`
`239`	`240`	`}`
`240`	`241`
`241`	`242`	`// prepare time step (+ initialize one-step-theta scheme correctly with`
`@@ -298,9 +299,9 @@ void ScaTra::ScaTraAlgorithm::set_velocity_field()`
`298`	`299`	`case Inpar::FLUID::timeint_npgenalpha:`
`299`	`300`	`case Inpar::FLUID::timeint_afgenalpha:`
`300`	`301`	`{`
	`302`	`+ scatra_field()->set_acceleration_field(*fluid_to_scatra(fluid_field()->accam()));`
`301`	`303`	`scatra_field()->set_convective_velocity(*fluid_to_scatra(fluid_field()->velaf()));`
`302`		`- scatra_field()->set_velocity_field(`
`303`		`- fluid_to_scatra(fluid_field()->accam()), fluid_to_scatra(fluid_field()->velaf()));`
	`304`	`+ scatra_field()->set_velocity_field(*fluid_to_scatra(fluid_field()->velaf()));`
`304`	`305`	`if (scatra_field()->fine_scale_velocity_field_required() and`
`305`	`306`	`fluid_field()->fs_vel() != nullptr)`
`306`	`307`	`{`
`@@ -312,9 +313,9 @@ void ScaTra::ScaTraAlgorithm::set_velocity_field()`
`312`	`313`	`case Inpar::FLUID::timeint_bdf2:`
`313`	`314`	`case Inpar::FLUID::timeint_stationary:`
`314`	`315`	`{`
	`316`	`+ scatra_field()->set_acceleration_field(*fluid_to_scatra(fluid_field()->hist()));`
`315`	`317`	`scatra_field()->set_convective_velocity(*fluid_to_scatra(fluid_field()->velnp()));`
`316`		`- scatra_field()->set_velocity_field(`
`317`		`- fluid_to_scatra(fluid_field()->hist()), fluid_to_scatra(fluid_field()->velnp()));`
	`318`	`+ scatra_field()->set_velocity_field(*fluid_to_scatra(fluid_field()->velnp()));`
`318`	`319`	`if (scatra_field()->fine_scale_velocity_field_required() and`
`319`	`320`	`fluid_field()->fs_vel() != nullptr)`
`320`	`321`	`{`
Original file line number	Diff line number	Diff line change
`@@ -142,7 +142,7 @@ void SSI::SSICouplingMatchingVolume::set_velocity_fields(`
`142`	`142`	`std::shared_ptr<const Core::LinAlg::Vector<double>> vel)`
`143`	`143`	`{`
`144`	`144`	`scatra->scatra_field()->set_convective_velocity(*convvel);`
`145`		`- scatra->scatra_field()->set_velocity_field(nullptr, vel);`
	`145`	`+ scatra->scatra_field()->set_velocity_field(*vel);`
`146`	`146`	`}`
`147`	`147`
`148`	`148`	`/----------------------------------------------------------------------/`
`@@ -290,7 +290,7 @@ void SSI::SSICouplingNonMatchingBoundary::set_velocity_fields(`
`290`	`290`	`scatra->scatra_field()->set_convective_velocity(`
`291`	`291`	`adaptermeshtying_->master_to_slave(extractor_->extract_cond_vector(*convvel)));`
`292`	`292`	`scatra->scatra_field()->set_velocity_field(`
`293`		`- nullptr, adaptermeshtying_->master_to_slave(extractor_->extract_cond_vector(vel)));`
	`293`	`+ adaptermeshtying_->master_to_slave(extractor_->extract_cond_vector(*vel)));`
`294`	`294`	`}`
`295`	`295`
`296`	`296`	`/----------------------------------------------------------------------/`
`@@ -410,7 +410,7 @@ void SSI::SSICouplingNonMatchingVolume::set_velocity_fields(`
`410`	`410`	`scatra->scatra_field()->set_convective_velocity(`
`411`	`411`	`volcoupl_structurescatra_->apply_vector_mapping21(convvel));`
`412`	`412`	`scatra->scatra_field()->set_velocity_field(`
`413`		`- nullptr, volcoupl_structurescatra_->apply_vector_mapping21(*vel));`
	`413`	`+ volcoupl_structurescatra_->apply_vector_mapping21(vel));`
`414`	`414`	`}`
`415`	`415`
`416`	`416`	`/----------------------------------------------------------------------/`
`@@ -606,7 +606,7 @@ void SSI::SSICouplingMatchingVolumeAndBoundary::set_velocity_fields(`
`606`	`606`	`std::shared_ptr<const Core::LinAlg::Vector<double>> vel)`
`607`	`607`	`{`
`608`	`608`	`scatra->scatra_field()->set_convective_velocity(*convvel);`
`609`		`- scatra->scatra_field()->set_velocity_field(nullptr, vel);`
	`609`	`+ scatra->scatra_field()->set_velocity_field(*vel);`
`610`	`610`	`}`
`611`	`611`
`612`	`612`	`/----------------------------------------------------------------------/`