fix buoyancy command

Syntax

fix ID group-ID buoyancy keyword value ...
  • ID, group-ID are documented in fix command
  • buoyancy = style name of this fix command
  • zero or more keyword/value pairs may be appended to args; the dim keyword must be used
  • keyword = density or dim or level or region
density value = density-value
  density-value = fluid density in mass/length^3 units
dim value = x or y or z
  x, y, z = define the water level as x=const, y=const, z=const plane
level value = lev
  lev = water level in length units, along the axis defined by dim
region value = region-ID
  region-ID = ID of region atoms must be in to have added buoyancy force

Examples

fix bu all buoyancy level 0.06 dim z density 1000

Description

Add a buoyancy force for each atom in the group. The water level is assumed to be a x=const, y=const or z=const plane, where the axis is defined via the dim keyword and const is defined via the level keyword. E.g. dim = x and level = 0.1 would define the water level as x = 0.1.

The buoyancy force is equivalent to the weight of the displaced fluid. Thus, if the particle is fully submerged, the buoyancy force is equal to particle volume * fluid density. If the particle is not submerged, the buoyancy force is 0. If the particle is partially submerged, the fix will calculate the force based on the submerged particle volume. The fluid density is defined by keyword density.

As buoyancy is triggered by static pressure difference in a fluid usually caused by gravity, this fix requires to use a fix gravity. The gravity vector has to be specified in the same axis as in this fix (via the dim keyword).

If the region keyword is used, the atom must also be in the specified geometric region in order to have buoyancy force added to it.

Restart, fix_modify, output, run start/stop, minimize info

No information about this fix is written to binary restart files. No fix_modify option applies to this fix. This fix computes a a global 3-vector of the total buoyancy force, which can be accessed by various output commands. No parameter of this fix can be used with the start/stop keywords of the run command.

The forces due to this fix are imposed during an energy minimization, invoked by the minimize command.

Restrictions

This fix requires to use a fix gravity. The gravity vector has to be specified in the same axis as in this fix (via the dim keyword).

Default

level = 0, density = 0