Experimental and Computational Investigation of a Fractal Grid Wake A. Fuchs W. Medjroubi H. Hochstein G. G ulker and J. Peinke Institute of Physics and ForWind University of Oldenburg K upkersweg 70 26129 Oldenburg Germany

2025-04-26 0 0 4.37MB 23 页 10玖币
侵权投诉
Experimental and Computational Investigation of a Fractal Grid Wake
A. Fuchs, W. Medjroubi, H. Hochstein, G. ulker, and J. Peinke
Institute of Physics and ForWind, University of Oldenburg, K¨upkersweg 70, 26129 Oldenburg, Germany
(Dated: October 13, 2022)
Fractal grids generate turbulence by exciting many length scales of different sizes simultaneously
rather than using the nonlinear cascade mechanism to obtain multi-scale structures, as it is the
case for regular grids. The interest in these grids has been further building up since the surprising
findings stemming from the experimental and computational studies conducted on these grids. This
work presents experimental wind tunnel and computational fluid dynamics (CFD) studies of the
turbulent flow generated by a space-filling fractal square grid. The experimental work includes
Particle Image Velocimetry (PIV) and hot-wire measurements. In addition, Delayed Detached Eddy
Simulations (DDES) with a Spalart-Allmaras background turbulence model are conducted using
the open-source package OpenFOAM. This is the first time DDES simulations are used to simulate
and characterize the turbulent flow generated by a fractal grid. Finally, this article reports on the
extensive statistical study and the direct comparison between the experimentally and numerically
acquired time series to investigate and compare one-point- and two-point statistics. Our goal is to
validate our computational results and provide enhanced insight into the complexity of the multi-scale
generation of turbulence using a fractal grid with a low number of fractal iterations. In particular,
we investigate the different turbulent structures and their complex interaction in the near-grid region
or the production regime of the fractal grid flow.
I. INTRODUCTION
Grid-generated turbulence has been investigated for
more than 70 years, being the most common method to
create turbulence in wind tunnels under controlled condi-
tions. The turbulence generated by the so-called regular
grids is nearly homogeneous and isotropic downstream
from the grid [
1
,
2
]. Nevertheless, the Reynolds numbers
and the turbulent length scales obtained are quite low
(
Reλ
=
u0λ/ν <
100, where
Reλ
is the Reynolds number
based on the Taylor microscale
λ
, the root mean square
of the velocity fluctuations
u0
=
pheu2i
and the kinematic
viscosity
ν
). This constitutes a handicap when it is
desired to investigate features of high Reynolds number
turbulence. Several other grids were proposed to overcome
this problem, including the use of mechanically moving
parts called active grids [
3
5
]. In [
6
], the use of an active
grid resulted in turbulence with higher Reynolds number
(
Reλ
= 390), higher turbulence intensities, and larger
integral length scales. Moreover, the resulting energy
spectrum had a wider inertial subrange, about two orders
of magnitude in wavenumber, compared to a regular grid.
Hurst and Vassilicos [
7
] proposed another method to
generate higher Reynolds numbers and high turbulence
intensities, which does not involve mechanically moving
grid components and is easier to use and construct. In
[
7
], the authors characterized the turbulence generated
by three types of fractal grids (cross, I and square shaped
grids), using 21 different grids in wind tunnel experiments.
The main idea behind fractal grids is to use grids that con-
sist of a multi-scale collection of obstacles and openings
based on a specific pattern, which is repeated in increas-
ingly numerous copies at smaller scales. Fractal therefore
emphasizes that these grids have multiple scales that are
self-similar. This property helps reduce the parameters
necessary to specify the design of these multi-scale grids.
Besides the easy way to generate high Re-number turbu-
lence, the topic of fractal grid turbulence has attracted
considerable attention due to the specific features of the
wake flow.
The concept behind fractal grids is to generate turbu-
lence by directly exciting a wide range of length scales at
once, rather than using the nonlinear cascade mechanism
to obtain multi-scale structures, as it is the case for regular
grids. These flow structures of different scales influence
each other and show very different properties compared
to all previously documented turbulent flows. Interesting
results were obtained using Fractal Square Grids (FSG)
with low blockage ratios (around 25%). These results
include obtaining higher Reynolds numbers, even in small
and conventional-size wind tunnels, compared to regular
grids. In addition, it was shown that there exists an ex-
tended region between the fractal grid and a downstream
distance, where the turbulence intensity increases progres-
sively to reach a maximum value at a distance
x
=
xpeak
[
7
]. This region is called the production region. The
turbulence intensity then decays downstream from
xpeak
,
following what was thought first to be an exponential
law [
7
,
8
] but was shown later on to be a fast power law
[
9
]. These findings contradict what is typically reported
for regular grid turbulence, namely that the turbulence
decays directly at the lee of the grid following a power-law
[
10
]. Another interesting characteristic of fractal grid tur-
bulence is that longitudinal and lateral turbulent integral
length scales and Taylor microscales were reported to vary
slowly downstream of the grid [7].
A more detailed study of the decay region mentioned
in [
7
] was conducted by Seoud and Vassilicos [
8
], which
confirmed the findings in terms of turbulent length scale
behavior and turbulence intensity decay law. Moreover, it
was shown that the ratio of the longitudinal length scale
and Taylor microscale remains approximatively constant
arXiv:2210.06208v1 [physics.flu-dyn] 12 Oct 2022
2
in the decay region. In [
8
,
9
,
11
,
12
] this striking be-
havior was set into the context of a non constant energy
dissipation rate
Cε
. Contrary to what is predicted by
the Richardson-Kolmogorov phenomenology, the authors
could collapse one-dimensional energy spectra at differ-
ent wake positions using only one length scale and
u0
.
They also presented assessments of the homogeneity and
isotropy of the turbulence produced by FSG [7].
Vassilicos and Mazelier [
11
] related
xpeak
to the wake-
interaction length scale
x
, which is the ratio of the length
of the biggest square bar and its thickness.
x
is found to
be the appropriate length scale to characterize first and
second-order statistics of the turbulent flow generated
by FSG. The authors also showed that the turbulence
statistics are non-homogeneous and non-Gaussian in the
production region but become homogeneous and Gaussian
in the decay region. They suggested that the FSG can
be used for flow control and to enhance turbulent mixing
properties by understanding and determination of how
xpeak
and the generated turbulence intensities depend on
fractal grid geometry.
Several Direct Numerical Simulations (DNS) were con-
ducted to investigate the flow generated by fractal grids
[
13
17
]. Although these simulations qualitatively repro-
duced some of the underlying characteristics of fractal grid
turbulence, quantitative comparisons with experimental
data was not possible for mainly three reasons. First,
the simulations were conducted at very small Reynolds
numbers compared to the experiments. As an example,
in [
15
] the Reynolds number based on the effective mesh
size (see Section II) is
ReMef f
= 4430 for the simula-
tions and
ReMef f
= 20800 for the experiments, which
is a reduction with a factor of 5. Second, fractal grids
with a smaller iteration number were simulated, mainly 3
iterations, compared to the available experimental studies
where the iterations vary from 4 to 6. Finally, the extent
of the computational domain in the streamwise direction
was limited due to the high computational cost associated
with DNS. As all the flow details have to be resolved, the
meshes used in the different DNS simulations are huge
in terms of the number of nodes. Therefore, substantial
computational power is needed to keep the computational
time realistic. In the following, we summarize two impor-
tant DNS reference papers.
Laizet and Vassilicos [
13
] conducted DNS simulations
using the numerical code Incompact3d based on sixth-
order compact schemes for spatial discretization and
second-order Adams-Bashforth schemes for temporal dis-
cretization. They modeled a regular and a fractal grid
(with the same effective mesh size) using Immersed Bound-
ary Method (IBM) and their total numerical mesh of 765
million points. They could recover the production and
decay regions as well as the decay behavior of the turbu-
lence intensity, although they did not test the simulation
results to determine the nature of the decay law. More-
over, the turbulence recovered in the wake of the fractal
grid was found to be not homogeneous, which contradicts
the findings in the experiments [
7
,
11
]. They related this
observation to the small iteration number of their fractal
grid (
N
= 3) and the domain size chosen to be insufficient
for reaching homogeneity.
In a more recent paper, Laizet and Vassilicos [
18
] used
a new version of Incompact3d, which scales better when
used with large numbers of computational cores. They
simulated one regular and three square fractal grids (with
different aspect ratios). The simulations required the
use of 3456 computational cores. The DNS results show
that the turbulence produced by the fractal grids is more
intermittent than the one produced by regular ones. Nev-
ertheless, the extent of the computational domain and
number of fractal iterations (
N
= 3) were still insufficient
to obtain homogeneous turbulence in the decay region.
Therefore, they emphasize the need for more simulations
with extended computational domain and bigger itera-
tions (
N >
3) and experiments with
N
= 3 to be able to
compare the results with DNS. The authors also showed
that the value of
xpeak
introduced in [
11
] does not apply
for their grids, which have a smaller N, and introduced a
new xpeak formula based on N and the blockage ratio σ.
IIn the present paper, the wake of the flow through an
N
= 3 fractal square grid is investigated numerically and
experimentally. The experimental data was acquired from
PIV and hot-wire anemometry. Our main purpose is to
validate our numerical simulations against the experimen-
tal results and to characterize the turbulence generated
by fractal grids in terms of turbulent kinetic energy decay
law and one-point statistical quantities such as turbulence
intensity, third and fourth central moment and probability
density functions of the velocity fluctuations. Further-
more, the characterization of such turbulent fields will
assess the capabilities of DDES in reproducing the ex-
perimental results. A further goal of this investigation
is to compare experimental and numerical data of the
wake of a fractal grid with the experimental examination
of a regular grid with the same mesh size and blockage
ratio. Our simulations and experiments also involve an ex-
tended domain to obtain more homogeneous and isotropic
turbulence than DNS obtained so far.
This paper is organized as follows. Section 2 introduces
the geometrical aspects of the fractal square grid and of
a regular grid having the same mesh size and blockage
ratio. The regular grid has been investigated only experi-
mentally and is used here for comparison purposes. The
experimental setup and details of the PIV and hot-wire
anemometry used are presented in Section 3. Section
4 is dedicated to the numerical setup, including a brief
overview of the turbulence modeling and the numerical
solver. Section 5 discusses the low-order and higher-order
statistics results. Finally, in Section 6 the conclusions are
presented with an outlook for further investigations.
3
II. GRID GEOMETRICAL PARAMETER
Figure 1shows the grids used in the present work. They
are placed at the inlet of the test section of a wind tunnel
in the experiments and considered when implementing
the corresponding numerical simulations (see also Figure
2b).
(a) (b)
FIG. 1. Illustration of (a) the used
N
= 3 space-filling square
fractal grid (FSG) geometry and (b) the regular grid used for
comparison.
In general, fractal grids are constructed from a multi-
scale collection of obstacles based on a single pattern
repeated in increasingly numerous copies with different
scales. As presented in [
7
], fractal grids can be constructed
using different geometrical patterns (see Figure 1 in [
15
]).
We selected a frequently used pattern for our fractal grid
[
8
,
9
,
11
], to be able to set our results in the context of
other research activities. The pattern of our fractal grid
is based on a square shape with
N
= 3 fractal iterations.
The fractal iteration parameter is the number the square
shape that is repeated at different scales. At each iteration
(
j
= 0
, ..., N
1), the number of squares is four times
higher than in the iteration
j
1. Each scale iteration
j
is defined by a length
Lj
and a thickness
tj
of the
square bars constituting the grid. The thickness of the
square bars in the streamwise direction is kept constant.
The dimensions of the square patterns are related by the
ratio of the length of subsequent iterations
RL
=
Lj
Lj1
and by the ratio of the thickness of subsequent iterations
Rt
=
tj
tj1
; respectively. The geometry of the fractal
grid used in this paper is completely characterized by
two further parameters, namely the ratio of the length
of the first iteration to the last one
Lr
=
L0
LN1
and the
ratio of the thickness of the first iteration to the last one
tr
=
t0
tN1
. Moreover, a fractal grid is said to be space-
filling, when its fractal dimension
Df
= 2 (see [
7
,
19
], for
the definition of Df), which is the case when RL= 0.5.
Contrary to regular grids, fractal grids, especially frac-
tal square grids, do not have a well-defined mesh size.
However, an equivalent effective mesh size was defined in
[7] as
Meff =4T2
P1σ, (1)
where
T2
is the cross-section of the wind tun-
nel/simulation domain,
P
is the perimeter of the fractal
grid and
σ
is the blockage ratio, which can be determined
by a calculation using the following formula
σ=A
T2=
L0t0
N1
P
j=0
4j+1Rj
LRj
tt2
0
N1
P
j=1
22j+1R2j1
t
T2.
(2)
A
is the sum of the areas
Aj
covered by each fractal
iteration
j
and corrected by the areas covered by two
iterations.
A
is therefore, the fractal grid’s total area.
It should be noted that the local blockage ratio varies
strongly for different parts of the grids. A complete
quantitative description of the fractal grid we used in this
study is shown in table I.
Another important length, which seems to collapse
several turbulent statistical quantities of the streamwise
fluctuating velocity, is the position of maximum turbu-
lence intensity
xpeak
, which also defines the production
(
x < xpeak
) and decay regions (
x > xpeak
). It was intro-
duced by Mazellier and Vassilicos [
11
] and defined using
the empirical formula
xpeak 0.45L2
0
t0
,(3)
which relies upon
xpeak
to the geometrical properties of
the fractal grid. This means that the turbulence generated
by fractal grids can be “manufactured”.
To account for the distance at which the different wakes
generated by the fractal grid bars interact, [
11
] introduced
the wake-interaction length scale
x
. It is defined with
the formula
x
=
L2
0
t0
and is in our case
x
= 953 mm
(see definitions of
L0
and
t0
in table I). The definitions
of the fractal grid effective mesh size eq.
(1)
and the
blockage ratio eq.
(2)
has the advantage of returning the
effective mesh size when applied to a fractal grid so that
a comparison can be made between regular and fractal
grids based on
Meff
. In this respect, we experimentally
investigated a regular grid (shown in Figure 1b) with the
same mesh size (
M
=
Meff
) and blockage ratio as the
used fractal grid for comparison.
III. EXPERIMENTAL SETUP
The experiments were conducted in a closed loop wind
tunnel with test section dimensions of 200 cm x 25 cm
x 25 cm (length x width x height) at the University of
Oldenburg. The wind tunnel has a background turbu-
lence intensity along the centerline of the complete test
section of approximately 2% for
U
10 m/s. The inlet
velocity was set to 10 m/s, corresponding to a Reynolds
4
TABLE I. Geometrical properties of the fractal grid used in this study.
Nσ/% L0/mm t0/mm RLRtLrtrMeff /mm T/mm
3 38.2 138.4 20.1 0.52 0.36 3.7 7.7 24 250
number related to the biggest grid bar length
L0
of about
ReL0=UL0= 83800,
where
ν
is the kinematic vis-
cosity. Optical access is provided through the wind tunnel
side walls. The experimental data were acquired from hot-
wire anemometry and Particle Image Velocimetry (PIV)
measurements.
Constant temperature anemometry measurements
of the velocity were performed using (Dantec 55P01
platinum-plated tungsten wire) single-hot-wire with a
wire sensing length of about
lw
= 2
.
0
±
0
.
1 mm and a
diameter of
dw
= 5
µm
which corresponds to a length-
to-diameter ratio of
lw/dw
400. A StreamLine mea-
surement system by Dantec in combination with CTA
Modules 90C10 and the StreamWare version 3.50.0.9 was
used for the measurements. The hot-wire was calibrated
with Dantec Dynamics Hot-Wire Calibrator. The over-
heat ratio was set to 0.8. In the streamwise direction,
measurements were performed on the centerline between
5 cm
x
176 cm distance to the grid. The data was
sampled with the frequency
fs
= 60 kHz with a NI PXI
1042 AD-converter and a total of 3.6 million samples were
collected per measurement point, representing 60 seconds
of measurement data. To satisfy the Nyquist condition,
the data was low-pass filtered at frequency fl= 30 kHz.
In addition to the hot-wire measurements, the flow
velocity is measured with two-component (2C) two-
dimensional (2D) Particle Image Velocimetry (PIV). A
double-pulsed Nd:YAG laser system illuminates the flow
with a maximum energy output of 2x380 mJ pulse
1
(
λ
= 532 nm, Quantel Brilliant B Twin). The laser system
is set at a sampling frequency of 10 Hz and a laser pulse
delay of 63
µs
(time between two image frames) in order
to record statistically independent images. To illuminate
the desired field of view, a laser light sheet is formed with
an optical arrangement (spherical lens
(f=100 mm)
and cylindrical lens (
f
= 200 mm)) to converge the sheet
into a minimum thickness of around 2 mm. The light
sheet plane is in the vertical mid-plane of the wind tun-
nel. The flow was seeded with Di-Ethyl-Hexyl-Sebacat
(DEHS) droplet aerosols generated by a Palas AGF 5
aerosol generator for atomizing liquids. The diameter of
the mist of DEHS droplets used in the present study is
about 0,3
µm
to 1
µm
. At each measuring station (there
are 15 stations in total) 3000 double images are recorded
by a cooled digital 14 bit CCD Camera (PCO1600) with
1600x1200 pixel resolution. The camera looks perpendic-
ular to the light sheet from the side and is synchronized
with the laser sheet pulse at 10 Hz. The camera is fit-
ted with a Nikon mikro Nikkor 55 mm objective lens
set with an aperture setting of f/8. A minimal effect of
peak locking was found for these experimental conditions.
The PIV analysis is done using commercial PIV software
(
PIV View
version 3.5.7) using an interrogation window
of 2 mm
2
. The vector fields are obtained by processing
the PIV images using a recursive cross-correlation engine
with a final interrogation window of size 16 x 16 pixels
with 75% overlap.
IV. NUMERICAL SETUP
A. Turbulence modelling
The three-dimensional, incompressible Navier-Stokes
equations describe the flow generated by a fractal grid.
The equations are discretized and solved using a turbu-
lence model. In this investigation, the Delayed Detached
Eddy Simulation (DDES) [
20
] with a Spalart-Allmaras
background turbulence model [
21
], commonly referred to
as SA-DDES is used. DDES is a hybrid method stem-
ming from the Detached Eddy Simulation method (DES)
[
20
], which involves the use of Reynolds Averaged Navier-
Stokes Simulation (RANS) at the wall and Large Eddy
Simulation (LES) away from it. This method combines
the simplicity of the RANS formulation and the accuracy
of LES, with the advantage of being less expensive in
terms of computational time when compared with pure
LES. DDES is an improvement of the original DES for-
mulation, where the so-called ”modeled stress depletion”
(or MSD), is treated [
22
,
23
]. MSD occurs when the LES
mode of the DES method is operating in the boundary-
layer [
23
]. In DDES, a shield function is derived, which
guarantees that the boundary layer is solved by the RANS
mode of the DES model. In the following, a brief descrip-
tion of SA-DDES is given. For more information refer to
[20,24].
In the framework of eddy viscosity models, the Reynolds
stresses are expressed in terms of the eddy viscosity
νt
and the strain-rate tensor
Sij
given by the formulation
huiuji
= 2
νtSij
. We consider the Spalart-Allmaras
model, which is a RANS model based on the eddy viscosity
approach [
21
], where a transport equation is defined for
the modified turbulent viscosity ˜νas follows
˜ν
t +uj
˜ν
xj
=cb1(1 ft2)˜
S˜νhcw1fwcb1
κ2ft2i˜ν
d2
+1
σ
xj(ν+ ˜ν)˜ν
xj+cb2
˜ν
xi
˜ν
xi.
(4)
˜ν
is related to the turbulent eddy viscosity by the equation
νt= ˜νfν1,(5)
5
and the function fν1is defined as
fν1=χ3
χ3+c3
ν1
,(6)
where
χ
=
˜ν
ν
,
ν
being the kinematic viscosity.
˜
S
is defined
as
˜
S=fν3S+˜ν
κ2+d2fν2,(7)
where
S
is the magnitude of the vorticity, and the func-
tions fν2and fν3are defined as follows
fν2=1
(1 + χ/cν2)3,
fν3=(1 + χfν1)(1 fν2)
χ,
(8)
where,
d
is the distance to the closest wall. Finally,
fω
is a function, which also contains the distance
d
and is
defined as
fω=g1 + c6
ω3
g6+c6
ω31
6
,(9)
where
g
and
r
are defined as
g
=
r
+
cω2
(
r6r
) and
r
=
˜ν
˜
Sκ2d2
. The remaining model constants are defined in
[21].
The DES model used in this contribution was obtained
by replacing the distance variable
d
with the distance
˜
d
defined by the expression
˜
d=min(lRAN S , lLES ) = min (d, CDES ∆) ,(10)
where,
CDES
= 0
.
65 is a constant and
∆ = max(∆x,y,z)
represents the grid size (the
filter size). One of the well-known shortcomings of DES
is the so-called Grid Induced Separation (GIS) which can
lead to the reduction of the RANS Reynolds stresses.
This problem occurs when the model operates in LES
mode in the boundary-layer [
23
], which results in the
separation point moving upstream (thus the name GIS).
This results from the fact that in refined regions of the
mesh,
lLES
becomes smaller and the switch from RANS
to LES occurs prematurely. To overcome this problem,
the ”detached” version of the DES method was proposed
in [
22
,
23
]. In DDES, the switching from RANS to LES
mode is performed in a more complex manner, where an
empirical shielding function
fd
is introduced. This results
in the following expression for the DDES length scale
˜
lDDES =lRANS fdmax(0, lRANS lLES ),(11)
where
fd
tends towards 0 inside the boundary-layer region
and towards 1 away from the boundary-layer. Moreover,
fd
is a continuous function that enables a smooth transi-
tion between the RANS and the LES regions.
B. Mesh and Numerical Method
The numerical simulation was set up analogous to the
experiments in order to compare the results consistently.
The open source code OpenFOAM [
25
] was used to solve
the incompressible Navier-Stokes equations. OpenFOAM
is based on a collocated formulation of the finite volume
method, and it consists of a collection of libraries writ-
ten in C++, which can be used to simulate a large class
of flow problems. For more information, see the official
documentation [
25
]. The solver used in this investigation
is the transient solver pimpleFoam, which is a merging
between the PISO (Pressure implicit with splitting of op-
erator) and SIMPLE (Semi-Implicit Method for Pressure
Linked Equation) algorithms. A second-order central-
differencing scheme is used for spatial discretization, and
a backward, second-order time-advancing schemes were
used. The solver is parallelized using the Message-Passing
Interface (MPI), which is necessary for problems of this
size. The SA-DDES model used in this work is already
implemented in the official release of OpenFOAM.
The numerical mesh was generated using the built-
in OpenFOAM meshing tools blockMesh and snappy-
HexMesh [
25
]. First, a background mesh constituted of
hexahedral cells filling the entire computational domain
is generated using blockMesh. The snappyHexMesh tool
is then applied on the generated hexahedral mesh, which
includes the fractal mesh. As a result, an unstructured
mesh is obtained, where regions of interest in the wake
are refined, as shown in Fig. 2a. Locally refining the
mesh required the use of the parallel option of the snap-
pyHexMesh tool. The mesh obtained is composed of a
total of 24 million cells.
The fractal grid is simulated in a domain with similar
dimensions as the real wind tunnel. The domain begins 2
m upstream of the fractal grid and covers a distance of 2
m downstream (see Fig. 2b). The flow-parallel boundaries
are treated as frictionless walls, where the slip boundary
condition was applied for all flow variables. At the inflow
boundary, Neumann boundary condition was used for the
pressure and Dirichlet condition for the velocity. At the
outflow boundary, the pressure was set to be equal to the
static pressure and a Neumann boundary condition was
used for the velocity. On the fractal grid, a wall function
is used for the modified viscosity
˜ν
, with the size of the
first cell of the mesh in terms of the dimensionless wall
distance is y+200 [26].
For each simulation, 480 processors were used, and for
each time step 4.5
GB
of data was collected for post-
processing. The data sampling frequency was 60
kHz
,
chosen to match the experimental one. It took approx-
imately 72 hours to simulate one second of data and a
total of 16 seconds of numerical data was collected. The
data was collected in the same positions as for the experi-
mental study. The numerical simulations were conducted
on the computer cluster of the ForWind Group [27].
摘要:

ExperimentalandComputationalInvestigationofaFractalGridWakeA.Fuchs,W.Medjroubi,H.Hochstein,G.Gulker,andJ.PeinkeInstituteofPhysicsandForWind,UniversityofOldenburg,Kupkersweg70,26129Oldenburg,Germany(Dated:October13,2022)Fractalgridsgenerateturbulencebyexcitingmanylengthscalesofdi erentsizessimultan...

展开>> 收起<<
Experimental and Computational Investigation of a Fractal Grid Wake A. Fuchs W. Medjroubi H. Hochstein G. G ulker and J. Peinke Institute of Physics and ForWind University of Oldenburg K upkersweg 70 26129 Oldenburg Germany.pdf

共23页,预览5页

还剩页未读, 继续阅读

声明:本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。玖贝云文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知玖贝云文库,我们立即给予删除!
分类:图书资源 价格:10玖币 属性:23 页 大小:4.37MB 格式:PDF 时间:2025-04-26

开通VIP享超值会员特权

  • 多端同步记录
  • 高速下载文档
  • 免费文档工具
  • 分享文档赚钱
  • 每日登录抽奖
  • 优质衍生服务
/ 23
客服
关注