FLUENT uses solid-body kinematics if the motion is prescribed based on the position and orientation of the center of gravity of a moving object. This is applicable to both cell and face zones.
The motion of the solid-body can be specified either as a profile or as a user-defined function (UDF). A profile may be defined by the following profile fields:
By default FLUENT assumes that the motion is specified in the inertial coordinate system. However, it is also possible to prescribe the motion relative to the coordinate system by selecting the Motion Type from the drop-down list for the respective fluid zone defined in Define Boundary Conditions.
Thus the motion may be prescribed relative to a Rotating Reference Frame, Moving Mesh, or as Stationary.
For in-cylinder simulations, the velocity profiles for valves can be expressed as a function of crank angle instead of time. In addition, transient boundary condition profiles can also be expressed as a function of crank angle instead of time. For more information about transient profiles, see Section 7.1.9.
Below are two examples of a profile format:
((movement_linear 3 point) (time 0 1 2 ) (x 2 3 4 ) (v_y 0 -5 0 ) )
((movement_angular 3 point) (time 0 1 2 ) (omega_x 2 3 4 ) )
For in-cylinder flows, crank angles can be included in transient tables as well as transient profiles, in a similar fashion to time. An example of a transient table using (crank) angle is as follows:
example 2 3 1 angle temperature 0 300 180 500 360 300
An example of a transient profile using (crank) angle is as follows:
((example transient 3 1) (angle 0.000000e+00 1.800000e+02 3.600000e+02) (temperature 3.000000e+02 5.000000e+02 3.000000e+02) )
In addition to the motion description, you must also specify the starting location of the center of gravity and orientation of the solid body. In 2D (and 3D non-6DOF), FLUENT automatically updates the center of gravity position and orientation at every time step such that
where and are the position and orientation of the center of gravity, and are the linear and angular velocities of the center of gravity. 3D, 6DOF cases use a more complex form of Equation 11.3-18 when updating .
Typically, is chosen to be an appropriate set of Euler angles. In this case, the solid-body motion must be specified using a user-defined function ( DEFINE_CG_MOTION).
The position vectors on the solid body are updated based on rotation about the instantaneous angular velocity vector . For a finite rotation angle = , the final position of a vector on the solid body with respect to can be expressed as (See Figure 11.3.21)
where can be shown to be
The unit vectors and are defined as
If the solid body is also translating with , the position vector on the solid body can be expressed as
where is given by Equation 11.3-19.