Force and Moment Coefficients#
Besides the non-dimensional output fields, there are also many coefficients commonly used in the computational fluid dynamics community, e.g. lift coefficient (\(C_L\)), drag coefficient (\(C_D\)), etc. Flow360 also exports the above coefficients as tabulated data.
Attention
The reference velocity for the force and moment coefficients is \(U_\text{ref}\), which is the reference velocity set in the operating condition and can be accessed through case.params.reference_velocity.
It is important to note that \(U_\text{ref}\) used here is different from the reference velocity scaling \(U_\text{scale}\) used for nondimensionalizing other variables (such as velocity fields, heat flux, angular speeds).
We want to emphasize that the non-dimensionalization of force and moment for BET Disk , Actuator Disk and Porous Media are different from the coefficients shown in this page.
The force coefficients and moment coefficients exported by Flow360 are listed in Table 3.8.5.
These coefficients can be obtained from total_forces_v2.csv, surface_forces_v2.csv files on WebUI. They can also be fetched by the case.results.total_forces and case.results.surface_forces Python API calls.
Property |
Definition |
|---|---|
CL |
\(\text{Lift}/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref}\) |
CD |
\(\text{Drag}/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref}\) |
CFx |
\(\text{Force}_x/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref}\) |
CFy |
\(\text{Force}_y/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref}\) |
CFz |
\(\text{Force}_z/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref}\) |
CMx |
\(\text{Moment}_x/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref} L_\text{ref}\left[0\right]\) |
CMy |
\(\text{Moment}_y/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref} L_\text{ref}\left[1\right]\) |
CMz |
\(\text{Moment}_z/\frac{1}{2}\rho_\infty U_\text{ref}^2 A_\text{ref} L_\text{ref}\left[2\right]\) |
Note
In the above table, all reference values can be accessed through Python API as shown in the Reference Variable Table.
The typical output options in the WebUI and various csv files are outlined in Table 3.8.6. The conventions assume z-axis upwards, y-axis spanwise (+ towards starboard side) and x-axis in the axial direction (+ in the freestream direction) for the global axes, as shown in Fig. 3.8.10.
Fig. 3.8.10 Axis conventions demonstrated using the CRM geometry#
Output |
|
|---|---|
\(CFx, CFy, CFz\) |
Force coefficients (global axes) |
\(CMx, CMy, CMz\) |
Moment coefficients (global axes) |
\(CL = CFz\cdot cos(\alpha) - CFx\cdot sin(\alpha)\) |
Lift coefficient (wind axes) |
\(CD = CFx\cdot cos(\alpha) cos(\beta) - CFy\cdot sin(\beta) +\) \(CFz\cdot sin(\alpha) cos(\beta)\) |
Drag coefficient (wind axes) |
\(CFxPressure, CFyPressure, CFzPressure\) |
Pressure contributions to force coefficients (global axes) |
\(CMxPressure, CMyPressure, CMzPressure\) |
Pressure contributions to moment coefficients (global axes) |
\(CFxSkinFriction, CFySkinFriction,\) \(CFzSkinFriction\) |
Skin friction contributions to force coefficients (global axes) |
\(CMxSkinFriction, CMySkinFriction,\) \(CMzSkinFriction\) |
Skin friction contributions to moment coefficients (global axes) |
\(CLPressure, CDPressure\) |
Pressure contributions to the lift and drag coefficients |
\(CLSkinFriction, CDSkinFriction\) |
Skin friction contributions to the lift and drag coefficients |
\(CS = CFx\cdot cos(\alpha) sin(\beta) + CFy\cdot cos(\beta) +\) \(CFz\cdot sin(\alpha) sin(\beta)\) |
Side force coefficient (wind axes). Not available in the csv file and can be calculated from body forces. |
Example: Lift Force and Pitching Moment#
1# The output values are averaged over the last 10% steps
2total_forces = case.results.total_forces.averages
3
4density = case.params.operating_condition.thermal_state.density
5A_ref = case.params.reference_geometry.area
6U_ref = case.params.reference_velocity
7moment_length = case.params.reference_geometry.moment_length
8
9# Compute Lift Force
10CL = total_forces["CL"]
11Lift = CL * 0.5 * density * U_ref**2 * A_ref
12
13# Compute Pitching Moment
14CMy = total_forces["CMy"]
15Pitching_Moment = CMy * 0.5 * density * U_ref**2 * A_ref * moment_length[1]
See also
Calculate Dimensional Forces — A comprehensive Python snippet for converting force and moment coefficients to dimensional values.