Acoustophoretic Radiation Force
The Acoustophoretic Radiation Force node is only available when Newtonian or Newtonian, first order is selected as the Formulation on the Settings window for The Particle Tracing for Fluid Flow Interface.
Use the Acoustophoretic Radiation Force node to exert forces on small particles due to acoustic radiation. To be subjected to this force, particles must be in a region where the acoustic pressure or acoustic velocity is not spatially uniform.
In addition, the Acoustophoretic Radiation Force is only applied if the acoustic pressure or acoustic velocity has already been computed in a Frequency Domain study type. This is because the magnitude of the acoustophoretic force depends on the angular frequency.
Radiation Force Model
This section determines how the equation for the acoustophoretic radiation force is defined. It also determines what inputs will be shown in subsequent sections.
Select an option from the Particle type list: Solid particle (the default) or Liquid droplet.
Select an option from the Thermodynamic loss model: Ideal (the default), Viscous, or Thermoviscous. The Thermoviscous loss model is the most detailed of the three but requires the most in-depth knowledge of the particle and fluid material properties.
Acoustic Fields
Enter a value or expression for the following:
Pressure p (SI unit: Pa). This should be computed from another physics interface using a Frequency Domain study type.
Acoustic velocity u (SI unit: m/s). Similarly, this should be solved for in another physics interface such as the Pressure Acoustics, Frequency Domain interface.
Particle Material Properties
Enter the material properties of the solid particles or liquid droplet in this section. Some additional particle properties, such as density, are instead specified in the settings for the Particle Properties node.
For all thermodynamic loss models, if the Particle type is Solid particle, enter values or expressions for the following:
Pressure-wave speed cp,p (SI unit: m/s, default 2,400 m/s), and
Shear-wave speed cs,p (SI unit: m/s, default 1150 m/s).
If instead the Particle type is Liquid droplet, just enter the Adiabatic speed of sound cp (SI unit: m/s, default 1445 m/s).
If the Thermodynamic loss model is Viscous or Thermoviscous, and the Particle type is Liquid droplet, also enter the Particle dynamic viscosity μp (SI unit: Pa·s, default 4.2 × 10-3 Pa·s).
If the Thermodynamic loss model is Thermoviscous and the Particle type is Liquid droplet, also enter the Bulk viscosity particle μB,p (SI unit: Pa·s, default 0).
If the Thermodynamic loss model is Thermoviscous, for either Particle Type enter the:
Heat capacity at constant pressure Cp,p (SI unit: J/(kg·K),
default 2.1 × 103 J/(kg·K)),
Isobaric coefficient of thermal expansion αp,p (SI unit: 1/K, default 7.05 1/K,)
Ratio of specific heats γp (dimensionless, default 1.15), and
Thermal conductivity kp (SI unit: W/(m K), default 0.17 W/(m K)). The thermal conductivity is always assumed to be a scalar for the purpose of computing the acoustophoretic radiation force.
Surrounding Fluid Properties
All of the inputs in this section are domain material properties that can either be taken From material or given a User defined expression. The defaults given below are shown by first selecting User defined.
For any Thermodynamic loss model, enter the:
Density ρ (SI unit: kg/m3, default 995 kg/m3), and
Speed of sound c (SI unit: m/s, default 1500 m/s).
If the Thermodynamic loss model is Viscous or Thermoviscous, also enter the Dynamic viscosity μ (SI unit: Pa·s, default 8.4 × 10-4 Pa·s).
If the Thermodynamic loss model is Thermoviscous, also enter the:
Bulk viscosity μB (SI unit: Pa·s, default 2.4 × 10-3 Pa·s),
Heat capacity at constant pressure Cp (SI unit: J/(kg·K),
default 4.18 × 103 J/(kg·K)),
Isobaric coefficient of thermal expansion αp (SI unit: 1/K, default 2.75 × 10-4 1/K),
Ratio of specific heats γ (dimensionless, default 1.01), and
Thermal conductivity k (SI unit: W/(m·K), default 0.61 W/(m·K)). The thermal conductivity is always assumed to be a scalar for the purpose of computing the acoustophoretic radiation force.
Advanced Settings
Select the Use piecewise polynomial recovery on field check box to smooth the radiation pressure using piecewise polynomial recovery. This can give a much more accurate representation of the radiation pressure because it uses information on adjacent mesh elements to reconstruct the field. If a coarse mesh is used to compute the field then this option can be especially useful.
Note on Backward Compatibility
In previous versions, the Acoustophoretic Radiation Force was called the Acoustophoretic Force, and it was possible to specify the particle bulk modulus Kp (SI unit: Pa) or compressibility βp (SI unit: 1/Pa) directly, instead of the speed of sound. However, now the particle speed of sound is always required, and these other inputs are not shown.
If a model built in a previous version is opened in the current version, then the speed of sound will be defined as follows to preserve backward compatibility, where <name> is the unique identifier for the physics interface (probably fpt):
If the old expression for the particle bulk modulus was <expr>, the new expression for the compressional speed of sound is sqrt(<expr>/<name>.rhop).
If the old expression for the particle compressibility was <expr>, the new expression for the compressional speed of sound is sqrt(1/(<expr>*<name>.rhop)).
Studies and Solvers and Frequency Domain in the COMSOL Multiphysics Reference Manual
Acoustophoretic Radiation Force in the theory section.