2D超临界二氧化碳压裂裂缝扩展带裂缝扩展动图展示的python代码

以下是一个用Python绘制2D超临界二氧化碳压裂裂缝扩展带的动画展示的代码示例：

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation

# 创建图形对象
fig = plt.figure()
ax = fig.add_subplot(111)

# 定义网格大小和时间步长
nx = 100  # X方向网格数
ny = 50  # Y方向网格数
nt = 100  # 时间步数

# 初始化裂缝扩展带
fracture = np.zeros((ny, nx))  # 初始时所有网格值为0

# 更新函数，用于更新裂缝扩展带状态
def update(frame):
    global fracture

    # 随机选择一个网格进行扩展
    x = np.random.randint(0, nx)
    y = np.random.randint(0, ny)

    # 扩展裂缝
    fracture[y, x] += 1

    # 绘制裂缝扩展带
    ax.clear()
    ax.imshow(fracture, cmap='hot', interpolation='nearest')

# 创建动画对象
ani = animation.FuncAnimation(fig, update, frames=nt, interval=200)

# 显示动画
plt.show()

这段代码使用了`numpy`和`matplotlib`库来生成并展示2D超临界二氧化碳压裂裂缝扩展带的动画。在每个时间步长，随机选择一个网格进行扩展，并使用热图展示裂缝扩展带的状态。动画通过`animation.FuncAnimation`函数创建，并使用`plt.show()`显示出来。