lattice_spheres#

Generates a cubic packing of spheres in a specified lattice arrangement.

import matplotlib.pyplot as plt

import porespy as ps
import inspect

b = inspect.signature(ps.generators.lattice_spheres)
print(b)
(shape: List, r: int = 5, spacing: int = None, offset: int = None, smooth: bool = True, lattice: Literal['sc', 'tri', 'fcc', 'bcc'] = 'sc')

radius#

Controls the size of the spheres.

fig, ax = plt.subplots(1, 2, figsize=[8, 4])

shape = [200, 200]
r = 10
im1 = ps.generators.lattice_spheres(shape=shape, r=r)
ax[0].imshow(im1, interpolation='none')
ax[0].axis(False)
ax[0].set_title(f'radius = {r}')

r = 20
im2 = ps.generators.lattice_spheres(shape=shape, r=r)
ax[1].imshow(im2, interpolation='none')
ax[1].axis(False)
ax[1].set_title(f'radius = {r}');
../../../_images/4f0aeff5f6ba87b4b70d06366024c0a6181789c9c715e452e10d24e86f258a58.png

spacing#

The center-to-center spacing between the spheres. If this value is less than the sphere diamter then the spheres will overlap.

fig, ax = plt.subplots(1, 2, figsize=[8, 4])

s = 35
im1 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s)
ax[0].imshow(im1, interpolation='none')
ax[0].axis(False)
ax[0].set_title(f'spacing = {s}')


s = 50
im2 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s)
ax[1].imshow(im2, interpolation='none')
ax[1].axis(False)
ax[1].set_title(f'spacing = {s}');
../../../_images/6c0be68a3fbf3fee5d6057baa7ab714fae8239e3fefa66536089a2511b87339b.png

offset#

Controls how far away from the edge the first sphere is located. The default is the sphere radius but it can be more or less depending on the desired effect:

fig, ax = plt.subplots(1, 2, figsize=[8, 4])

o = 0
im1 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o)
ax[0].imshow(im1, interpolation='none')
ax[0].axis(False)
ax[0].set_title(f'offset = {o}')

o = 25
im2 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o)
ax[1].imshow(im2, interpolation='none')
ax[1].axis(False)
ax[1].set_title(f'offset = {o}');
../../../_images/cf2e85af4a13bbc79bf07e48ae6ee24dc4cc624cd12576fdf9fc064d85a52f93.png

lattice#

Controls the arrange of spheres. In 2D the options are simple cubic (‘sc’) and triangular (‘tri’). Note that the offset and spacing apply to the outer spheres when lattice='tri'.

fig, ax = plt.subplots(1, 2, figsize=[8, 4])

L = 'sc'
r = 10
im1 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o, lattice=L)

ax[0].imshow(im1, interpolation='none')
ax[0].axis(False)
ax[0].set_title(f'lattice = {L}')

L = 'tri'
im2 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o, lattice=L)

ax[1].imshow(im2, interpolation='none')
ax[1].axis(False)
ax[1].set_title(f'lattice = {L}');
../../../_images/45163c170a027bb0a0513483939ff87e16ce7cca5b9adba255a37a8202dd9d3b.png

In 3D the options are simple cubic (‘sc’), face centered cubic (‘fcc’), and body centered cubic (‘bcc’). It’s more difficult to visualize in 3D but PoreSpy has a basic function called “show_3D” that works if the image is small:

fig, ax = plt.subplots(1, 3, figsize=[9, 3])

r = 10
s = 25
shape = [100, 100, 100]
L = 'sc'
im1 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o, lattice=L)
ax[0].imshow(ps.visualization.show_3D(im1))
ax[0].axis(False)
ax[0].set_title(f'lattice = {L}')

L = 'fcc'
im2 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o, lattice=L)
ax[1].imshow(ps.visualization.show_3D(im2))
ax[1].axis(False)
ax[1].set_title(f'lattice = {L}')

L = 'bcc'
im3 = ps.generators.lattice_spheres(shape=shape, r=r, spacing=s, offset=o, lattice=L)
ax[2].imshow(ps.visualization.show_3D(im3))
ax[2].axis(False)
ax[2].set_title(f'lattice = {L}');
../../../_images/a4610418679b51ccafb7af6591f79b3ebe3facadf38e544edf84575cff856b6b.png