Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
1/6
Organization (S): EDF/MTI/MFTT, CS IF
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and the plates
Document: V2.03.103
SDLS103 - Coques coaxial under flow
annular: inertial coupling between modes
Summary
One considers a hardware configuration made up of two coaxial cylindrical hulls, in interaction
with a fluid running out in annular space separating the hulls.
The goal of the case-test is to validate the model of coupling fluid-structure developed in the operator
CALC_FLUI_STRU for this type of configuration.
One is interested here more particularly in the taking into account of the inertial coupling between modes, obtained out of water
at rest (mean velocity of flow null). The reference solution is provided by a calculation carried out with
Code_Aster implementing operator CALC_MATR_AJOU.
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
2/6
1
Problem of reference
1.1 Geometry
The studied configuration consists of two 4 meters height coaxial cylindrical hulls:
2m
4m
0.1m
The internal hull has an average radius of 1 meter and a thickness of 1 centimetre.
The external hull has an average radius of 1,10 meter and a thickness of 1 centimetre.
1.2 Properties
material
The material constituting the two hulls is steel. Its physical characteristics are:
= 7800 kg/m3
E = 2. 1011 Pa
= 0,3
1.3
Boundary conditions and loadings
The conditions of self-supporting quality are the same ones for the two hulls: partly embedded ends
lower (Z = 0) and free partly higher (Z = 4 m).
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
3/6
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is provided by a calculation carried out by means of Code_Aster implementing
operator CALC_MATR_AJOU.
With this intention, one uses a grid on which the elements of the structure (hulls internal are defined
and external), voluminal elements of fluid (annular space) and elements of interface
fluid-structure. This method is described completely in [bib1].
On the elements of interface, the condition normal speed limit null is imposed, translating
condition of nonpenetration of the fluid in the hulls.
On the sections of input and output of annular space, the limiting condition of null potential is
imposed, translating the condition of null disturbed pressure at the ends.
One carries out the first modal calculation of the structure in air. Operator CALC_MATR_AJOU allows
then to calculate the matrix of mass added by the fluid, projected on the basis of modal
structure in air. One can then recombine this matrix of mass added with the matrix of mass
generalized of the structure in air, then to solve a new modal problem which leads to
characteristics of the water system at rest. Results obtained for the Eigen frequencies
system out of water at rest constitute the reference solution.
Characteristics of the grid 17 nodes on a vertical generator
:
60 nodes on a crown
960 meshs QUAD4 on each hull
=> 1920 meshs for the interface fluid-structure
2880 meshs HEXA8 for the fluid field
180 meshs QUAD4 for the section of input of annular space
180 meshs QUAD4 for the section of output of annular space
2.2
Results of reference
N°
mode
1 2 3 4 5 6 7 8 9 10 11
12
Freq
(Hz) 5,65 5,65 6,48 6,48 9,34 9,34 20,82 20,82 28,22 28,22 31,48 31,4
8
2.3
Uncertainty on the solution
In the studied case, the half-thicknesses of hull account for 10% of the dimension of the annular play.
In the calculation of reference realized by operator CALC_MATR_AJOU, the half-thicknesses are
neglected for the definition of the fluid field. One thus awaits a systematic error about
5% on the frequencies, because of this approximation.
Moreover, differences of modeling and resolution between the numerical method implemented
in operator CALC_MATR_AJOU and the analytical model developed in the operator
CALC_FLUI_STRU induce an additional variation.
2.4 References
bibliographical
[1]
G. ROUSSEAU: “Schedule of conditions and principle of realization of the calculation of stiffness and
of damping added in Code_Aster ", HP-51/96/005/B.
[2]
L. PEROTIN, “Note of principle of model MOCCA_COQUE”, HT-32/95/021/A.
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
4/6
3 Modeling
With
3.1
Characteristics of modeling
The geometry of the structures and the characteristics of material constituting the hulls were
presented before.
Concerning the absolute roughness of wall of the structures, one takes a value of 105 meter.
The surrounding fluid is water. Values taken for the density and viscosity
kinematics are respectively
= 106 m2/S.
F = 1000 kg/m3 and
F
One considers the flow not confined upstream and downstream from the structures.
The problem of coupling fluid-structure is solved by analytical model MOCCA_COQUE [bib2]
integrated in Code_Aster (operator CALC_FLUI_STRU), for a mean velocity of flow
null: one thus obtains the modal characteristics of the water system at rest, by catch in
count effects of added mass.
3.2 Characteristics
grid
Compared to the calculation of reference, the characteristics of the grid are similar, with the difference
close the fluid field is not represented any more. In order to make the two hulls interdependent one of
the other, one adds a group of meshs connecting the nodes to the level of embedding.
One a: 17 nodes on a vertical generator,
60 nodes on a crown,
960 meshs QUAD4 on each hull,
60 meshs QUAD4 to solidarize the two hulls (bases embedded).
3.3
Stages of calculation
The definition of the characteristics of a hardware configuration made up of two hulls
cylindrical coaxial for a calculation of coupling fluid-structure is given via
operator DEFI_FLUI_STRU key word factor COQUE_COAX.
The resolution of the coupling fluid-structure for a configuration of the type “coaxial hulls” and it
calculation of the modal parameters (reduced frequencies and depreciation) and deformed modal out of water
at rest are realized with operator CALC_FLUI_STRU.
3.4 Functionalities
tested
Commands
Key word factor
Key word
DEFI_FLUI_STRU COQUE_COAX
MASS_AJOU
“OUI”
CALC_FLUI_STRU BASE_MODALE
MODE_MECA
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
5/6
4
Results of modeling A
4.1 Values
tested
The comparisons relate to the water frequencies at rest of the first 12 modes of the system.
Numbers of the modes
CALC_FLUI_STRU
CALC_MATR_AJOU
variation
(H = 9cm)
(H = 10cm)
1 and 2
5,30 Hz
5,65 Hz
6,2%
3 and 4
6,12 Hz
6,48 Hz
5,5%
5 and 6
8,69 Hz
9,34 Hz
6,9%
7 and 8
21,96 Hz
20,82 Hz
5,5%
9 and 10
29,34 Hz
28,22 Hz
3,9%
11 and 12
33,75 Hz
31,48 Hz
7,2%
One gives, for information, the values of the frequencies of these modes in air:
Numbers of
Hull in mvt
Command of hull
Command of beam
Frequency
modes
1 and 2
external
3
1
25,15 Hz
3 and 4
intern
3
1
26,12 Hz
5 and 6
external
4
1
31,91 Hz
7 and 8
intern
2
1
36,85 Hz
9 and 10
intern
4
1
37,42 Hz
11 and 12
external
2
1
39,49 Hz
4.2 Remarks
The results are in conformity so that one could wait. One observes indeed:
· a systematic error of about 5%, because of not taken into account thicknesses
of hull for the definition of the fluid field in the calculation of reference;
· a residual variation due to the differences of modeling and resolution between the two
operators CALC_MATR_AJOU and CALC_FLUI_STRU.
Modes 1 to 6 are modes for which the structure is strongly coupled with the fluid. In
practical, these modes correspond to movements of the hulls internal and external overall in
opposition of phase. For these modes, the terms of added mass are theoretically proportional
3
F
R
with
, where R indicates the average radius and H the thickness of annular space.
H
For these the first six modes, the results provided by CALC_FLUI_STRU lead to frequencies
clean lower than those calculated by CALC_MATR_AJOU. Indeed, H being weaker, them
terms of added mass are larger.
Modes 7 to 12 are modes for which the structure is slightly coupled with the fluid. In
practical, these modes correspond to movements of the hulls internal and external almost in
phase. Thus, the terms of added mass are theoretically proportional to the water mass
involved, i.e. with
F
RH.
In this case, it is normal that the results provided by CALC_FLUI_STRU lead to
Eigen frequencies higher than those calculated by CALC_MATR_AJOU. Indeed, H being more
weak, the terms of added mass are smaller.
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Code_Aster ®
Version
5.0
Titrate:
SDLS103 - Coques coaxial under annular flow
Date:
09/10/01
Author (S):
A. ADOBES, Key Mr. LAINET
:
V2.03.103-A Page:
6/6
4.3 Parameters
of execution
Version: NEW 5.02.20
Machine: SGI ORIGIN 2000
Obstruction memory: 64 Mo
Time CPU To use: 50,63 seconds
5
Summary of the results
The comparison of the water frequencies at rest calculated with operators CALC_FLUI_STRU and
CALC_MATR_AJOU is satisfactory. The variations met between these two operators are explained by
the fact that they use a different modeling.
Handbook of Validation
V2.03 booklet: Linear dynamics of the hulls and plates
HI-86/01/021/A
Outline document