Acoustofluidics 2022: Minimal Example Plotting Scattering Fields

This is the minimal example on how to plot the scattered pressure field around a polystyrene particle suspended in water and subjected to an acoustic standing wave from the abstract book from the Acoustofluidics 2022 conference.

We import osaft and initialise the model from Yosioka & Kawasima.

from matplotlib import pyplot as plt

import osaft

yosioka = osaft.yosioka1955.ScatteringField(
    f=1e6,
    R_0=1e-6,
    rho_s=1020,
    c_s=2350,
    rho_f=997,
    c_f=1498,
    p_0=1e5,
    wave_type=osaft.WaveType.STANDING,
    position=osaft.pi / 4,
)

The model is passed to the plotting class and the pressure field is plotted. The option incident=False makes sure that only the scattered field is plotted but not the incident field.

plot = osaft.FluidScatteringPlot(yosioka, r_max=4e-6)

plot.plot_pressure(incident=False)

plt.show()

Estimated memory usage: 29 MB