Mapping Between Fluid Flow and Acoustics Mesh
When the aeroacoustic and the Fluid Flow (CFD) model are solved on different computational mesh, careful mapping of the CFD solution from the CFD mesh onto the acoustics mesh should be done. This step is crucial in order not to introduce nonphysical numerical noise into the acoustic solution (Ref. 20). If the solution is not properly mapped terms containing gradients of the background mean flow variables (especially the reactive terms) become very noisy and are typically the source of the numerical error. This mapping is handled by the multiphysics coupling Background Fluid Flow Coupling and the dedicated Mapping study.
If the acoustic problem is solved on the same mesh as the CFD then the mapping problem is less of an issue, especially if both physics use the same shape function order (they are both linear per default). Using the Background Fluid Flow Coupling is also necessary, if perfectly matched layers (PMLs) are present in the model, in order to get the flow field inside the PML regions. In any case, the coupling can still be useful as the CFD solution will be smoothed by the mapping (the degree of smoothing can be controlled).
The multiphysics coupling of turbulent flow and acoustics is described in the application library model Helmholtz Resonator with Flow: Interaction of Flow and Acoustics. The Application Library path Acoustics_Module/Aeroacoustics_and_Noise/helmholtz_resnoator_with_flow
Manual Mapping
The multiphysics coupling Background Fluid Flow Coupling can couple all Single-Phase Flow and Nonisothermal Flow interfaces including the high Mach Number Flow interfaces. If such a coupling and mapping is necessary, manual mapping can be set up using the procedure described here.
The manual mapping is defined by a study that solves an additional set of equations which maps and smooths the background flow variables onto the acoustics mesh. The mapping equations can be set up using the Weak Form PDE interface from the Mathematics branch. In the most general case the mean background flow pressure p0, velocity field u0, density ρ0, temperature T0, and turbulent viscosity μT variables should be mapped onto corresponding variables on the acoustics mesh, for example, p0,aco, u0,aco, ρ0,aco, T0.aco, and μT,aco. These new variables should then be used as the model inputs in the aeroacoustics model.
A simple but efficient mapping and smoothing is achieved by solving
where the term on the right-hand side adds smoothing using isotropic diffusion. The amount of diffusion is controlled by the parameter δ (a constant that can be tuned with a typical value of 0.01) and the mesh size squared h2. The term corresponds to so-called source term stabilization as known from CFD.
In the Weak Form PDE interface define as many dependent variables as necessary and give them the same shape order as the order used for the acoustics (typically all linear). Assume that the source CFD variable for the x-velocity component is u and the new destination variable is U0, then the above mapping is achieved with the following Weak Expression, entered in the user interface:
(U0-withsol('sol1',u,setval(Ma,0.1)))*test(U0)+delta*h^2*(U0x*test(U0x)+U0y*test(U0y)+U0z*test(U0z)
For the pressure map the total pressure spf.pA to P0 using:
(P0-withsol('sol1',spf.pA,setval(Ma,0.1)))*test(P0)+delta*h^2*(P0x*test(P0x)+P0y*test(P0y)+P0z*test(P0z))
Notice the use of the withsol() operator, which is an extrusion-like coupling operator that can refer directly to a solution object and parameter value. In this example it is used to fetch data from the CFD mesh, in the solution generated by solver sol1 and for the parameter value Ma = 0.1. The setval() statement is optional. It is important to use such an operator such that the CFD solution is mapped and interpolated correctly to the integration (Gauss) points on the acoustics mesh.
In the COMSOL Multiphysics Reference Manual:
Operators, Functions, and Constants and Built-In Operators for details about the withsol() operator
Numerical Stabilization (general introduction)
In some cases it can be necessary to extend the mapping equation with boundary conditions for the mapped variables. This should be done at symmetry planes for the mapped velocity field. This is also the cases when flow details near walls are important, for example, resolving the acoustic boundary layer. In this case add a no-slip condition on u0,aco. In general constraining the value of the mapped variables to the value in the CFD model values at boundaries can be necessary.