Models - Cfd.fluid Damper
Models - Cfd.fluid Damper
Models - Cfd.fluid Damper
This model is licensed under the COMSOL Software License Agreement 6.2.
All trademarks are the property of their respective owners. See www.comsol.com/trademarks.
Introduction
Fluid dampers are used in military devices for shock isolation and in civil structures for
suppressing earthquake-induced shaking and wind-induced vibrations, among many other
applications. Fluid dampers work by dissipating the mechanical energy into heat (Ref. 1).
This example shows the phenomenon of viscous heating and consequent temperature
increase in a fluid damper. Viscous heating is also important in microflow devices, where
a small cross-sectional area and large length of the device can generate significant heating
and affect the fluid flow consequently (Ref. 2).
Model Definition
The structural elements of a fluid damper are relatively few. Figure 1 depicts a schematic
of the fluid damper modeled herein with its main components: damper cylinder housing,
piston rod, piston head, and viscous fluid in the chamber. There is a small annular space
between the piston head and the inside wall of the cylinder housing. This acts as an
effective channel for the fluid. As the piston head moves back and forth inside the damper
cylinder, fluid is forced to pass through the annular channel with large shear rate, which
leads to significant heat generation. The heat is transferred in both the axial and radial
directions. In the radial direction, the heat is conducted through the cylinder house wall
and convected to the air outside the damper, which is modeled using the Newton’s
convective cooling law.
Cylinder Fluid
P(t)
Chamber 1 Chamber 2
Piston head
You make use of the axially symmetric nature of the fluid damper and model it in a 2D-
axisymmetric geometry as shown in Figure 2. The geometric dimensions and other
parameters of the damper are taken according to Ref. 1 to represent the smaller, 15 kip
damper experimentally studied therein. Thus, the piston head has a diameter of 8.37 cm,
the piston rod diameter is 2.83 cm, and the gap thickness is about 1/100 of the piston
Figure 2: Geometry and mesh. The domains (from left to right) represent: piston rod, piston
head and damping fluid space, and the damper outer wall.
FLUID FLOW
The fluid flow in the fluid damper is described by the incompressible Navier–Stokes
equations, solving for the velocity field u u w and the pressure p:
u T
+ u u = – pI + u + u
t
u = 0
= 0 – T – T0 (1)
T T
C p + C p u T + q = Q + u + u : u
t
where the second term on the right-hand side represent the heat source from viscous
dissipation. Hence, the problem is a fully coupled fluid-thermal interaction problem.
In the solid domain of the cylinder house wall, this equation reduces to conductive heat
transfer equation without any heating source.
The heat flux boundary condition based on the Newton’s cooling law is applied on the
outside boundaries of the cylinder house wall. The temperature field is continuous
between the fluid and solid domains. The ends of the damper connected to the structures
outside are kept at constant temperature.
The piston head movement is provided as harmonic oscillations with given amplitude and
frequency, z a0sin2ft, the piston head starting the simulation in its middle position.
The motion is modeled using the arbitrary Lagrangian–Eulerian (ALE) deformed mesh.
The ALE method handles the dynamics of the deforming geometry and the moving
boundaries with a moving grid. The Navier–Stokes equations for fluid flow and heat
equations for temperature variation are formulated in these moving coordinates.
Note that the simulation results for the temperature are presented in degrees Fahrenheit
for the sake of easier comparison with the experimental measurements.
Figure 4 shows the temperature of the inner wall of the damper at the end-of-stroke
position z U0. This corresponds to the internal probe position under experiments
performed in Ref. 1. The simulation results show very good agreement with the
experimental measurements (see Fig. 9 in Ref. 1).
Figure 5 shows the temperature variation along the inner wall of the damper after 10 s and
40 s of loading. It clearly shows that the temperature at the probe position does not
represent the maximum temperature within the damper. This supports the conclusion
drawn in Ref. 1, where the choice of the probe positioning was limited by the construction
of the outer shell of the damper. Figure 5 also shows that the temperature near the center
of the damper increases by about 100 degrees already after few loading cycles.
The steel material needed for the damper solid parts is available in the built-in material
library. You create a user-defined material for the silicone oil. Such damping fluids are
typically characterized by the density, kinematic viscosity at the temperature 25°C, and so-
called viscosity temperature coefficient, VTC 1viscosity at 98.9°C)/(viscosity at
37.8°C). Using this parameter, you create the linear correlation for the dynamic viscosity
given by Equation 1.
2. G.L. Morini, “Viscous Heating in Liquid Flows in Micro-Channels,” Int. J. Heat Mass
Transfer, vol. 48, pp. 3637–3647, 2005.
Modeling Instructions
From the File menu, choose New.
NEW
In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D Axisymmetric.
2 In the Select Physics tree, select Heat Transfer>Conjugate Heat Transfer>Laminar Flow.
3 Click Add.
4 Click Study.
5 In the Select Study tree, select General Studies>Time Dependent.
6 Click Done.
GLOBAL DEFINITIONS
Parameters 1
1 In the Model Builder window, under Global Definitions click Parameters 1.
2 In the Settings window for Parameters, locate the Parameters section.
3 Click Load from File.
4 Browse to the model’s Application Libraries folder and double-click the file
fluid_damper_parameters.txt.
Piston Displacement
1 In the Home toolbar, click Functions and choose Global>Analytic.
Argument Unit
t s
8 Click Plot.
9 Click the Zoom Extents button in the Graphics toolbar.
10 In the Function name text field, type zp.
11 In the Label text field, type Piston Displacement.
DEFINITIONS
Variables 1
1 In the Home toolbar, click Variables and choose Local Variables.
2 In the Settings window for Variables, locate the Variables section.
3 Click Load from File.
4 Browse to the model’s Application Libraries folder and double-click the file
fluid_damper_variables.txt.
GEOMETRY 1
Rectangle 1 (r1)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type Dr/2.
4 In the Height text field, type 2*Ld.
5 Locate the Position section. In the z text field, type -Ld.
Rectangle 3 (r3)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type Dd/2-Hw.
4 In the Height text field, type 2*Ld.
5 Locate the Position section. In the z text field, type -Ld.
Rectangle 4 (r4)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type Dd/2.
4 In the Height text field, type 2*Ld.
5 Locate the Position section. In the z text field, type -Ld.
Rectangle 5 (r5)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type Dd/2.
4 In the Height text field, type 2*Lp.
5 Locate the Position section. In the z text field, type -Lp.
6 In the Geometry toolbar, click Build All.
7 Click the Zoom Extents button in the Graphics toolbar.
The model geometry is now complete.
Fluid 1
1 In the Model Builder window, under Component 1 (comp1)>
Heat Transfer in Solids and Fluids (ht) click Fluid 1.
2 Select Domains 4 and 6–9 only.
3 In the Settings window for Fluid, locate the Thermodynamics, Fluid section.
4 From the Fluid type list, choose Gas/Liquid.
5 From the list, choose User defined.
ADD MATERIAL
1 In the Home toolbar, click Add Material to open the Add Material window.
2 Go to the Add Material window.
3 In the tree, select Built-in>Steel AISI 4340.
4 Click Add to Component in the window toolbar.
5 In the Home toolbar, click Add Material to close the Add Material window.
MATERIALS
In the following steps, you create a new material for the damping fluid, Silicone Oil.
Silicone Oil
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, type Silicone Oil in the Label text field.
3 Select Domains 4 and 6–9 only.
4 In the Model Builder window, expand the Component 1 (comp1)>Materials>
Silicone Oil (mat2) node, then click Basic (def).
5 In the Settings window for Basic, locate the Model Inputs section.
6 Click Select Quantity.
7 In the Physical Quantity dialog box, type temperature in the text field.
8 Click Filter.
9 In the tree, select General>Temperature (K).
10 Click OK.
11 In the Settings window for Basic, locate the Local Properties section.
Initial Values 1
1 In the Model Builder window, under Component 1 (comp1)>
Heat Transfer in Solids and Fluids (ht) click Initial Values 1.
2 In the Settings window for Initial Values, locate the Initial Values section.
3 In the T text field, type T0.
Temperature 1
1 In the Physics toolbar, click Boundaries and choose Temperature.
2 Select Boundaries 2, 7, 9, 14, 16, 21, 23, and 28 only. These are the upper and lower
boundaries of the cylinder.
3 In the Settings window for Temperature, locate the Temperature section.
4 In the T0 text field, type T0.
Because the damper is a closed container, you need to pinpoint the pressure level within.
To achieve that, use the point constraint as follows.
Wall 2
1 In the Physics toolbar, click Boundaries and choose Wall.
2 Select Boundaries 22, 24, and 26 only.
3 In the Settings window for Wall, click to expand the Wall Movement section.
4 From the Translational velocity list, choose Manual.
COMPONENT 1 (COMP1)
Prescribed Deformation 1
1 In the Physics toolbar, click Moving Mesh and choose Prescribed Deformation.
2 In the Settings window for Prescribed Deformation, locate the Prescribed Deformation
section.
3 Specify the dx vector as
0 R
zp(t) Z
0 R
zlin1 Z
Prescribed Deformation 3
1 In the Moving Mesh toolbar, click Prescribed Deformation.
2 In the Settings window for Prescribed Deformation, locate the Prescribed Deformation
section.
3 Specify the dx vector as
0 R
zlin2 Z
MESH 1
Mapped 1
In the Mesh toolbar, click Mapped.
Distribution 1
1 Right-click Mapped 1 and choose Distribution.
2 Select Boundaries 23, 25, 27, and 28 only.
3 In the Settings window for Distribution, locate the Distribution section.
4 From the Distribution type list, choose Predefined.
5 In the Number of elements text field, type 4.
6 In the Element ratio text field, type 4.
7 From the Growth rate list, choose Exponential.
8 Select the Reverse direction check box.
Distribution 2
1 In the Model Builder window, right-click Mapped 1 and choose Distribution.
Distribution 3
1 Right-click Mapped 1 and choose Distribution.
2 Select Boundaries 9, 11, 13, and 14 only.
3 In the Settings window for Distribution, locate the Distribution section.
4 From the Distribution type list, choose Predefined.
5 In the Number of elements text field, type 30.
6 In the Element ratio text field, type 10.
7 From the Growth rate list, choose Exponential.
8 Select the Symmetric distribution check box.
Distribution 4
1 Right-click Mapped 1 and choose Distribution.
2 Select Boundaries 16, 18, 20, and 21 only.
3 In the Settings window for Distribution, locate the Distribution section.
4 In the Number of elements text field, type 8.
Distribution 5
1 Right-click Mapped 1 and choose Distribution.
2 Select Boundaries 3, 10, 17, 24, and 30 only.
3 In the Settings window for Distribution, locate the Distribution section.
4 In the Number of elements text field, type 32.
Distribution 6
1 Right-click Mapped 1 and choose Distribution.
2 Select Boundaries 2, 4, 6, and 7 only.
3 In the Model Builder window, right-click Mesh 1 and choose Build All.
The mesh is now complete. It should look similar to that shown in Figure 2.
STUDY 1
1 In the Model Builder window, under Study 1 click Step 1: Time Dependent.
2 In the Settings window for Time Dependent, locate the Study Settings section.
3 In the Output times text field, type range(0,tstep,(ncycle-1)/f)
range((ncycle-1)/f,tstep/2,tmax).
4 Click to expand the Results While Solving section. From the Update at list, choose
Times stored in output.
5 In the Home toolbar, click Compute.
RESULTS
Temperature (ht)
Change the unit of the temperature results to degrees Fahrenheit for the sake of easier
comparison with the experimental measurements.
Fluid Temperature
1 In the Model Builder window, expand the Temperature and Fluid Flow (nitf1) node, then
click Fluid Temperature.
2 In the Settings window for Surface, locate the Expression section.
3 From the Unit list, choose degF.
Solid Temperature
1 In the Model Builder window, click Solid Temperature.
2 In the Settings window for Surface, locate the Expression section.
3 From the Unit list, choose degF.
To produce a 3D plot of the temperature field and the flow streamlines within the damper,
use one of the predefined plots. The plot should appear similar to that shown in Figure 3.
Start with creating additional Cut Plane datasets for plotting the streamlines.
Cut Plane 1
1 In the Results toolbar, click Cut Plane.
2 In the Settings window for Cut Plane, locate the Plane Data section.
3 From the Plane type list, choose General.
4 From the Plane entry method list, choose Point and normal vector.
5 Find the Normal vector subsection. In the x text field, type 1.
6 In the z text field, type 0.
7 Click Plot to visualize the orientation of the cut plane.
8 Right-click Cut Plane 1 and choose Duplicate.
Cut Plane 2
1 In the Model Builder window, click Cut Plane 2.
2 In the Settings window for Cut Plane, locate the Plane Data section.
3 Find the Normal vector subsection. In the y text field, type 1.
4 Click Plot.
1 In the Results toolbar, click Add Predefined Plot to open the Add Predefined Plot
window.
2 Go to the Add Predefined Plot window.
3 In the tree, select Study 1/Solution 1 (sol1)>Heat Transfer in Solids and Fluids>
Temperature (ht).
4 Click Add Plot in the window toolbar.
5 In the Results toolbar, click Add Predefined Plot to close the Add Predefined Plot
window.
RESULTS
Streamline Surface 1
In the Temperature and Velocity Streamlines toolbar, click More Plots and choose
Streamline Surface.
Volume 1
1 In the Model Builder window, click Volume 1.
2 In the Settings window for Volume, locate the Expression section.
3 From the Unit list, choose degF.
Streamline Surface 1
1 In the Model Builder window, click Streamline Surface 1.
2 In the Settings window for Streamline Surface, locate the Data section.
3 From the Dataset list, choose Cut Plane 1.
4 From the Solution parameters list, choose From parent.
5 Click Replace Expression in the upper-right corner of the Expression section. From the
menu, choose Component 1 (comp1)>Laminar Flow>Velocity and pressure>u,v,w -
Velocity field (spatial frame).
6 Locate the Streamline Positioning section. From the Positioning list, choose
Uniform density.
7 In the Separating distance text field, type 0.025.
Streamline Surface 2
1 In the Model Builder window, click Streamline Surface 2.
2 In the Settings window for Streamline Surface, locate the Data section.
3 From the Dataset list, choose Cut Plane 2.
4 Click to expand the Title section. From the Title type list, choose None.
5 In the Temperature and Velocity Streamlines toolbar, click Plot.
6 Click the Zoom Extents button in the Graphics toolbar.
3 Locate the Legend section. Clear the Show legends check box.
This plot shows the temperature variation at the probe position over the complete
loading time period, it should look similar to that shown in Figure 4.
1 In the Home toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, type Temperature Along Inner Wall in
the Label text field.
Line Graph 1
1 Right-click Temperature Along Inner Wall and choose Line Graph.
2 In the Settings window for Line Graph, locate the y-Axis Data section.
3 In the Expression text field, type z/U0.
4 Locate the x-Axis Data section. From the Parameter list, choose Expression.
5 From the Unit list, choose degF.
6 Select Boundaries 22, 24, and 26 only.
7 Click to expand the Legends section. Select the Show legends check box.
Surface Average 1
Finally, plot the average velocity within the damper over the last cycle.
8 Click Evaluate.