PIC "default" settings

MFIX-Exa Help
1

Is there a set of “default” (read: suggested starting point) PIC parameter settings for setting up a new case. I have something like this (below) but I’m not sure if this is a good default or if I’ve adjusted these out of whack for some problem years ago.

#----------------------------------------------------------------------
# PIC model settings
#----------------------------------------------------------------------
#dem pic.solve = None 
pic.solve = solid0

pic.verbose = 1

pic.max_iter = 3

pic.beta = 3.0
pic.pressure_coefficient = 100.0

pic.close_pack = 0.60
pic.parcels_per_cell_at_pack  =  30.

pic.small_number = 1.0e-8

pic.damping_factor = 0.40
pic.damping_factor_wall_tangent = 0.9
pic.damping_factor_wall_normal  = 0.7

pic.velocity_reference_frame = 0.50

pic.advance_vel_p = 0.5

pic.initial_step_type = nth_eps
2

What you posted is a great starting point for PIC settings and should work well for most simulations.

  • pic.close_pack should match whatever corresponds to the solids you are modeling.

  • pic.parcels_per_cell_at_pack can be thought of in terms of resolution with 24 being an approximate lower bound. This input controls how much effective volume each parcel represents. If you use a small value (e.g., 4), each parcel represents a lot of solids volume. Therefore, small changes in parcel positions will cause large swings in the solids volume fraction, and in turn, the solids stress.

  • I hardly ever use different values for beta an pressure_coefficient. The same can be said about the damping_factor settings.

As for the rest – I do not recommend changing them unless you have a good understanding of how they affect the PIC update algorithm. Arbitrarily setting initial_step_type, advance_vel_p, or velocity_reference_frame to different values most likely will not end well.

3

you can’t go as low as 4. you put a limit at 12. I tried to force it so I could get a stat wt of 1 on a dx = dp and got swatted at haha.

4

you can’t go as low as 4. you put a limit at 12.

I regret nothing.