40 using VectorType = dealii::LinearAlgebra::distributed::Vector<Number>;
49 parameters = parameters_in;
50 acoustic_solver = acoustic_solver_in;
51 fluid_solver = fluid_solver_in;
53 acoustic_solver_in->pde_operator->initialize_dof_vector_pressure(source_term_acoustic);
54 fluid_solver_in->pde_operator->initialize_vector_pressure(source_term_fluid);
57 if(parameters.fluid_to_acoustic_coupling_strategy ==
58 FluidToAcousticCouplingStrategy::ConservativeInterpolation)
60 non_nested_grid_transfer.reinit(fluid_solver_in->pde_operator->get_dof_handler_p(),
61 acoustic_solver_in->pde_operator->get_dof_handler_p(),
62 *fluid_solver_in->pde_operator->get_mapping(),
63 *acoustic_solver_in->pde_operator->get_mapping());
67 AssertThrow(
false, dealii::ExcMessage(
"FluidToAcousticCouplingStrategy not implemented."));
72 data.dof_index_pressure = fluid_solver_in->pde_operator->get_dof_index_pressure();
73 data.dof_index_velocity = fluid_solver_in->pde_operator->get_dof_index_velocity();
74 data.quad_index = fluid_solver_in->pde_operator->get_quad_index_pressure();
75 data.density = parameters_in.density;
76 data.consider_convection = parameters_in.source_term_with_convection;
77 data.blend_in = parameters.blend_in_source_term;
78 data.blend_in_function = field_functions_in->source_term_blend_in;
80 source_term_calculator.setup(fluid_solver_in->pde_operator->get_matrix_free(), data);
86 if(parameters.fluid_to_acoustic_coupling_strategy ==
87 FluidToAcousticCouplingStrategy::ConservativeInterpolation)
89 source_term_calculator.evaluate_integrate(source_term_fluid,
90 fluid_solver->time_integrator->get_velocity(),
91 fluid_solver->time_integrator->get_pressure(),
92 fluid_solver->get_pressure_time_derivative(),
93 fluid_solver->time_integrator->get_time());
95 non_nested_grid_transfer.restrict_and_add(source_term_acoustic, source_term_fluid);
99 AssertThrow(
false, dealii::ExcMessage(
"FluidToAcousticCouplingStrategy not implemented."));
102 acoustic_solver->pde_operator->set_aero_acoustic_source_term(source_term_acoustic);
109 std::shared_ptr<SolverAcoustic<dim, Number>> acoustic_solver;
110 std::shared_ptr<SolverFluid<dim, Number>> fluid_solver;
113 std::shared_ptr<FieldFunctions<dim>> field_functions;
116 dealii::MGTwoLevelTransferNonNested<dim, VectorType> non_nested_grid_transfer;
122 VectorType source_term_acoustic;
125 VectorType source_term_fluid;