Fluid–Structure Interaction Simulation of a Coriolis Mass Flowmeter Using a Lattice Boltzmann Method

: In this paper, we use a fluid–structure interaction (FSI) approach to simulate a Coriolis mass flowmeter (CMF). The fluid dynamics is calculated by the open-source framework OpenLB, based on the lattice Boltzmann method (LBM). For the structural dynamics we employ the open-source software Elmer, an implementation of the finite element method (FEM). A staggered coupling approach between the two software packages is presented. The finite element mesh is created by the mesh generator Gmsh to ensure a complete open source workflow. The Eigenmodes of the CMF, which are calculated by modal analysis, are compared with measurement data. Using the estimated excitation frequency, a fully coupled, partitioned, FSI simulation is applied to simulate the phase shift of the investigated CMF design. The calculated phase shift values are in good agreement to the measurement data and verify the suitability of the model to numerically describe the working principle of a CMF.


Introduction
The exact measurement of mass flow of fluids is important in many branches of technology, for example chemical, oil and gas industry.It is needed to control processes and ensure safety, filling batches, inventory and others.The Coriolis mass flow meter (CMF) is an accurate instrument, which is becoming increasingly important in various applications [1].It consists of one or multiple measuring tubes that are stimulated to vibrate by an electromagnetic pulse generator.The fluid to be investigated is directed through the tubes.Due to inertia, the Coriolis force causes a phase shift of the vibration, which is detected by sensors on both ends of the system.As the mass flow of the conveyed fluid is proportional to the Coriolis force, it can be determined directly.
CMFs have been widely described by analytical and structural models [2,3,4,5,6,7].These models have helped to understand the fundamental principle of CMF devices.Nevertheless, the influence of the fluid was greatly simplified and the practical operation could not be described completely.Therefore, fluid-structure interaction (FSI) models were developed to realize the operating principle, which means that the fluid motion is affected by the measuring pipe oscillation and the pipe motion in turn by the hydrodynamic forces.In recent years, iterative two-way FSI models, which consist of a separated computational fluid dynamics (CFD) solver and a computational structural mechanics (CSM) solver, were applied to simulate CMF.Bobovnik et al. [8] used two different solvers to simulate a straight tube.Commercially available finite volume code for three dimensional turbulent fluid flow and finite element code for a shell structure were coupled.Five different tube lengths were studied simulating free tube vibration.The results for phase shift and frequency were similar to an analytical Flügge shell and potential flow model.In 2008, Mole et al. [9] extended the three dimensional numerical model of Bobovnik et al. [8] to deal with forced vibration.The study comprises the investigation of meter sensitivity at different Reynolds numbers.A maximum decrease of 0.4 % was observed for the lowest Reynolds number.This deviation is known as the low Reynolds number effect.The same numerical model was used by Bobovnik et al. [10] to study the influence of the design parameters on the installation effects of a CMF.Installation effects are measured as change of meter sensitivities from fully developed to disturbed fluid flow.Considering a single straight tube the errors vary with sensor positions and decrease with increasing tube length.In contrast, Kumar [11] claimed that a CMF is not sensitive to flow profiles.The FSI model of ANSYS-CFX was used to consider a straight single tube.The results were quite similar for the shorter tube lengths in comparison to previous studies [8].In contrast, the longer tubes showed a higher deviation, which was attributed to the different resolution.By changing the viscosity, the Reynolds number was varied and the deviation in meter sensitivity could be captured.It was found that at low Reynolds numbers the oscillating viscous fluid forces become relatively strong and interact with the oscillating Coriolis force, which changes the measurement results.To further investigate the effect of the Reynolds numbers, Kumar and Anklin [12] investigated a curved double tube CMF with an FSI simulation.The meter deviation at low Reynolds numbers were found in good agreement to measurement data.The low Reynolds number effect was indicated as correctable, if the viscosity of the examined fluid is known.Also Rongmo and Jian [13] used the ANSYS-CFX FSI module to study the low Reynolds number effect in a U-tube CMF.They assumed that arising deviations may be due to those different damping factors.Damping influences the natural frequency of the tube and was expected to change the meter sensitivity.
The aforementioned studies employ traditional discretization methods like the finite volume method (FVM) for the fluid solver.Meanwhile, alternative approaches, such as the lattice Boltzmann method (LBM), have received increasing attention.Its highly efficient parallel algorithm [14,15] and the applicability to a wide range of flow phenomena, e.g.flows in complex geometry [16,17] or turbulent flows [18,19] offer a high potential.
One of the first approaches that couple LBM to a structural solver can be found by Scholz et al. [20].They propose an anisotropic p-adaptive method for elastodynamic problems and show a higher convergence rate in comparison to a uniform p-version.Especially, the load transfer between the fluid and structural mesh were discussed.Geller et al. [21] used a partitioned approach to address the famous two-dimensional FSI benchmark case proposed by Turek and Hron [22].The proposed coupling approach by Geller et al. [21] leads to consistent quantitative result.A further study to validate an LBM solver coupled to a p-FEM solver with the Turek and Hron [22] benchmark was published by Kollmannsberger et al. [23].The staggered coupling was shown to be sufficient for simulating the reference case due to the weaker impact of the additional mass effect at small time steps.In contrast, Li et al. [24] claimed that the added mass effect has a major influence on accuracy and stability.They shown that the use of a non-staggered coupling approach based on subiterations reduce the effect of artificially added mass.Based on the previously mentioned studies [23,21,20], Geller et al. [25] extended the developed FSI approach to address three dimensional benchmark problems.
In contrast, this paper aims to demonstrate the feasibility of a complete open source FSI workflow to simulate a CMF.Therefore, OpenLB [26,27], an open source implementation of LBM, is coupled to the open source FEM framework Elmer [28].The implemented coupling procedure uses a staggered approach.A modal analysis of the CMF geometry is executed to extract the excitation frequency.The obtained excitation frequency is applied in a frequency response test to evaluate the transient structural setting.The partitioned FSI approach is used to calculate the phase shift.Both the Eigenfrequencies and the phase shift values are compared to measurement data.The evaluation and validation of a complex engineering problem with a partitioned FSI approach using LBM is a novelty.As a further highlight the new FSI workflow is built on open source frameworks to ensure additional adaptions in the coupling interface.
The paper is structured as follows, Section 2 introduces the applied FSI approach covering the fluid and structural models.In Section 3 the CMF test case is depicted in detail.The related modal analysis and the subsequent phase shift calculation results, using the FSI approach, are presented and compared to the measurement data in Section 4. Finally, Section 5 summarizes the findings and draws a conclusion.

Methodology
Firstly, the governing equations for the fluid domain presented.Afterwards the LBM and the moving boundary approach is introduced.Next, the Navier-Cauchy equation and the applied solution routine for the structural domain are shown.Finally, the FSI approach to coupling the structural an the fluid domain is presented, including the coupling conditions and implementation details.

Navier-Stokes Equations
The incompressible Navier-Stokes equations consist of the continuity equation and the momentum equation, which reads where Greek indices obey the Einstein notation, u f α is the fluid velocity, p is the pressure field, T f αβ is the fluid stress tensor and ρ f is the fluid density.Assuming a Newtonian fluid, the fluid stress tensor is given by where ν f is the kinematic viscosity.

(4)
The choice of the collision operator is justified by the higher computation performance and the lower memory demand in the used LBM implementation.The violation of the rotational invariance [29] in comparison to D3Q27 can be neglected in the laminar flow regime.
The lattice Boltzmann equation without external forces is given by where f i is the particle distribution function at discrete lattice position x LB and time step t LB .The collision operator Ω i is implemented by a single-relaxation time model proposed by Bhatnagar, Gross and Krook [30].It can be defined as where τ is the relaxation time towards the discrete particle distribution function at equilibrium state f eq i , ρ LB is the lattice density and u LB the velocity field.Hence, the collision operator conserve mass and momentum.The particle distribution function equilibrium f eq i is described by a low Mach number truncated Maxwell-Boltzmann distribution where ω i are the lattice weights obtained by the Gauss-Hermite quadrature [31,32], c s = 1/ √ 3 is the speed of sound of the lattice and δ αβ is the Kronecker delta.
The discrete moments of the particle distribution functions f i result in macroscopic flow quantities.The density ρ LB , the momentum ρ LB u LB and the momentum flux Π are respectively obtained by the zeroth, first and second moments, which are given by The relaxation time τ is coupled with the lattice kinematic viscosity ν LB through Taking a simplified isothermal equation of state into account, the lattice pressure is related to the lattice density by The lattice Mach number Ma LB is written as where u LB char is the characteristic lattice velocity.In the incompressible limit (Ma LB → 0), the incompressible Navier-Stokes equations (see Eqs. ( 1) and ( 2)) are recovered.
Finally, the lattice Boltzmann algorithm is parted into two steps: local collision step and subsequent streaming step.The local collision step is represented by the right-hand side of Eq. ( 5) and the streaming step is associated with the left-hand side of Eq. ( 5).

Moving Boundary Methods
A fluid-solid interface is required for the simulation of FSI.The lattice Boltzmann method typically uses three groups of approaches to describe this type of interface, namely the partially saturated methods [33,34], the immersed boundary methods [35,36,37] and the moving boundary methods [38,39,40].In the present study, the latter type of interface description is applied.
A moving fluid-solid interface inside the fluid domain can be described by the position of its boundary, which changes over time.The current boundary position indicates fluid and solid nodes.If a former solid becomes a fluid node, a refill algorithm is applied to reconstruct the unknown particle distribution functions.Hence, moving boundary methods are conceptually parted into a velocity boundary formulation and a refill algorithm.For a better comprehension, the introduced index conventions are displayed in Fig. 1.The index b is related to the boundary node positioned inside the solid domain.In direction c i the wall is intersected at x LB w .The locations x LB f , x ff and x fff denote the corresponding fluid nodes in this direction.The distance between x w and x f is given by the normalized distance q, which is calculated by Position x LB new indicates nodes, where a refill algorithm is necessary.The nodes in discrete normal direction c n are identified by x LB n , x LB nn and x LB nnn .This subscript convention is also used for velocity u LB and density ρ LB .
The present work uses the curved boundary condition proposed by Bouzidi et al. [41], which represents an extension of a half-way bounce back boundary scheme.Thereby a linear interpolation is utilized to take the distance to the boundary into account for increased accuracy.The unknown populations fī(x LB f , t LB + 1) after the streaming step are calculated by where index ī denotes a quantity in the opposite direction of the one with index i.The half-way bounce back condition is recovered for q = 1/2.
For the refill algorithm, a second order extrapolation scheme can be found in [40] Hereby the particle distributions f i (x LB new , t LB ) are extrapolated by the particle distribution functions in discrete normal direction c n .
A further necessary step for FSI is the calculation of the hydrodynamic forces that act on the interface.Therefore, a Galilean invariant momentum exchange approach [42] is used.The boundary force that acts on a solid node x LB b can be calculated by where L is the set of fluid-solid links.This formulation is suitable for the precise description of the boundary force of moving fluid-solid interfaces and avoids the disadvantages of a conventional momentum exchange calculation [43,42].

Navier-Cauchy Equation
The present work uses the Navier-Cauchy equation to describe the structural motion.Therefore, the structral moition is assumed to be linear elastic.The equation of motion for a linear elastic structure in differential form reads where Φ β is the structural displacement, ρ s is the solid density, F s β is the body-force acting on the structure.Thereby, the Cauchy stress tensor T s αβ can be written as where µ s and λ s represents the first and second Lamé constants.Both Lamé constants can be defined by Young's modulus E and the Poisson's ratio ν s as and

Direct Methods
This linear Navier-Cauchy equation can be solved by a direct method.Therefore, the Elmer solver module is used, which provides the LAPACK collection to address band matrices.Direct methods are known for their robustness, but their scaling of order n 3 leads to a high memory demand.Nevertheless, the present work uses a direct method solution procedure due the stability advantages.Further information can be found, e.g. in the book of Larson and Bengzon [44].

Coupling Conditions
The FSI problem has to fulfill certain coupling conditions on the interface I(t), based on physical principles [45].

Kinematic condition
The kinematic condition describes the continuity of the velocities on the interface, i.e.
The use of Lagrangian and Eulerian coordinate systems for the different solvers requires a mapping procedure.

Dynamic condition
The dynamic condition ensures that the forces that act on the interface are balanced due to Newton's third law 'Actio est Reactio'.The coupling condition reads Due to the different coordinate systems, a mapping procedure is also required here.

Geometric condition
The third coupling condition is the geometric condition.The condition ensures that the domain is continuous at the interface, i.e.
x f (t) = x s (t) on I(t).
Hence, the fluid and solid domains cannot overlap or separate at the interface.

Segregated Approaches
There are several mathematical and technical problems involved in the analysis of FSI.In most cases, two different subsystems are used for the governing equations.Even for one subsystem the uniqueness of the solution can be shown only locally in time.One fundamental problem is the different nature of the partial differential equations.Incompressible Navier-Stokes are of parabolic type, but the structural equation is of hyperbolic type.Therefore, the different coupling conditions are difficult to ensure on the interface.Segregated approaches are typically used to address FSI.The idea is to combine two different solvers, where each solver is specialized to solve either a fluid or a structural problem.The coupling is then fulfilled by an outer control instance.Due the maturity of each solver, this approach is often a quick possibility to treat complex application problems.A common segregated two-way coupled FSI workflow is illustrated in Fig. 2. The CFD solver on the left side of the sketch, numerically

Implementation
The FSI process which uses Elmer and OpenLB is depicted in Fig. 3.Note that a data based workflow is used to exchange information between the applications.Currently, the interface allows parallel execution of OpenLB, while Elmer is running in serial mode.A detailed explanation of each step in the workflow is given in the following procedure, which is executed for every coupling step.
1.The OpenLB instance calculates the hydrodynamic forces acting on the boundary for each solid node according to Eq. ( 17).
Figure 3: Fluid structure interaction workflow using OpenLB and Elmer.
2. The hydrodynamic forces are communicated and collected from each worker to the master process.
3. The master process maps the collected boundary forces to the finite element grid by integrating the force on each finite element mesh point.
4. The mapped boundary forces are written into an Elmer input deck file (.sif).

5.
Elmer is restarted by the master process using the input deck file (.sif) and a related restart file (.dat).
6.The Elmer instance is closed after the displacement velocity and the deformed mesh is written to disk as an unstructured mesh file (.vtu) and a new Elmer restart file (.dat) is created.
7. The master process reads the mesh file (.vtu) and uses the built-in OpenLB voxelizer, which decides whether a point is outside or inside the fluid domain and allows the later distance calculation.
8. The master process maps the displacement velocity of the FEM grid to the LBM link intersection points x w by a linear interpolation procedure and distributes the information to each worker process.
9. The OpenLB instance reconstructs the particle distribution functions for the fresh nodes by using the extrapolation refill algorithm (see Eq. ( 16)).
11.After the streaming step is executed, the unknown particle distribution function are calculated by the curved boundary approach using the mapped displacement velocity (see Eq. ( 15)).

Setup of the Coriolis Mass Flowmeter Test Case
The investigated CMF geometry is depicted in Figure 4.The CMF geometry consists of a flow divider that distributes the incoming mass flow in two U-shaped measuring pipes.After the flow passed both measuring pipes, a flow combiner unite the streams.The oscillation of the measuring pipe is initialized by an electromagnetic exciter at the top of both measuring pipes.The resulting oscillation signal is captured at sensor position 1 and 2. In addition, two node plates are used to damp the oscillation at the end of the pipes.The structural and fluid properties used in the Figure 4: Geometry representation and description of the investigated CMF without outer housing.
simulations are listed in Table 1 unless otherwise specified.The structural properties correspond to steel.The fluid density ρ f is chosen according to the density of water, but the dynamic viscosity η f is greatly increased to ensure a laminar flow.

Structural Domain
For the structural simulation setting, a zero displacement condition at the flanges is used, i.e.
Figure 5 indicates the flange faces in green, where this boundary condition is set.At the sensor exciter position an excitation load is applied where F x,max = 100 N and f exc is the excitation frequency.The force is only acting in the first period to excite the Eigenmode.

Fluid Domain
A Dirichlet velocity condition is applied as a boundary condition for the fluid domain at the inlet where R is the inlet radius and u y,max is the maximum velocity calculated according the used mass flow.This Poiseuille profile assumes a fully developed laminar pipe flow which is justified by an inlet Reynolds number of Re inlet < 337.The pressure on the outlet is set by a Dirichlet condition to As FSI is known to be time-consuming it is recommended to start the simulation with a converged flow field.Therefore, the fluid geometry is simulated without the structural simulation to initialize the flow field.On the measuring pipe walls no-slip boundary conditions are set.The flow velocity at the inflow is increased at the inlet for 5 seconds until the desired mass flow is reached.This initialization procedure prevents occurring pressure waves due to high gradients and adjusts the non-equilibrium parts of the particle distribution functions.

Coupling Conditions
On each grid point at the interface I(t) of fluid and solid, the mapped time dependent hydrodynamic force is applied via Thus the dynamic coupling condition is fulfilled.The velocity on the fluid structure interface is given by u f (t) = u s (t + ∆t c ) on I(t + ∆t c ), (30) where ∆t c is the coupling period.The occurring time shift is related to the staggered coupling scheme (see Sec. 2.3.2).The geometric condition of the interface is also influenced by the coupling period and is written as 3.2 Mesh Generation

Structural Domain
The computational mesh for the structural simulation is generated by the open source FEM preprocessor Gmsh [46].The mesh element is chosen according to the geometry pre-processor of OpenLB, because the extraction of a triangulated surface mesh out of tetrahedral volume mesh is straightforward.The choice of other mesh element shapes would lead to an additional triangulation step in every coupling period.In Table 2, the characteristic length scales of the FEM elements in the according regions are shown.Regions, where simulation results are extracted or high gradients may occur, are refined.Therefore, the sensor positions and the measuring pipes require small mesh elements.The maximal mesh element size is chosen with respect to the largest mesh element size that is used for the housing ∆x s = 0.035 m.The generated volume mesh in clip representation is depicted in Figure 6.The mesh contains 52624 nodes and 163164 elements.The five predefined regions support different refinement layers.This geometry adaptive mesh allows to reduce the amount of mesh points by consistent accuracy of the displacement signal at the sensor positions.The quality of the mesh was evaluated with the mesh criteria of Gmsh.Furthermore, the connection of critical mesh regions were checked, see Figure 7. Particularly in locations where elements are perpendicular to each other, defects may occur.

Fluid Domain
The meshing procedure for LBM is straightforward due to the equidistant uniform Cartesian mesh.The used discretization parameters for the two desired mass flows 20 000 kg h and 40 000 kg h are shown in Table 3.The resulting grid consists of 1.286 million grid cells.Figure 8 shows the voxelized measuring tubes at the beginning of the simulation.Additional two layers of solid cells cover the measuring pipes to allow the fluid-solid interface motion due to the pipe oscillation.

Results of the Coriolis Mass Flowmeter Test Case
After the mesh generation is completed, the Eigenfrequencies for the FEM mesh are calculated.The detection of the excitation frequency is a preliminary for the later phase shift calculation.Therefore, a modal analysis is performed with the structural solver Elmer.

Analysis
The first modal analysis describes the condition for the measuring pipes filled with resting air.The structural parameters of steel are listed in Table 1.Due to the low density of air compared to steel, the additional mass of air can be neglected.
Using the zero displacement boundary condition (see Eq. ( 25)), the first ten Eigenfrequencies of the FEM grid are calculated.The resulting values are shown in Table 4.A closer look to each Eigenfrequency reveals the physical meaning.The searched excitation mode is found at mode number 2 and the Coriolis twist mode corresponds to mode number 8. The excitation mode is related to a parallel movement of the pipes towards and away from each other.On the contrary, the  Coriolis twist introduces an additional twist of the pipes.For a better illustration both modes are displayed in a front and top view in Figure 9.The next step is the test of the transient structural simulation.Two major aspects are investigated: on the one hand the stability of the transient settings are estimated and on the other hand the resonant behavior are tested.The used structural boundary conditions are described in Section 3.1.1.In the first case, an excitation frequency different from the Eigenfrequency is selected to f exc = 50 Hz.In Figure 10, the structural response over time is plotted.It can be seen that the amplitude is strongly decreasing after the first period and no resonance is observable.This behavior was expected, because the excitation frequency and the Eigenfrequency are mismatched.Nevertheless, the transient simulation is stable over the entire simulation time.In the second configuration the excitation frequency is chosen with the Eigenfrequency to f exc = 104.28Hz.The displacement signal is depicted in Figure 11.The resonance is now clearly visible, which indicates that the results of the modal analysis are reliable and the transient simulation is also stable in the resonant case.Additionally a further modal analysis is examined for water conveying tubes, which are used in the FSI case.Hereby, the additional mass of water cannot be neglected.The water filled tubes are approximated by a fictitious density of the tubes ρ s f ictitious = 12 319 kg m 3 , which is calculated by the total mass of the measuring pipes divided by the volume of the structural pipe domain.The results are summarized and compared to the measurement data in Table 5.The excitation frequencies for air and water are in good agreement to the measurement data (error ≈ 3%).The errors for the Coriolis frequency seems to be squared due to the higher mode.

Phase Shift Calculation
After the modal analysis has determined the Eigenfrequency of the pipes filled with water, the transient fluid structure simulation is used to extract the phase shift.Firstly, the fluid field is initialized according to Section 3.1.2.The simulation procedure, which is described in Section 2.   be seen that the signal is almost sinusoidal in the first 5 cycles and the amplitude slowly decays over time.The last depicted periods show irregularities and differ from the expected pure sinusoidal course of the displacement signal.Furthermore a frequency analysis is performed to estimate resonance frequency.The results can be seen in Figure 13.The highest peak at 84.12 Hz in the frequency analysis is in good agreement with the estimated excitation frequency.A discrete Hilbert transformation is applied on the displacement signal to calculate the phase shift, see Figure 14.The stability of the FSI system is given for the first 8 cycles of the simulation.The symmetry condition, which should be fulfilled due to a axial-symmetric geometry, is only slightly violated.The error of the averaged phase shift value φ sim with respect to the experimental data φ ref is smaller than 5%, which is shown in Table 6.The relative errors for a mass flow of 20 000 and 40 000 kg h are less than 5%.Numerical experiments with a lower amount of coupling steps are diverging in the first period, which indicates that the reduction of coupling steps does not lead to satisfactory results.Consequentially, 202 coupling periods are necessary to stabilize the simulation.The simulation runtime was evaluated on a single node which consists of two deca-core Intel Xeon E5-2660 v3 processors.The comparison of the runtime to other numerical FSI simulations is depicted in Table 7.It can be seen that both computation time and calculated periods of the present  study are comparable to literature values.The computation runtime is estimated to 65 hours and over 3000 coupling steps are performed.It is noticeable that the computing time has hardly changed over the years.This is a consequence of the segregated approach, if two solvers are involved in the FSI approach.The partitioning of the fluid and the solid domain differs due to the numerical method and geometrical constraints.This implies that the exchanged information are collected and communicated between the solvers, which is a time consuming step that is very difficult to parallelize.

Conclusion and Outlook
An FSI approach was presented for the simulation of a CMF.Thereby, the open source framework OpenLB and Elmer were used to create a segregated approach.The target equations of the structural and fluid domain were described.In addition, the coupling conditions and the implementation were outlined in detail.The FEM mesh generation process utilized the open source meshing tool Gmsh to ensure a complete open source workflow.A modal analysis was performed to extract the excitation frequency of water and air conveying pipes.The found excitation frequency was in good agreement to experimental measurements (error ≈ 3%).Afterwards, the FSI simulation, which uses the determined excitation frequency, was executed.The FSI simulation was stable for several cycles and allows to extract the phase shift with a sufficient precision (error ≈ 5%).Therefore, the presented FSI approach for CMF is able to describe the operating principle of a CMF.Furthermore, the runtime time of the created FSI coupling were comparable to literature approaches using commercial software.Nevertheless, certain issues should be addressed in future studies.The FSI simulation becomes unstable after several periods.The reasons for this upcoming instability could be diverse.Firstly, the coupling time step could be decreased to reduce the time shift problem of the staggered coupling approach.Unfortunately this leads to an extended calculation time.Another possibility is the introduction of a subiteration scheme [47] that reduces the added mass effect due to the time shift.Further improvements can be made by the calculation of the hydrodynamic force, because momentum exchange based approaches suffer from inaccuracy, if too few points are used for integration.Therefore, a stress based calculation proposed in Geller et al. [21] may be an alternative.Furthermore, the applied linear mapping method between the uniform Cartesian LBM grid and the unstructured FEM grid can be improved by using more complex mapping methods [25].

Figure 1 :
Figure 1: Illustration of the used indexing convention.

Figure 6 :
Figure 6: Clip representation of the volume mesh.

Figure 7 :
Figure 7: Connection of the critical mesh regions.

3 . 3 ,
is executed in every coupling step ∆t c .The coupling period is chosen to the fluid time step ∆t f to minimize the time shift problem of the staggered approach.The simulation takes a total of 15 cycles which corresponds to approx.0.18 s at the Eigenfrequency.Every cycle consists of 202 coupling steps.The displacement signal is extracted at the sensor positions S1_plus, S1_minus, S2_plus and S2_minus, where plus and minus indicate the left and right measuring pipe, respectively.The written data files are post processed to extract the phase shift and the frequency of the displacement signals.The displacement signals of sensor S1_plus and S2_plus are depicted in Figure12.It can (a) Mode 2, front view (b) Mode 2, top view (c) Mode 8, front view (d) Mode 8, top view

Figure 13 :
Figure 13: Discrete Fourier analysis of the displacement signal.

Table 1 :
Structural and fluid properties.

Table 2 :
Averaged characteristic length scales of the finite element mesh regions.

Table 4 :
First ten Eigenfrequencies of the modal analysis and their physical meaning.Mode ω 2 in Hz 2 f in Hz

Table 5 :
Excitation and Coriolis twist frequency for water and air filled tubes in comparison to measurement data.

Table 6 :
Phaseshift errors for different mass flows and coupling steps per period.Mass flow in kg h φ sim in mrad φ ref in mrad Error in % Coupling steps

Table 7 :
Comparison of computation time between the present approach to literature values.