Rocky Module Manual
Rocky Module Manual
ANSYS, Ansys Workbench, AUTODYN, CFX, FLUENT and any and all ANSYS, Inc. brand, product, service and feature
names, logos and slogans are registered trademarks or trademarks of ANSYS, Inc. or its subsidiaries located in the
United States or other countries. ICEM CFD is a trademark used by ANSYS, Inc. under license. CFX is a trademark
of Sony Corporation in Japan. All other brand, product, service and feature names or trademarks are the property
of their respective owners. FLEXlm and FLEXnet are trademarks of Flexera Software LLC.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFID-
ENTIAL AND PROPRIETARY PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products
and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates under a software license agreement
that contains provisions concerning non-disclosure, copying, length and nature of use, compliance with exporting
laws, warranties, disclaimers, limitations of liability, and remedies, and other provisions. The software products
and documentation may be used, disclosed, transferred, or copied only in accordance with the terms and conditions
of that software license agreement.
ANSYS, Inc. and ANSYS Europe, Ltd. are UL registered ISO 9001: 2015 companies.
For U.S. Government users, except as specifically granted by the ANSYS, Inc. software license agreement, the use,
duplication, or disclosure by the United States Government is subject to restrictions stated in the ANSYS, Inc.
software license agreement and FAR 12.212 (for non-DOD licenses).
Third-Party Software
See the legal information in the product help files for the complete Legal Notice for ANSYS proprietary software
and third-party software. If you are unable to access the Legal Notice, contact ANSYS, Inc.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. iii
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
iv of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. v
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
vi of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. vii
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
viii of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. ix
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
x of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xi
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xii of ANSYS, Inc. and its subsidiaries and affiliates.
List of Figures
1.1. Rocky tab in Adams. ............................................................................................................................. 21
1.2. Rocky Coupling dialog in Adams. ......................................................................................................... 22
1.3. Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module. ................................................. 23
1.4. Standard Adams Coupling Curves for a geometry that was exported out of Adams as an Output. ............ 25
1.5. Additional Adams Coupling Curves for a geometry that was exported out of Adams as an Input. ............. 26
2.1. Standard Rocky Motion Frames vs. Ansys Motion Coupling Cabilities. ..................................................... 29
2.2. Ansys Motion Coupling FMU file lifecycle. .............................................................................................. 30
2.3. Co-Simulator tab in Ansys Motion. ......................................................................................................... 32
2.4. Rocky Coupling dialog in Ansys Motion. ................................................................................................ 33
2.5. Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module. ................................................. 34
2.6. Standard Ansys Motion Coupling Curves for a geometry that was exported out of Ansys Motion as an
Output. ...................................................................................................................................................... 37
2.7. Additional Ansys Motion Coupling Curves for a geometry that was exported out of Ansys Motion as an
Input. ......................................................................................................................................................... 38
2.8. "Flexible Bodies" message in Ansys Motion. ........................................................................................... 40
2.9. Ansys License Manager Error in Rocky. ................................................................................................... 42
3.1. Example application of a drilling process using the Bond Model module. ............................................... 43
3.2. Physics | Momentum tab options when the Bond Model module is enabled. ............................................. 44
3.3. Contacts tab options. ............................................................................................................................ 45
3.4. Bond Model module options. ................................................................................................................. 46
3.5. Bond Model module materials interactions options. ............................................................................... 48
3.6. Bond Model module inputs options. ..................................................................................................... 50
3.7. Additional Contact Properties created by the Bond Model module. ......................................................... 51
3.8. Additional Particle Properties created by the Bond Model module. .......................................................... 52
3.9. Geometry of the bond between two spherical particles of different sizes. ............................................... 54
3.10. Forces and moments exerted by a bond on a particle. .......................................................................... 56
4.1. Boundary Outlet module options for an individual geometry component that you have imported. .......... 60
5.1. Particle combustion over time using the CFD-Coupled Particle Combustion module. ............................. 64
5.2. CFD-Coupled Particle Combustion module options under the Species tab. ............................................. 65
5.3. CFD-Coupled Particle Combustion module under the Reactions Tab. ..................................................... 66
5.4. CFD-Coupled Particle Combustion module options under the Advanced Tab. ......................................... 66
5.5. Options for a Continuous Injection input ............................................................................................... 69
5.6. Options for a Volumetric Inlet input. ...................................................................................................... 70
5.7. Options for a Custom Inlet input. ........................................................................................................... 71
5.8. Particles properties when the CFD-Coupled Particle Combustion is enabled. .......................................... 73
5.9. Setting the mixture species in Fluent. .................................................................................................... 74
6.1. Particle drying rate over time as a function of the moisture content. ...................................................... 79
6.2. CFD-Coupled Particle Drying module input parameters. .......................................................................... 81
6.3. Options for a Particle Inlet Input (left) and a Volumetric Inlet Input (right). ................................................ 83
6.4. Options for a Particle Custom Inlet. ......................................................................................................... 85
6.5. Particles Properties when the CFD-Coupled Particle Drying module is enabled. ........................................ 87
6.6. Fluent species model configuration. ..................................................................................................... 89
6.7. Fluent species mixture-template configuration. .................................................................................... 89
6.8. Particle Drying - problem description. ................................................................................................... 90
6.9. Damping factor as a function of the moisture content for different damping function B coefficients with
nc = nv = 1.0. .............................................................................................................................................. 92
6.10. Damping factor as a function of the moisture content for different damping function nc coefficients
with B = 0.501. ............................................................................................................................................ 93
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xiii
Rocky Module Manual
6.11. Damping factor as a function of the moisture content for different damping function nv coefficients
with B = 0.501. ............................................................................................................................................ 94
7.1. Example application of grains packing using the Coordination Number module. .................................... 95
7.2. Example application of particle mixing using the Coordination Number module. ..................................... 96
7.3. Coordination Number module options. .................................................................................................. 97
7.4. Coordination Number module options. .................................................................................................. 98
8.1. Interactions tab options for the 1-Way Constant CFD Coupling mode when the Dhole, Chhabra & Eswaran
Drag Law module is enabled. ..................................................................................................................... 102
8.2. Dhole, Chhabra & Eswaran Drag Law module options. ........................................................................... 103
9.1. Snow accretion 1-way Rocky- Fluent coupling simulation using the Disable Collision Outside Region
module. Blue particles are enabled whereas light pink particles are removed from contact detection calcula-
tions. ........................................................................................................................................................ 108
9.2. Disable Collision Outside Region module options. .................................................................................. 109
9.3. Additional Particle Properties created by the Disable Collision Outside Region module. .......................... 110
10.1. Electromagnetic Coupling Module options under the Model tab. .......................................................... 114
10.2. Electromagnetic Coupling Module options under the Fields tab. ........................................................... 115
10.3. Electromagnetic Coupling Module options under the Advanced tab. .................................................... 115
10.4. Electromagnetic Coupling Module options under Materials. ................................................................. 118
10.5. Electromagnetic Coupling Module options under Material Interactions. ................................................ 119
10.6. Electromagnetic Coupling Module options under Geometries. ............................................................. 120
10.7. Electromagnetic Coupling Module options under Particle Groups. ........................................................ 121
10.8. Electromagnetic Coupling Module options under Particle Inlets and Volumetric Inlets. .......................... 122
10.9. Electromagnetic Coupling Module flowchart. ....................................................................................... 125
11.1. Eroded Mass Analyses module options under tab. ............................................................................... 132
11.2. Eroded Mass Analyses module options under tab. ............................................................................... 133
11.3. Eroded Mass Analyses module options under tab. ............................................................................... 133
11.4. Additional Geometry Properties created by the Eroded Mass Analyses module. ................................... 136
11.5. Additional Geometry Curves created by the Eroded Mass Analyses module. ......................................... 137
12.1. Options for the Frictional Heating module. ......................................................................................... 142
12.2. Frictional Heating module options. ..................................................................................................... 143
13.1. Options for the Heat Generation Region module. ................................................................................ 146
14.1. Tag column on the custom fiber definition file. .................................................................................. 149
14.2. Pairs of Tags table .............................................................................................................................. 150
14.3. Tags mapped to joints ....................................................................................................................... 150
14.4. Elastic Ratio and Joint Damping Ratio. ............................................................................................... 151
14.5. Element Damping Ratio. .................................................................................................................... 152
14.6. Options for the Heterogeneous Custom Fiber (Beta) module. ................................................................ 152
14.7. Options for the Heterogeneous Custom Fiber (Beta) module. ................................................................ 153
14.8. —Additional particle properties created by the Heterogeneous Custom Fiber (Beta) module. ................ 154
15.1. Options for the Leblicq Joint Model module. ....................................................................................... 157
15.2. Example paths to stiffness input files. ................................................................................................ 158
15.3. Joint Model options for a multi-element Fiber when the Leblicq Joint Model module is enabled. .......... 159
15.4. Example Loading Stiffness table. ....................................................................................................... 161
15.5. Example Unloading Stiffness table. .................................................................................................... 161
15.6. Folder location of the example .csv files. ............................................................................................ 162
16.1. Physics | Momentum tab options when the Leeds Contact Model module is enabled. .......................... 166
16.2. Options for the Leeds Contact Model module. ..................................................................................... 167
16.3. Materials Interactions options when the Leeds Contact Model module is enabled. ............................... 168
16.4. Schematic diagram of a typical normal force overlap response of the Leeds Contact Model. ............... 170
16.5. Normal force in Leeds Contact Model as a composition of repulsive and adhesive forces. .................... 172
17.1. Options for the Liquid Absorption module. ......................................................................................... 174
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xiv of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
17.2. Liquid Absorption module options for a Continuous Injection Input. ..................................................... 175
17.3. Particles Properties created by the Liquid Absorption module. ............................................................. 177
18.1. Example mixing process using the Liquid Bridge Model ..................................................................... 179
18.2. Liquid Bridge Model module options. .................................................................................................. 180
18.3. Liquid Bridge Model overriding the Adhesive Force setting. ................................................................... 181
18.4. Liquid Bridge Model Materials Interactions options. ............................................................................. 183
18.5. Options for a Volumetric Inlet (left) and a Particle Inlet (right). ............................................................ 184
18.6. Options for a Particle Custom Inlet. ..................................................................................................... 185
18.7. Options for a Particle Custom Inlet. ..................................................................................................... 186
18.8. Geometry of a liquid bridge between two particles. ........................................................................... 187
18.9. Geometry of a liquid bridge between a particle and a boundary. ........................................................ 188
19.1. Options for the Liquid Drying by Heat Flux module. ............................................................................. 194
19.2. Options for a Continuous Injection Input (left) and a Volume Fill Input (right). ..................................... 195
19.3. Options for a Custom Input. ................................................................................................................ 196
19.4. Particle Properties added by the Liquid Drying by Heat Flux module. .................................................... 198
20.1. Options for the Liquid Drying Region module. ..................................................................................... 202
20.2. Options for a Continuous Injection Input (on the left) and a Volume Fill Inputs (on the right). ............... 204
20.3. Options for a Custom Input. ................................................................................................................ 206
20.4. Particle Properties added by the Liquid Drying Region module. ............................................................ 208
21.1. Defining a Point Cloud in Rocky. ......................................................................................................... 213
21.2. Example text file defining a Point Cloud. ............................................................................................. 214
21.3. Example spheres of influence formed around the points defined in a Point Cloud. ............................... 214
21.4. Liquid Injection Sphere Clouds module options. ................................................................................... 215
21.5. Liquid Injection Sphere Clouds module options for a Continuous Injection Input. ................................... 217
21.6. Particles property created by Liquid Injection Sphere Clouds module. ................................................... 218
22.1. Options for the Liquid Injection Region module. .................................................................................. 222
22.2. Options for a Continuous Injection Input (on the left) and a Volume Fill Input (on the right). ............... 223
22.3. Options for a Custom Input. ................................................................................................................ 225
22.4. Particle Properties added by the Liquid Injection Region module. ......................................................... 227
23.1. Physics | Momentum tab options when the Thornton-Ning Normal Force module is enabled. ............. 229
23.2. Options for the Thornton-Ning Normal Force module. ........................................................................ 230
23.3. Materials Interactions options when the Thornton-Ning Normal Force module is enabled. .................. 231
24.1. Physics | Momentum tab options when the Nonlinear Hysteretic Normal Force module is enabled. ......... 236
24.2. Options on the Nonlinear Hysteretic Normal Force module. .................................................................. 237
24.3. Materials Interactions options when the Nonlinear Hysteretic Normal Force module is enabled. ............. 238
25.1. Example application of the SPH Non-Newtonian (Beta) module in a food modelling scenario. ............. 242
25.2. SPH Non-Newtonian Models (Beta) module options. ............................................................................ 243
25.3. Non-newtonian models flow curves. .................................................................................................. 244
26.1. Physics | Thermal tab options when the Nonlinear Thermal Contact Model module is enabled. .............. 248
26.2. Nonlinear Thermal Contact Model module options. ............................................................................. 248
26.3. Nonlinear Thermal Contact Model module options for an individual material component. .................... 250
26.4. Illustration of a piece-wise scheme. .................................................................................................... 252
27.1. Inter-group Collisions Statistics module Curves upon which the Normal and Tangential Energy Dissipation
module is based. ....................................................................................................................................... 255
27.2. Inter-particle Collision Statistics module Property upon which the Normal and Tangential Energy Dis-
sipation module is based. .......................................................................................................................... 256
27.3. Normal and Tangential Energy Dissipation module options. ................................................................. 257
27.4. Particles Curves created by the Normal and Tangential Energy Dissipation module. .............................. 258
27.5. Particles properties created by the Normal and Tangential Energy Dissipation module. ......................... 259
27.6. Normal Force and Tangential Force Contact Models. .......................................................................... 261
28.1. Defining a Point Cloud in Rocky. ......................................................................................................... 263
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xv
Rocky Module Manual
28.2. Example text file defining a fluid velocity field point cloud. ................................................................ 264
28.3. Options for the 1-Way CFD Field Coupling module. .............................................................................. 265
29.1.Breakage | Criteria tab options (left) and Breakage | Fragments tab options (right) for the Particle set when
the Overlap-Based Breakage Model module is enabled. ............................................................................... 267
29.2. Overlap-Based Breakage Model module options. ................................................................................. 268
29.3. Overlap-Based Breakage Model module options for a Particle set. ........................................................ 270
29.4. Particles Properties created by the Overlap-Based Breakage Model module. ......................................... 272
30.1.Interactions tab options for the 2-Way Fluent Unresolved Coupling mode when the Parametric Gidaspow,
Bezburuah & Ding module is enabled. ........................................................................................................ 278
30.2. Parametric Gidaspow, Bezburuah & Ding module options. .................................................................... 279
31.1.Interactions tab options for the 1-Way Constant CFD Coupling mode when the Parametric Ranz & Marshall
Law module is enabled. ............................................................................................................................. 284
31.2. Parametric Ranz & Marshall Law module options. ................................................................................ 285
32.1. Disk harrow simulation using the Particle Freeze Outside Region module. Blue particles are enabled
whereas light pink particles are frozen and removed from contact detection calculations. .......................... 289
32.2. Particle Freeze Outside Region module options. ................................................................................... 291
32.3. Additional Particle properties created by the Particle Freeze Outside Region module. ........................... 292
33.1. Filling of a silo using the Particle Freeze when Settled module. Grey particles are frozen whereas blue
particles are active. ................................................................................................................................... 296
33.2. Particle Freeze when Settled module options. ...................................................................................... 298
33.3. Additional Particle Properties created by the Particle Freeze when Settled module. ............................... 299
34.1. Defining a Point Cloud in Rocky. ......................................................................................................... 304
34.2. Example text file defining a fluid velocity field point cloud. ................................................................ 304
34.3. Example application of the SPH Point Cloud Air Drag (Beta) module in a gearbox. .............................. 305
34.4. SPH Point Cloud Air Drag (Beta) module options. ................................................................................. 306
34.5. Additional property enabled on the SPH entity when the SPH Point Cloud Air Drag (Beta) is used. ........ 308
34.6. Flowchart illustrating the steps involved in computing the drag force for each SPH element. .............. 310
34.7. Deformation of a droplet from the spherical shape with radius L to the disk-like shape with radius
......................................................................................................................................................... 312
34.8. Occlusion factor computation accounting for the angle between the direction of fluid flow, vi,?rel, and
the neighbor direction xi,?j. ....................................................................................................................... 313
35.1. Options for the Region Outlet module. ............................................................................................... 316
36.1. Example application of the SPH Air Drag (Beta) module in a fire hose jet scenario. .............................. 320
36.2. SPH Air Drag (Beta) module options. ................................................................................................... 321
36.3. Additional property enabled on the SPH entity when the SPH Air Drag (Beta) is used. .......................... 322
36.4. Flowchart illustrating the steps involved in computing the drag force for each SPH element. .............. 324
36.5. Deformation of a droplet from the spherical shape with radius L to the disk-like shape with radius
......................................................................................................................................................... 326
36.6. Occlusion factor computation accounting for the angle between the direction of fluid flow, vi,?rel, and
the neighbor direction xi,?j. ....................................................................................................................... 327
37.1. Simplified wheel simulation using the SPH Freeze Outside Region module. Blue SPH elements are enabled
whereas light red elements are frozen and removed from contact detection calculations. .......................... 329
37.2. Cube creation inside Region of Interest setup. ....................................................................................... 330
37.3. Cylinder creation inside Region of Interest setup. .................................................................................. 331
37.4. SPH Freeze Outside Region module options. ........................................................................................ 332
37.5. Additional SPH properties created by the SPH Freeze Outside Region module. ...................................... 333
38.1. Stress-Based Multi-Contact Normal Force module options. ................................................................. 336
38.2. Stress-Based Multi-Contact Normal Force module options for a particle group. ................................... 337
38.3. Illustration of a contact pair interaction with the average stress terms of particle i (green) and particle
j (red) from the neighboring particles contacts. ......................................................................................... 338
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xvi of ANSYS, Inc. and its subsidiaries and affiliates.
Rocky Module Manual
39.1.Breakage | Criteria tab options (left) and Breakage | Fragments tab options (right) for the Particle set when
the TBM Breakage Model module is enabled. .............................................................................................. 342
39.2. TBM Breakage Model module options. ................................................................................................ 343
39.3. TBM Breakage Model module options for a Particle set. ....................................................................... 344
39.4. Particles Property created by the TBM Breakage Model module. .......................................................... 346
39.5. Breakage specific energy versus particle size. ..................................................................................... 348
40.1. Heating particles inside a cylinder using the Thermal Radiation module.The upper surface is set to a
high temperature, which heats the particles at the top.The particles then exchange heat with each other in
the lower layers. ........................................................................................................................................ 351
40.2. Options for the Thermal Radiation module. ........................................................................................ 353
40.3. Thermal Radiation module options for an individual material component. .......................................... 355
40.4. Illustration of the calculation domain for radiation heat transfer according to the defined screening
distances. ................................................................................................................................................. 358
40.5. Particle-particle interaction. .............................................................................................................. 360
40.6. Triangle-particle interaction. .............................................................................................................. 361
41.1. Physics | Momentum tab options when the Velocity-Dependent Adhesive Force module is enabled. ....... 364
41.2. Velocity-Dependent Adhesive Force module options. ............................................................................ 364
41.3. Velocity-Dependent Adhesive Force module options for a pair of materials interactions. ........................ 365
41.4. Contacts Properties created by the Velocity-Dependent Adhesive Force module. ................................... 366
42.1. Velocity-Dependent Restitution Model module options. ........................................................................ 369
42.2. Velocity-Dependent Restitution Model module options for a pair of materials interactions. .................... 371
42.3. Contacts Property created by the Velocity-Dependent Restitution Model module. ................................. 373
42.4. Coefficient of restitution as a function of impact velocity, for different values of the exponent K. ......... 375
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xvii
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xviii of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 1: Adams Coupling
1.1. Introduction
Adams Coupling enables users to prescribe specific, complex motions for their equipment in Adams and
then have those same motions available to be used in a DEM simulation that is carried out by Rocky.
When compared to the standard method of defining Motion Frames in Rocky, movements defined with
Adams Coupling have the added benefit of simulating geometry-to-geometry interactions, stress and
deformation of geometry components, chained motions and more.
It is not a requirement for both Adams and Rocky to be opened at the same time. Adams needs a license
for each instance of it running. Therefore, in order to have an open instance of Adams while a coupled
simulation is running, more than one license is needed. However, in order to make use of the Adams
solver in Rocky, you must still have a valid Adams license and have it installed in the same machine
that is running the Rocky simulation.
The Adams Coupling module allows for each portion of the coupling setup to be accomplished by dif-
ferent users and have only the FMU file shared between them. Therefore, separate setup instructions
for each product will be provided in this document; follow one or both sets of instructions as needed.
– You have installed the Adams Coupling Module plugin. For information about installation, refer
to the Adams Coupling Installation Guide.
– You have a valid license for your version of Adams and you have it installed on the same machine
as Rocky.
1.1.2. Limitations
This version of Adams Coupling has the following limitations:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 19
Adams Coupling
• Gravity direction must match in both Adams and Rocky. Make sure that the gravity acceleration is
set in Adams in all directions (x, y, z). If a direction has no gravity force, make sure that the field is
filled with 0 (zero).
• The Simulation Duration in Adams must be greater than or equal to the one set in Rocky.
In addition, there are specific limitations for each product, as described below.
• Flexible bodies are supported on the Adams side of the coupled simulation; however, Rocky will
be unable to calculate particle interactions on those flexible bodies.
• Including flexible bodies in your coupled simulation might also increase the processing time.
• Motion Frames : A Rocky project with Adams Coupling cannot also have Motion Frames. Your
simulation will not start processing until all Motion Frames are removed.
2. Open an Existing Model or create a New Model, setting up your equipment and motions as you
normally would. Important:.
• Ensure that the Gravity direction in Adams is set to what you later want to use in Rocky. These
must match for the coupling process to work. Check also if the gravity acceleration is filled in
all directions (with a value set to zero if there is no gravity in that direction).
• Ensure that the simulation time (T) in Adams is equal or greater than the Simulation Duration
you later want to use in Rocky.
3. Ensure that the model that you want to export is correctly selected in the drop down menu.
4. From the Rocky tab select Generate FMU for Rocky Coupling (see Figure 1.1:Rocky tab in
Adams. (p. 21)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
20 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup procedures
5. From the Select parts for Rocky coupling dialog (Figure 1.2:Rocky Coupling dialog in
Adams. (p. 22)), do the following:
a. For each of the rigid Body parts without beam elements, do one or both of the following:
• Add to the Output box the parts you want motions to be shared with Rocky.
• Add to the Input box the parts you want forces and moments to be collected due to inter-
actions with particles. This data will then be used by the Adams solver during processing.
Tip: to add the parts to a box, right click an empty area of it. Then point to Part, then
Browse... and then select the parts to be added to the box.
b. Click the Export button. If no errors are reported, the operations were successful.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 21
Adams Coupling
At this point, the FMU file has been generated and is stored in the Working Directory defined in
Adams. It can now be shared with a Rocky user for them to complete the second part of the coupling
procedure, or you can continue with the next step yourself.
2. Begin a new simulation project. (From the File menu, click New Project.)
3. From the Data panel, under Modules, select Multibody Dynamics FMU Coupling.
4. From the Data Editors panel, on the Multibody Dynamics FMU Coupling tab, do the following
(Figure 2.5: Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module. (p. 34)):
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
22 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup procedures
Figure 1.3: Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module.
a. Click the FMU filename: button, and then from the Select file to import dialog, locate and
select the FMU file you want to use, and then click Open. Save the project now if you are
asked to do so.
b. From the Communication step size control list, select one of the following options:
• Automatic: For each Adams time step, Rocky computes 100 timesteps.
• Manual: Adams' time step is defined by the Communication step size value.
c. If Manual was selected for Communication step size control, then in the Manual commu-
nication step size field, enter what time step you want Adams to use.
5. From the Data panel, right-click Geometries, and then click Import Geometry.
6. From the Select file to import dialog, locate and select the CAD files of the geometry to be sim-
ulated.
Important:
The geometry imported in Rocky should be exactly the same used to export the
FMU from Adams. Therefore it is recommended that you use the same CAD files
to setup both Adams and Rocky simulations. If the geometry was created inside
Adams or you do not have access to the CAD files that were used to setup the
Adams simulation, the geometry can be exported from Adams and then imported
in Rocky. A CAD software may be required to convert CAD files from the formats
exported by Adams to formats that are readable in Rocky.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 23
Adams Coupling
7. From the Import File Info dialog, select the options you want, and then click OK. The Geometries
list is populated with the parts that were imported.
8. Set up the rest of your simulation as you normally would, being careful to ensure the following:
• You have defined the gravity direction in Rocky to match what was set in the Adams.
• You have defined the Simulation Duration in Rocky to be less than or equal to the duration
that was set in the Adams.
9. Process your simulation and analyze the results in Rocky as you normally would.
1. With the Rocky project open that contains information from the FMU file you want to update, do
the following:
a. Remove and then re-import the geometries added by the FMU file by doing the following:
i. From the Data panel, under Geometries, select or multi-select all the parts that were ori-
ginally imported by the FMU file.
ii. Right-click the selected geometries, and then click Remove Geometries.
Alternatively, if you made only small changes to the bodies in Adams, you can make those
changes to individual Geometries in Rocky without removing all of them and importing them
again. However, you must ensure that the changes you make exactly match the information
in the FMU file you import.
2. Set up, process, and analyze the results of your simulation as you normally would.
1.3. Post-processing
After processing your coupled simulation, you will have new Curves for your geometry components,
depending upon whether they were selected to be exported out of Adams as an Output body (Fig-
ure 1.4: Standard Adams Coupling Curves for a geometry that was exported out of Adams as an Out-
put. (p. 25)) – which has no interactions with Particles – or an Input body (Figure 1.5: Additional Adams
Coupling Curves for a geometry that was exported out of Adams as an Input. (p. 26)) – which does have
interactions with Particles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
24 of ANSYS, Inc. and its subsidiaries and affiliates.
Post-processing
Figure 1.4: Standard Adams Coupling Curves for a geometry that was exported out of Adams as
an Output.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 25
Adams Coupling
Figure 1.5: Additional Adams Coupling Curves for a geometry that was exported out of Adams as
an Input.
• Acceleration : Rotational : X: The amount of rotational (angular) acceleration for the geometry in
the X direction as measured from the Center of Mass.
• Acceleration : Rotational : Y: The amount of rotational (angular) acceleration for the geometry in
the Y direction as measured from the Center of Mass.
• Acceleration : Rotational : Z: The amount of rotational (angular) acceleration for the geometry in
the Z direction as measured from the Center of Mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
26 of ANSYS, Inc. and its subsidiaries and affiliates.
Post-processing
• Acceleration : Translational : X: The amount of translational acceleration for the geometry in the X
direction as measured from the Center of Mass.
• Acceleration : Translational : Y: The amount of translational acceleration for the geometry in the Y
direction as measured from the Center of Mass.
• Acceleration : Translational : Z: The amount of translational acceleration for the geometry in the Z
direction as measured from the Center of Mass.
• Force : X : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the force applied to the geometry along the X axis.
• Force : Y : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the force applied to the geometry along the Y axis.
• Force : Z : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the force applied to the geometry along the Z axis.
• Moment : X : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the moment (torque) applied to the geometry’s Center of Mass around the X axis.
• Moment : Y : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the moment (torque) applied to the geometry’s Center of Mass along the Y axis.
• Moment : Z : Instantaneous: Only for geometries that were exported out of Adams as an Input body,
this provides the moment (torque) applied to the geometry’s Center of Mass along the Z axis.
• Orientation Angle: The angle the geometry is rotated around the Orientation Vector.
• Orientation Vector X: The X component defining the vector around which the geometry is rotated.
• Orientation Vector Y: The Y component defining the vector around which the geometry is rotated.
• Orientation Vector Z: The Z component defining the vector around which the geometry is rotated.
• Position X: The global coordinate of the geometry’s Center of Mass as specified in the X direction.
• Position Y : The global coordinate of the geometry’s Center of Mass as specified in the Y direction.
• Position Z : The global coordinate of the geometry’s Center of Mass as specified in the Z direction.
• Velocity : Rotational : X : The amount of rotational (angular) velocity for the geometry in the X dir-
ection as measured from the Center of Mass.
• Velocity : Rotational : Y: The amount of rotational (angular) velocity for the geometry in the Y direction
as measured from the Center of Mass.
• Velocity : Rotational : Z: The amount of rotational (angular) velocity for the geometry in the Z direction
as measured from the Center of Mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 27
Adams Coupling
• Velocity : Translational : X: The amount of translational velocity for the geometry in the X direction
as measured from the Center of Mass.
• Velocity : Translational : Y: The amount of translational velocity for the geometry in the Y direction
as measured from the Center of Mass.
• Velocity : Translational : Z: The amount of translational velocity for the geometry in the Z direction
as measured from the Center of Mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
28 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 2: Ansys Motion Coupling Module
2.1. Introduction
Ansys Motion Coupling enables users to prescribe specific, complex motions for their equipment in Ansys
Motion and then have those same motions be used in a DEM simulation that is carried out by Rocky.
As illustrated in Figure 2.1: Standard Rocky Motion Frames vs. Ansys Motion Coupling Cabilities. (p. 29),
when compared to the standard method of defining Motion Frames in Rocky, movements defined with
Ansys Motion Coupling have the added benefit of simulating geometry-to-geometry interactions, stress
and deformation of geometry components, chained motions, and more.
Figure 2.1: Standard Rocky Motion Frames vs. Ansys Motion Coupling Cabilities.
The current implementation of Ansys Motion Coupling includes the installation of one module for the
Ansys Motion product, which enables Ansys Motion users to export all applicable geometries, motion
definitions, and motion solver files needed for coupling to a Functional Mock-Up Unit (FMU) file.
For Rocky, an internal module named Multibody Dynamics FMU Coupling is used to import that same
FMU file and have Rocky use it to simulate both the geometries and geometry motions during the DEM
simulation (Figure 2.2: Ansys Motion Coupling FMU file lifecycle. (p. 30)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 29
Ansys Motion Coupling Module
– You have installed one of the compatible versions of Ansys Motion 2025 R1
– You have installed the Ansys Motion Module: FMU Setup for Rocky component.
– You have a valid license for your version of Ansys Motion. It is also required that you have Ansys
Motion installed on the same machine as Rocky.
2.3. Limitations
This version of Ansys Motion Coupling has the following limitations:
• The Simulation Duration in Ansys Motion must be greater than or equal to the one set in Rocky.
In addition, there are specific limitations for each product, as described below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
30 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup procedures
• Flexible bodies are supported on the Ansys Motion side of the coupled simulation; however, Rocky
will be unable to calculate particle interactions on those flexible bodies.
• Including flexible bodies in your coupled simulation might also increase the processing time.
• Motion Frames: A Rocky project with Ansys Motion Coupling cannot also have Motion Frames.
Your simulation will not start processing until all Motion Frames are removed.
• Resume processing a stopped simulation: Once you click Stop on the Solver tab, you might be
unable to resume processing from the point you stopped.
• Extend the duration of a simulation: Once you click Stop on the Solver tab, you might be unable
to extend your processing time. Or if you are able to extend your processing time, you might have
errors resuming processing.
• Save a copy of the partially processed simulation for restart purposes: If you choose to Save
as a New Project for Restart your Ansys Motion Coupling project, you might get errors when
creating the file. Or, if you are able to create the file, you might be unable to restart it.
2. Open an existing subsystem file or create a new subsystem file, setting up your equipment and
motions as you normally would.
Tips:
• Ensure your Gravity direction in Ansys Motion is set to what you later want to use in Rocky.
These must match for the coupling process to work.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 31
Ansys Motion Coupling Module
• Ensure your simulation time (T) in Ansys Motion is greater than Simulation Duration the you
later want to use in Rocky.
Because you cannot resume or extend a stopped simulation in Rocky, it is recommended that
you make the simulation time longer than you think you need.
• It’s recommended to use the default Ansys Motion unit system (m, kg, N, s, V, A), using Ansys
Motion Pre or Ansys Motion inside Ansys Mechanical environment. If the user needs to work
with millimeters, change the Postprocessor units settings (inside Analysis Settings) to "Convert
to Motion-specific (mm kg N s rad)". Manually changing the import scale, when importing
geometries to Rocky, will not fix the movements unit scale on Ansys Motion.
3. From the Work Navigator panel, under Generate Information, select the .dfsub file representing
the geometries and motions you want to couple with Rocky (Figure 2.3: Co-Simulator tab in Ansys
Motion. (p. 32)).
4. From the Co-Simulator tab, click the Generate Information button on the Rocky Coupling
toolbar (Figure 2.3: Co-Simulator tab in Ansys Motion. (p. 32)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
32 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup procedures
5. From the Rocky Coupling dialog (Figure 2.4: Rocky Coupling dialog in Ansys Motion. (p. 33)), do
the following:
a. For each of the rigid Body components without beam elements listed, do one or both of the
following:
• Enable the Output checkbox for each component you want motions to be shared with
Rocky.
• Enable the Input checkbox for each component you want forces and moments to be col-
lected due to interactions with particles. This data will then be used by the Ansys Motion
solver during processing.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 33
Ansys Motion Coupling Module
6. From the Save As dialog, verify the file name and location, and then click Save. The progress bar
on the Rocky Coupling dialog indicates export status.
Important: To avoid errors during geometry import later, an FMU file should not be manually
renamed after it is generated.
At this point, you can share the FMU file with a Rocky user for them to complete the second part of
the coupling procedure, or continue with the next step yourself.
Sometimes, if inside the Ansys Motion set-up there are two or more geometries with the same name,
or maybe a geometry, a function expression, and a joint load properties named equally, Motion
automatically changes the name of the geometry with a prefix ID inside its post-processing. The user
can still export the FMU file with this configuration, but he will need to manually fix the geometry
name in Rocky’s side (remove the ID).
2. Begin a new simulation project. (From the File menu, click New Project.)
3. From the Data panel, under Modules, select Multibody Dynamics FMU Coupling.
4. From the Data Editors panel, on the Multibody Dynamics FMU Coupling tab, do the following
(Figure 2.5: Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module. (p. 34)):
Figure 2.5: Rocky Coupling dialog in Multybody Dynamics FMU Coupling Module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
34 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup procedures
a. Click the FMU filename: button, and then from the Select file to import dialog, locate and
select the FMU file you want to use, and then click Open. Save the project now if you are
asked to do so.
b. From the Communication step size control list, select one of the following options:
• Automatic: For each time step in the motion tool, Rocky computes 100 timesteps.
• Manual: The motion tool’s time step is defined by the Communication step size value.
• Rocky timestep: The time step in the motion tool is the same as in Rocky.
c. If Manual was selected for Communication step size control, then in the Manual commu-
nication step size field, enter what time step you want the motion tool to use.
5. From the Data panel, right-click Geometries, and then click Import Geometry.
6. From the Select file to import dialog, locate and select one or more files in STL format containing
the geometries used in the motion tool project.
Important:
It is critical that you import the exact same geometry parts required for the
coupling and with the same exact names. Otherwise, the motions may not work.
If you later need to re-import a newer version of this FMU file, ensure that you
also update the geometry files. See Reimport an updated FMU file into Rocky
procedure for more details.
7. From the Import File Info dialog, select the options you want, and then click OK. The Geometries
list is populated with the components that were exported out of the motion tool.
8. Set up the rest of your simulation as you normally would, being careful to ensure the following:
• You have defined the gravity direction in Rocky to match what was set in the motion tool.
• You have defined the Simulation Duration in Rocky to be less than or equal to the duration
that was set in the motion tool.
9. Process your simulation and analyze the results in Rocky as you normally would.
Once you have started processing the simulation, avoid stopping the simulation as you will be
unable to resume processing once it has started. Refer to the Limitations section for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 35
Ansys Motion Coupling Module
1. With the Rocky project open that contains information from the FMU file you want to update, do
the following:
a. Remove and then re-import the geometries added by the FMU file by doing the following:
i. From the Data panel, under Geometries, select or multi-select all the components that
were originally imported by the FMU file.
ii. Right-click the selected geometries, and then click Remove Geometries.
iii. Re-import the geometries via the updated FMU file as you normally would.
Alternatively, if you made only small changes to the bodies in Ansys Motion, you can make
those changes to individual Geometries in Rocky without removing all of them and importing
them again. However, you must ensure that the changes you make exactly match the inform-
ation in the FMU file you import into the Ansys Motion Coupling module.
b. Re-import the FMU file into the Ansys Motion Coupling module by following the Rocky setup
for coupling with Ansys Motion procedure.
2. Set up, process, and analyze the results of your simulation as you normally would.
2.5. Post-processing
After processing your coupled simulation, you will have new Curves for your geometry components,
depending upon whether they were selected to be exported out of Ansys Motion as an Output body
(Figure 2.6: Standard Ansys Motion Coupling Curves for a geometry that was exported out of Ansys Motion
as an Output. (p. 37)) – which has no interactions with Particles – or an Input body (Figure 2.7: Additional
Ansys Motion Coupling Curves for a geometry that was exported out of Ansys Motion as an Input. (p. 38))
– which does have interactions with Particles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
36 of ANSYS, Inc. and its subsidiaries and affiliates.
Post-processing
Figure 2.6: Standard Ansys Motion Coupling Curves for a geometry that was exported out of Ansys
Motion as an Output.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 37
Ansys Motion Coupling Module
Figure 2.7: Additional Ansys Motion Coupling Curves for a geometry that was exported out of
Ansys Motion as an Input.
• Acceleration : Rotational : X: The amount of rotational (angular) acceleration for the geometry in
the X direction as measured from the Center of Mass.
• Acceleration : Rotational : Y: The amount of rotational (angular) acceleration for the geometry in
the Y direction as measured from the Center of Mass.
• Acceleration : Rotational : Z: The amount of rotational (angular) acceleration for the geometry in
the Z direction as measured from the Center of Mass.
• Acceleration : Translational : X: The amount of translational acceleration for the geometry in the X
direction as measured from the Center of Mass.
• Acceleration : Translational : Y: The amount of translational acceleration for the geometry in the Y
direction as measured from the Center of Mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
38 of ANSYS, Inc. and its subsidiaries and affiliates.
Post-processing
• Acceleration : Translational : Z: The amount of translational acceleration for the geometry in the Z
direction as measured from the Center of Mass.
• Force : X : Instantaneous: Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the force applied to the geometry along the X axis.
• Force : Y : Instantaneous: Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the force applied to the geometry along the Y axis.
• Force : Z : Instantaneous: Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the force applied to the geometry along the Z axis.
• Moment : X: Instantaneous: Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the moment (torque) applied to the geometry’s Center of Mass around the X axis.
• Moment : Y: Instantaneous : Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the moment (torque) applied to the geometry’s Center of Mass along the Y axis.
• Moment : Z: Instantaneous: Only for geometries that were exported out of Ansys Motion as an Input
body, this provides the moment (torque) applied to the geometry’s Center of Mass along the Z axis.
• Orientation Angle: The angle the geometry is rotated around the Orientation Vector.
• Orientation Vector X: The X component defining the vector around which the geometry is rotated.
• Orientation Vector Y: The Y component defining the vector around which the geometry is rotated.
• Orientation Vector Z: The Z component defining the vector around which the geometry is rotated.
• Position X: The global coordinate of the geometry’s Center of Mass as specified in the X direction.
• Position Y: The global coordinate of the geometry’s Center of Mass as specified in the Y direction.
• Position Z: The global coordinate of the geometry’s Center of Mass as specified in the Z direction.
• Velocity : Rotational : X: The amount of rotational (angular) velocity for the geometry in the X dir-
ection as measured from the Center of Mass.
• Velocity : Rotational : Y: The amount of rotational (angular) velocity for the geometry in the Y direction
as measured from the Center of Mass.
• Velocity : Rotational : Z: The amount of rotational (angular) velocity for the geometry in the Z direction
as measured from the Center of Mass.
• Velocity : Translational : X: The amount of translational velocity for the geometry in the X direction
as measured from the Center of Mass.
• Velocity : Translational : Y: The amount of translational velocity for the geometry in the Y direction
as measured from the Center of Mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 39
Ansys Motion Coupling Module
• Velocity : Translational : Z: The amount of translational velocity for the geometry in the Z direction
as measured from the Center of Mass.
2.6. Troubleshooting
If you have trouble with either creating the FMU file in Ansys Motion or importing the FMU file in Rocky,
see the sections below for solutions.
If you have trouble finding the modules in either program, refer instead to the Ansys Motion Coupling
Installation Guide.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
40 of ANSYS, Inc. and its subsidiaries and affiliates.
Troubleshooting
If you require these flexible bodies for proper motion setup, you can click OK and continue with
your FMU file generation as usual. Just be aware that the flexible bodies will not be able to interact
with particles in Rocky when the coupled simulation is processed. Also know that a coupled simu-
lation with flexible bodies might take longer to process.
If you do not require these flexible bodies for proper motion setup, you can help reduce processing
time later by removing the flexible bodies from the Ansys Motion project, and then generating the
FMU file again. (Click OK on this message, cancel out of the FMU generation dialog, remove those
flexible bodies from your Ansys Motion project, and then click the Generate Information | Rocky
Coupling button again.) If you have properly removed the flexible bodies, then you should no
longer see the message shown in Figure 2.7: Additional Ansys Motion Coupling Curves for a geometry
that was exported out of Ansys Motion as an Input. (p. 38).
If you do not have this kind of license for Ansys, please contact your Ansys and/or Rocky represent-
ative for assistance.
If you have duplicated geometries in your Rocky project, then you might have re-imported geometries
from the FMU file before removing the old geometry components. Refer to the Re-import an updated
FMU file into Rocky section for specific steps.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 41
Ansys Motion Coupling Module
While you do not need to have a copy of Ansys Motion on the same machine you are running
Rocky, you do need to have an active license of Ansys in order to process the coupled simulation
in Rocky. Ensure your system meets the System Requirements, and then try processing your
coupled simulation again.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
42 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 3: Bond Model Module
3.1. Usage details
Figure 3.1: Example application of a drilling process using the Bond Model module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 43
Bond Model Module
where neighboring particles are bonded to represent the concrete. As the rotary drill bit forwards
against the wall, the bonds strength limits are exceeded and the hole is formed.
The bonds created by the Bond Model exert elastic and viscous forces and moments on the neighboring
particles as a reaction to deformations caused by their relative motion. Moreover, the bond model
includes also the possibility of adding particle-to-boundary bonds to the simulation.
Internally, the module associates every bond in the simulation to an adhesive contact and, as a con-
sequence, Rocky recognizes the module as an adhesive model. Therefore, even if there are no adhesion
effects enabled by this module, the iterative method being used takes the place of an adhesion
model and will therefore override what you have set for Adhesive Force on the Momentum tab of
the Physics entity (Figure 3.2:Physics | Momentum tab options when the Bond Model module is en-
abled. (p. 44)).
Figure 3.2: Physics | Momentum tab options when the Bond Model module is enabled.
Simulations using the Bond Model require the contacts data collection to be enabled by the Collect
Contacts Data checkbox from Contacts entity (Figure Figure 3.3:Contacts tab options. (p. 45)).
Moreover, these simulations can be resumed if you enable the option to save the adhesive contact
data. To do that, you must also enable the Include Adhesive Contacts checkbox for Contacts.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
44 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Simulations using the Bond Model module cannot be saved as project for Restart.
If you are running a simulation with both Bond Model and Particle Freeze Outside Region enabled, keep
in mind that any bond that leaves the region of interest will be lost. Thus, particles that were bonded
will no longer have bonds if they reenter the region of interest.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 45
Bond Model Module
where:
• Reference Time Step for non-Viscous Damping: This defines the reference amount of time
to dissipate energy through the non-viscous damping model.
Note:
This only needs attention if the non-viscous damping model will be used. Otherwise,
you can skip this setting.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
46 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
This only needs attention if you wants to customize the Bond radius. Otherwise,
you can skip this setting.
– Region Of Interest: The Region of Interest that you want to apply a radius multiplier to the
bond that are inside it when the bond is activated.
Range: [Any Region of Interest that has already been defined in the simulation]
– Radius Multiplier: This defines the radius multiplier that will be applied to the Bond radius.
There are ten additional parameters that you must set in the Materials Interactions entity, as shown
in Figure 3.5:Bond Model module materials interactions options. (p. 48).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 47
Bond Model Module
• Enable Bonds: This enables the bonds for the particles associated with the selected materials.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
48 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Approximate Maximum Elongation: The estimated value of the largest elongation of any bond
during a simulation. The value of this parameter must be chosen carefully because a bond might
be deactivated if its actual elongation exceeds that maximum value, even if the bond breakage
criteria are not satisfied. On the other hand, the simulation performance can be compromised if
too large values are prescribed for this parameter.
• Damping Ratio: This defines the damping ratio attributed to the bond.
• Distance Factor: This defines the factor for the bonds activation distance calculation.
• Maximum Activation Distance: This defines the maximum distance for the bonds activation.
• Normal Stiffness Per Area: This defines the normal stiffness per area, which is used to calculate
the normal component of the linear elastic force applied on the particles in each bonded pair.
• Reference non-Viscous Damping Coefficient: This defines the reference non-viscous damping
coefficient attributed to the bond.
Note:
This is the input parameter of the Non-viscous damping model. If this parameter is
defined as zero, the Non-viscous damping model will not have effect on the simulation.
• Shear Stress Limit: This defines the shear stress limit of the bonds.
• Tangential Stiffness Per Area: This defines the tangential stiffness per area, which is used to
calculate the tangential component of the linear elastic force applied on the particles in each
bonded pair.
• Tensile Stress Limit: This defines the tensile stress limit of the bonds.
Since there are no adhesion effects enabled by this module, the parameter Adhesive Distance will
not be used in any calculations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 49
Bond Model Module
For the Inputs step, a new Modules sub-tab will appear to help you define the properties of the
particles, as shown below.
Note:
This sub-tab will appear only for the Particle Inlet and Volume Fill inputs. For Particle
Custom Inlets, there is not a Modules sub-tab; rather, the Modules parameters will appear
on the main tab for Particle Custom Inlets.
• Bond Activation Time: Defines the specific time when the bonds activation happens.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
50 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
In addition, some extra contact properties are available as shown in Figure 3.7: Additional Contact
Properties created by the Bond Model module. (p. 51) and explained below.
Figure 3.7: Additional Contact Properties created by the Bond Model module.
Note:
Properties generated by this module will only appear if the Meshed Particles Upscaling is
disabled. Refer to Rocky User Manual (About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 51
Bond Model Module
• Bond Shear Stress: is the bond shear stress at the current time.
• Bond State: is a criterion that classifies whether the bond is active, inactive or broken. When equals
to 0, the bond is inactive; when equals to 1, the bond is active; and when equals to 2, the bond is
broken.
• Bond Tensile Stress: is the bond tensile stress at the current time.
Moreover, two additional properties are available for the particles as shown in Figure 3.8: Additional
Particle Properties created by the Bond Model module. (p. 52) and explained below.
Figure 3.8: Additional Particle Properties created by the Bond Model module.
• Bond Activation Time: is the time when the bond was activated.
• Particle Broken Bonds: is the number of broken bonds for the particle.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
52 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
3. From the Data Editors panel, enter the values you want for the bond model parameters.
5. Prior to processing your simulation, ensure that you do the following for each material that will
be assigned to particles that you want to have bonds in your project:
a. From the Data panel, select Materials, and from the Materials Interactions tab, select the
materials that you assigned to the particles that will have bonds, and then from the Data Ed-
itors panel do the following:
i. From the Materials Interactions tab, enable the Enable Bonds option.
ii. From the Materials Interactions tab, enter the values you want for the bond model
parameters.
b. From the Data panel, under Inlets and Outlets, select the inlet entity, and then from the Data
Editors panel do the following:
7. When you are ready to analyze your simulation results, you may choose to make use of any contact-
related property and may also make use of the additional module-related properties.
Bonds are created between pairs of particles (or a particle and boundary) satisfying the following two
conditions:
• Bonds are explicitly enabled for the material interaction associated to the pair of particles (or a particle
and boundary) by turning on the checkbox Enable Bonds in the Rocky UI.
• The closest points of the particles pair surfaces (or a particle and boundary) have a distance from
each other, which is lower/equal than the activation distance given by:
(3.1)
where and are the radii of the bonded particles, as shown in Figure 3.9: Geometry of the bond
between two spherical particles of different sizes. (p. 54), while and are input parameters
listed as Distance Factor and Maximum Activation Distance in the Rocky UI, respectively. The ac-
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 53
Bond Model Module
tivation of bonds happens only once during a simulation, at a specific time defined by the value of
the parameter listed as Bond Activation Time in the Rocky UI.
Note:
For non-spherical particles, the radius is replaced by one half of the sieve size in
equation Equation 3.1 (p. 53) and other equations including it.
Since Rocky calculates forces and solves motion equations at discrete times, the actual
activation of the bonds will happen at the discrete time nearest to the specified activ-
ation time.
Figure 3.9: Geometry of the bond between two spherical particles of different sizes.
The search algorithm has an increment to compensate the bond tensile deformations during the simu-
lation. This will guarantee the Bond remains inside the search algorithm and avoids deactivating a Bond
wrongly:
(3.2)
The radius of the bond between two particles (or a particle and boundary) depends only on the radii
of those particles, by means of the following equation:
(3.3)
Note:
If you skipped the settings of Radius Multiplier Regions on Rocky UI, a default value of
1 will be used.
At the time of its activation, a bond is in an undeformed state. From that time on, any relative motion
of the bonded particles will cause a deformation on the bond, to which it will react exerting forces and
moments opposing to that relative motion. The bond deformation can be linear or angular, being the
former caused by both the translational and rotational motion of the bonded particles, while the latter
is exclusive consequence of their rotational motion. In the module, both types of deformation are cal-
culated incrementally, starting from the time of bond activation, . Thus, the linear deformation of a
bond at any given time is given by:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
54 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(3.4)
where is the instantaneous relative velocity at the center point of the bond, , while is the sim-
ulation timestep. Internally, the module associates every bond in the simulation to an adhesive contact.
Therefore, the center point of a bond will coincide with the point where Rocky locates the associated
adhesive contact. That point is located in the middle of the shortest segment between the two bonded
particles, as depicted in Figure 3.9: Geometry of the bond between two spherical particles of different
sizes. (p. 54).The relative velocity is calculated as:
(3.5)
where:
• and are the rotational velocities of the bonded particles. For particle-to-boundary bonds,
and will be the translational and rotational velocity of the boundary.
• is the vector that joins the centroid of the particle to the bond's center point . The definition
is equivalent for the vector .
The bond in the considered model will oppose to the linear deformation by exerting an elastic force
on both bonded particles. The normal and tangential components of this force will be given, respectively,
by:
(3.6)
(3.7)
where:
• and are the normal and tangential stiffnesses per unit area, respectively. These are input
parameters which are listed as Normal Stiffness Per Area and Tangential Stiffness Per Area, re-
spectively, in the Rocky UI.
• and are, respectively, the normal and tangential components of the bond's linear deformation
.
The decomposition of the vector is carried out by using the following simple expressions:
(3.8)
(3.9)
in which is the normal unit vector parallel to bond's axis, as depicted schematically in Figure 3.10: Forces
and moments exerted by a bond on a particle. (p. 56). At any given time, the shortest distance between
the two bonded particles will be in the direction of this unit vector. It is worth noting that depending
on the relative motion of the particles, can be positive or negative. Therefore, the normal force
can be attractive or repulsive.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 55
Bond Model Module
As previously stated, a bond will also react to angular deformations caused by the relative rotation of
the bonded particles, by exerting a moment opposing to that deformation. As in the linear deformation
case, the angular deformation of a bond, , is also calculated incrementally. The expression employed
in this case is:
(3.10)
The normal component of the angular deformation vector, , represents a torsional deformation of
the bond, while the tangential component, , characterizes the bond's bending. The corresponding
torsional and bending components of the moment that oppose to those angular deformations are cal-
culated, respectively, as:
(3.12)
(3.13)
where:
• and are the normal and tangential stiffnesses per unit area, respectively. These are input
parameters which are listed as Normal Stiffness Per Area and Tangential Stiffness Per Area, re-
spectively, in the Rocky UI.
In order to obtain the normal and tangential components of the bonds angular deformation, this vector
is decomposed using expressions equivalent to Equation 3.8 (p. 55) and Equation 3.9 (p. 55).
In addition to the elastic forces and moments described above, in this module a bond may also exert
viscous forces and moments intended for reducing oscillations that may arise between bonded particles.
The normal and tangential components of those forces and moments are calculated using the following
expressions:
(3.14)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
56 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(3.15)
(3.16)
(3.17)
where:
• is the damping ratio attributed to the bond. The value of this parameter can be specified by users
and is listed in the Rocky UI as Damping Ratio. For this parameter is allowed any value in the interval
between 0 and 1, the larger the value, the faster oscillations will be damped out.
• is the equivalent mass of the pair of bonded particles, whose masses are and . For
particle-to-boundary bonds, this equivalent mass will be equal to the mass of the bonded particle.
• and are the normal and tangential components, respectively, of the relative velocity vector
defined by Equation 3.5 (p. 55). Expressions equivalent to equations Equation 3.8 (p. 55) and Equa-
tion 3.9 (p. 55) are used to decompose that vector.
• and are the normal and tangential components, respectively, of the relative rotational velocity
vector defined by Equation 3.11 (p. 56). Expressions equivalent to equations Equation 3.8 (p. 55) and
Equation 3.9 (p. 55) are used to decompose that vector.
This implementation can also include a non-viscous damping model to the bond, which is defined as
a relaxation factor. The bond forces and moments are calculated in a time incremental algorithm and
the relaxation factor is applied to the previous forces and moments during the calculation [2] (p. 377).
The description for the time incremental and next value are written below for the normal bond force,
but the same incremental method is applied to the other forces and moments:
(3.18)
(3.19)
where is the normal bond force from previous timestep and is the non-viscous damping coeffi-
cient. If , the previous timestep force value is fully used to the next timestep and no force/energy
is dissipated, and if , a dissipation occurs since the previous timestep value is only partially carried
to the next timestep.
The non-viscous damping coefficient is simulation timestep independent and defined as:
(3.20)
where:
• is the Reference non-Viscous Damping Coefficient parameter available in the Rocky UI.
• is the Reference Time Step for non-Viscous Damping parameter available in the Rocky UI.
Another distinctive feature of the bond model implemented in this module is the ability to break bonds
whose strength capacity has been exceeded, and remove them from the simulation. The breakage cri-
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 57
Bond Model Module
terion is based on the maximum values of the tensile and shear stresses acting on the periphery of the
bond. Those values are given, respectively, by [3] (p. 377):
(3.21)
(3.22)
During a simulation the values of the maximum tensile and shear stresses are constantly monitored for
all active bonds. When exceeds the tensile strength limit or exceeds the shear strength limit,
the bond breaks and, consequently, is deactivated immediately. After such a breakage event, a bond
cannot be reactivated again during the simulation. The values of the tensile and shear stress limits can
be specified through the parameters Tensile Stress Limit and Shear Stress Limit, respectively, in the
Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
58 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 4: Boundary Outlet Module
4.1. Usage details
Selecting this module will not override any models or settings in the Rocky UI.
For the Geometries setup step, a new Modules sub-tab will appear to help you define the properties
of the geometry, as shown in Figure 4.1:Boundary Outlet module options for an individual geometry
component that you have imported. (p. 60) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 59
Boundary Outlet Module
Figure 4.1: Boundary Outlet module options for an individual geometry component that you
have imported.
Boundary Outlet
• Treat Boundary as an Outlet for Particles: When enabled, this boundary will be used as an
outlet for particles; any particles that come in contact with this boundary will be removed from
the simulation. When cleared, this boundary will not be used as an outlet for particles.
• Start Time:This defines the time you want the boundary to begin acting as an outlet for particles.
• Stop Time: This defines the time you want the boundary to stop acting as an outlet for particles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
60 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Boundary Outlet checkbox is enabled.)
3. From the Data panel, under Geometries, select (or multi-select) the individual geometry (or
geometries) that you want to define as an Boundary Outlet and then from the Data Editors panel,
on the Geometry | Modules tab (or the ...Conveyor | Modules tab), do the following:
b. Define the Start Time and Stop Time values you want
4. Set up and process the simulation as you normally would. During the simulation, any particle that
comes into contact with a boundary defined as an outlet for particles will be immediately removed
from the simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 61
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
62 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 5: CFD Particle Combustion Module (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However,
your feedback will help us improve the overall quality of the product. We will not
guarantee that the projects using this beta feature will run successfully when the feature
is finally released so you may, therefore, need to modify the projects.
In this problem, conduction and convection are considered in the heat transfer modeling, as the
local temperature affects the combustion rates. Radiation effects can be included through the Thermal
Radiation module. Because of these abilities, it is designed to work with Rocky's 2-Way Fluent coupling
method, and Thermal Model.
Note:
For more information about the CFD Coupling mode and the Thermal Model options
included in Rocky, refer to the Rocky User Manual.
Two reactions occur in parallel. In the first one, the oxygen from the air oxidizes the particle's carbon
content. This leads to the oxygen and carbon consumption, carbon dioxide generation and the con-
sequent particle size reduction. Similarly, in the second one, the carbon dioxide from the air reacts
with the particle's carbon content to release carbon monoxide, reducing the particle size as well.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 63
CFD Particle Combustion Module (Beta)
Figure 5.1: Particle combustion over time using the CFD-Coupled Particle Combustion
module.
As the reactions occur, the burn-off (quantity of reactive, i.e. carbon, mass consumed) is updated,
depending on how the particle mass decreases, also reflecting in its volume. This is done, preserving
the particle density constant. On the fluid side, the mass consumption/production of species is sent
to Fluent in the form of mass sink/source terms.
The heat of reaction for both reactions are considered to affect only the particle energy equation,
and therefore, no corresponding energy source terms are sent to Fluent.
• The module only works for oxygen, nitrogen, carbon dioxide and carbon monoxide, so the fluid
species must be named "o2", "n2", "co2" and "co", respectively.
Important:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
64 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 5.2: CFD-Coupled Particle Combustion module options under the Species tab.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 65
CFD Particle Combustion Module (Beta)
Figure 5.3: CFD-Coupled Particle Combustion module under the Reactions Tab.
Figure 5.4: CFD-Coupled Particle Combustion module options under the Advanced Tab.
Oxygen:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
66 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
Refer to the Technical details (p. 75) section for more information.
Carbon Dioxide:
Note:
Refer to the Technical details (p. 75) section for more information.
Carbon Monoxide:
Note:
Refer to the Technical details (p. 75) section for more information.
R1:
Note:
Refer to the Technical details (p. 75) section for more information.
• R1 activation temperature [K]: Activation temperature for the carbon-oxygen reaction (R1).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 67
CFD Particle Combustion Module (Beta)
Note:
Refer to the Technical details (p. 75) section for more information.
R2:
• R2 pre-exponential factor [m/s.K] : Pre-exponential factor for the carbon-carbon dioxide re-
action (R2).
Note:
Refer to the Technical details (p. 75) section for more information.
• R2 activation temperature [K]: Activation temperature for the carbon-carbon dioxide reaction
(R2).
Note:
Refer to the Technical details (p. 75) section for more information.
Advanced:
Important:
Rocky Solver SDK cannot get the reference pressure from Fluent, so this is mistake-prone.
• Maximum allowed burn-off: Maximum burn-off values. Above this value, the particle is removed.
• Combustion cut-off particle size ratio: Minimum ratio between the carbon core size and the
particle total size. Below this value, the particle is no longer reactive, but is not removed from the
simulation.
Range: [0,1]
• Combustion Start Time [s]: Defines the time that combustion starts to take place.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
68 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 69
CFD Particle Combustion Module (Beta)
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab.
• Carbon content: The mass fraction of carbon that a particle initially has.
Range: [0,1]
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
70 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
Refer also to the Graphing (Plot or Histogram) a Data Set Within Rocky topic in the Rocky
User Manual.
Note:
Refer also to the View Geometries, Particles, Points, and Fluids in 3D topic in the Rocky
User Manual.
• Carbon content: The mass fraction of carbon that a particle initially has.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 71
CFD Particle Combustion Module (Beta)
• Burn-off: The instantaneous burn-off, i.e., the percentage of reactive mass (carbon) already con-
sumed.
• R1 Ash Layer Resistance: the ash layer resistance for the carbon-oxygen reaction.
• R1 Mass Diffusion Resistance: the mass diffusion resistance for the carbon-oxygen reaction.
• R1 Surface Reaction Resistance: the surface reaction resistance for the carbon-oxygen reaction.
• R2 Ash Layer Resistance: the ash layer resistance for the carbon-carbon dioxide reaction.
• R2 Mass Diffusion Resistance: the mass diffusion resistance for the carbon-carbon dioxide reaction.
• R2 Reaction Rate: the rate in which the carbon-carbon dioxide reaction occurs.
• R2 Surface Reaction Resistance: the surface reaction resistance for the carbon-carbon dioxide
reaction.
• Unreacted Core Size Ratio: the ratio between the carbon core and the particle size.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
72 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 5.8: Particles properties when the CFD-Coupled Particle Combustion is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 73
CFD Particle Combustion Module (Beta)
Step 2: Ensure that the module is enabled (From the Data panel, select Modules and then from the
Data Editors panel, ensure the CFD-Coupled Particle Combustion checkbox is enabled).
Step 3: From the Data panel, under Modules, select the new entry for CFD-Coupled Particle Combus-
tion.
Step 4: From the Data Editors panel, on the CFD-Coupled Particle Combustion tab, enter the values
you want.
Note:
Step 5: Continue setting up the simulation as you normally would, with the following exceptions:
1. When setting up your Physics parameters, ensure that the Thermal Model is enabled.
Note:
Refer to the Rocky Enable Thermal Modeling Calculations topic in the Rocky User Manual
for details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
74 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
2. When setting up your Physics parameters, ensure that the Thermal Model is enabled.
Note:
Refer to the Rocky User Manual for setup details and/or refer to the Set or Modify Fluid
and/or Air Flow Properties topic in the Rocky User Manual for details.
Note:
Step 7: When you are ready to post-process your simulation results, you can make use of the new
parameters on the Properties tab for the main Particles entity.
Note:
The carbon from the particle reacts with the oxygen, producing carbon dioxide. In parallel to it, the
carbon from the particle reacts with carbon dioxide, producing carbon monoxide. The two reactions
can be represented by the following equations:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 75
CFD Particle Combustion Module (Beta)
where and are the reaction rates (mol/s) for and , respectively.
As in Zhang et al. (2013), the reaction rates are generically given as:
where:
• is the stoichiometric coefficient. As it is one for both reactions, it is suppressed for simplicity;
where:
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
76 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
The ash layer resistance is given as in Wang and Shen (2021) by:
where:
Assuming that the particle has a constant density, the unreacted core diameter ratio is:
where:
and imply the consumption and the generation of species, which reflects on how the particle
mass changes over time:
Conversely, this also happens for each fluid species, and the mass source terms are sent to Fluent.
For the energy conservation, the first law of thermodynamics is given by:
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 77
CFD Particle Combustion Module (Beta)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
78 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 6: CFD Particle Drying Module
6.1. Usage details
In many applications, wet granular material is allowed to dry by the action of a drying media, such
as hot air, whose main role is to remove moisture and carry it downstream. This module enables you
to specify how much moisture is in a particle at the start of the simulation, and then define both how
it dries and over what time period.
Because of these abilities, it is designed to work in combination with several Rocky CFD Coupling
modes (specifically, the 1-Way Constant, 1-Way Fluent Steady State, and 2-Way Fluent methods), the
Rocky Thermal Model, and the Coarse-Grain Model (CGM).
Figure 6.1: Particle drying rate over time as a function of the moisture content.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 79
CFD Particle Drying Module
As drying takes place, diffusion of the water vapor across the air-moisture interface is the rate-con-
trolling mechanism. Towards the end, moisture has to be transported from inside the particle to the
surface by capillary forces and the drying rate may still be constant until the critical moisture content
is reached.
At this point, the surface film of moisture has been so reduced by evaporation that further drying
causes dry spots to appear upon the surface. This stage is called the falling rate period and it proceeds
until the surface film of the liquid is entirely evaporated.
In order to keep on drying, the ability with which moisture can move through the solid as a result of
concentration gradients between the deeper parts of the particle and its surface controls the rate of
evaporation. This is known as the 2nd stage of the falling rate period. Overall, the drier the particle,
the harder it is to keep drying. When the condition to which the particle is subjected is no longer
capable of removing any additional moisture, the equilibrium moisture content is reached.
Species transfer between Rocky and Fluent is only considered for the water vapor when this module
is used with the 2-Way Fluent coupling mode. Therefore, when used with the other supported coupling
modes (1-Way Constant and 1-Way Fluent Steady State), humidity is considered constant and local
humidity changes are not taken into account.
• Either the 1-Way Constant, 1-Way Fluent Steady State, or 2- Way Fluent CFD Coupling Mode
must be enabled in Rocky. Also:
– In the case of 1-Way Fluent Steady State coupling, a temperature field must be available in the
Fluent case.
– In the case of 2-Way Fluent coupling, the solution for the Energy equation must be activated
in the Fluent case.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
80 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 81
CFD Particle Drying Module
• Solvent Molar Mass: The solvent molar mass, which should be provided in the SI unit.
• Solvent Mass Diffusivity (Drying Media): The solvent mass diffusivity, which should be provided
in the SI unit.
• Solute Mass Fraction *: This defines the ratio mass of solute to the liquid solution mass. As the
liquid evaporates, the equivalent solute mass will be added to the particle mass and the particle
volume will be updated accordingly.
Range: [0,1]
• Critical Moisture Content: The particle moisture content below which evaporation transitions from
a constant to a falling rate regime.
Range: [0,1]
• Equilibrium Moisture Content: The particle moisture content at which evaporation no longer
takes place.
Range: [0,1]
• Drying Air Relative Humidity: The relative humidity of the drying air. This will be treated as a
constant scalar field.
Range: [0,1]
• Drying Start Time *: Defines the time that drying starts to take place.
• Drying Stop Time *: Defines the time that drying stops taking place.
• Drying Duration: Defines for how long drying will take place.
Range: [Positive values smaller than Drying Stop Time - Drying Start Time]
• Drying Period: Drying cycle period. This parameter combined with Drying Start Time, Drying
Stop Time, and Drying Duration, allows the user to define periodic drying cycles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
82 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Switch to Absorption?: Allows the module to remove moisture from the surrounding media.
For Particle Inlet Inputs and Volumetric Inlet Inputs, a new Modules sub-tab will appear to help you
define the initial moisture content of the particles, as shown in Figure 6.3: Options for a Particle Inlet
Input (left) and a Volumetric Inlet Input (right). (p. 83).
Figure 6.3: Options for a Particle Inlet Input (left) and a Volumetric Inlet Input (right).
Within the table on the Modules sub-tab, one parameter (column) can be set per Particle set (row),
as defined below:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 83
CFD Particle Drying Module
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab.
• Moisture Content: This defines the initial moisture content that will be applied to each particle
within the selected Particle set when it enters into the simulation.
• Mass Increment: This defines the initial increment mass-in relation to the original mass of each
individual whole particle or fragment-that will be applied to each particle within the selected
Particle set when it enters into the simulation.
• Volume Increment: This defines the initial increment volume-in relation to the original volume
of each individual whole particle or fragment-that will be applied to each particle within the se-
lected Particle set when it enters into the simulation.
For Particle Custom Inlet, a new section will appear to help you define the moisture content input
properties of the particles, as shown in Figure 6.4: Options for a Particle Custom Inlet. (p. 85). The
parameters in this section are defined exactly the same as indicated in the preceding section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
84 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 85
CFD Particle Drying Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
86 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 6.5: Particles Properties when the CFD-Coupled Particle Drying module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 87
CFD Particle Drying Module
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Or, you can choose to display this property graphically in a 3D View window.
• Evaporation Heat Transfer Rate: The instantaneous evaporation heat transfer rate being experi-
enced by the particle.
• Mass Increment: The increment mass in relation to the original mass of each individual whole
particle or fragment. This can help you to visualize the amount of mass that was absorbed by each
particle over time.
• Rate of Evaporation: The instantaneous rate of evaporation being experienced by the particle.
• Volume Increment: The increment volume in relation to the original volume of each individual
whole particle or fragment. This can help you to visualize the amount of volume that was absorbed
by each particle over time.
• If you plan to use this module with a 2-Way Fluent simulation, ensure the Fluent case has the
Energy equation activated.
• If you plan to use this module with a 1-Way Fluent Steady State simulation, ensure the Fluent
case has a temperature field available.
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the CFD-Coupled Particle Drying checkbox is enabled.)
3. From the Data panel, under Modules, select the new CFD-Coupled Particle Drying entry.
4. From the Data Editors panel, on the CFD-Coupled Particle Drying tab, enter the values you
want.
5. Continue setting up the simulation as you normally would, with the following exceptions:
a. When setting up your Physics parameters, ensure that the Thermal Model is enabled
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
88 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
b. When setting up your Inputs, be sure to define the parameters on the Modules tab (for Particle
Inlet Inputs and Volumetric Inlet Inputs) or in the CFD-Coupled Particle Drying section (for
Particle Custom Inlet).
c. When setting up your CFD Coupling options, ensure you select either 1-Way Constant, 1-
Way Fluent Steady State, or 2-Way Fluent for Coupling Mode.
7. When you are ready to post-process your simulation results, you topic in the Rocky User Manual.
can make use of the new parameters on the Properties tab for the main Particles entity.
For the Fluent configuration, it is important to have the Species Transport model turned on, as shown
below:
Another step required for the module to work is setting the mixture-template to be formed by air +
H2O, ensuring that air is the Last Species:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 89
CFD Particle Drying Module
If we neglect the contribution of radiation, we can see that the only missing quantity is the heat transfer
rate due to evaporation. Since this is a phase-change process, we can write it as:
(6.2)
where is rate of evaporation, and is the latent heat of evaporation at , which can be expressed
as:
(6.3)
where:
is the ratio between the particle temperature and the critical temperature.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
90 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(6.4)
with the moisture content Y defined as the ratio between the moisture mass and the particle bone-
dry mass . The rate of evaporation can be obtained from a general mass transfer
correlation:
(6.5)
where:
and are the mass density of water vapor at the particle surface and free stream, respectively.
(6.6)
The approach adopted in this module relies on empirical correlations based on the Sherwood number,
from which can be obtained by solving:
(6.7)
In equation [eq-sherwood] (p. 91), represents the mass diffusivity coefficient and should not be
mistaken with the particle diameter . itself is obtained from:
(6.8)
where and are the local Reynolds and Schmidt numbers, respectively.
In order to mimic the falling rate regime, we define a damping function , such that:
(6.9)
(6.10)
(6.11)
The damping function is applied to the rate of evaporation and therefore, if the particle moisture content
approaches the equilibrium moisture content, it will no longer decrease. On the other hand, if the
moisture content approaches the critical moisture content, the particle will lose moisture at a constant
rate. The CFD-Coupled Particle Drying module uses a simple correlation to obtain , as defined below.
(6.12)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 91
CFD Particle Drying Module
(6.13)
(6.14)
(6.15)
(6.16)
Figure [f-damping_b] (p. 92) presents the damping factor as a function of the moisture content for
different coefficients when . As the coefficient is increased, the damping factor behaves
in a more linear way.
Figures [f-damping_nc] (p. 93) and [f-damping_nv] (p. 94) present the damping factor as a function of
the moisture content for different and coefficients when , respectively. Higher values of
the and coefficients will yield higher moisture content for the same damping factor.
Figure 6.9: Damping factor as a function of the moisture content for different damping
function B coefficients with nc = nv = 1.0.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
92 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Figure 6.10: Damping factor as a function of the moisture content for different damping
function nc coefficients with B = 0.501.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 93
CFD Particle Drying Module
Figure 6.11: Damping factor as a function of the moisture content for different damping
function nv coefficients with B = 0.501.
Finally, the rate of evaporation considering the damping function can be rewritten as:
(6.17)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
94 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 7: Coordination Number Module
7.1. Usage details
Figure 7.1: Example application of grains packing using the Coordination Number module.
Figure 7.2: Example application of particle mixing using the Coordination Number module. (p. 96)
shows another example of this module application. Initially, particles of two different particle groups
are placed separately inside the bowl (left), when there is low interaction between the particles and
consequently low coordination numbers. After some time of mixing, we can see higher coordination
numbers for most of the particles, indicating a mixing homogeneity (right).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 95
Coordination Number Module
Figure 7.2: Example application of particle mixing using the Coordination Number module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
96 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Calculate inter-group CRN: When disabled, the module will calculate the instantaneous coordination
number for each particle present in the simulation and the average value for each particle. When
enabled, the module will also calculate instantaneous coordination numbers for all the combinations
of particle group pairs present in the simulation.
In addition, some extra particle properties are available if you enabled the module to process your
simulation, as shown in Figure 7.4:Coordination Number module options. (p. 98) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 97
Coordination Number Module
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• CRN | All Particle Groups: This provides the average coordination number considering all particle
groups. The contacts between the particles from all particle groups are considered for the coordin-
ation number calculation.
• CRN Average | All Particle Groups: This provides the time average coordination number consid-
ering all particle groups.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
98 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• CRN | Particle <01> X Particle <01>: This provides the coordination number considering only the
particle group Particle <01>. In other words, only the contacts between the particles from the
particle group Particle <01> are considered for the coordination number calculation.
Note:
This property's name composition carries the same name(s) you have defined for
the particle group(s).
This property will be available only if the option Calculate inter-group CRN is en-
abled.
A new property is created for each particle group pair, so the coordination number
is calculated for all the particle group combinations.
• CRN | Particle <02> X Particle <01>: This provides the coordination number considering the
particle groups Particle <02> and Particle <01>. In other words, only the contacts between the
particles from the particle group Particle <02> and particles from the particle group Particle
<01>are considered for the coordination number calculation.
Note:
This property's name composition carries the same name(s) you have defined for
the particle group(s).
This property will be available only if the option Calculate inter-group CRN is en-
abled.
A new property is created for each particle group pair, so the coordination number
is calculated for all the particle group combinations.
1. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Coordination Number checkbox is enabled.)
2. According to your needs, enable or not the module option Calculate inter-group CRN.
4. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related property and may also make use of the additional module-related property.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 99
Coordination Number Module
(7.1)
where:
Additional calculations are performed by the module to show the instantaneous coordination numbers
related to each particle group pairs, so you can know the number of contact points of each particle
from a particle group with particles from other specific particle group.
(7.2)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
100 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 8: Dhole, Chhabra & Eswaran Drag Law
Module
8.1. Usage details
Non-Newtonian fluids can be found in multiple industries such as chemical, oil and gas, food and
beverage, pharmaceutical, paints, and cosmetics. Therefore, having a model that allows the modeling
of particle flow in such fluids is crucial.
• 2-Way Fluent
• 1-Way Constant
This module does not work with the 2-Way Fluent Semi-Resolved nor the 1-Way LBM CFD Coupling
methods. This correlation is best suited for spherical particles in dilute flows.
Note:
The correlation implemented in this module includes a factor correction that accounts
for the non-Newtonian behavior of the fluid in the drag force, it does not account for
the particle shape nor the local volume fraction.
In addition, because drag laws must be specified per particle group in Rocky, after enabling this
module you will therefore need to select it when defining your CFD Coupling Particle Interactions
(Figure 8.1:Interactions tab options for the 1-Way Constant CFD Coupling mode when the Dhole, Ch-
habra & Eswaran Drag Law module is enabled. (p. 102)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 101
Dhole, Chhabra & Eswaran Drag Law Module
Figure 8.1: Interactions tab options for the 1-Way Constant CFD Coupling mode when the Dhole,
Chhabra & Eswaran Drag Law module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
102 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 8.2: Dhole, Chhabra & Eswaran Drag Law module options.
• Flow Behavior Index: This defines the power-law behavior index of the non-Newtonian drag law.
• Flow Consistency Index: This defines the power-law consistency index of the non-Newtonian drag
law.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 103
Dhole, Chhabra & Eswaran Drag Law Module
Tip: If you want to collect particle-fluid statistics, ensure you also enable the CFD Coupling Particle
Statistics module, and choose what kind of statistics you want to collect. (Refer to the About Fluid-
Related Statistics for Particles topic in the Rocky User Manual for details.)
2. From the Data panel, under Modules, select the new Dhole, Chhabra & Eswaran Drag Law
entry.
3. From the Data Editors panel, on the Dhole, Chhabra & Eswaran Drag Law tab, enter the values
you want.
4. Continue setting up the simulation as you normally would, with the following exceptions:
• You have selected a supported CFD Coupling mode. (Either 1-Way Constant, 1-Way Fluent
Steady State, or 2-Way Fluent).
• For the CFD Coupling mode you have enabled, on the Interactions tab for at least one
Particle, you have defined the Drag Law as Dhole, Chhabra & Eswaran Drag Law.
6. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related Property. Or if you chose to enable the CFD Coupling Particle Statistics module,
you may also make use of those additional fluid-particle properties.
The Dhole, Chhabra & Eswaran Drag Law module must be applied in simulations including non-Newto-
nian fluids whose behavior is adequately described by the power law, also known as Ostwald-de Waele
relationship. In this relationship, the shear stress, , relates to the shear rate, , according to [7] (p. 377):
(8.1)
where:
• is the power-law behavior index. This dimensionless parameter is defined by the user in the Dhole,
Chhabra & Eswaran Drag Law module.
• is the power-law consistency index, measured in , which is defined by the user in the Dhole,
Chhabra & Eswaran Drag Law module.
According to the power-law behavior index, fluids that obey are classified into three types:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
104 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• Newtonian fluids: Power-law fluids with n=1, where the shear stress is directly proportional to the
shear rate. The majority of the common fluids have this constant viscosity across all shear rates, in-
cluding water, oils, and gases.
• Pseudoplastic, or shear-thinning fluids: Power-law fluids with n<1, whose apparent viscosity decreases
with increasing shear rates. This behavior is commonly encountered in solutions or suspensions,
particularly when polymer solutions are used. Examples of pseudoplastic fluids are blood, milk, paints,
and drilling fluids.
• Dilatant, or shear-thickening fluids: Power-law fluids with n>1, and which shows a proportional increase
in the shear stress as a function of shear rate. One example is the mixture of cornstarch and water,
which hardens upon application of high enough forces, allowing people to run across large pools
filled with the mixture. Potential applications of dilatant fluids are liquid armor and smart speed
bumps.
When considering these types of fluids, the Reynolds number definition has to be modified to account
for the power-law index and the consistency index. The modified particle Reynolds number is
defined as follows [8] (p. 377):
(8.2)
where:
• is the fluid density, measured in [kg/m3], which is defined by the user on the CFD Coupling setup.
• is the fluid velocity vector, measured in [m/s], which is defined by the user on the CFD Coupling
setup when CFD Coupling mode is 1-Way Constant.
• is the equivalent diameter of the particle, measured in [m], which is defined as the diameter of a
sphere whose volume is equal to the particle's volume.
Based on numerical simulations of the flow of a power-law fluid past a sphere, Dhole, Chhabra & Eswaran
[6] (p. 377) obtained the following correlation for the drag coefficient, :
(8.3)
This correlation for the drag coefficient is valid for the following ranges:
(8.4)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 105
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
106 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 9: Disable Collision Outside Region Module
9.1. Usage details
The idea of the module is to apply a hybrid approach: outside the Regions of Interests particles are
tracked as body forces (such as gravity, eletrostatic and magnectic) and fluid forces (such as drag)
are computed, but no contact detection calculations involving these particles are performed, whereas
inside the ROIs the DEM model is turned on and collisions are solved.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 107
Disable Collision Outside Region Module
Figure 9.1: Snow accretion 1-way Rocky- Fluent coupling simulation using the Disable
Collision Outside Region module. Blue particles are enabled whereas light pink particles
are removed from contact detection calculations.
In this example, the performance gain is reached by disabling the contact detection for particles that
are located outside the blue box that encompasses the car. Outside the defined ROI, the snow particles
continue to move as body and fluid forces are computed but the contact detection calculations are
not performed. Near the car, inside the blue box, contact detection calculations are performed and
contact forces are accounted for when tracking the particles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
108 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Start Time: Defines the simulation time at which the module can begin disabling particles from
contact detection calculations.
• Contact Detection Zones: Defines the regions in which particles will be enabled for contact detec-
tion calculations.
In addition, one extra particle property is available as shown in Figure 9.3: Additional Particle Properties
created by the Disable Collision Outside Region module. (p. 110) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 109
Disable Collision Outside Region Module
Figure 9.3: Additional Particle Properties created by the Disable Collision Outside Region module.
• Disable Criterion: Criterion that classifies whether the particle is enabled or disabled from contact
calculations. If equals to 0.0, the particle is enabled and included in contact calculations performed
by the DEM solver. If equals to 1.0, the particle is disabled and excluded from such calculations.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
2. Ensure that the module is enabled. (From the Data panel, for setup details. select Modules and
then from the Data Editors panel, ensure the Disable Collision Outside Region checkbox is
enabled.)
3. Define the Start Time and the Contact Detection Zones. (From the Data panel, select Disable
Collision Outside Region, and then from the Data Editors panel, input the required values.)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
110 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
5. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related Property and may also make use of the additional module-related property.
When using this module, if the particle is outside of all the Regions of Interest, its Disable criterion, ,
receives a value of 1.0, eliminating the particle from the computation of contacts. The particle is auto-
matically re-enabled once it enters one of the selected Regions of Interest.
(9.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 111
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
112 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 10: Electromagnetic Coupling Module (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However,
your feedback will help us improve the overall quality of the product. We will not
guarantee that the projects using this beta feature will run successfully when the feature
is finally released so you may, therefore, need to modify the projects.
• Charging and discharging effects on particles due contact with boundary surfaces;
• Interaction force between point charge and a charged boundary (triangle point charge) or a con-
ductive boundary;
• Interaction force between dipole moment (either electric or magnetic) and a point cloud (repres-
enting an electric or magnetic field);
• Interaction torque between dipole moment (either electric or magnetic) and a point cloud (repres-
enting an electric or magnetic field).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 113
Electromagnetic Coupling Module (Beta)
• Short-Range Effect: When calculating particle-particle and particle-wall forces, only the nearby
particles and triangles surrounding a given particle will be considered. For more information, refer
to the Screening Distance Factor in Module options (p. 114)
• Point Charge: Charged particles and triangles are modeled as point charges concentrated at the
centroid of each particle or triangle.
• Dipole Moment: Dipole particles are represented by a single dipole moment, which has a constant
magnitude and orientation defined relative to the particle's reference frame. Currently, dipole
particles can interact only with electromagnetic fields, and interactions with charged particles or
other dipoles are not considered.
• 1-way approach: In this approach, the electromagnetic fields are unaffected by charged particles
or dipoles. Each particle interacts only with the closest point to its centroid. This formulation is best
suited for scenarios where the electromagnetic field can be considered nearly constant over the
particle span.
• Charging mechanism: Only collisions between particles and walls are considered for triboelectri-
fication. The charge of a particle is not influenced by neighboring charged particles, and the charge
density on the triangles remains constant throughout the simulation.
Important:
Figure 10.1: Electromagnetic Coupling Module options under the Model tab.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
114 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 10.2: Electromagnetic Coupling Module options under the Fields tab.
Figure 10.3: Electromagnetic Coupling Module options under the Advanced tab.
Model tab:
• Medium:
– Relative Permittivity: the ratio of the medium permittivity to the vacuum permittivity.
– Relative Permeability: the ratio of the medium permeability to the vacuum permeability.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 115
Electromagnetic Coupling Module (Beta)
– Dielectric Strength: the maximum electric field that a material can withstand without breaking
down or becoming conductive. The default value is set for air.
• Numerics:
– Screening Distance Factor: the ratio of the interaction distance to the particle size. It defines the
maximum distance that each particle will search for interactions.
Fields tab:
• EM Fields:
– EM Field: selects one of the pre-imported point clouds for use in the simulation. For more details
on the Point Cloud requirements, refer to About the Point Cloud field (p. 117).
– Start Time: defines when the selected point cloud will begin to influence particle interactions in
the simulation.
– Stop Time: defines when the selected point cloud will stop to influence particle interactions in
the simulation.
Advanced tab:
• 0-way Interactions:
– Coulomb Force Limiter: enables/disables charged particle-wall interaction limiter. Refer to Coulomb
Force Limiter (p. 127) for details. For more details, refer to Surface Triboelectrification (p. 127).
– Surface Triboelectrification: enables/disables particle charging on wall collisions. For more details,
refer to Surface Triboelectrification (p. 127).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
116 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• 1-way Interactions:
– Lorentz Force: enables/disables charged particle-point cloud interaction. For more details, refer
to Lorentz Force (p. 128).
– Electric Dipole: enables/disables electric dipole-point cloud interaction. For more details, refer to
Dipole on a Field (p. 129).
– Magnetic Dipole: enables/disables magnetic dipole-point cloud interaction. For more details,
refer to Dipole on a Field (p. 129).
• Electric field:
– Point cloud header:X Y Z Ex Ey Ez dEx_dx dEx_dy dEx_dz dEy_dx dEy_dy dEy_dz dEz_dx dEz_dy
dEz_dz
• Magnetic field:
– Point cloud header:X Y Z Bx By Bz dBx_dx dBx_dy dBx_dz dBy_dx dBy_dy dBy_dz dBz_dx
dBz_dy dBz_dz
Note:
the Field Gradient is required only for calculating the force on dipoles; it can be omitted
in simpler cases. For more details, refer to Dipole on a Field (p. 129).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 117
Electromagnetic Coupling Module (Beta)
10.1.4.2.1. Materials
A new sub-group called Electromagnetic Coupling has been added for the particle Material, as
shown in the figure below. The following section provides an explanation of the available para-
meter.
• Electromagnetic Coupling:
– Pauthenier Constant: is the material constant used to compute the saturation level, during
particle triboelectrification. For more details, refer to Surface Triboelectrification (p. 127).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
118 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Electromagnetic Coupling:
– Charge Dissipation Constant: the dissipation constant , defines the discharging rate when
in contact with the boundary. When , there is no charge dissipation. For more details,
refer to Surface Triboelectrification (p. 127).
– Charge Generation Constant: the generation constant , defines the charging rate when in
contact with the boundary with relative motion. A value of indicates that the particle
will be positively charged, whereas results in a negative charge. For more details, refer
to Surface Triboelectrification (p. 127).
10.1.4.2.3. Geometries
A new Modules sub-tab has been added for the Walls and Conveyor Templates, as shown below.
The following section provides an explanation of the available parameters.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 119
Electromagnetic Coupling Module (Beta)
• Electromagnetic Coupling:
– Conducting Wall: when enabled, charged particles will account for interaction forces due to
the image charge effect. For more details, refer to Image Charge Force (p. 127).
– Charge density: defines the triangle charge density used to compute particle-wall electrostatic
forces. For more details, refer to Coulomb Force Particle-Wall (p. 126).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
120 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Electromagnetic Coupling:
– Electric Dipole Moment Orientation : X: defines the x-component of the particle electric
dipole moment orientation vector .
– Electric Dipole Moment Orientation : Y: defines the y-component of the particle electric
dipole moment orientation vector .
– Electric Dipole Moment Orientation : Z: defines the z-component of the particle electric
dipole moment orientation vector .
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 121
Electromagnetic Coupling Module (Beta)
Figure 10.8: Electromagnetic Coupling Module options under Particle Inlets and Volumetric
Inlets.
• Electric Dipole Moment : Magnitude: defines the electric dipole moment magnitude.
• Electric Charge
• Remanence Magnitude
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
122 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
The Custom Inlet Input file must include each parameter name in the header, exactly as specified
above, to correctly identify the desired parameters. For more details, refer to Particle Inlets and
Volumetric Inlets (p. 122).
Note:
For additional details, you may refer to the Graphing (Plot or Histogram) a Data Set Within
Rocky topic in the Rocky User Manual and also to the View Geometries, Particles, Points,
and Fluids in 3D topic in the Rocky User Manual.
The list of available properties will depend on the active models, as described later in Technical de-
tails (p. 124).
• Wall Properties
– Charge Density: the amount of electric charge per unit area within a defined triangle
– Exchanged Charge: the total amount of electric charge that has been transferred across the triangle
during a time interval due to surface triboelectrification
– Force : Electric : Magnitude: the magnitude of the total electric force exerted on a triangle by a
particle charge
– Force : Electric : X/Y/Z: the components of the total electric force exerted on a triangle by a
particle charge, along the global frame of reference axis
• Particle Properties
– Electric Dipole Moment : Magnitude: the magnitude of the electric dipole moment for a particle
– Electric Dipole Moment : Orientation : X/Y/Z: the components of the particle electric dipole mo-
ment orientation vector, along the global frame of reference axis
– Force : Electric : Magnitude: the magnitude of the total electric force exerted on a particle
– Force : Electric : X/Y/Z: the components of the total electric force exerted on a particle, along the
global frame of reference axis
– Force : Magnetic : Magnitude: the magnitude of the total magnetic force exerted on a particle
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 123
Electromagnetic Coupling Module (Beta)
– Force : Magnetic : X/Y/Z: the components of the total magnetic force exerted on a particle, along
the global frame of reference axis
– Magnetic Dipole Moment : Magnitude: the magnitude of the magnetic dipole moment for a
particle
– Magnetic Dipole Moment : Orientation : X/Y/Z: the components of the particle magnetic dipole
moment orientation vector, in the global frame of reference
– Remanence : Magnitude: the residual magnetic flux density in the particle after an external
magnetic field has been removed
– Torque : Electric : Magnitude: the magnitude of the total electric torque exerted on a particle
– Torque : Electric : X/Y/Z: the components of the total electric torque exerted on a particle, along
the global frame of reference axis
– Torque : Magnetic : Magnitude: the magnitude of the total magnetic torque exerted on a particle
– Torque : Magnetic : X/Y/Z: the components of the total magnetic torque exerted on a particle,
along the global frame of reference axis
• Step 2: From the Data panel, under Modules, select the new entry for Electromagnetic Coupling for
the Module options, and continue setting up the simulation as you normally would.
Note:
• Step 4: When you are ready to post-process your simulation results, you can make use of the new
parameters on the Properties tab for the main Particles entity.
Note:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
124 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
outlines the requirements for enabling each of the available models. Each individual model is described
in detail in the following sections.
During the simulation initialization, the module will loop through the defined particle groups and
store the maximum particle size . The absolute screening distance is then computed from:
(10.1)
where:
For particle-particle interactions, the screening distance is compared against the magnitude of the
vector connecting their centroids, whilst for particle-triangle interactions, connects the particle
centroid to the triangle centroid.
The electrostatic force as defined by Coulomb's law is expressed by the following equations:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 125
Electromagnetic Coupling Module (Beta)
(10.2)
where:
where:
The electrostatic force as defined by Coulomb's law is expressed by the following equations:
(10.4)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
126 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
To mitigate this overestimation, a correction factor can be applied, based on the analytical solution
for the electric field produced by a charged disk. Assuming a point charge is positioned along the
axis of the disk at a distance from the disk's centroid, with the disk having a radius , the resulting
electric field is given by [58] (p. 380):
(10.5)
From :
(10.6)
(10.7)
(10.8)
where is the normal distance vector conecting the particle centroid to the closest triangle.
The sign of the particle's charge depends on the material properties of both the particle and the
surface, including factors such as the work function. This model encapsulates the direction of electri-
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 127
Electromagnetic Coupling Module (Beta)
fication in the charge generation constant, . A value of indicates that the particle will be posit-
ively charged, whereas results in a negative charge. The model also considers discharge currents
that lead to charge decay, represented by the charge dissipation constant, . When , there is no
source of charge leakage.
The expression for the charge transferred during a timestep is defined as follows:
(10.9)
(10.10)
In this model, the magnitude of the particle's charge can increase until it reaches the saturation level,
denoted as . This saturation level depends on the particle's size, material, and the surrounding
environment. The value of can be estimated using the Pauthenier equation [62] (p. 380):
(10.11)
(10.12)
where:
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
128 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
The dipole forces and torques are expressed by the following equations:
(10.16)
(10.17)
(10.18)
(10.19)
(10.20)
(10.21)
where:
The dipole moments are described as a function of the dipole moment orientation vector:
(10.22)
(10.23)
where
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 129
Electromagnetic Coupling Module (Beta)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
130 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 11: Erosion Model Module
11.1. Usage details
Unlike the surface wear modification model implemented by default in Rocky, this module does not
change the surface of the geometry. Rather, it calculates the cumulative eroded mass associating it
to a boundary-triangle property that can later be post-processed.
Selecting this module will not override any models or settings in the Rocky UI.
Important:
Even though this module supports CGM, special care is suggested when using this
feature with the Eroded Mass Analyses module. The use of CGM requires that the
general dynamics of the simulation is not changed with/without CGM. However, this
module has rather sensitive parameters to the dynamics of the simulation (especially
contact angles), which makes it harder to keep the same dynamics with the use of
CGM.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 131
Erosion Model Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
132 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Start Time: This defines the time you want to begin calculating the cumulative eroded mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 133
Erosion Model Module
• Finnie Model: When enabled, the eroded mass will be calculated using the Finnie erosion model.
– Angle of Maximum Erosion: This defines the impact angle of maximum erosion.
• Oka Model: When enabled, the eroded mass will be calculated using the Oka erosion model.
– Wall Material Vickers Hardness: This defines the wall material Vickers hardness.
• McLaury Model: When enabled, the eroded mass will be calculated using the McLaury erosion
model.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
134 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 135
Erosion Model Module
Figure 11.4: Additional Geometry Properties created by the Eroded Mass Analyses module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
These new Properties can help you visualize the evolution of the removed mass from each boundary
triangle due to collisions with particles, and are explained below:
• Surface Eroded Mass: Finnie Model: When Finnie Model is enabled in the Eroded Mass Analyses
module, this provides the erosion mass as calculated by the Finnie model.
• Surface Eroded Mass: McLaury Model: When McLaury Model is enabled in the Eroded Mass
Analyses module, this provides the erosion mass as calculated by the McLaury model.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
136 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Surface Eroded Mass: Oka Model: When Oka Model is enabled in the Eroded Mass Analyses
module, this provides the erosion mass as calculated by the Oka model.
• Surface Erosion Rate: Finnie Model: When Finnie Model is enabled in the Eroded Mass Analyses
module, this provides the erosion rate as calculated by the Finnie model.
• Surface Erosion Rate: McLaury Model: When McLaury Model is enabled in the Eroded Mass
Analyses module, this provides the erosion rate as calculated by the McLaury model.
• Surface Erosion Rate: Oka Model: When Oka Model is enabled in the Eroded Mass Analyses
module, this provides the erosion rate as calculated by the Oka model.
In addition to the properties, a new curve will also be created for each empirical model enabled in
the module, that can help you to evaluate the amount of eroded mass for each geometry present in
the simulation. After processing is complete, these new curves are available on the Curves tab when
an individual geometry component is selected on the Data panel under Geometries as shown in
Figure 5 Figure 11.5: Additional Geometry Curves created by the Eroded Mass Analyses module. (p. 137).
Figure 11.5: Additional Geometry Curves created by the Eroded Mass Analyses module.
2. From the Data panel, under Modules select Eroded Mass Analyses.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 137
Erosion Model Module
a. From the Finnie Model tab, enable the checkbox if you want to consider this empirical model
for the calculations and enter the proper values for the model parameters.
b. From the Oka Model tab, enable the checkbox if you want to consider this empirical model
for the calculations and enter the proper values for the model parameters.
c. From the McLaury Model tab, enable the checkbox if you want to consider this empirical
model for the calculations and enter the proper values for the model parameters.
d. Enable the checkbox for each empirical model you want calculated in the simulation.
4. Set up and process your simulation as you normally would, with the following exception:
5. After processing is complete, you may post-process the erosion data you collected by doing the
following:
a. From the Data panel, under Geometries, select the geometry component for which you would
like to analyze erosion.
b. From the Data Editors panel, select the Properties tab, and then under Transient, display in
a 3D View window or plot the Erosion... property you are interested in.
(11.2)
Where:
is the dimensionless mass, which can be defined as the ratio between the eroded mass and the
mass of the particle that collided with the boundary triangle.
is the velocity exponent, generally in the range 2.3 to 2.5 [rad] for metals.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
138 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(11.4)
Where:
is the dimensionless mass, which can be defined as the ratio between the eroded mass and the
mass of the particle that collided with the boundary triangle.
is calculated based on the combination of the two equations for the calculation of the angle function,
is automatically generated when ??=??lim.
(11.5)
(11.6)
Where:
is the dimensionless mass, which is the ratio between the eroded mass and the mass of the particle
that collided with the boundary triangle.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 139
Erosion Model Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
140 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 12: Frictional Heating Module
12.1. Usage details
Selecting this module will not override any models or settings in the Rocky UI.
In addition, this model works only in conjunction with the Rocky Thermal Model. Ensure that the
Thermal Model checkbox is enabled on the Thermal tab of the Physics entity and that all thermal-
related properties have been defined for your geometries, materials, and particle inputs before pro-
cessing your simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 141
Frictional Heating Module
• Start Time: This defines the time you want particles to begin being influenced by the heat friction
effects.
• Stop Time: This defines the time you want particles to stop being influenced by the heat friction
effects.
• Heat Generation Coefficient: This defines the fraction of the dissipated energy due to friction that
is transformed into heat.
In addition, an extra particle property is available if you enabled the module to process your simulation,
as shown in Figure 12.2:Frictional Heating module options. (p. 143) and explained below.
• Frictional Generated Heat: This provides the amount of heat generated on the particles due to
the collisions with other particles and walls.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
142 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
The property generated by this module will only appear if the Upscaling is disabled. Refer
to Rocky User Manual (About Meshed Particles Upscaling) for more details.
a. The Thermal Model checkbox is enabled on the Thermal tab of the Physics entity.
b. The module is enabled. (From the Data panel, select Modules and then from the Data Editors
panel, ensure the Frictional Heating checkbox is enabled.)
3. From the Data Editors panel, define the parameters as you want.
4. Set up your simulation as you normally would, with the following exceptions:
a. When setting up your Geometries, ensure that no active geometries have Motion Frames as-
signed.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 143
Frictional Heating Module
b. Ensure that you define all thermal-related properties for your geometries, materials, and particle
inputs.
6. When you are ready to analyze your simulation results, you may choose to make use of the fol-
lowing thermal-related Property:
Where:
is the heat, measured in [W], that will be added to the particle per unit time.
is the heat generation coefficient [-], which is an option defined by the user in the Frictional Heating
module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
144 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 13: Heat Generation Region Module
13.1. Usage details
Selecting this module will not override any models or settings in the Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 145
Heat Generation Region Module
• Start Time: This defines the time you want particles to begin being influenced by the heat gener-
ation region.
• Stop Time: This defines the time you want particles to stop being influenced by the heat generation
region.
• Region of Interest: The Region of Interest that you want to apply a heat source to the particles
that are inside it.
Range: [Any Region of Interest that has already been defined in the simulation]
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
146 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
a. The Thermal Model checkbox is enabled on the Thermal tab of the Physics entity.
b. The module is enabled. (From the Data panel, select Modules and then from the Data Editors
panel, ensure the Heat Generation Region checkbox is enabled.)
2. From the Data panel, under Modules, select Heat Generation Region.
b. Under Heat Generation Regions, click the Add button (green plus) to add as many rows
to the table as you want separate regions of interest to be heat generation regions.
c. For each separate row in the table, define the Heat Generation and the related Region
of Interest.
Set up your simulation as you normally would, with the following exceptions:
Tips:
• To insert a new row beneath the selected row, click the Insert button.
• To remove all rows in the table, click the Remove All button.
4. Set up your simulation as you normally would, with the following exceptions:
a. Ensure that you define all thermal-related properties for your geometries, materials, and particle
inputs.
6. When you are ready to analyze your simulation results, you may choose to make use of the fol-
lowing thermal-related Property:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 147
Heat Generation Region Module
Where:
is the heat transferred, measured in [W], per unit time to the particle while it is inside the defined
region.
is the projected area of the particle, measured in [m2], which is computed considering the
particle is a sphere.
is the heat source, measured in [W/m2], which is an option defined by the user in the Heat Gen-
eration Region module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
148 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 14: Heterogeneous Fiber Module (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However, your
feedback will help us improve the overall quality of the product. We will not guarantee that
the projects using this beta feature will run successfully when the feature is finally released
so you may, therefore, need to modify the projects.
Then, the properties for each pair of tag values are defined on the module (Figure 14.2: Pairs of Tags
table (p. 150)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 149
Heterogeneous Fiber Module (Beta)
The information of the custom fiber definition file and the pairs of tags table is mapped to the joints
of custom fiber, as shown in Figure 14.3: Tags mapped to joints (p. 150). Thus, the elastic, plastic and
discrete breakage models can be applied per joint.
The module uses the Bilinear Elastoplastic model for elastic and plastic behavior of the joints and the
Tensile or Shear Stress for joint breakage model. Both already present in Rocky for homogeneous fiber
modeling.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
150 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
If the element damping is enabled in the simulation, although the Rocky’s original Elastic Ratio is
not used in the calculations of joints forces, moments and breakage criterion calculations, this para-
meter is used in the fundamental frequency calculation, which is, in its turn, used on the element
damping coefficient calculation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 151
Heterogeneous Fiber Module (Beta)
Figure 14.6: Options for the Heterogeneous Custom Fiber (Beta) module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
152 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 14.7: Options for the Heterogeneous Custom Fiber (Beta) module.
• Show tags on particles: Defines whether tags given in the Custom Fiber Definition File will be saved
as particle properties for post-processing.
• First Element Tag: Sets the first tag value that defines a pair of tags.
• Second Element Tag: Sets the second tag value that defines a pair of tags.
• Elastic Ratio: Defines the Elastic Ratio of a joint formed by elements with tag values defined by
the First Element Tag and Second Element Tag values.
• Damping Ratio: Defines the Joint Damping Ratio of a joint that is formed by elements with tag
values defined by the First Element Tag and Second Element Tag values. This value is not used
for element damping.
• Plastic Ratio: Defines the Plastic Ratio of a joint formed by elements with tag values defined by
the First Element Tag and Second Element Tag values.
• Torsion Angle Limit: Defines the Torsion Angle Limit of a joint formed by elements with tag values
defined by the First Element Tag and Second Element Tag values.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 153
Heterogeneous Fiber Module (Beta)
• Bending Angle Limit: Defines the Bending Angle Limit of a joint formed by elements with tag values
defined by the First Element Tag and Second Element Tag values.
• Tensile Stress Limit: Defines the Tensile Stress Limit of a joint formed by elements with tag values
defined by the First Element Tag and Second Element Tag values.
• Shear Stress Limit: Defines the Shear Stress Limit of a joint formed by elements with tag values
defined by the First Element Tag and Second Element Tag values.
In addition, one extra particle property is available as shown in Figure 14.8: —Additional particle
properties created by the Heterogeneous Custom Fiber (Beta) module. (p. 154) and explained below.
Figure 14.8: —Additional particle properties created by the Heterogeneous Custom Fiber (Beta)
module.
• Tag: The tag defined in the Custom Fiber Definition File. This property will only appear if the Show
tags on particles module property is enabled and if the Meshed Particles Upscaling is disabled.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
154 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Heterogeneous Custom Fiber (Beta) checkbox is enabled.)
3. From the Data panel, under Modules, select Heterogeneous Custom Fiber (Beta).
4. From the Data Editors panel, within the Element Pair Properties list, click on the Add button
to add all the possible pairs of tags combinations, setting up all the parameters needed.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 155
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
156 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 15: Leblicq Joint Model Module
15.1. Usage details
In this version of the Leblicq Joint Model module (Figure 15.1: Options for the Leblicq Joint Model
module. (p. 157)), you must include your own Loading Stiffness and Unloading Stiffness values as
separate input files (Figure 15.2: Example paths to stiffness input files. (p. 158)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 157
Leblicq Joint Model Module
Enabling this module adds a new Leblicq Joint Model option to the Joint Model list on the Particle
| Composition tab when a multi-element particle shape is defined (Figure 15.3: Joint Model options
for a multi-element Fiber when the Leblicq Joint Model module is enabled. (p. 159)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
158 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 15.3: Joint Model options for a multi-element Fiber when the Leblicq Joint Model module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 159
Leblicq Joint Model Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
160 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
In brief, the values in the .csv files represent the lower triangular part of a square matrix, which define
the stiffness as a function of the deformation angle and the maximum deformation angle.
In order to obtain those values, you must first perform three-point bending tests on your fiber, do
force and deformation angle measurements, and then use the procedure in the referenced paper to
obtain the resulting stiffness tables.
This experimental work is required before using this module as a set of tables is specific to a fiber of
a given shape and material.
For illustration purposes, example loading and unloading stiffness tables were provided with this
module. As shown in Figure 15.6: Folder location of the example .csv files. (p. 162), these two .csv files
are located in the same folder as the compiled Leblicq Joint Model module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 161
Leblicq Joint Model Module
• Loading Stiffness Input File: The full path to the .csv file that defines the loading stiffness
values.
Range: [.csv files formatted according to the Loading Stiffness procedure defined in the refer-
enced paper.]
• Unloading Stiffness Input File: The full path to the .csv file that defines the unloading stiffness
values.
Range: [.csv files formatted according to the Unloading Stiffness procedure defined in the ref-
erenced paper.]
Once the module is enabled and your input files specified, you are able to select the new Leblicq
Joint Model option from the Joint Model list on the Particle | Composition tab.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
162 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical Details
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Leblicq Joint Model checkbox is enabled.)
3. From the Data panel, under Modules, select the new Leblicq Joint Model entry, and then from
the Data Editors panel, enter the two separate paths to the .csv files you want to use.
When setting up your Particles, ensure that for any Particle sets to which you want the Leblicq
Joint Modelto apply, you do all of the following:
b. From the Composition sub-tab, enable Multiple Elements from the Composition list, and
then from the Joint Model list, select Leblicq Joint Model.
If you choose to enable this module, you must have at least one Particle set that makes use
of it. Otherwise, you should disable the module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 163
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
164 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 16: Leeds Contact Model Module
16.1. Usage details
The Normal Force is used to calculate the normal components of the contact forces. The Adhesive
Force defines how adhesion between materials is calculated, and will therefore affect Material-to-
Material interactions on the Materials Interactions step.
Selecting this module overrides what you have set for both Normal Force and Adhesive Force in
the Physics | Momentum tabFigure 16.1: Physics | Momentum tab options when the Leeds Contact
Model module is enabled. (p. 166)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 165
Leeds Contact Model Module
Figure 16.1: Physics | Momentum tab options when the Leeds Contact Model module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
166 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
However, once the module is enabled, it does provide additional options in the Materials Interactions
section of the UI, as explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 167
Leeds Contact Model Module
Figure 16.3: Materials Interactions options when the Leeds Contact Model module is enabled.
• Damping Ratio: This defines the amount of extra viscous dissipation added in reloading
periods, in order to suppress small oscillations when a particle is near to a still position, as
shown in Equation 16.7 (p. 171).
Range: [0, 1]
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
168 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• Surface Energy: This defines the work per area unit expended in separating two contacting
surfaces. In the Leeds model, it is used to calculate the pull-off force, as shown in Equa-
tion 16.2 (p. 170).
2. Set up your simulation as you normally would, with the following exceptions:
• This model overrides both the Normal Force and Adhesive Force parameters on the
Physics | Momentum tab so you will be unable to set those parameters when this module
is enabled.
• This model is incompatible with Coarse Grain Modeling. (From the Physics | Coarse-Graining
tab, ensure that the Coarse Grain Modeling checkbox is cleared.)
b. When setting up Materials Interactions, ensure that for each material-material pair you are
making use of in the simulation, you define what you want for Damping Ratio and Surface
Energy.
Figure 16.4: Schematic diagram of a typical normal force overlap response of the Leeds Contact Mod-
el. (p. 170) shows schematically a typical normal force overlap response of the Leeds Contact Model,
which is determined by the piecewise function of the normal overlap ,
(16.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 169
Leeds Contact Model Module
Figure 16.4: Schematic diagram of a typical normal force overlap response of the Leeds Contact
Model.
In Figure 16.1: Physics | Momentum tab options when the Leeds Contact Model module is enabled. (p. 166),
and are the plastic and elastic stiffnesses, as they are called in the original formulation of the
model [14] (p. 377). In Rocky, and are taken to be equal to the loading stiffness, , and the
unloading stiffness, , respectively, defined for the hysteretic linear spring model.
Moreover, in Equation 16.1 (p. 169) is the elastic pull-off force predicted by the JKR theory, which is
given by:
(16.2)
where is the effective or equivalent radius, whereas is the module parameter listed as Surface
Energy in the Rocky UI.
The overlap values that delimit the different stages in Equation 16.1 (p. 169) are given by:
(16.3)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
170 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(16.4)
(16.5)
where is the maximum overlap reached during the loading, as can be seen in Figure 16.4: Schem-
atic diagram of a typical normal force overlap response of the Leeds Contact Model. (p. 170).
Additionally, the horizontal coordinate of point A, at which the adhesive force begins to exist, is given
by:
(16.6)
On the other hand, the contact is considered to break at point F, when the adhesive force reduces to
, where is the adhesive force at point E, corresponding to the overlap given by Equa-
tion 16.5 (p. 171). If before reaching point F the decreasing in the overlap is reversed and a reloading
process begins, the force will follow the same unloading path, but in the opposite direction, as indicated
by the arrows in Figure 16.4: Schematic diagram of a typical normal force overlap response of the Leeds
Contact Model. (p. 170). Moreover, if in that reloading process the overlap increases beyond , starting
at point D, the force will follow an extension of the original loading path C-D, with slope given by .
The Leeds Contact Model is closely related to the JKR model, in fact, the initial part of the force response
in the Leeds Contact Model is designed to be a linearized version of JKR, with the same pull-off force.
However, as the Leeds Contact Model is based on the hysteretic linear spring model, the main mechanism
of energy dissipation is plastic rather than viscous as in the JKR model implemented in Rocky. Still, in
some cases it can be desirable to add an extra viscous dissipation mechanism, in order to suppress low-
amplitude oscillations when particles are in a near-static state. With that purpose, a viscous force is
added in Rocky, which is activated only when a reloading period begins. This extra viscous force is added
also if a unloading begins before reaching point C. Afterwards, if the overlap gets bigger than , the
viscous force is deactivated. This force is given by:
(16.7)
where is the Damping Ratio, which is a user input in the Rocky UI. Additionally, is the equivalent
or effective mass, whereas is the time derivative of the overlap.
The normal contact force defined in the Leeds Contact Model can be considered as composed by a re-
pulsive force force plus an adhesive one. Internally in Rocky, the data structures associated to those
types of forces are different, therefore, they must be calculated separately. This separation is useful also
for the computation of some statistics in which repulsive forces and adhesive forces are accounted
separately. Then, as shown schematically in Figure 16.5: Normal force in Leeds Contact Model as a
composition of repulsive and adhesive forces. (p. 172), the normal force for the Leeds model in Rocky
is considered to be the sum of a repulsive force (in purple), which coincides with the force predicted
by the hysteretic linear spring model, and an adhesive force (in blue).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 171
Leeds Contact Model Module
Figure 16.5: Normal force in Leeds Contact Model as a composition of repulsive and adhesive
forces.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
172 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 17: Liquid Absorption Model Module
17.1. Usage details
Selecting this module will not override any models or settings in the Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 173
Liquid Absorption Model Module
• Activate Mass Increment: When enabled, the liquid mass removed from the liquid film will be
added to the particle mass. When cleared, the liquid mass will still be removed but the removed
mass will not be added to the particle mass.
• Absorption Rate Coefficient: This defines the absorption rate coefficient at which it will define
the liquid mass absorption rate of the liquid film.
• Start Time: This defines the time you want particles to begin being influenced by the liquid adsorp-
tion model. Because the liquid in the film will either be removed or absorbed by the particle over
time, the effective Stop Time will be different for each particle.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
174 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
For the Inputs step, a new Modules sub-tab will appear to help you define the properties of the
particles, as explained below and shown in Figure 17.2:Liquid Absorption module options for a
Continuous Injection Input. (p. 175).
Note:
This sub-tab will appear only for the Continuous Injection and Volume Fill inputs. For
Custom Inputs, there is not a Modules sub tab; rather, the Module s parameters will
appear on the main tab for Custom Inputs.
Figure 17.2: Liquid Absorption module options for a Continuous Injection Input.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 175
Liquid Absorption Model Module
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab for the Input.
• Liquid Mass: This defines the mass of the liquid film that will be applied to each particle within
the selected Particle set when it enters into the simulation. Important: Every particle in the set,
regardless of the particle size or volume, will contain the same amount of liquid mass and mass
increment when injected.
• Mass Increment: This defines the mass increment that will be applied to each particle within
the selected Particle set.
Normally, this value should be set to zero, since it will be the Liquid Absorption module itself that
will calculate the mass increment during the simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
176 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• Liquid Mass: This provides the mass of the liquid film around each individual whole particle or
fragment. This can help you to visualize the evolution of the particle's liquid content, as the liquid
mass will tend to decrease over time due the absorption process being simulated.
• Mass Increment: This provides the increment mass in relation to the original mass of each individual
whole particle or fragment. This can help you to visualize the amount of mass that was absorbed
by each particle over time.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 177
Liquid Absorption Model Module
3. From the Data Editors panel, define the parameters as you want
4. Set up your simulation as you normally would, with the following exception:
• For the Inputs step, ensure that for each Input and Particle set you define, you also define
the Liquid Mass and Mass Increment values you want.
6. When you are ready to analyze your simulation results, do the following:
b. From the Data Editors panel, select the Properties tab, and then display in a 3D View window
or plot any of the following Properties:
• Liquid Mass
• Mass Increment
Where:
• is the absorption rate coefficient [1/s], which is a user-defined option in the Liquid Absorption
module.
• is the instantaneous mass of the liquid film around a particle, as measured in [kg].
If the Activate Mass Increment checkbox is enabled in the Liquid Absorption module, the mass of the
liquid film will be reduced by and will be added to the particle s mass. If the checkbox remains
cleared, the mass of the liquid film will still be reduced by but that will not be added to the particle's
mass.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
178 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 18: Liquid Bridge Module
18.1. Usage details
Figure 18.1: Example mixing process using the Liquid Bridge Model
If you want to include the effect of liquid films around particles in the simulation, you can choose to
enable the Liquid Bridge Model module (Figure 18.2:Liquid Bridge Model module options. (p. 180)) prior
to processing your simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 179
Liquid Bridge Module
Enabling the Liquid Bridge Model module overrides the Adhesive Force model that is defined on the
Physics | Momentum tab, as shown in Figure 18.3:Liquid Bridge Model overriding the Adhesive Force
setting. (p. 181).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
180 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 18.3: Liquid Bridge Model overriding the Adhesive Force setting.
Only particles in Rocky are considered to have liquid films associated with them. However, it is still
possible for a particle to create a liquid bridge with a boundary. In these cases, only the particle
contributes liquid to the bridge's formation, and liquid bridge forces are considered when computing
the interaction between particles and boundaries.
Note: The liquid mass is not included in the particle mass when solving motion equations for the
particle. For example, if a dry particle has 0.1 kg and the liquid film around it has 0.001 kg, the motion
equations for the particle will be solved considering only a mass of 0.1 kg.
Additionally, if the particle groups of the simulation have an aspect ratio smaller than 0.9 or bigger
than 1.1, the forces calculated by this module may be inaccurate.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 181
Liquid Bridge Module
• Liquid Viscosity: This defines the viscosity of the liquid film around particles.
• Liquid Density: This defines the density of the liquid film around particles.
• Surface Tension: This defines the surface tension of the liquid film around particles.
The default value of 0.072 N/m is equal to the surface tension of water.
• Materials Interactions
The additional options on the Materials Interaction tab (Figure 18.4:Liquid Bridge Model Materials
Interactions options. (p. 183)) help you define how the liquid bridge will act when it is formed
between particles, or between particles and boundaries.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
182 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Bridge Volume Fraction: This defines the fraction of the liquid bridge volume that the
particle pair will each contribute to when they are close enough for their liquid films to form
a bridge. Each particle will contribute to the bridge with a liquid volume equal to the Bridge
Volume Fraction times the volume of the film around it. For particle-to-boundary interaction
pairs, only the particle liquid contributes to the formation of the bridge.
Range: [0-1]
• Contact Angle: This defines the angle between the surface of the liquid bridge and the
particle (and/or boundary) surfaces between which it formed. It is used to compute the
critical rupture distance for the liquid bridge between particles, and is also used to calculate
the capillary force.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 183
Liquid Bridge Module
• Minimum Separation Ratio: This defines the smallest (minimum) separation distance below
which the capillary and viscous liquid bridge forces will remain with a constant value, which
is computed using the minimum distance. The minimum separation distance is obtained by
multiplying the Min. Separation Ratio and the radius of the bigger particle of the pair (or
just the radius of the particle when it is a particle-boundary interaction). Above this separation
distance, the liquid bridge force is computed using the real separation values. This numerical
artifact was implemented to avoid extremely high values when the separation is close to
zero.
Figure 18.5: Options for a Volumetric Inlet (left) and a Particle Inlet (right).
Within the table on the Modules sub-tab, one parameter (column) can be set per Particle Group
(row), as defined below:
• Particle: For each row, lists the Particle Group that is defined on the Particles sub-tab.
• Liquid Mass: This defines the liquid mass that will be applied to each particle within the
selected Particle Group when it enters into the simulation. Important: Every particle in the
set, regardless of the particle size or volume, will contain the same amount of liquid mass
when it enters into the simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
184 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
For Particle Custom Inlets, a new Liquid Bridge Model section will appear to help you define the
liquid film input properties of the particles, as shown in Figure 18.6: Options for a Particle Custom
Inlet. (p. 185). The only parameter in this section is Liquid Mass, which is defined exactly the same
as indicated in the preceding section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 185
Liquid Bridge Module
film around a particle. This can help you to visualize the evolution of the particle's liquid content, as
the liquid can be transferred from particle to particle during each collision.
Or, you can choose to display this property graphically in a 3D View window.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
2. From the Data panel, under Modules, select the new Liquid Bridge Model entry.
3. From the Data Editors panel, on the Liquid Bridge Model tab, enter the values you want.
4. Continue setting up the simulation as you normally would, with the following exceptions:
a. When setting up your Materials Interactions, be sure to define the Liquid Bridge Model
values you want for each material interaction pair.
b. When setting up your Inlets and Outlets, be sure to define the parameters on the Modules
tab (for Particle Inlet and Volumetric Inlet) or in the Liquid Bridge Model section (for Particle
Custom Inlets).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
186 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical Details
6. When you are ready to post-process your simulation results, you can make use of the new para-
meter on the Properties tab for the main Particles entity
A liquid bridge exerts an attractive capillary force over the particles, which is calculated in Rocky according
to a correlation derived by Mikami et al.[15] (p. 377). For a bridge between two particles, the magnitude
of this force is given by the expression:
(18.1)
In this expression, is the surface tension, is the geometric mean of the radio of the two
bridged particles, and is the dimensionless separation distance between the particles. Additionally,
A, B and C are coefficients that depend on the dimensionless volume of the liquid bridge, and the
contact angle, , through the following expressions:
(18.2)
(18.3)
(18.4)
The dimensionless parameters included in the expressions above are defined as follows:
(18.5)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 187
Liquid Bridge Module
(18.6)
where is the separation distance between the bridged particles, as depicted in Figure 18.9: Geometry
of a liquid bridge between a particle and a boundary. (p. 188), and is the volume of the liquid contained
in the bridge. Moreover, is defined as:
(18.7)
where and are the radii of the two bridged particles. For non-spherical particles, the radius of each
particle is taken one half of its sieve size.
(18.10)
where and are the liquid masses and and are the densities of the liquid films around the
two particles, respectively.
Additionally, is a user-defined parameter listed as Bridge Volume Fraction in the Materials Inter-
actions section of the Rocky UI.
Mikami et al.[15] (p. 377) also provide expressions for the coefficients in Equation 18.1 (p. 187) correspond-
ing to a liquid bridge between a particle and a boundary, as the one depicted in Figure 18.9: Geometry
of a liquid bridge between a particle and a boundary. (p. 188). Those expressions are the following:
(18.11)
(18.12)
(18.13)
In this specific case, both and will be equal to the radius of the particle interacting with the
boundary, . In the current Liquid Bridge Model in Rocky, boundaries do not have associated a liquid
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
188 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical Details
film. Therefore, in particle-boundary bridges only particles will contribute with liquid to the formation
of the bridge. This means that the bridge volume in those cases will be given by:
(18.14)
where is the volume of the liquid film around the particle in a particle-boundary bridge interaction.
In Rocky, the amount of liquid carried by each particle is tracked though a special particle property
called Liquid Mass. This property can be visualized and post-processed after a simulation has been
completed. Additionally, users can specify the value of the initial mass of liquid carried by the particles
when they enter to the domain. At any given time during the simulation, the volumes of the liquid
films that appear on equations Equation 18.7 (p. 188) and Equation 18.14 (p. 189) will be computed using
the liquid mass of each particle and the liquid density, which is a property that users must specify.
The mass of the liquid films attached to particles can vary throughout a simulation because it can be
redistributed when a bridge interaction between particles ends. The liquid redistribution scheme imple-
mented in Rocky is quite simple. When a liquid bridge between two particles is broken, half of the mass
of liquid that was contained in the bridge is added to the liquid mass associated to each particle. On
the other hand, when a particle-boundary bridge is ruptured, all the liquid mass goes to the particle,
since boundaries are not able to retain liquid films in the current model implemented in Rocky.
The criterion for the rupture of a liquid bridge is based only on the value of the separation distance, .
When that distance exceeds a critical value , the liquid bridge will break and, therefore, the bridge
interaction will finish. Rocky uses the following expressions for the dimensionless critical distance,
provided by Mikami et al.[15] (p. 377):
(18.15)
for particle-to-particle liquid bridges, and
(18.16)
for particle-to-boundary liquid bridges.
In addition to capillary forces, the Liquid Bridge Model in Rocky includes viscous forces that oppose the
relative motion between particles or between a particle and a boundary, joined by a bridge. The expres-
sions used for those viscous forces, derived from the lubrication theory, are taken from Nase et al.
[16] (p. 377). Unlike the capillary force that acts always in the normal direction, the viscous forces have
a normal component and a tangential component. The expressions for those components are, respect-
ively:
(18.17)
(18.18)
where:
Values of fluid properties, such as viscosity and density, as well as other properties such as surface
tension and contact angle, must be specified by users through the Rocky UI.
(18.19)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 189
Liquid Bridge Module
• and are the normal and tangential components of the relative velocity, in relation to the mid
plane between the two particles or the particle and the boundary.
The minus sign in equations Equation 18.17 (p. 189) and Equation 18.18 (p. 189) indicates that the viscous
forces always oppose to the relative motion, as already mentioned. As a consequence, the normal viscous
force will behave as a repulsive force when particles approach each other, and as an attractive force
when they move away.
As the separation distance decreases to zero, both components of the viscous force tend to infinity.
In order to avoid stability issues caused by very large forces, Rocky uses a cut-off value for the separation
distance. This cut-off value is determined by an input parameter listed in the Rocky UI as Minimum
Separation Ratio. So, if is the value of this parameter, the cut-off value for the separation distance
will be given by:
(18.21)
Whenever a separation distance in a liquid bridge falls below , this latter value will be used in
equations Equation 18.19 (p. 189) and Equation 18.5 (p. 187) instead of the actual value, .
(18.22)
Where is the CGM scale factor, the original contact force, and is the CGM scaled force.
A liquid bridge exerts an attractive capillary force over the particles, which is calculated in Rocky according
to a correlation derived by Mikami et al. For a bridge between two particles, the magnitude of this force
is given by the expression:
(18.23)
In this expression, is the surface tension, the geometric meanof the radii of the two bridged
particles, and is the dimensionless separationdistance between the particles. Additionally,A,B and C
are coefficients thatdepend on the dimensionless volume of the liquid bridge, and the contactangle,
through the following expressions:
(18.24)
The dimensionless parameters included in the expressions above are defined as follows:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
190 of ANSYS, Inc. and its subsidiaries and affiliates.
Coarse Grain Modeling (CGM) for Liquid Bridge
(18.25)
Where is the volume of the liquid contained in the bridge. Moreover, is defined as:
(18.26)
Where and are the radii of the two bridged particles. For non-spherical particles, the radius of each
particle is taken one-half of its sieve size.
Once the CGM radius is scaled by , if one scales the by the constant A,B, and C are
equal to the original particle without scaling.
Finally, if one multiple the surface tension by the , one gets the CGM Force as:
(18.27)
Therefore, in order to make the Liquid Bridge Model compatible with CGM one just needs to cubic scale
the volume of liquid bridge , and scale the surface tension by the factor.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 191
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
192 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 19: Liquid Drying by Heat Flux Module
19.1. Usage details
• The liquid contains a certain amount of solid solute that you specify, uniformly dissolved. When
the solvent evaporates, the solute is deposited on the surface of the particle.
• The liquid film evaporation occurs when the particle temperature is above the liquid boiling
point and/or receives a certain amount of heat.
• Both the particle-particle and particle-boundary conductive heat transfer and fluid-particle
convective heat transfer are taken into account.
• The liquid film has always the same temperature as the particle.
• The liquid film mass is small and does not add to the particle mass (no impact on the particle
motion).
• No energy is required to increase the liquid film temperature (no sensible heat exchange).
• Energy is only consumed in the evaporation of the liquid film (latent heat exchange).
• When a certain liquid film mass evaporates, the corresponding solute mass is uniformly depos-
ited around the particle.
• The particle's volume and mass are increased according to the amount of solute deposited.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 193
Liquid Drying by Heat Flux Module
Figure 19.1: Options for the Liquid Drying by Heat Flux module.
• Boiling Point: This defines the boiling temperature of the liquid film.
The default values of Boiling Point and Heat of Evaporation correspond to water
• Heat of Evaporation: This defines the amount of heat necessary to evaporate a unit mass of
liquid.
• Solute Density *: This defines the density of the solute dissolved in the liquid.
Asterisks (*) following a parameter name indicate that the parameter and value entered are
shared with other modules or UI entities that also make use of this parameter. Shared paramet-
ers may have only one value stored for it; the value you enter most recently will override any
shared values set earlier in other locations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
194 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Solute Density is shared with the Liquid Drying Zone module, if it is active in the project.
• Solute Mass Fraction *: This defines the ratio of the mass of solute to the mass of the liquid
solution. In turn, the mass of the liquid solution is the sum of the masses of the solute and
the solvent.
Solute Mass Fraction is shared with the Liquid Drying Zone module, if it is active in the project.
Figure 19.2: Options for a Continuous Injection Input (left) and a Volume Fill Input (right).
The meaning of the parameters appearing in the columns of the table inside the Modules tab is
the following:
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab.
• Liquid Mass: This defines the liquid mass that will be applied to each particle within the
selected Particle set when it enters into the simulation. Important: Every particle in the set,
regardless of the particle size or volume, will contain the same amount of liquid mass when
it enters into the simulation.
• Mass Increment: The initial value of the mass of solute deposited on the surface of particle.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 195
Liquid Drying by Heat Flux Module
Since usually one of the objectives of using a Liquid Drying module is to determine the
amount of solute deposited on the particles during a simulation, the initial values of Mass
Increment and Volume Increment normally should be set to zero.
• Volume Increment: The initial value of the volume of solute deposited on the surface of
particle.
The definition of the three parameters in this section is exactly the same specified above in the
section about the other two kinds of particle inputs.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
196 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Liquid Mass: This provides the mass of the liquid film around each individual whole particle or
fragment. This can help you to visualize the evolution of the particle's liquid content during a
simulation, as the drying process will tend to reduce this content over time.
• Mass Increment: This provides the mass of the solute accumulated over the particle's surface, due
to the gradual evaporation of the liquid film. During a simulation, both the Particle Mass and the
Particle Volume properties will be increased over time by the values of the Mass Increment and the
Volume Increment, respectively, in relation to their intrinsic values.
• Volume Increment: This provides the volume of the solute accumulated over the particle s surface,
due to the gradual evaporation of the liquid film.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 197
Liquid Drying by Heat Flux Module
Figure 19.4: Particle Properties added by the Liquid Drying by Heat Flux module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
2. From the Data panel, under Modules, select Liquid Drying by Heat Flux.
3. From the Data Editors panel, define the values you want for the parameters displayed.
4. Set up the simulation as you normally would, with the following exception:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
198 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• When setting up your Inputs, be sure to define the parameters on the Modules tab (for Con-
tinuous Injection and Volume Fill Inputs) or in the Liquid Drying by Heat Flux section (for Custom
Inputs).
6. When you are ready to post-process your simulation, you can make use of the new parameters
on the Properties tab for the main post-processing details. Particles entity.
The evaporation of the solvent occurs when the particle's temperature is above or equal to the the
boiling temperature of the solvent. The mass of the solvent evaporated during a simulation timestep
will be given by:
(19.1)
where is the heat transferred from the particle to the liquid film during the timestep, while
is the latent heat of evaporation of the solvent, whose value is specified through the module parameter
Heat of Evaporation.
The heat that the particle is able to transfer to the liquid film is composed of two sources, as expressed
in this equation:
(19.2)
where:
is the heat transfer rate to the particle . This includes both conductive heat transfer through contacts
or joints, and convective heat transfer from CFD coupling.
and are the mass and the specific heat of the particle, respectively.
and are, respectively, the instantaneous temperature of the particle and the boiling temper-
ature of the solvent in the liquid film. The value of this last temperature is specified through the module
parameter Boiling Point. The model does not consider any thermal resistance between the particle
and the liquid film. Therefore, the particle may reduce its temperature up to and transfer the
corresponding amount of sensible heat to the liquid film instantaneously, as expressed in equation
Equation 19.2 (p. 199).
When the mass of solvent given by Equation 19.1 (p. 199) evaporates, the corresponding mass of solute
that will be deposited on the particle's surface will be given by:
(19.3)
where is the value specified through the module parameter Solute Mass Fraction. In order to
take into account this added mass, its value will be summed to the current value of the particle's
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 199
Liquid Drying by Heat Flux Module
property Mass Increment. At any given time, the mass of a particle will be the sum of its original mass
(the mass the particle had at the time it entered into the simulation) and the instantaneous value of
Mass Increment.
Additionally, the corresponding volume of the solute deposited on the particle's surface will be calculated
as:
(19.4)
where is the specified value of the Solute Density. This value is added to the current value of
the particle's property Volume Increment. Rocky will use the values of both the mass and volume in-
crements for updating the instantaneous values of the particle's mass, volume and size.
Finally, the mass of the liquid solution that will be removed from the liquid film as a consequence of
the solvent evaporation and the solute deposition is given by:
(19.5)
The process of evaporation described in this section will cease if the mass of the liquid film reduces to
zero or the temperature of the particle falls below .
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
200 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 20: Liquid Drying Region Module
20.1. Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 201
Liquid Drying Region Module
• Drying Liquid Mass Rate: This defines the constant rate at which the mass liquid will be re-
moved from the particles inside the drying region.
• Solute Density: This defines the density of the solid dissolved in the liquid around particles.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
202 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Solute Mass Fraction: This defines the ratio of the mass of solute to the mass of the liquid
solution. In turn, the mass of the liquid solution is the sum of the masses of the solute and
the solvent.
• Start Time: This defines the time you want particles in the regions of interest to begin losing
liquid film.
• Stop Time: This defines the time you want particles in the regions of interest to stop losing
liquid film.
• Region of Interest: The Region of Interest in which you want liquid film evaporation to occur.
Range: [Any Region of Interest that has already been defined in the simulation]
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 203
Liquid Drying Region Module
Figure 20.2: Options for a Continuous Injection Input (on the left) and a Volume Fill Inputs
(on the right).
Within the table on the Modules sub tab, three parameters (columns) can be set per Particle set
(row), as defined below:
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab.
• Liquid Mass: This defines the liquid mass that will be applied to each particle within the
selected Particle set when it enters into the simulation. Important: Every particle in the set,
regardless of the particle size or volume, will contain the same amount of liquid mass when
it enters into the simulation
• Mass Increment: The initial value of the mass of solute deposited on the surface of particle.
Since usually one of the objectives of using a Liquid Drying Region module is to determine
the amount of solute deposited on the particles during a simulation, the initial values of
Mass Increment and Volume Increment normally should be set to zero.
• Volume Increment: The initial value of the volume of solute deposited on the surface of
particle.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
204 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
For Custom Inputs, a new Liquid Drying Region section will appear to help you define the liquid
film input properties of the particles, as shown in Figure 20.3: Options for a Custom Input. (p. 206).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 205
Liquid Drying Region Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
206 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
The definition of the three parameters in this section is exactly the same specified above in the
section about the other two kinds of particle inputs.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 207
Liquid Drying Region Module
Figure 20.4: Particle Properties added by the Liquid Drying Region module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
208 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 209
Liquid Drying Region Module
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• Liquid Mass: This provides the mass of the liquid film around each individual whole particle
or fragment. This can help you to visualize the evolution of the particle's liquid content during
a simulation, as the drying process will tend to reduce this content over time.
Asterisks (*) following a parameter name, Figure 20.3: Options for a Custom Input. (p. 206), indicate
that the parameter and value entered are shared with other modules or UI entities that also
make use of this parameter. Shared parameters may have only one value stored for it; the
value you enter most recently will override any shared values set earlier in other locations.
• Mass Increment: This provides the mass of the solute accumulated over the particle's surface,
due to the gradual evaporation of the liquid film.
During a simulation, both the Particle Mass and the Particle Volume properties will be increased
over time by the values of the Mass Increment and the Volume Increment, respectively, in relation
to their intrinsic values. In addition, the Particle Size property will also be increased as a result
of the volume increment.
• Volume Increment: This provides the volume of the solute accumulated over the particle's
surface, due to the gradual evaporation of the liquid film.
2. From the Data panel, under Modules, select Liquid Drying Region.
b. Under Liquid Drying Regions, click the Add button (green plus) to include as many separate
rows for the number of regions as you want, and then for each separate row, define the Region
of Interest.
You may add as many rows as regions you want in the simulation, and you may also remove
rows, if necessary. These operations must be performed using the buttons on the top of the
table.
4. Set up the simulation as you normally would with the following exception:
• When setting up your Inputs, be sure to define the parameters on the Modules tab (for Con-
tinuous Injection and Volume Fill Inputs) or in the Liquid Drying Region section (for Custom In-
puts).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
210 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
6. When you are ready to post-process your simulation, you can make use of the new parameters
on the Properties tab for the main Particles entity.
A fraction of the liquid solution mass given by the solvent that evaporated by an unspecified mechanism.
The remaining mass will correspond to the solute dissolved into the solution. This mass will be given
by:
(20.2)
where is the value specified for the Solute Mass Fraction parameter. It is assumed that this
mass of solute will be deposited uniformly on the particle's surface when the solvent evaporates. In
order to take into account this added mass, its value will be summed to the current value of the particle's
property Mass Increment. Note that at any given time, the mass of a particle will be the sum of its
original mass (the mass the particle had at the time it entered into the simulation) and the instantaneous
value of Mass Increment.
Finally, the corresponding volume of the solute deposited on the particle s surface will be calculated
as:
(20.3)
where is the specified value of the Solute Density. This value is added to the current value of
the particle's property . Rocky will use the values of both the mass and volume increments for updating
the instantaneous values of the particle's mass, volume and size.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 211
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
212 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 21: Liquid Sphere Injection Module
21.1. Usage details
To use this module, you first need to import through a Rocky Point Cloud (Figure 21.1: Defining a Point
Cloud in Rocky. (p. 213)) a text file containing one or more points that define a region of influence
(Figure 21.2: Example text file defining a Point Cloud. (p. 214)). The text file must be generated in ex-
ternal software that produces a list containing the spatial coordinates of a set of points.
The points within the Point Cloud will then be used by the module to define spheres of influence
(Figure 21.3: Example spheres of influence formed around the points defined in a Point Cloud. (p. 214)),
inside which any particle that enters will receive a liquid injection.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 213
Liquid Sphere Injection Module
Figure 21.3: Example spheres of influence formed around the points defined in a Point Cloud.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
214 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Depending upon the particle path and the location of the spheres of influence, a single particle can
pass through several spheres and therefore receive several injections of liquid during the simulation.
• Influence Ratio: This parameter is defined as the ratio between the radius of each sphere of influence
that is centered at each point within the Point Cloud, and the particle diameter.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 215
Liquid Sphere Injection Module
• Injection Clouds:
• Point Cloud: This defines the pre-imported Point Cloud upon which this set of spheres of influence
will be created.
• Liquid Injection Flux Density: This defines the liquid injection per unit area and unit time that
you want applied to any particles that enter a sphere of influence based upon the selected Point
Cloud.
• Start Time: This defines the time you want particles to begin receiving liquid injection from the
spheres of influence for the selected Point Cloud.
• Stop Time: This defines the time you want particles to stop receiving liquid injection from the
spheres of influence for the selected Point Cloud.
For the Inputs step, a new Modules sub-tab will appear to help you define a module-related
property for the particles upon injection, as shown in Figure 21.5:Liquid Injection Sphere Clouds
module options for a Continuous Injection Input. (p. 217) and defined below.
Note:
This sub-tab will appear only for the Continuous Injection and Volume Fill inputs. For
Custom Inputs, there is not a Modules sub-tab; rather, the Module's parameters will
appear on the main tab for Custom Inputs.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
216 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 21.5: Liquid Injection Sphere Clouds module options for a Continuous Injection Input.
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab for the Input.
• Liquid Mass: This defines the liquid mass that will be applied to each particle within the selected
Particle set when it enters into the simulation. Important: Every particle in the set, regardless of
the particle size or volume, will contain the same amount of liquid mass when injected.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 217
Liquid Sphere Injection Module
Figure 21.6: Particles property created by Liquid Injection Sphere Clouds module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
218 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• Liquid Mass: This provides the mass of the liquid film around each individual whole particle
or fragment. This can help you to visualize the evolution of the particle's liquid content as the
particles move in and out of the spheres of influence.
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Liquid Injection Sphere Clouds checkbox is enabled.)
3. From the Data panel, under Modules, select Liquid Injection Sphere Clouds.
i. Click the Add button (green plus) to include a new Injection Cloud row.
ii. Within the new row, select a Point Cloud you imported in #1, and then define the Liquid
Injection Flux Density, Start Time, and Stop Time parameters for that Point Cloud set.
c. Repeat step #4b for each individual Point Cloud you want considered by the module.
5. Set up your simulation as you normally would, with the following exception:
• For the Inputs step, ensure that for each Input and Particle set you define, you also define
the Liquid Mass value you want.
b. From the Data Editors panel, select the Properties tab, and then display in a 3D View window
or plot the Liquid Mass property.
Where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 219
Liquid Sphere Injection Module
is the liquid mass, measured in , that is added to the particle during a simulation timestep
when it is inside a sphere of influence.
is the liquid injection flux density, , which is an option defined by the user in the Liquid
Injection Sphere Clouds module.
is the projected area of the particle, measured in [m2], which no matter the shape of the
particle, is computed assuming a sphere with equivalent volume.
is the distance between the centroid of a particle and the closest point in the associated Point Cloud,
measured in [m].
is the particle diameter, measured in [m], which no matter the shape of the particle, is always equal
to the particle sieve size.
is the influence ratio, which is a dimensionless parameter defined by the user in the Liquid Injection
Sphere Clouds module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
220 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 22: Liquid Injection Region Module
22.1. Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 221
Liquid Injection Region Module
• Start Time: This defines the time you want particles in the regions of interest to begin receiving
liquid film.
• Stop Time: This defines the time you want particles in the regions of interest to stop receiving
liquid film.
• Region of Interest: The Region of Interest in which you want liquid film injection application
to occur.
Range: [Any Region of Interest that has already been defined in the simulation]
• Liquid Injection Mass Rate: This defines the constant rate at which mass of liquid film is applied
to the particles in this region.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
222 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 22.2: Options for a Continuous Injection Input (on the left) and a Volume Fill Input (on
the right).
Within the table on the Modules sub-tab, one parameter (column) can be set per Particle set (row),
as defined below:
• Particle: For each row, lists the Particle set that is defined on the Particles sub-tab.
• Liquid Mass: This defines the liquid mass that will be applied to each particle within the
selected Particle set when it enters into the simulation. Important: Every particle in the set,
regardless of the particle size or volume, will contain the same amount of liquid mass when
it enters into the simulation.
For Custom Inputs, a new Liquid Injection Region section will appear to help you define the liquid
film input properties of the particles, as shown in Figure 22.3: Options for a Custom Input. (p. 225).
The only parameter in this section is Liquid Mass, which is defined exactly the same as indicated
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 223
Liquid Injection Region Module
in the preceding section. Asterisks (*) following a parameter name indicate that the parameter and
value entered are shared with other modules or UI entities that also make use of this parameter.
Shared parameters may have only one value stored for it; the value you enter most recently will
override any shared values set earlier in other locations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
224 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 225
Liquid Injection Region Module
• Liquid Mass: This provides the mass of the liquid film around each individual whole particle
or fragment. This can help you to visualize the evolution of the particle's liquid content during
a simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
226 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 22.4: Particle Properties added by the Liquid Injection Region module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
2. From the Data panel, under Modules, select Liquid Injection Region.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 227
Liquid Injection Region Module
b. Under Liquid Injection Regions, click the Add button (green plus) to include as many separate
rows for the number of regions as you want, and then for each separate row, select the Region
of Interest and define the Liquid Injection Mass Rate. You may add as many rows as you
want in the simulation, and you may also remove rows, if necessary. These operations must
be performed using the buttons on the top of the table.
4. Set up the simulation as you normally would, with the following exception:
• When setting up your Inputs, be sure you define the parameters on the Modules tab (for
Continuous Injection Inputs and Volume Fill Inputs) or in the Liquid Injection Region section
(for Custom Inputs).
6. When you are ready to post-process your simulation, you can make use of the new parameters
on the Properties tab for the main Particles entity.
where:
is the value of Liquid Injection Mass Rate defined for the region, in the table associated to Liquid
Injection Regions.
If the regions defined in the project have overlapping portions, then particles located inside those
overlaps will receive liquid simultaneously from more than one region.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
228 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 23: Thornton-Ning Normal Force Module
23.1. Usage details
Figure 23.1: Physics | Momentum tab options when the Thornton-Ning Normal Force module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 229
Thornton-Ning Normal Force Module
However, once the module is enabled, it does provide additional options in the Materials Interactions
section of the UI, as explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
230 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 23.3: Materials Interactions options when the Thornton-Ning Normal Force module is
enabled.
where:
• Yield Ratio: This defines the ratio between the yield contact overlap and equivalent diameter,
which is used to determine the transition between the elastic and plastic regimes in contact
loading/unloading phases.
Range: [0,1]
2. Set up the simulation as you normally would, with the following exceptions:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 231
Thornton-Ning Normal Force Module
a. When setting up your Physics parameters, note this model overrides the Normal Force setting on
the Momentum sub-tab, so you will be unable to set that parameter when this module is enabled.
b. When setting up your Materials | Materials Interactions, ensure that for each material-material
pair you are making use of in the simulation, you define what you want for the parameter in the
Thornton-Ning Normal Force section.
(23.1)
where:
The yield normal contact overlap is defined thourgh the Yield Ratio input parameter [65] (p. 380):
(23.2)
where:
(23.4)
(23.5)
(23.6)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
232 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
The plastic behavior is described by a force-displacement relationship defined by the yield normal
contact overlap and the critical yield normal force :
(23.7)
(23.8)
(23.9)
(23.10)
where:
During the unloading phase, the decompression force can be purely elastic (Hertz model) or elastic
with plastic stiffness, depending on the maximum overlap reached by the contact and the yield normal
contact overlap [64] (p. 380):
(23.11)
where:
For unloadings which the maximum overlap is under the yield point, the contact force is purely elastic
and follows the same equation as defined for loading. When the maximum overlap is beyond the yield
point, the contact force is defined by a modified Hertz model with a plastic stiffness :
(23.12)
(23.13)
where:
The plastic effective radius is derived based on the maximum loading force [64] (p. 380):
(23.14)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 233
Thornton-Ning Normal Force Module
where:
The residual plastic overlap is calculated based on Herz force model using the plastic effective
radius for the maximum compression force and the maximum overlap [64] (p. 380):
(23.15)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
234 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 24: Nonlinear Hysteretic Normal Force
Module
24.1. Usage details
In Rocky, the Normal Force is used to calculate the normal components of the contact forces, and will
therefore affect particle-particle and particle-boundary interactions.
Selecting this module overrides what you have set for Normal Force in the Physics | Momentum
tab, as shown in Figure 24.1:Physics | Momentum tab options when the Nonlinear Hysteretic Normal
Force module is enabled. (p. 236).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 235
Nonlinear Hysteretic Normal Force Module
Figure 24.1: Physics | Momentum tab options when the Nonlinear Hysteretic Normal Force module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
236 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
However, once the module is enabled, it does provide additional options in the Materials Interactions
setup section of the UI, as explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 237
Nonlinear Hysteretic Normal Force Module
Figure 24.3: Materials Interactions options when the Nonlinear Hysteretic Normal Force module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
238 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
2. Set up the simulation as you normally would, with the following exceptions:
a. When setting up your Physics parameters, note this model overrides the Normal Force setting
on the Momentum sub-tab, so you will be unable to set that parameter when this module is
enabled.
b. When setting up your Materials Interactions, ensure that for each material-material pair you
are making use of in the simulation, you define what you want for the three parameters in
the Nonlinear Hysteretic Normal Force section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 239
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
240 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 25: SPH Non-Newtonian Models (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However,
your feedback will help us improve the overall quality of the product. We will not
guarantee that the projects using this beta feature will run successfully when the feature
is finally released so you may, therefore, need to modify the projects.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 241
SPH Non-Newtonian Models (Beta)
Figure 25.1: Example application of the SPH Non-Newtonian (Beta) module in a food modelling
scenario.
• The current module (like any other fluid model available in the literature) is not able to simulate a
real solid behavior below the yield stress, so mathematical manipulations are required to increase
the viscosity of the fluid to model a solid-like behavior.
• The implemented models do not consider more complex non-newtonian characteristics such as
elasticity and time-dependency of the properties.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
242 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Non-Newtonian Model: List of the three available non-Newtonian viscosity models: Power Law,
Bingham Plastic and Herschel-Bulkley.
• Power Law and Herschel-Bulkley Parameters: The required input parameters for the models
Power Law and Herschel-Bulkley.
– Consistency Index: This defines the consistency index used in the consitutive equation.
– Behavior Index: This defines the behavior index used in the consitutive equation.
• Bingham Plastic Parameters: The required input parameter for the Bingham Plastic model.
– Plastic Viscosity: This defines the plastic viscosity used in the consitutive equation.
• Bingham Plastic and Herschel-Bulkley Parameters: The required input parameters for the models
Bingham Plastic and Herschel-Bulkley.
– Yield Stress: This defines the yield stress limit used in the consitutive equation.
– Viscosity Factor: This defines a viscosity factor used in the consitutive equation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 243
SPH Non-Newtonian Models (Beta)
2. Select SPH Non-Newtonian Models (Beta) from the Data panel under Modules.
a. Choose one of the three available models in the Non-Newtonian Model list.
b. Depending on the model you have chosen, set the desired parameters in the Power Law and
Herschel-Bulkley Parameters, Bingham Plastic Parameters and Bingham Plastic and Her-
schel-Bulkley Parameters sections.
These models can also be written as a generalized Newtonian fluid, which is defined as a fluid for which
the shear stress depends on the shear rate at a given time, but not depends on the deformation history
as in other non-Newtonian fluids. The constitutive law for this type of fluids is usually written as:
(25.1)
where is the shear stress tensor, is the shear rate tensor and is the so-called apparent or effective
viscosity, which is a function of the shear rate. The symbol in Equation 25.1 (p. 244) denotes the
Frobenius norm, defined for a second-order tensor as:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
244 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(25.2)
The shear rate is a symmetric tensor related to the velocity gradient by means of:
(25.3)
The specific constitutive equation for the Power Law model is given by:
(25.4)
where and are model parameters listed in the module UI as Consistency Index and Behavior Index
, respectively. The behavior index is a dimensionless parameter, while the S.I. units of the consistency
index are . However, due to the impossibility of representing these units within the module UI,
both parameters will appear as dimensionless there.
The constitutive equation for the Bingham Plastic model is usually formulated as:
(25.5)
in which is the Yield Stress value, below which the modeled material behaves like a solid. On the
other hand, , usually known as Plastic Viscosity, is the slope of the stress-shear rate line beyond the
yield point, at which the material begins to flow. Instead of using Equation 25.5 (p. 245), the following
approximate alternative representation suggested in [53] (p. 380) is used in the module:
(25.6)
where is an additional parameter listed in the module UI as Viscosity Factor. Instead of treating the
material as a solid if the stress is below the yield point, which is very difficult using a SPH formulation,
it is represented alternatively as very viscous fluid, with a viscosity value given by . The shear rate
at the yield point for this approach is given by:
(25.7)
The Herschel-Bulkley model is a generalization of the Bingham Plastic model, that introduces an exponent
as an additional parameter for representing the stress-shear rate relation beyond the yield point in a
similar way to the Power Law. In the same way as in the Bingham Plastic model, the behavior of the
material below the yield point is modeled as a very viscous fluid in the module. Thus, the expression
used for the Herschel-Bulkley model used the module has the following form:
(25.8)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 245
SPH Non-Newtonian Models (Beta)
where and are the Consistency Index and the Behavior Index , respectively, while is the Yield
Stress and is the Viscosity Factor . All those parameters have the same meaning already mentioned
when describing the preceding models. For the Herschel-Bulkley model, the value of the shear rate at
the yield point that guarantees the continuity of the stress-shear rate curve at that point is:
(25.9)
The module calculates the viscous stresses defined according to one of the preceding models, the one
selected on the Non-Newtonian Model dropdown list. Those stresses cause accelerations on the SPH
elements, that are also calculated by the module. The acceleration values are then passed to the SPH
solver, which applies them to the elements in the simulation. The acceleration of an element caused
by the viscous stresses is approximated by means of the expression:
(25.10)
where the sum involves all neighbor elements on the kernel support region of particle , is the
mass of the elements, is their density, and is the gradient of the kernel function. In this expres-
sion, and are the stress tensor values calculated considering the selected non-Newtonian model,
for elements and , respectively.
All the numerical approximations involved in the calculation of viscous accelerations are explicit,
therefore, the whole simulation will be subjected to a stability constraint related to those approximations.
This constraint can be more severe than the one found for Newtonian flow simulations, due to the fact
that the apparent viscosity can achieve very large values. More specifically, the SPH timestep must sat-
isfy the following condition for the simulation to be stable:
(25.11)
where is the fluid density value, is the initial spacing between SPH elements. With an estimation
of the upper bound for the apparent viscosity during a simulation it would be possible to determine
an stable timestep using Equation 25.11 (p. 246). During initialization, the module will use a coarse es-
timation of such upper bound and will adjust the simulation timestep if the value obtained with Equa-
tion 25.11 (p. 246) is smaller than the value estimated internally by the solver. However, since only in
the case of the Bingham Plastic model it is possible to estimate with confidence an upper bound for
the apparent viscosity, instabilities may still arise when using the other models, depending on the pre-
scribed parameter values. In such cases, users may need to reduce the simulation timestep manually.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
246 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 26: Heat Conduction Model Module
26.1. Usage details
For each material that will be used in the simulation, you will need to fill both lookup tables (Specific
Heat and Thermal Conductivity) with one or more rows.
In addition, because this module relies on the same code structure used by the Heat Conduction
Model and the Thermal Integration Model included in Rocky, enabling this module will therefore
override these selections on the Physics | Thermal tab (Figure 26.1:Physics | Thermal tab options
when the Nonlinear Thermal Contact Model module is enabled. (p. 248)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 247
Heat Conduction Model Module
Figure 26.1: Physics | Thermal tab options when the Nonlinear Thermal Contact Model module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
248 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
When you select a material component under Materials, two new tables will appear to help you
define the specific heat and thermal conductivity properties of the material at different temperatures,
as shown in Figure 26.3:Nonlinear Thermal Contact Model module options for an individual material
component. (p. 250). You must add at least one row per table. If you choose to add 2 or more rows,
the values that you enter can be in any order you wish.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 249
Heat Conduction Model Module
Figure 26.3: Nonlinear Thermal Contact Model module options for an individual material
component.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
250 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
– Temperature: This defines the temperature of reference related to the specific heat property.
Range: [Positive values]
– Specific Heat: This defines the specific heat property to be interpolated. Range: [Positive values]
– Temperature: This defines the temperature of reference related to the thermal conductivity property.
Range: [Positive values]
a. The Thermal Model checkbox is enabled on the Thermal tab of the Physics entity
b. The module is enabled. (From the Data panel, select Modules and then from the Data Editors
panel, ensure the Nonlinear Thermal Contact Model checkbox is enabled.)
2. Set up your simulation as you normally would, with the following exceptions:
• For the Materials step, for each material that you want to use in your simulation, select it, and
then from the Data Editors panel, do the following:
a. For both the Specific Heat Table and the Thermal Conductivity Table, click the Add button
(green plus) to add as many rows to the tables as you want to have points for interpolation.
b. For each separate row in the table, define the Temperature and your related property value.
Tips:
– To insert another row at the bottom of the table, click the Add button.
– To remove all rows in the table, click the Remove All button.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 251
Heat Conduction Model Module
• Ensure that you define all thermal-related properties for your geometries and particle inputs.
4. When you are ready to analyze your simulation results, you may choose to make use of the fol-
lowing thermal-related Property:
The piece-wise scheme will search for the backward and forward points based on the particle's current
temperature; both specific heat and thermal conductivity properties will be calculated following a linear
ramp between the two consecutive points:
(26.1)
where:
• is the interpolated value of the particle s specific heat or thermal conductivity corresponding
to the temperature , measured in [J/kg.K] or [W/m.K], depending on the property.
• are the property values (specific heat or thermal conductivity) at the points that
limit the interval inside which the temperature is located, measured in [J/kg.K] or [W/m.K], depending
on the property.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
252 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• are the temperatures that limit the interval inside which the temperature is located, measured
in [K].
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 253
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
254 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 27: Normal and Tangential Energy
Dissipation Module
27.1. Usage details
• The same Particle Curve calculations you get in Rocky when the Energy Dissipation option is en-
abled on the Inter-group Collision Statistics module (Figure 27.1: Inter-group Collisions Statistics
module Curves upon which the Normal and Tangential Energy Dissipation module is based. (p. 255))
• Similar Particle Property calculations you get in Rocky when the Power option is enabled on the
Inter-particle Collision Statistics module, only instead of power, it calculates energy (Fig-
ure 27.2: Inter-particle Collision Statistics module Property upon which the Normal and Tangential
Energy Dissipation module is based. (p. 256)).
Figure 27.1: Inter-group Collisions Statistics module Curves upon which the Normal and
Tangential Energy Dissipation module is based.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 255
Normal and Tangential Energy Dissipation Module
Figure 27.2: Inter-particle Collision Statistics module Property upon which the Normal and
Tangential Energy Dissipation module is based.
Neither of these two embedded Collisions Statistics modules need to be enabled in order for the
Normal and Tangential Energy Dissipation module to work. Note: An embedded Module comes by
default installed with the Rocky product. This kind of Module is documented in the User Manual and
Technical manuals.
Once the Normal and Tangential Energy Dissipation module is enabled, the data for Inter-particle and
Inter-group energy dissipation will be collected during simulation processing.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
256 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer
to Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 257
Normal and Tangential Energy Dissipation Module
Figure 27.4: Particles Curves created by the Normal and Tangential Energy Dissipation module.
• Energy: Dissipation: Normal: Provides the total energy dissipated by the normal component of
contact forces for all collisions recorded for the indicated pair (particle-particle or particle-geometry)
that occurred during an interval between two consecutive output times.
• Energy: Dissipation: Tangential: Provides the total energy dissipated by the tangential component
of contact forces for all collisions recorded for the indicated pair (particle-particle or particle-
geometry) that occurred during an interval between two consecutive output times.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
258 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 27.5: Particles properties created by the Normal and Tangential Energy Dissipation
module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 259
Normal and Tangential Energy Dissipation Module
• Energy : Dissipation : Normal: Provides the total energy dissipated by the normal component
of contact forces for all collisions recorded for each individual whole particle or fragment that oc-
curred during an interval between two consecutive output times.
• Energy : Dissipation : Tangential: Provides the total energy dissipated by the tangential com-
ponent of contact forces for all collisions recorded for each individual whole particle or fragment
that occurred during an interval between two consecutive output times.
2. Set up your simulation as you normally would, keeping in mind the following:
• When setting up your Physics parameters, remember that this module will not take into account
the effects from Adhesive Forces and Rolling Resistance Torque. (In other words, if on the
Physics | Momentum tab, you choose to define an Adhesive Force and/or a Rolling Resistance
Model, this module will ignore those factors when generating its results.)
4. When you are ready to analyze your simulation results, do the following:
• From the Data Editors panel, select the Properties tab, and then display in a 3D View
window or plot the and/or the Energy : Dissipation : Normal and/or the Energy : Dissipa-
tion : Tangential property.
• From the Data Editors panel, select the Curves tab, and then plot the Energy : Dissipation
: Normal and/or the Energy : Dissipation : Tangential Curves for the particleparticle or
particle-group pair you want to analyze.
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
260 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(27.2)
where:
• and are, respectively, the normal and tangential components of the contact force, which are
calculated forboth the Normal Force and Tangential Force contact models that you select on the
Physics | Momentum tab (Figure 27.6: Normal Force and Tangential Force Contact Models. (p. 261))
for your Rocky project.
Combining equations [eq-work] (p. 260) and [eq-contact_force] (p. 261) yields:
(27.3)
This way, the dissipation work done by the normal and tangential components of the contact force are,
respectively:
(27.4)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 261
Normal and Tangential Energy Dissipation Module
(27.5)
Equations [eq:work_normal] (p. 261) and [eq:work_tangential] (p. 262) calculate the works of the normal
and tangential components of the contact force in a collision between two particles. In order to calculate
energy dissipated by each particle, the dissipated work of the collision must be split between particles
using a splitting factor, which is inversely proportional to the particle/geometry stiffness. Thus, the works
and split between particle, A, and particle/geometry, B, are given by the following expressions:
(27.6)
(27.7)
(27.8)
(27.9)
(27.10)
where and are the stiffnesses attributed to particle A and particle/geometry B, respectively.
The rationale behind this practice is that the energy dissipated during a collision by each one of the
colliding entities will be proportional to its own deformation, which, in turn, will be inversely propor-
tional to its stiffness.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
262 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 28: 1-Way CFD Field Coupling Module
28.1. Usage details
To use this module, you first need to import a fluid velocity field through a Rocky Point Cloud (Fig-
ure 28.1: Defining a Point Cloud in Rocky. (p. 263)). The structure of this text file is explained in the
next section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 263
1-Way CFD Field Coupling Module
Figure 28.2: Example text file defining a fluid velocity field point cloud.
• X, Y, Z point coordinates.
This text file must then be imported during the process of creating a Point Cloud.
This module works only with an imported Point Cloud that meets the exact specifications defined in
the About the Fluid Velocity Field Point Cloud (p. 264) section.
Although it is possible to load multiple Point Clouds into Rocky, this module can make use of only
one Point Cloud.
The external imported fluid velocity field will be active throughout the whole simulation time.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
264 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 28.3: Options for the 1-Way CFD Field Coupling module.
• Fluid Density: This defines the fluid density that you want to apply to the imported fluid velocity
field.
• Fluid Dynamic Viscosity: This defines the fluid dynamic viscosity that you want to apply to the
imported fluid velocity field.
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the 1-Way CFD Field Coupling checkbox is enabled.)
3. From the Data panel, under Modules, select 1-Way CFD Field Coupling.
c. Under Point Cloud, select the Point Cloud you imported earlier in step #1.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 265
1-Way CFD Field Coupling Module
where:
In this module, the drag coefficient, , is calculated using the modified version of the Schiller & Naumann
correlation, which is valid for spherical particles and is given by the following equation
[24] (p. 378)[22] (p. 378):
(28.2)
where the relative particle Reynolds number, , is defined using the particle diameter and the relative
particle-fluid velocity according to [10] (p. 377):
(28.3)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
266 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 29: Overlap-Based Breakage Model Module
29.1. Usage details
These models are especially useful for simulating machines that fragment particles due to compression,
like jaw, gyratory, or roll crushers. The contact overlap is used as a measurement of the particle de-
formation.
Once the module is enabled, you must use it when setting up your Particle sets to define both the
Model parameters (defined from the Particle's Breakage | Critera tab) and the Distribution model
parameters (defined from the Particle's Breakage | Fragments tab) as shown in Figure 29.1:Breakage
| Criteria tab options (left) and Breakage | Fragments tab options (right) for the Particle set when the
Overlap-Based Breakage Model module is enabled. (p. 267).
Figure 29.1: Breakage | Criteria tab options (left) and Breakage | Fragments tab options (right)
for the Particle set when the Overlap-Based Breakage Model module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 267
Overlap-Based Breakage Model Module
• A custom instantaneous Breakage Model that analyzes particle's deformation in order to determine
if a particle must break.
In addition, because this module relies on the same code structure used by the Breakage Model and
by the Fragments Distribution Model included by default in Rocky, after enabling this module you will
therefore need to select it when defining your particle breakage criteria.
However, there are additional module-related parameters that you must set for each Particle set. Refer
to the next section for details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
268 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 269
Overlap-Based Breakage Model Module
Figure 29.3: Overlap-Based Breakage Model module options for a Particle set.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
270 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• [Breakage] Kf0 - Kf4: These five settings are coefficients of a 4th order polynomial that describes
the breakage probability as a function of the particle's relative deformation, both expressed as a
percentage.
Range
• [d0] Ks0 - Ks1: These two settings define the coefficients of a linear polynomial that gives the
maximum relative size of the fragments as a function of the particle's relative deformation at the
moment of breakage, for relative deformations smaller than 5%.
• [d50] Km0 - Km1: These two settings define the coefficients of a linear polynomial that specifies
the relative median size of the fragment distribution as a function of the particl's relative deformation
at the moment of breakage.
• [sigma] Ksl0 - Ksl1: These two settings define the coefficients of a linear polynomial that defines
the standard deviation of the fragment size distribution as a function of the particle's relative de-
formation at the breakage moment, when this relative deformation is larger than 5%.
• [sigma] Kss0 - Kss2: These three settings define the coefficients of a quadratic polynomial that
defines the standard deviation of the fragment size distribution as a function of the particle's relative
deformation at the breakage moment, when this relative deformation is smaller than 5%.
Besides that, new Properties are available for the main Particles entity, as shown in Figure 29.4: Particles
Properties created by the Overlap-Based Breakage Model module. (p. 272) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 271
Overlap-Based Breakage Model Module
Figure 29.4: Particles Properties created by the Overlap-Based Breakage Model module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
272 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• Maximum Overlap With Particles: For a given particle, this provides the maximum overlap value
between all particle-particle contacts in which the particle is involved at the output time. It is used
for model calculation purposes.
• Maximum Overlap With <Geometry Component>: For a given particle, this provides the maximum
overlap value between the particle and the indicated geometry, at the current output time. It is
used for model calculation purposes.
There will be a separate Property entry for each particle-geometry pair in the simulation.
2. Set up the simulation as you normally would, with the following exception:
• On the Particle | Breakage | Criteria tab for at least one Particle set, you have defined the
Model as Overlap-Based Breakage Model.
• On the Particle | Breakage | Fragments tab for at least one Particle set, you have defined
the Distribution model as Overlap-Based Breakage Model.
• On the Particle | Modules tab, you have entered the values you want.
4. When you are ready to analyze your simulation results, do one or more of the following:
• Make use of the default breakage Curves and/or Properties for Particles.
• Make use of the additional module-related Properties for Particles listed below:
(29.1)
(29.2)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 273
Overlap-Based Breakage Model Module
• are coefficients for the breakage probability polynomial. These dimensionless parameters
are defined by the user for Particle sets when the Overlap-Based Breakage Model module is enabled.
• is the equivalent diameter of the particle, which is defined as the diameter of a sphere whose
volume is equal to the particle's volume.
• and are the maximum normal deformations of a particle with regards to the particle-
geometry and particle-particle contacts, respectively. These maximum normal deformations are the
highest values of normal deformation among all particle-geometry and particle-particle contacts that
a particle is part of at a given time.
(29.3)
where:
• and are the normal deformations of a particle with regards to the particle-geometry and
particle-particle contacts, respectively.
• and are the stiffnesses of the particle and its contacting particle, respectively, on a particle-
particle contact.
A particle will break if at any moment the value of computed with equation [eq:eq1] (p. 273) is
larger than the strength of the particle. When that condition is met and the particle breaks, the fragments
are generated following a truncated logarithmic normal distribution of the particle diameter, :
(29.4)
(29.5)
(29.6)
(29.7)
(29.8)
(29.9)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
274 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• are coefficients for the relative maximum fragment size polynomial. These dimensionless
parameters are defined by the user for Particle sets when the Overlap-Based Breakage Model module
is enabled.
• are coefficients for the relative median fragment size polynomial. These dimensionless
parameters are defined by the user for Particle sets when the Overlap-Based Breakage Model module
is enabled.
• are coefficients for the upper part of the standard deviation polynomial. These dimensionless
parameters are defined by the user for Particle sets when the Overlap-Based Breakage Model module
is enabled.
• are coefficients for the lower part of the standard deviation polynomial. These dimensionless
parameters are defined by the user for Particle sets when the Overlap-Based Breakage Model module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 275
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
276 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 30: Parametric Gidaspow, Bezburuah & Ding
Module
30.1. Usage details
These three parameters can be adjusted in order to allow a fine tuning of the drag coefficient. It can
help you to correlate the simulation results with experimental data you have available.
While using this module you can also enable the CFD Coupling Particle Statistics module to evaluate
the drag force for each particle present in your simulation.
After enabling this module, the parametric Gidaspow, Bezburuah & Ding Drag Law will be available
for selection among the other drag laws embedded in Rocky-Fluent coupling. Therefore, for using
this model to calculate the drag law for a particular group of particles, you need to select it when
defining your CFD Coupling Particle Interactions for this group (Figure 30.1:Interactions tab options
for the 2-Way Fluent Unresolved Coupling mode when the Parametric Gidaspow, Bezburuah & Ding
module is enabled. (p. 278)).
It is important to emphasize that by enabling this module, you must define the parametric Gidaspow,
Bezburuah & Ding Drag Law for at least one particle group in your project.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 277
Parametric Gidaspow, Bezburuah & Ding Module
Figure 30.1: Interactions tab options for the 2-Way Fluent Unresolved Coupling mode when the
Parametric Gidaspow, Bezburuah & Ding module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
278 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• parameter 1: This parameter defines the fluid volume fraction power coefficient of the parametric
Gidaspow, Bezburuah & Ding Drag Law for solid volume fractions less than 20%.
• parameter 2: This parameter defines the viscous energy losses coefficient of the parametric Gi-
daspow, Bezburuah & Ding Drag Law for solid volume fractions greater than 20%.
• parameter 3: This parameter defines the kinetic energy losses coefficient of the parametric Gidaspow,
Bezburuah & Ding Drag Law for solid volume fractions greater than 20%.
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Parametric Gidaspow, Bezburuah & Ding checkbox is enabled.)
3. From the Data panel, under Modules, select the new Parametric Gidaspow, Bezburuah & Ding
entry.
4. From the Data Editors panel, on the Parametric Gidaspow, Bezburuah & Ding tab, enter the
values you want for the three parameters.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 279
Parametric Gidaspow, Bezburuah & Ding Module
• On the Interactions tab for at least one Particle, you have defined the Drag Law as Parametric
Gidaspow, Bezburuah & Ding.
where is the relative velocity between particle and fluid and is the projected particle area in
the flow direction.
Various drag correlations based on particle shape (spherical and non-spherical) and particle concentration
(dilute or dense flows) are available within the Rocky package for the calculation of the drag coefficient.
For a dense flow of particles, Gidaspow, Bezburuah & Ding[33] (p. 378) have developed a simple connec-
tion between other two drag correlations (Wen & Yu [34] (p. 378) and Ergun[31] (p. 378)) to represent
the complete range of fluid volume fraction in a single drag law, by simply applying each law
over its valid range. The Gidaspow, Bezburuah & Ding correlation is then given by:
(30.2)
where:
The particle's sphericity is defined as the ratio between the actual particle surface area and the surface
area of a sphere having the same volume as the particle.
The parameterization provided by this module consists of the introduction of three parameters that
can be modified by the user to achieve a better correlation with experimental data. The parameterized
form of the Gidaspow, Bezburuah & Ding correlation is then given by:
(30.3)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
280 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
where:
• P1 is the parameter to define the fluid volume fraction power coefficient for solid volume fractions
less than 20%.
• P2 is the parameter to define the viscous energy losses coefficient for solid volume fractions greater
than 20%.
• P3 is the parameter to define the kinetic energy losses coefficient for solid volume fractions greater
than 20%.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 281
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
282 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 31: Parametric Ranz & Marshall Law Module
31.1. Usage details
• The project is CFD-coupled with one of the following modes: 1-Way Fluent Steady State, 2-Way
Fluent, or 1-Way Constant.
• Note: While you are able to use this module with any particle type, the correlation enabled by this
module is valid mainly for spherical particles.
This module does not work with the 2-Way Fluent Semi-Resolved nor the 1-Way LBM CFD Coupling
methods.
In addition, because this module relies on the same code structure used by the Convective Heat
Transfer Law included by default in Rocky, after enabling this module you will therefore need to select
it when defining your CFD Coupling Particle Interactions (Figure 31.1:Interactions tab options for the
1-Way Constant CFD Coupling mode when the Parametric Ranz & Marshall Law module is en-
abled. (p. 284)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 283
Parametric Ranz & Marshall Law Module
Figure 31.1: Interactions tab options for the 1-Way Constant CFD Coupling mode when the
Parametric Ranz & Marshall Law module is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
284 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• a - b: These two settings define the coefficients of the parametric Ranz & Marshall Convective Heat
Transfer Law.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 285
Parametric Ranz & Marshall Law Module
• If you plan to use this module with a 2-Way Fluent simulation, ensure the Fluent case has the
Energy equation activated.
• If you plan to use this module with a 1-Way Fluent Steady State simulation, ensure the Fluent
case has a temperature field available.
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the Parametric Ranz & Marshall Law checkbox is enabled.)
3. From the Data panel, under Modules, select the new Parametric Ranz & Marshall Law entry.
4. From the Data Editors panel, on the Parametric Ranz & Marshall Law tab, enter the values you
want.
5. Continue setting up the simulation as you normally would, with the following exceptions:
a. When setting up your Physics parameters, ensure that the Thermal Model is enabled
• You have selected a supported CFD Coupling mode. (Either 1-Way Constant, 1-Way Fluent
Steady State, or 2-Way Fluent).
• For the CFD Coupling mode you have enabled, on the Interactions tab for at least one
Particle, you have defined the Convective Heat Transfer Law as Parametric Ranz & Mar-
shall Law.
7. When you are ready to analyze your simulation results, you may choose to make use of the fol-
lowing thermal-related Property:
This model is made parametric by adding two free parameters to the calculation of the Nusselt number,
, as described by the following equations:
(31.1)
(31.2)
(31.3)
where:
• a and b are the coefficients of the Parametric Ranz & Marshall Convective Heat Transfer Law. These
dimensionless parameters are defined by the user in the Parametric Ranz & Marshall Law module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
286 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• is the relative particle Reynolds number, which is defined using the equivalent diameter of the
particle and the relative particle-fluid velocity.
• is the fluid density, measured in [kg/m3], which is defined by the user on the CFD Coupling setup.
• is the fluid velocity vector, measured in [m/s], which is defined by the user on the CFD Coupling
setup when CFD Coupling mode is 1-Way Constant.
• is the equivalent diameter of the particle, measured in [m], which is defined as the diameter of a
sphere whose volume is equal to the particle's volume.
• is the fluid dynamic viscosity, measured in [kg/m.s], which is defined by the user on the CFD
Coupling setup.
• is the fluid specific heat, measured in [J/kg.K], which is defined by the user on the CFD Coupling
setup.
• is the fluid thermal conductivity, measured in [W/m.K], which is defined by the user on the CFD
Coupling setup.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 287
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
288 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 32: Particle Freeze Outside Region Module
32.1. Usage details
Figure 32.1: Disk harrow simulation using the Particle Freeze Outside Region module. Blue
particles are enabled whereas light pink particles are frozen and removed from contact
detection calculations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 289
Particle Freeze Outside Region Module
The main assumption in this model is that the particles which are far from the disk harrow are steady,
and therefore freezing them will not affect the simulation results. In this case, the region of interest
is delimited by the blue cube, and translates according to the motion of the disk harrow.
The performance gain is reached by freezing and disabling the particles that are outside the regions
of interest from contact detection calculations.
The position of the Cubes and Cylinders used to delimit the region of interest is updated at every
timestep according to the translation of the geometries to which are linked to.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
290 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Start Time: Defines the simulation time at which the module can begin freezing and disabling
particles from contact detection calculations.
In the Active Regions, it is possible to select each Cubes and Cylinders that were created for the
Regions of interesting, in which Particle elements inside it will be frozen and disabled for contact
detection calculations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 291
Particle Freeze Outside Region Module
Figure 32.3: Additional Particle properties created by the Particle Freeze Outside Region module.
• Freeze Criterion: Criterion that classifies whether the particle is frozen and disabled from contact
calculations or not. If equals to 0.0, the particle is free to move and included in contact calculations
performed by the solver. If equals to 1.0, the particle is frozen and excluded from such calculations.
1. Set up all the desired Cubes and Cylinders whose motion will be used to translate the Regions
of Interest.
2. Set up all the required Motion Frames that will be assigned to the Cubes and Cylinders.
4. Ensure that the module is enabled (from the Data panel, select Modules and then from the Data
Editors panel, ensure the Particle Freeze Outside Region checkbox is enabled).
5. Define the Start Time and the Regions of Interest (from the Data panel, select Particle Freeze
Outside Region, and then from the Data Editors panel, input the required values).
7. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related property and may also make use of the additional module-related property.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
292 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
When using this module, if the particle is outside all the regions of interest, its Freeze Criterion, , receives
a value of 1.0, and the particle is freezed and eliminated from the contacts computation. The particle
is automatically unfrozen and re-enabled once it enters one of the regions of interested selected.
(32.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 293
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
294 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 33: Particle Freeze when Settled Module
33.1. Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 295
Particle Freeze when Settled Module
Figure 33.1: Filling of a silo using the Particle Freeze when Settled module. Grey particles are
frozen whereas blue particles are active.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
296 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
The Particle Freeze when Settled module works by first determining which particles in the simulation
have settled by checking the particles and their neighbors velocities. Then, all particles defined as
settled are frozen (grey particles), which prevents them from moving, and disabled, which removes
them from the contact evaluation step performed by the DEM solver. This way, although the number
of particles increases continuously as more particles are injected, the number of particles that are
actually involved in the collision detection and calculation steps is relative small, keeping the simulation
pace to a reasonable value.
In this velocity checking step, all particles/triangles that are in contact with the particle i are considered
neighbors. For simulations that do not have adhesive models, these are the frictional contacts. For
simulations that involve the addition of an adhesion model, the neighbors are both the frictional and
adhesive contacts. Therefore, the number of neighbors is expected to be higher is an adhesive model
is used.
Once a particle is frozen, it cannot move anymore during the remaining simulation. Therefore the
Particle Freeze when Settled module should be used only for simulations in which it is expected that
particles remain steady after settling down.
The contact detection is disabled for a pair of frozen particles, but will continue to be performed for
a frozen-unfrozen pair of particles. This avoids issues for moving particles that reach frozen particles
in this situation, the contact is detected, the collision is solved, and the moving particle trajectory
will account for the contact forces due to this collision. However, as the moving step is disabled for
a frozen particle, the previously frozen particle will not move after the collision.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 297
Particle Freeze when Settled Module
• Freezing Velocity: Defines the critical velocity that is used to determine which particles should be
frozen.
• Start Time: Defines the simulation time at which the module can begin freezing/disabling particles.
Note: It is important to set a value above zero when Volume Fill is used as the input method.
Otherwise, particles may be frozen as soon as they enter the simulation since Volume Fill particles
have an initial momentum equal to zero.
In addition, three additional particle properties are available as shown in Figure 33.3: Additional Particle
Properties created by the Particle Freeze when Settled module. (p. 299) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
298 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 33.3: Additional Particle Properties created by the Particle Freeze when Settled module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• Freeze Criterion: Criterion that classifies whether the particle is frozen. If it is equal to 0.0, the
particle is not frozen (can move as usual) and is included in contact calculations performed by the
DEM solver. If is equal to 1.0, the particle is frozen and cannot move, and is excluded from contact
calculations.
• Number of Neighbors: Provides the total number of neighboring particles and triangles surrounding
this particle.
• Neighbors Average Relative Velocity: Provides the average relative velocity between the particle
and its neighbors.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 299
Particle Freeze when Settled Module
2. Define the module parameters. (From the Data panel, select Particle Freeze when Settled, and
then from the Data Editors panel, input the required values.)
4. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related Property and may also make use of the additional module-related properties.
The crucial step of the Particle Freeze when Settled module is to decide whether one particle can be
considered settled or not. For this purpose, the Particle Freeze when Settled module checks not only the
particle translational velocity but also its neighbors velocities. In this velocity checking, all particles/tri-
angles that are in contact with the particle i are considered neighbors. For simulations that do not have
adhesive models, these are the frictional contacts. For simulations that involve the addition of an adhesion
model, the neighbors are both the frictional and adhesive contacts.
For every particle i, the magnitude of the average relative velocity between this particle and its neighbors,
, is computed as follows:
(33.1)
where:
When using this module, the particle is assumed as settled if its velocity magnitude Vi and the average
relative velocity of its neighbors are both smaller than the Freezing Velocity magnitude, Vc. Once a
particle is settled, it is frozen and its Freeze criterion, (Fi) receives a value of 1.0:
(33.2)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
300 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
If a particle is frozen , the particle is also eliminated from the computation of contacts.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 301
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
302 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 34: SPH Point Cloud Air Drag (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However,
your feedback will help us improve the overall quality of the product. We will not
guarantee that the projects using this beta feature will run successfully when the feature
is finally released so you may, therefore, need to modify the projects.
The SPH Point Cloud Air Drag (Beta) module estimates and applies the drag force caused by surround-
ing air flowing at velocity predefined in a point cloud onto a liquid flow simulated with SPH.
To use this module, you first need to import a fluid velocity field through a Rocky Point Cloud (Fig-
ure 34.1: Defining a Point Cloud in Rocky. (p. 304)). The structure of this text file is explained in the
next section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 303
SPH Point Cloud Air Drag (Beta)
Figure 34.2: Example text file defining a fluid velocity field point cloud.
• X, Y, Z point coordinates.
• x-velocity, y-velocity, z-velocity values defining the fluid velocity field vector, .
This text file must then be imported during the process of creating a Point Cloud.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
304 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
By employing this one-way coupled model, it becomes possible to evaluate the impact of the winding
effects and directions on the fluid distribution inside the gearbox.
Figure 34.3: Example application of the SPH Point Cloud Air Drag (Beta) module in a gearbox.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 305
SPH Point Cloud Air Drag (Beta)
• Inability to reproduce entrapped air effects: the module cannot accurately reproduce effects associated
with entrapped air, such as rising air bubbles. These phenomena are not accounted for in the
simulation due to the absence of a fully simulated air phase.
• One-Way coupling and limited fluid-to-air interactions: since the module operates with a one-way
coupling approach, interactions where the fluid influences the air cannot be accurately modeled.
The module focuses primarily on the effect of air on the fluid, without considering the reciprocal
influence.
• Lack of turbulent behavior and lift effects: as the module does not simulate the air phase, it cannot
capture turbulent behavior or lift effects. These features, which are essential for simulating scenarios
like a flying airplane, are not accounted for in the current implementation.
Figure 34.4: SPH Point Cloud Air Drag (Beta) module options.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
306 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
• Activate Output: This defines the outputs of the module. Depending on the selection of the items
below, additional properties are created for each SPH element and can be used for post-processing.
– Drag Force: When enabled, the drag force exerted by the air will be added as an additional
property of the SPH elements. When cleared, the drag force due to the air phase will not be
available for post-processing.
• Air Velocity Point Cloud: This field is for the selection of the Point Cloud for air velocities and the
time range it will be or not activated.
– Point Cloud: Defines the Point Cloud that will be used for the air velocity field.
– Start Time: Defines the time that drag forces calculations due to this module start/take into ac-
count the air velocity defined in the Point Cloud.
– Stop Time: Defines the time that drag forces calculations due to this module stop.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 307
SPH Point Cloud Air Drag (Beta)
Figure 34.5: Additional property enabled on the SPH entity when the SPH Point Cloud Air Drag
(Beta) is used.
You can choose to analyze this property in a plot or histogram window, or display this property
graphically in a 3D View window.
• Drag Force: This provides the drag force acting on each SPH element due to the interaction with
the surrounding air phase, which can help you to visualize how the liquid phase is affected by the
air phase.
2. Ensure that the module is enabled. Go to the Data panel, select Modules, and make sure the SPH
Point Cloud Air Drag (Beta) checkbox is enabled in the Data Editors panel.
3. Select SPH Point Cloud Air Drag (Beta) from the Data panel under Modules.
a. Specify the density and viscosity of the air in the Air Properties section.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
308 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
b. Choose whether to activate the module-related output properties in the Activate Output
section.
c. Select the Point Cloud you just imported and define the time range for its activation.
6. When you are ready to analyze your simulation results, you may choose to make use of the Drag
Force property.
The drag force acting on each SPH element is computed based on the methodology developed by
[38] (p. 379). This model not only estimates the exposed surface area of each SPH element but also ap-
proximates its deformation. By considering deformation, the model enhances the accuracy of drag force
estimation by taking into account its influence on the drag coefficient and cross-sectional area.
Figure 34.6: Flowchart illustrating the steps involved in computing the drag force for each SPH ele-
ment. (p. 310) provides an overview of the steps involved in calculating the drag force for each SPH
element. At each time step, the deformation of each SPH element is computed as described in section
Element deformation (p. 310). Based on this deformation, the drag coefficient is obtained as described
in section Drag coefficient (p. 311). Subsequently, the cross-area of the deformed SPH element and the
occlusion factor are used to compute the exposed area of each SPH element to the flow, as outlined
in section SPH exposed area (p. 312).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 309
SPH Point Cloud Air Drag (Beta)
Figure 34.6: Flowchart illustrating the steps involved in computing the drag force for each SPH
element.
This model assumes that the external force, the spring force and the damping force are dependent
on the phase of the droplet and the surrounding air phase. An additional constant is used to
convert the displacement to a dimensionless deformation , which ranges from 0 when the droplet
is a sphere to 1 when the droplet is a disk. The deformation evolves over time according to a
second-order differential equation (ODE):
(34.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
310 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Here CF, Ck, Cd, and Cb are coefficients, is the relative velocity between the air and SPH element,
and are the rest densities of the air and the liquid, is the surface tension coefficient, L is the
droplet radius and is the dynamic viscosity of the liquid.
The simplified approach used by the SPH Point Cloud Air Drag (Beta) module assumes that the SPH
element is instantaneously and completely deformed based on the current relative velocity. Con-
sequently, Equation 34.1 (p. 310) simplifies to:
(34.2)
By combining the independent variables, the approximate deformation of the SPH element is computed
as:
(34.3)
The radius of the droplet, L, is determined as the radius of a sphere corresponding to the volume h3
of an SPH element:
(34.4)
Since the SPH elements cannot break into smaller elements, deformations greater than 1 are not al-
lowed:
(34.5)
(34.7)
(34.8)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 311
SPH Point Cloud Air Drag (Beta)
To ensure that a fluid element that is part of a larger cluster of fluid elements has a drag coefficient
of 1, the final drag coefficient is interpolated between Equation 34.6 (p. 311) and 1 based on the
number of SPH neighbors n:
(34.9)
Here, represents the number of neighbors when the fluid element is fully surrounded, which
depends on the kernel type.
(34.10)
Figure 34.7: Deformation of a droplet from the spherical shape with radius L to the
disk-like shape with radius
To compute the correct cross-sectional area of the SPH elements, the SPH Point Cloud Air Drag
(Beta) module utilizes the approximate deformation to estimate the increase in radius of
the deformed droplet, denoted as , based on the work of [40] (p. 379):
(34.11)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
312 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Similar to the drag coefficient calculation, the SPH area is linearly interpolated between the square
area of the SPH element with neighbors and the area when there are no neighbors.
This is expressed by the equation:
(34.13)
An occlusion factor , which represents the fraction of the SPH element that is occluded, ranging
between 0 and 1, is used to calculate the occluded area as follows:
(34.14)
The SPH Point Cloud Air Drag (Beta) module employs a surface detection algorithm to compute
the occlusion factor . This algorithm checks the cosine of the angle between the direction of
fluid flow, , and the neighbor direction , as illustrated inFigure 34.7: Deformation of a droplet
from the spherical shape with radius L to the disk-like shape with radius (p. 312).
Figure 34.8: Occlusion factor computation accounting for the angle between the direction of
fluid flow, vi,?rel, and the neighbor direction xi,?j.
The occlusion value is computed according to Equation 34.15 (p. 313), and is clamped to obtain
a value of 1 if the SPH element is not occluded at all, and 0 if it is completely occluded.
(34.15)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 313
SPH Point Cloud Air Drag (Beta)
Here, is a vector pointing in the direction of the relative velocity difference between the air ve-
locity and the SPH element velocity .
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
314 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 35: Region Outlet Module
35.1. Usage details
Selecting this module will not override any models or settings in the Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 315
Region Outlet Module
• Start Time: This defines the time you want the removal region to begin acting as an outlet for
particles.
• Stop Time: This defines the time you want the removal region to stop acting as an outlet for
particles.
• Removal Regions:
– Regions Of Interest: This allows the selection of multiple Regions of Interest previously created
that will act as an outlet for particles.
Range: [Any Region of Interest that has already been defined in the simulation]
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
316 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
b. Under Removal Regions, click the Add button (green plus) to add as many regions as you
want separate Regions of Interest.
4. Set up and process the simulation as you normally would. During the simulation, any particles
entering a removal region will be immediately removed from the simulation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 317
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
318 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 36: SPH Air Drag Module (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However, your
feedback will help us improve the overall quality of the product. We will not guarantee that
the projects using this beta feature will run successfully when the feature is finally released
so you may, therefore, need to modify the projects.
The SPH Air Drag (Beta) module estimates and applies the drag force caused by surrounding air
flowing at constant velocity onto a liquid flow simulated with SPH.
By employing this one-way coupled model, it becomes possible to evaluate the impact of varying
wind velocities and directions on the cooling pattern of the surrounding container.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 319
SPH Air Drag Module (Beta)
Figure 36.1: Example application of the SPH Air Drag (Beta) module in a fire hose jet scenario.
• Uniform air velocity throughout the domain: the air velocity is constant and uniform throughout the
domain and its effects may be sensed by all SPH elements in the entire domain.
• Inability to reproduce entrapped air effects: the module cannot accurately reproduce effects associated
with entrapped air, such as rising air bubbles. These phenomena are not accounted for in the
simulation due to the absence of a fully simulated air phase.
• One-Way coupling and limited fluid-to-air interactions: since the module operates with a one-way
coupling approach, interactions where the fluid influences the air cannot be accurately modeled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
320 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
The module focuses primarily on the effect of air on the fluid, without considering the reciprocal
influence.
• Lack of turbulent behavior and lift effects: as the module does not simulate the air phase, it cannot
capture turbulent behavior or lift effects. These features, which are essential for simulating scenarios
like a flying airplane, are not accounted for in the current implementation.
• Air Velocity: The constant velocity vector of the air field, the magnitude of which is made up of
the x, y, and z components respectively.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 321
SPH Air Drag Module (Beta)
• Activate Output: This defines the outputs of the module. Depending on the selection of the items
below, additional properties are created for each SPH element and can be used for post-processing.
– Drag Force: When enabled, the drag force exerted by the air will be added as an additional
property of the SPH elements. When cleared, the drag force due to the air phase will not be
available for post-processing.
Figure 36.3: Additional property enabled on the SPH entity when the SPH Air Drag (Beta) is
used.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
322 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
You can choose to analyze this property in a plot or histogram window., or display this property
graphically in a 3D View window.
• Drag Force: This provides the drag force acting on each SPH element due to the interaction with
the surrounding air phase, which can help you to visualize how the liquid phase is affected by the
air phase.
2. Select SPH Air Drag (Beta) from the Data panel under Modules.
a. Specify the density and viscosity of the air in the Air Properties section.
b. Set the velocity vector of the air in the Air Velocity section.
c. Choose whether to activate the module-related output properties in the Activate Output
section.
5. When you are ready to analyze your simulation results, you may choose to make use of the Drag
Force property.
The drag force acting on each SPH element is computed based on the methodology developed by
[38] (p. 379). This model not only estimates the exposed surface area of each SPH element but also ap-
proximates its deformation. By considering deformation, the model enhances the accuracy of drag force
estimation by taking into account its influence on the drag coefficient and cross-sectional area.
Figure 36.4: Flowchart illustrating the steps involved in computing the drag force for each SPH ele-
ment. (p. 324) provides an overview of the steps involved in calculating the drag force for each SPH
element. At each time step, the deformation of each SPH element is computed as described in section
Element deformation (p. 324). Based on this deformation, the drag coefficient is obtained as described
in section Drag coefficient (p. 325). Subsequently, the cross-area of the deformed SPH element and the
occlusion factor are used to compute the exposed area of each SPH element to the flow, as outlined
in section SPH exposed area (p. 326).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 323
SPH Air Drag Module (Beta)
Figure 36.4: Flowchart illustrating the steps involved in computing the drag force for each SPH
element.
This model assumes that the external force, the spring force and the damping force are dependent
on the phase of the droplet and the surrounding air phase. An additional constant is used to
convert the displacement to a dimensionless deformation , which ranges from 0 when the droplet
is a sphere to 1 when the droplet is a disk. The deformation evolves over time according to a
second-order differential equation (ODE):
(36.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
324 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Here CF, Ck, Cd, and Cb are coefficients, is the relative velocity of the SPH element, and are
the rest densities of the air and the liquid, is the surface tension coefficient and is the dynamic
viscosity of the liquid.
The simplified approach used by the SPH Air Drag (Beta) module assumes that the SPH element is
instantaneously and completely deformed based on the current relative velocity. Consequently,
Equation Equation 36.1 (p. 324) simplifies to:
(36.2)
By combining the independent variables, the approximate deformation of the SPH element is computed
as:
(36.3)
The radius of the droplet, L, is determined as the radius of a sphere corresponding to the volume h3
of an SPH element:
(36.4)
Since the SPH elements cannot break into smaller elements, deformations greater than 1 are not al-
lowed:
(36.5)
(36.7)
(36.8)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 325
SPH Air Drag Module (Beta)
To ensure that a fluid element that is part of a larger cluster of fluid elements has a drag coefficient
of 1, the final drag coefficient is interpolated between Equation [eq:drag_1] (p. 325) and 1 based on
the number of SPH neighbors n:
(36.9)
Here, represents the number of neighbors when the fluid element is fully surrounded, which
depends on the kernel type.
(36.10)
Figure 36.5: Deformation of a droplet from the spherical shape with radius L to the
disk-like shape with radius
To compute the correct cross-sectional area of the SPH elements, the SPH Air Drag (Beta) module
utilizes the approximate deformation to estimate the increase in radius of the deformed
droplet, denoted as , based on the work of [40] (p. 379):
(36.11)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
326 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Similar to the drag coefficient calculation, the SPH area is linearly interpolated between the square
area of the SPH element with neighbors and the area when there are no neighbors.
This is expressed by the equation:
(36.13)
An occlusion factor , which represents the fraction of the SPH element that is occluded, ranging
between 0 and 1, is used to calculate the occluded area as follows:
(36.14)
The SPH Air Drag (Beta) module employs a surface detection algorithm to compute the occlusion
factor . This algorithm checks the cosine of the angle between the direction of fluid flow, ,
and the neighbor direction , as illustrated inFigure 36.6: Occlusion factor computation accounting
for the angle between the direction of fluid flow, vi,?rel, and the neighbor direction xi,?j. (p. 327).
Figure 36.6: Occlusion factor computation accounting for the angle between the direction of
fluid flow, vi,?rel, and the neighbor direction xi,?j.
The occlusion value is computed according to Equation [eq:occlusion2] (p. 327), and is clamped
to obtain a value of 1 if the SPH element is not occluded at all, and 0 if it is completely occluded.
(36.15)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 327
SPH Air Drag Module (Beta)
Here, is a vector pointing in the direction of the relative velocity difference between the air ve-
locity and the SPH element velocity .
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
328 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 37: SPH Freeze Outside Region (Beta)
Important:
Note that beta features have not been fully tested and validated. Ansys, Inc. makes no
commitment to resolve defects reported against these prototype features. However,
your feedback will help us improve the overall quality of the product. We will not
guarantee that the projects using this beta feature will run successfully when the feature
is finally released so you may, therefore, need to modify the projects.
Figure 37.1: Simplified wheel simulation using the SPH Freeze Outside Region module. Blue SPH
elements are enabled whereas light red elements are frozen and removed from contact detection
calculations.
In this module, a SPH element is enabled once it enters one of the Cubes and Cylinders that define the
regions of interest and therefore can move while collisions involving this SPH elements are solved.
Otherwise, if the SPH element is outside the regions of interest, it is frozen and removed from the
contacts computation. Also, the motion equations for the frozen elements are not resolved.
The position of the Cubes and Cylinders used to delimit the region of interest is updated at every
timestep according to the translation of the geometries to which are linked to.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 329
SPH Freeze Outside Region (Beta)
Each of the Cubes and Cylinders used to delimit the regions of interest must be linked to a geometry's
motion frame and its position is updated according to the geometry translation. If the Cube and the
Cylinder should not move, it does not be linked to any motion frame.
After enabling the SPH Freeze Outside Region module, there are some parameters you need to define.
These parameters are explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
330 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 331
SPH Freeze Outside Region (Beta)
• Start Time: Defines the simulation time at which the module can begin freezing and disabling SPH
elements from contact detection calculations.
In the Active Regions, it is possible to select each Cubes and Cylinders that were created for the Regions
of interesting, in which SPH elements inside it will be frozen and disabled for contact detection calcula-
tions.
After processing your simulation, you can analyze the results using the same SPH-related properties
that are available by default in Rocky.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
332 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 37.5: Additional SPH properties created by the SPH Freeze Outside Region module.
• Freeze Criterion: Criterion that classifies whether the SPH element is frozen and disabled from contact
calculations or not. If equals to 0.0, the element is free to move and included in contact calculations
performed by the SPH solver. If equals to 1.0, the element is frozen and excluded from such calcula-
tions.
1. Set up all the desired Cubes and Cylinders whose motion will be used to translate the Regions of
Interest.
2. Set up all the required Motion Frames that will be assigned to the Cubes and Cylinders.
4. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure the SPH Freeze Outside Region checkbox is enabled.)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 333
SPH Freeze Outside Region (Beta)
5. Define the Start Time and the Regions of Interest. (From the Data panel, select SPH Freeze Outside
Region, and then from the Data Editors panel, input the required values.)
7. When you are ready to analyze your simulation results, you may choose to make use of any SPH-
related property and may also make use of the additional module-related property.
When using this module, if the SPH element is outside all the regions of interest, its Freeze Criterion, ,
receives a value of 1.0, and the SPH element is freezed and eliminated from the contacts computation.
The element is automatically unfrozen and re-enabled once it enters one of the regions of interested
selected.
(37.1)
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
334 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 38: Stress-Based Multi-Contact Normal Force
Module
38.1. Usage details
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 335
Stress-Based Multi-Contact Normal Force Module
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
336 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 38.2: Stress-Based Multi-Contact Normal Force module options for a particle group.
• Empirical Geometric Prefactor: defines the empirical geometric prefactor (dimensionless input
parameter that accounts for the particle geometry).
2. Set up your simulation as you normally would, with the following exceptions:
• From the Data Editors panel, select each particle group that you want to use in your
simulation and then do the following:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 337
Stress-Based Multi-Contact Normal Force Module
The particle stress tensor computes the stress level on a particle exerted by its neighboring particles.
It can be calculated by averaging the sum of all contact stresses within the particle volume, making use
of a dyadic product, and it is given by:
(38.1)
where:
• Ic is the branch vector, that connects the particle centroid to the contact point.
The trace of stresses is considered by adding a new pressure term into the normal force calculation
equation. This pressure term ( ) refers to the isotropic component of the stress and it is given by:
(38.2)
illustrates the additional term, dependent on the average stress between two particles in contact, which
is calculated by equation [e:eq2] (p. 338) and added to the conventional normal force model.
Figure 38.3: Illustration of a contact pair interaction with the average stress terms of particle i
(green) and particle j (red) from the neighboring particles contacts.
Finally, the normal force is calculated using the Hertzian spring-dashpot formulation. Regarding the
trace of stress tensor term, the equation for the normal force calculation is then given by:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
338 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
(38.3)
where:
• is the overlap.
• ? is the empirical geometric prefactor (dimensionless input parameter that accounts for the particle
geometry).
The empirical geometric prefactor = 1.71 (default value) was calibrated with experimental data
from an uniaxial compression test of elastic rubber spheres between two rigid plates (Tatara,
2021[43] (p. 379)). The calibration is performed in a way to capture the maximum force exerted on
the top plate and also to follow the path during compression. The parameter depends on the ma-
terial properties, especially on the stiffness. For a given material, the calibration should be done once
and then be used for other test configurations without the need of re-calibration.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 339
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
340 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 39: TBM Breakage Model Module
39.1. Usage details
This breakage specific impact energy limit is determined by Rocky based upon the parameter values
you enter for the model.
This model also requires use of the Fragment distribution type TBM Breakage Model, which is a
uniform size distribution. This means that all broken fragments generated from the TBM Breakage
Model will be the same size. In addition, the minimum number of fragments in a breakage event is
2, with a maximum number you can define via the Maximum number of fragments parameter.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 341
TBM Breakage Model Module
Figure 39.1: Breakage | Criteria tab options (left) and Breakage | Fragments tab options (right)
for the Particle set when the TBM Breakage Model module is enabled.
• A custom instantaneous Breakage Model where calculations will be applied according to the model
developed by TAKRAF, which uses statistical data from experimental breakage tests.
In addition, because this module relies on the same code structure used by the Breakage Model and
by the Fragments Distribution Model included by default in Rocky, after enabling this module you will
therefore need to select it when defining your particle breakage criteria.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
342 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
However, there are additional module-related parameters that you must set for each Particle set. Refer
to the next section for details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 343
TBM Breakage Model Module
Figure 39.3: TBM Breakage Model module options for a Particle set.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
344 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Besides that, a new Property is available for the main Particles entity, as shown in and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 345
TBM Breakage Model Module
Figure 39.4: Particles Property created by the TBM Breakage Model module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
• Breakage Specific Energy: For a given particle, this provides the value which is compared with
the breakage specific energy limit to determine whether the particle will break.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
346 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
2. Set up the simulation as you normally would, with the following exception:
• On the Particle | Breakage | Criteria tab for at least one Particle set, you have defined the
Model as TBM Breakage Model.
• On the Particle | Breakage | Fragments tab for at least one Particle set, you have defined
the Distribution model as TBM Breakage Model.
• On the Particle | Modules tab, you have entered the values you want.
4. When you are ready to analyze your simulation results, do one or more of the following:
• Make use of the default breakage Curves and/or Properties for Particles.
The breakage specific energy limit, , in the TBM model is determined using statistical parameters
from experimental tests. depicts a typical log-log plot of the breakage specific energy versus the
particle size, obtained in a sequence of experimental tests. The equation of a line that best fits the
experimental points in the log-log plane can be determined by means of linear regression. This line
defines a mean value of breakage specific energy, , for a given particle size, . Then, the rela-
tionship between those two variables can be expressed simply as:
(39.1)
where and are the coordinates of a point on the fitted line, while is its slope in the log-
log plot. In the Rocky UI, and are listed as Reference Size and Reference Specific Energy,
respectively, while re is listed as Slope Coefficient.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 347
TBM Breakage Model Module
Assuming that the experimental points obey a log-normal distribution, lower and upper limits for a
distribution of breakage probability can be established as:
(39.2)
(39.3)
where is the standard deviation of the experimental data set in relation to the fitted line, given
by:
(39.4)
in which and , for , are measured values of particle size and breakage specific
energy, respectively. For further details, please refer to Schmidt (2011) [44] (p. 379).
In addition, in Equation 39.2 (p. 348) and Equation 39.3 (p. 348) is a parameter that determines the
amplitude of the band limited by and . For example, for , the band defined will
include 95% of the measured values. In the Rocky UI, and are listed as Standard Deviation and
Scatter Band Parameter, respectively.
In Rocky, a random strength value is assigned to every particle in a simulation during the
particle generation process. The corresponding breakage specific energy, assuming a cumulative log-
normal distribution, will be given by:
(39.5)
However, in the TBM model, breakage is only allowed when . On the other hand, if ,
the particle will always break. This effect is achieved by defining the breakage specific energy as:
(39.6)
Thus, a particle will break if at any given moment the instantaneous specific impact energy in the
particle is larger than the specific breakage energy, , given by Equation 39.6 (p. 348). At a breakage
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
348 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
event, the size of the resulting fragments will be determined by the described in the following section,
which is implemented within the same module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 349
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
350 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 40: Thermal Radiation Module
40.1. Usage details
Note:
Geometries, however, are not able to receive heat transfer from particles. See the Limitations
and other considerations (p. 353) section for more details.
Typically, the heat transfer due to radiation becomes significant at higher temperatures, such as 700
K or greater. Therefore, the effects of thermal radiation can play an important role in some high-heat
applications, such as industrial heating, cooling, and drying processes; as well as combustion processes
and solar radiation.
Figure 40.1: Heating particles inside a cylinder using the Thermal Radiation module. The upper surface
is set to a high temperature, which heats the particles at the top. The particles then exchange heat
with each other in the lower layers. (p. 351) shows an example of how this module can be applied. In
a cylinder and particle simulation set up with Thermal Modeling, the lid geometry at the top of the
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 351
Thermal Radiation Module
cylinder is defined—via its assigned material—as an emissive geometry, and a high prescribed tem-
perature is applied to the lid. The group of particles within the cylinder is defined—via its assigned
material—as emissive, which enables the particles within the group to receive heat transfer via radi-
ation.
A screening distance defines the area around the lid as the radiation heat source for the affected
particles. Below the lid, a second screening distance defines the area encompassing the particles as
the location within which radiation heat can be exchanged.
As the temperatures of the top layer of particles increase due to radiation effects from the lid, those
top particles exchange heat with each other inside the cylinder by conduction and radiation. Thus,
the heat spreads from the particles at the top of the cylinder through to the particles at the bottom.
For each emissive particle-particle or triangle-particle pair, a radiation view factor that takes into account
their respective sizes, angular positions and the distances between each other is automatically calcu-
lated. This radiation view factor is then used to consider the fraction of the radiation leaving the
source surface that is intercepted by the target surface.
It is important to remember that by enabling the Thermal Model, heat transfer will also occur due to
conduction and if/when CFD Coupling is enabled convection in addition to the radiation effects enabled
by this module.
Depending upon what kind of thermal effects you want to simulate in Rocky, you may need to set
up your simulation differently. Below is a summary of how you may apply the three different types
of thermal modes:
• To simulate both conduction and convection, enable the Thermal Model and then couple the
simulation with a fluid case (CFD) that has energy calculations.
• To simulate both conduction and radiation, enable the Thermal Model and the Thermal Radiation
module.
• To simulate the three types of heat transfer conduction, convection, and radiation enable the
Thermal Model and the Thermal Radiation module, and then couple the simulation with a fluid
case (CFD) that has energy calculations.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
352 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
This module does not take into account the occlusion effects when computing the radiation view
factors.
It is also important to note that in Rocky, geometries do not receive any heat transfer from particles.
For geometries, temperature is defined as a constant value and therefore, they may only participate
as heat sources inside the simulation. This is true for radiation effects as well as for other types of
heat transfer (conduction and convection).
When using the Feed Conveyor template to create an emissive geometry, be aware that the belt element
will not participate in the radiation calculations.
• Particles Screening Distance: Defines the maximum influence distance in which the particle-particle
interactions will be accounted. Emissive particles within this location will exchange radiation heat
between them.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 353
Thermal Radiation Module
• Geometries Screening Distance: Defines the maximum influence distance in which the triangle-
particle interactions will be accounted. Each emissive geometry present in the simulation within
this location will participate as a radiation heat source for the respective particles.
• Start Time: Defines the time you want particles within the screening distance to begin being influ-
enced by the thermal radiation.
• Stop Time: Defines the time you want particles within the screening distance to stop being influ-
enced by the thermal radiation.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
354 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 40.3: Thermal Radiation module options for an individual material component.
• Material Emissivity: This parameter defines the thermal emissivity that will be applied to the
selected material. When a material with emissivity greater than zero is assigned to a geometry,
that geometry becomes emissive and can act as a radiation heat source. When a material with
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 355
Thermal Radiation Module
emissivity greater than zero is assigned to a particle group, those particles become emissive
and can receive and exchange radiation heat transfer.
Note:
Even though you are able to define this parameter for any material, only those
materials assigned to geometries or particle groups within the Geometries
Screening Distance and/or Particles Screening Distance will participate in the ra-
diation heat transfer process.
Any materials with this parameter set to 0 will have their assigned geometries
and/or particle groups ignored by radiation effects.
Range: [0 - 1]
• Define the Geometries Screening Distance no lower than 2.5 times the largest triangle size of
any emissive geometry present in the simulation.
• For all emissive geometries, define the geometry Triangle Size as equal to or up to 2 times bigger
than the equivalent diameter of the smallest particle present in the simulation.
• For your geometries, avoid defining Triangle Size values smaller than the equivalent diameter of
the smallest particle present in the simulation. Triangles that are too small can reduce the solver
performance, increasing the computational cost with no substantial enhancement in the results
accuracy.
• For simulations with high particle density, avoid using a Particles Sreening Distance larger than
2 times the equivalent dimeter of the biggest particle present in simulation. Otherwise, it can
overestimate the heat transfer due to occlusions neglect, reduce the solver performance, and increase
the computational cost.
• If you have geometries or particle groups in your simulation that you do not want to be considered
as a source for thermal radiation, ensure that you assign a material to these geometries and/or
particle groups that has a Material Emissivity equal to 0. By doing so, you not only deactivate the
radiation heat transfer calculation for these items but also improve the solver performance.
By following these recommendations, the Thermal Radiation model will deliver confident results to-
gether with good solver performance.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
356 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
2. Ensure that the module is enabled. (From the Data panel, select Modules and then from the Data
Editors panel, ensure that the Thermal Radiation checkbox is enabled.)
6. Prior to processing your simulation, ensure that you do the following for each emissive geometry
in your project:
a. From the Data panel, under Materials, select a material that you want to assign to an emissive
geometry, and then from the Data Editors panel define the Material Emissivity.
b. From the Data panel, under Geometries, select a geometry that you want to consider as an
emissive geometry, and then from the Data Editors panel do the following:
i. From the Geometry tab, select from the Material list the emissive material previously
defined in step 6a.
ii. From the Geometry | Geometry tab, enable the Prescribed Temperature option for
Thermal Boundary Type, and then define the emissive geometry's Temperature.
7. Moreover, ensure that you do the following for each emissive particle group in your project:
a. From the Data panel, under Materials, select a material that you want to assign to an emissive
particle group, and then from the Data Editors panel define the Material Emissivity.
b. From the Data panel, under Particles, select a particle group that you want to consider as an
emissive particle group, and then from the Data Editors panel do the following:
i. From the Particle tab, select from the Material list the emissive material previously defined
in step 7a.
8. Ensure that you defined initial temperature for the emissive particles in your project:
a. From the Data panel, under Inputs, select the particle input that you want to define the initial
temperature, and then from the Data Editors panel do the following:
i. If you selected a Continuous Injection or a Volume Fill, from the Particles tab, define the
temperature for each particle group.
ii. If you selected a Custom Input, define the Default Temperature value or define the tem-
perature values directly in the csv file you want to import.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 357
Thermal Radiation Module
During the simulation, each emissive geometry present in the simulation within the defined Geometries
Screening Distance will participate as a radiation heat source for the respective particles. Also, emissive
particles within the defined Particles Screening Distance will exchange radiation heat between each
other.
Figure 40.4: Illustration of the calculation domain for radiation heat transfer according to the
defined screening distances.
In the example, the green particles and green line represent the items experiencing radiation effects,
as defined by the screening distance. The "i" in the center of this area and the wavy arrows emitting
from it represent the radiation interactions happening between participating particles and geometries.
• Gray bodies are physical entities that partially emit and absorb radiation. The level of emission and
absorption is determined by the emissivity parameters defined in the module. Emissivity is a surface
radiation property defined as the ratio of the radiation emitted by the surface to the radiation
emitted by a black body at the same temperature. A Black body, in turn, represents an idealized,
perfect absorber and emitter of radiation.
• Bodies with diffuse surfaces emit and absorb radiation equally in all directions.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
358 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
40.2.2. Calculations
In the Thermal Radiation module, the particle-particle and particle-triangle heat exchange is calculated
using the thermal radiation equation for heat transfer between two gray, diffuse and opaque surfaces,
defined as follows:[45] (p. 379)[48] (p. 379)[47] (p. 379)
(40.1)
where:
• is the source surface index, which can be either a geometry triangle or particle within the defined
screening distance.
• is the target surface index, which can be a particle within the defined screening distance.
• is the radiation heat transfer rate between the source surface, , and the target surface, .
Note:
Note:
• is the radiation view factor, which is defined as the fraction of the radiation leaving the source
surface, , that is intercepted by the target surface, .
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 359
Thermal Radiation Module
The particle-particle view factor, , for any sphere sizes and spacings is defined as follows:[49] (p. 379)
(40.2)
where:
The radiation view factor for each triangle-particle interaction is also independent of the target particle
shape. This factor is calculated considering the approximation of a small surface to a sphere with
equivalent diameter, as shown in Figure 40.6: Triangle-particle interaction. (p. 361). The view factor for
a small surface to a sphere is defined as follows:
(40.3)
where:
• is the distance between the target particle centroid and the geometry triangle centroid.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
360 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• is the angle between the triangle surface normal vector and the particle position vector.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 361
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
362 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 41: Velocity-Dependent Adhesive Force
Module
41.1. Usage details
Note also that enabling this module will override what you have set for Adhesive Force on the Mo-
mentum tab of the Physics entity (Figure 41.1:Physics | Momentum tab options when the Velocity-
Dependent Adhesive Force module is enabled. (p. 364)).
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 363
Velocity-Dependent Adhesive Force Module
Figure 41.1: Physics | Momentum tab options when the Velocity-Dependent Adhesive Force module
is enabled.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
364 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
When you select a pair of materials interactions, two new parameters will appear as shown in Fig-
ure 41.3:Velocity-Dependent Adhesive Force module options for a pair of materials interactions. (p. 365)
and explained below.
Figure 41.3: Velocity-Dependent Adhesive Force module options for a pair of materials
interactions.
• Diameter Exponent: This defines the diameter exponent of the parametric adhesive force
model.
• Force Coefficient: This defines the force coefficient of the parametric adhesive force model.
Note:
This parameter is defined with the unit [kg/m(1+b)] (with b being the Diameter
Exponent) but as Rocky does not recognize this unit, the parameter will be shown
as dimensionless in the Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 365
Velocity-Dependent Adhesive Force Module
In addition, if you chose to enable the Collect Contacts Data checkbox on the Contacts entity, a
new Property is available for Contacts, as shown in Figure 41.4: Contacts Properties created by the
Velocity-Dependent Adhesive Force module. (p. 366) and explained below.
Figure 41.4: Contacts Properties created by the Velocity-Dependent Adhesive Force module.
• Adhesive Force: For a given contact, this provides the adhesive force calculated by the module.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
2. Set up your simulation as you normally would, with the following exceptions:
• On the Materials Interactions step, from the Data Editors panel, select each pair of materials
interactions that you want to use in your simulation, and then under Velocity-Dependent
Adhesive Force, define the parameters you want.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
366 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
• (Optional) On the Contacts step, from the Data Editors panel, enable the Collect Contacts
Data checkbox. Collecting contacts will give you access to a new Contact property generated
by this module, which is used for calculation purposes.
4. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related Property. Or if you chose to collect contacts data, you may also make use of that
additional module-related property.
(41.2)
where:
• is the force coefficient, measured in [kg/m(1+b)], which is defined by the user for each pair of mater-
ials on the Materials Interactions step when the Velocity-Dependent Adhesive Force module is enabled.
• is the diameter exponent. This dimensionless parameter is defined by the user for each pair of
materials on the Materials Interactions step when the Velocity-Dependent Adhesive Force module is
enabled.
• and are the equivalent diameters of the two contacting particles, measured in [m]. The
equivalent diameter is defined as the diameter of a sphere whose volume is equal to the particle?s
volume.
• and are the magnitudes of the normal and tangential components of the relative velocity at
the beginning of the contact, respectively, measured in [m/s].
• is the normal restitution coefficient. This dimensionless parameter is defined by the user for each
pair of materials on the Materials Interactions setup.
Note:
Both Constant and Velocity Dependent restitution model are compatible with this
module.
The tangential restitution coefficient ( ) used in equation Equation 41.1 (p. 367) is computed as:
[37] (p. 379)
(41.3)
where:
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 367
Velocity-Dependent Adhesive Force Module
• is the dynamic friction coefficient. This dimensionless parameter is defined by the user for each
pair of materials on the Materials Interactions setup.
• is a parameter that characterizes the mass distribution within the particles. For homogeneous
spheres, this parameter takes the value 2/7, which is the value used for the implementation of this
module.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
368 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 42: Velocity-Dependent Restitution Model
42.1. Usage details
Also note that enabling this module will override what you have set for Restitution Coefficient in
the Materials Interactions entity.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 369
Velocity-Dependent Restitution Model
When you select a pair of materials interactions, 4 new parameters will appear as shown in Fig-
ure 42.2:Velocity-Dependent Restitution Model module options for a pair of materials interac-
tions. (p. 371) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
370 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 42.2: Velocity-Dependent Restitution Model module options for a pair of materials
interactions.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 371
Velocity-Dependent Restitution Model
• Enable for This Materials Interaction: This enables or disables the velocity dependent restitution
model.
Note:
If Enable for This Materials Interaction is disabled, the restitution coefficient will
be constant, with value equal to Restitution Coefficient Base Value.
• Restitution Coefficient Base Value: This defines the base value for the restitution coefficient.
• Velocity Exponent: This defines the velocity exponent of the expression that calculates the ve-
locity dependent restitution.
• Velocity Limit: This defines a maximum velocity value, below which the coefficient of restitution
is considered constant.
In addition, if you chose to enable the Collect Contacts Data checkbox on the Contacts entity, a
new Property will be available for Contacts, as shown in Figure 42.3: Contacts Property created by the
Velocity-Dependent Restitution Model module. (p. 373) and explained below.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
372 of ANSYS, Inc. and its subsidiaries and affiliates.
Usage details
Figure 42.3: Contacts Property created by the Velocity-Dependent Restitution Model module.
• Restitution Coefficient: Provides the restitution coefficient calculated by the module to a given
contact.
Note:
Properties generated by this module will only appear if the Upscaling is disabled.Refer to
Rocky User Manual (5.1.7. About Meshed Particles Upscaling) for more details.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 373
Velocity-Dependent Restitution Model
2. Set up your simulation as you normally would, with the following exceptions:
• On the Materials Interactions step, from the Data Editors panel, select each pair of materials
interactions that you want to use in your simulation, and then under Velocity-Dependent
Restitution Model, define the parameters you want.
• (Optional) On the Contacts step, from the Data Editors panel, enable the Collect Contacts
Data checkbox. Collecting contacts will give you access to a new Contact property generated
by this module, which is used for calculation purposes.
4. When you are ready to analyze your simulation results, you may choose to make use of any
particle-related Property. Or if you chose to collect contacts data, you may also make use of that
additional module-related property.
(42.1)
where:
• is an input parameter labeled as Restitution Coefficient Base Value in the Rocky UI.
• is a limit velocity, below which the coefficient of restitution is considered constant. This is also
an input parameter, which is listed in the Rocky UI as Velocity Limit.
• is a exponent whose value specifies how quickly the coefficient of restitution decreases when the
impact velocity grows. Figure 42.4: Coefficient of restitution as a function of impact velocity, for dif-
ferent values of the exponent K. (p. 375) shows four curves obtained from Equation 42.1 (p. 374) for
different values of the exponent . This parameter is listed as Velocity Exponent in the Rocky UI.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
374 of ANSYS, Inc. and its subsidiaries and affiliates.
Technical details
Figure 42.4: Coefficient of restitution as a function of impact velocity, for different values of the
exponent K.
The impact velocity considered in Equation 42.1 (p. 374) is defined as the normal component of the rel-
ative velocity at the contact point, , at the time when a collision starts. That is:
(42.2)
where in the contact's normal unit vector and is the time when the collision starts, that is, the
time when the physical contact between the colliding entities begins. In turn, the relative velocity
is defined as:
(42.3)
where, for a particle-to-particle collision:
• is the vector that joins the center of mass of the particle to the instantaneous contact point
. The definition is equivalent for the vector .
In the case of a particle-to-boundary collision, and will be, respectively, the translational and ro-
tational velocities of the boundary, whereas will be the vector that joins the rotation center of the
boundary to the contact point.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 375
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
376 of ANSYS, Inc. and its subsidiaries and affiliates.
Bibliography
[1] Akinci et al.. Versatile surface tension and adhesion for SPH fluids.. ACM Trans. Graph. 32, 6, Article
182. 8 pages. 2013.
[2] Gimenez, C. S. and Finke, B. and Nowak, C. and Schilde, C. and Kwade, A. Structural and mechanical
characterization of lithium-ion battery electrodes via {DEM} simulations. Advanced Powder Tech-
nology. 29. 2312-2321. 2018.
[3] Potyondy, D. O. and Cundall, P. A.. A bonded-particle model for rock. International Journal of Rock
Mechanics & Mining Sciences. 41. 1329-1364. 2004.
[4] Smith, IWi. The combustion rates of coal chars: a review. Symposium (International) on combustion.
19. 1045-1065. 1982.
[5] Nicoud, Lucre ce and Licordari, Filippo and Myerson, Allan S. Estimation of the solubility of metastable
polymorphs: a critical review. Crystal Growth & Design. 18. 7228-7237. 2018.
[6] Dhole, S.D. and Chhabra, R.P. and Eswaran, V.. Flow of power-law fluids past a sphere at intermediate
Reynolds numbers. Industrial & Engineering Chemistry Research. . . 2006.
[7] Bird, R.B. and Armstrong, R.C. and Hassager, O.. Dynamics of polymeric liquids. John Wiley & Sons. .
. 1987.
[8] Dazhi, G. and Tanner, R.I.. The drag on a sphere in a power-law fluid. Journal of Non-Newtonian Fluid
Mechanics. . . 1985.
[9] Wang, Zekun and Liu, Moubin and Yang, Xi. A four-way coupled CFD-DEM modeling framework for
charged particles under electrical field with applications to gas insulated switchgears. Powder
Technology. 373. 433-445 . 2020.
[10] . ANSYS\textsuperscript{\tiny\textregistered} Fluent, Release 2019 R3, Help System, Theory Guide. ANSYS,
Inc.. 29=019.
[11] Svoboda, J. Extraction of ferromagnetic particles by a suspended magnet. Magnetic and Electrical
Separation. 4. 223-239. 1993.
[12] Tom Leblicq and Bart Smeets and Simon Vanmaercke and Herman Ramon and Wouter Saeys. A
discrete element approach for modelling bendable crop stems. Computers and Electronics in
Agriculture. 124. 141-149. 2016.
[13] Pasha, M. and Dogbe, S. and Hare, C. and Hassanpour, A. and Ghadiri, M.. A linear model of elasto-
plastic and adhesive contact deformation. Granular Matter. 16. 151-162. 2014.
[14] Massih Pasha. Modelling of flowability measurement of cohesive powders using small quantities.
PhD Thesis. University of Leeds. 2013.
[15] Mikami, T. and Kamiya, H. and Horio, M.. Numerical simulation of cohesive powder behavior in a
fluidized bed. Chemical Engineering Science. 53. 1927-1940. 1997.
[16] Nase, S. T. and Vargas, W. L. and Abatan, A. A. and McCarthy, J. J.. Discrete characterization tools
for cohesive granular material. Powder Technology. 116. 214-223. 2001.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 377
Bibliography
[17] Vasilakis, Andrew. The Lorentz force law and Maxwell's equations in symmetric form. International
Journal of Mathematical Education in Science and Technology. 67. 443-446. 1994.
[19] Carlo A. Gonano*, Riccardo E. Zich, and Marco Mussetta. Definition for Polarization P and Magnet-
ization M Fully Consistent with Maxwell’s Equations. Progress In Electromagnetics Research. 64.
83–101. 2015.
[20] . Contact radius and curvature corrections to the nonlocal contact formulation accounting for
multi-particle interactions in elastic confined granular systems. International Journal of Engineering
Science. 133. 26–46. 2018.
[21] Y. Bommireddy, A. Agarwal, V. Yettella, V. Tomar, and M. Gonzalez. Loading-unloading contact law
for micro-crystalline cellulose particles under large deformations. arXiv:1812.06173. arXiv preprint.
. 2018.
[22] Flemmer, R. L. C. and Banks, C. L.. On the Drag Coefficient of a Shape. Powder Technology. 48. 217-
221. 1986.
[23] Pritchard, Philip J.. Fox and McDonald's Introduction to Fluid Mechanics. John Wiley & Sons. 2011.
[24] Crowe, C.T. and Schwarzkopf, J.D. and Sommerfeld, M. and Tsuji, Y.. Multiphase Flows with Droplets
and Particles, Second Edition. Taylor & Francis. 2011.
[25] . Ansys\textsuperscript{\tiny\textregistered} Fluent, Release 2019 R3, Help System, User's Guide. Ansys,
Inc.. 2019.
[27] . Ansys\textsuperscript{\tiny\textregistered} CFX, Release 2019 R3, Help System, Theory Guide. Ansys,
Inc.. 2019.
[28] Unland, Georg and Szczelina, Piotr. Coarse crushing of brittle rocks by compression. International
Journal of Mineral Processing. 74. . 2004.
[29] Huilin, Lu and Gidaspow, Dimitri. Hydrodynamics of binary Fluidization in a riser: CFD simulation
using two granular temperatures. Chemical Engineering Science. 58. 3777-3792. 2003.
[30] Gidaspow, D.. Multiphase Flow and Fluidization. Academic Press. 1994.
[31] Ergun, Sabri. Fluid flow through packed columns. Chem. Eng. Prog.. 48. 89-94. 1952.
[32] Gidaspow, D.. Multiphase Flow and Fluidization: Continuum and Kinetic Theory Descriptions. Elsevier
Science. 2012.
[33] Gidaspow, Dimitri and Bezburuah, Rukmini and Ding, J.. Hydrodynamics of Circulating Fluidized
Beds: Kinetic Theory Approach. Engineering foundation conference on fluidization. 7. 75-82. 1993.
[34] Wen, C.Y. and Yu, Y.H.. Mechanics of Fluidization. The Chemical Engineering Progress Symposium
Series. 162. 100-111. 1966.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
378 of ANSYS, Inc. and its subsidiaries and affiliates.
[35] Bergman, T. and Incropera, F. and Dewitt, D. and Lavine, A.. Fundamentals of Heat and Mass
Transfer. Wiley. . . 2011.
[36] E., M.. Particles, Bubbles & Drops: Their Motion, Heat and Mass Transfer. World Scientific Publishing
Company. . . 2006.
[37] Schwager, T. and Becker, V. and Pöschel, T.. Coefficient of tangential restitution for viscoelastic
spheres. The European Physical Journal E. . . 2008.
[38] Christoph Gissler and Stefan Band and Andreas Peer and Markus Ihmsen and Matthias Teschner.
Generalized drag force for particle-based simulations. Computers & Graphics. 69. 1-11. 2017.
[39] O'Rourke, P J and Amsden, A A. The TAB method for numerical calculation of spray droplet breakup.
. . . 1987.
[40] Liu, Alex B. and Mather, Daniel and Reitz, Rolf D.. Modeling the Effects of Drop Drag and Breakup
on Fuel Sprays. International Congress & Exposition. SAE International. . . 1993.
[41] O'Rourke, Peter J. and Amsden, Anthony A.. The Tab Method for Numerical Calculation of Spray
Droplet Breakup. 1987 SAE International Fall Fuels and Lubricants Meeting and Exhibition. . . 1987.
[42] Giannis, K. and Schilde, C. and Finke, J. and Kwade, A. and Caligueta, M. and Taghizadeh, K. and
Luding, S.. Stress based multi-contact model for discrete-element simulations. Granular Matter.
. . 2021.
[43] Tatara, Y.. Extensive theory of force-approach relations of elastic spheres in compression and in
impact. J. Eng. Mater. Technol.. 111. 163–168. 2021.
[44] Schmidt, M.. Modellierung der Zerkleinerung in Profilwalzenbrechern. Technischen Universität Ber-
gakademie Freiberg. 2011.
[45] Bergman, T. L. and Lavine, A. S. and Incropera, F. P. and Dewitt, D. P.. Radiation Exchange between
Surfaces. . . 861-931. .
[46] Bergman, T. L. and Lavine, A. S. and Incropera, F. P. and Dewitt, D. P.. Fundamentals of Heat and
Mass Transfer. Wiley & Sons. 2011.
[47] Hou, Q. and Gan, J. and Zhou, Z. and Yu, A.. Particle Scale Study of Heat Transfer in Packed and Fluidized
Beds. Advances in Chemical Engineering. 46. 2015.
[48] Chen, J. C. and Grace, J.R. and Golriz, M. R.. Heat transfer in fluidized beds: design method. Powder
Technology. 150. 123-132. 2005.
[49] Felske, J. D.. Approximate Radiation Shape Factors between Two Spheres. Journal of Heat Transfer.
100. 547-548. 1978.
[50] Hauptmann, E.G.. Angle Factors between a Small Flat Plate and a Diffusely Radiating Sphere.
American Institute of Aeronautics and Astronautics (AIAA) Journal. 6. 938-939. 1967.
[51] Hoghe, M. D. and Calle, C. I. and Curry, D. R.and Weitzman P. S.. Discrete element modeling ({DEM})
of triboelectrically charged particles: Revised experiments. Journal of Electrostatics. 67. 691-694.
2009.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 379
Bibliography
[52] Wu, C. and Li, L. and Thornton, C.. Rebound behaviour of spheres for plastic impacts. International
Journal of Impact Engineering. 28. 929–946. 2003.
[53] S. M. Hosseini, M. T. Manzari, and S. K. Hannani. A fully explicit three-step SPH algorithm for simu-
lation of non-newtonian fluid flow.. International Journal of Numerical Methods for Heat and Fluid
Flow. 17. 715-735. 2006.
[54] M. Ihmsen, J. Cornelis, B. Solenthaler, C. Horvath, and M. Teschner. Implicit incompressible SPH..
IEEE Transactions on Visualization and Computer Graphics. 20. 426-435. 2014.
[55] Q. Zhang, L. Dor, A. K. Biswas, W. Yang and W. Blasiak. Modeling of steam plasma gasification for
municipal solid waste.. Fuel Processing Technology. 106. 546-554. 2013.
[56] S. Wang and Y. Shen. CFD-DEM modelling of raceway dynamics and coke combustion in an iron-
making blast furnace.. Fuel. 302. 2021.
[57] Pei, C., Wu, C.-Y, England, D., Byard, S., Berchtold, H. and Adams, M.. DEM-CFD modeling of particle
systems with long-range electrostatic interactions. AIChE. 61 . 1792-1803 . 2015.
[58] Walker, J., Halliday, D., Resnick, R.. Fundamentals of Physics, 10th edition. 2 . 2014.
[59] Tinkle, M. and Barlow, S.. Image Charge Forces Inside Conducting Boundaries. Non-Neutral Plasma
Physics IV. 525-530 . 2002.
[60] Alobaid, F., Almohammed, N., Farid, M. M., May, J., Rössger, P., Richter, A., Epple, B.. Progress in CFD
Simulations of Fluidized Beds for Chemical and Energy Process Engineering. Progress in Energy
and Combustion Science. 2021.
[61] Greason, W.. Investigation of a test methodology for triboelectrification. Journal of Electrostatics.
245-256 . 2000.
[63] Kholmetskii, A., Missevitch, O.. Electric/magnetic dipole in an electromagnetic field: force, torque
and energy. European Physical Journal Plus. 1-13. 2014.
[64] Thornton, C., Ning, Z.. A theoretical model for the stick/bounce behaviour of adhesive, elastic-plastic
spheres. Powder Technology. 1998.
[65] Giménez, C.S., Finke, B., Shilde, C., Froböse, L., Kwade, A.. Numerical simulation of the behavior of
lithium-ion battery electrodes during the calendaring process via the discrete element method.
Powder Technology. 2019.
Release 2025 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
380 of ANSYS, Inc. and its subsidiaries and affiliates.