燃烧仿真.燃烧数值模拟方法：大涡模拟（LES）：高级燃烧仿真技术.pdfVIP

燃烧仿真.燃烧数值模拟方法：大涡模拟（LES）：高级燃烧

仿真技术

1燃烧仿真基础

1.1燃烧理论简介

燃烧是一种复杂的物理化学过程，涉及到燃料与氧化剂的化学反应、热量

的产生与传递、以及流体动力学的相互作用。燃烧理论主要研究燃烧的机理、

燃烧的条件、燃烧过程中的能量转换和物质转化。燃烧可以分为几个阶段：燃

料的蒸发或分解、燃料与氧化剂的混合、化学反应的发生、以及燃烧产物的冷

却和扩散。

1.1.1燃烧的类型

扩散燃烧：燃料和氧化剂在燃烧前是分开的，燃烧发生在它们混

合的界面。

预混燃烧：燃料和氧化剂在燃烧前已经充分混合，燃烧速度由化

学反应速率决定。

层流燃烧：在低雷诺数条件下，燃烧过程是层流的，火焰传播速

度稳定。

湍流燃烧：在高雷诺数条件下，燃烧过程受到湍流的影响，火焰

传播速度和燃烧效率会显著提高。

1.1.2燃烧的化学反应

燃烧的化学反应通常涉及多个步骤，包括链引发、链传递和链终止。例如，

甲烷（CH4）在氧气（O2）中的燃烧可以简化为以下反应：

CH4+2O2-CO2+2H2O

但实际上，这个过程涉及多个中间反应，如自由基的生成和消耗。

1.2数值模拟方法概述

数值模拟是通过计算机程序来解决数学模型的方法，广泛应用于燃烧仿真

中。它能够处理复杂的非线性方程组，提供燃烧过程的详细信息，如温度、压

力、浓度分布等。数值模拟方法主要包括：

有限差分法：将连续的偏微分方程离散化，用差分方程近似求解。

有限体积法：基于控制体积原理，将计算域划分为多个控制体积，

然后在每个控制体积上应用守恒定律。

有限元法：将计算域划分为多个单元，用单元内的插值函数来逼

1

近解。

1.2.1大涡模拟（LES）

大涡模拟（LargeEddySimulation,LES）是一种用于模拟湍流燃烧的高级数

值方法。LES的目标是直接模拟大尺度涡旋，而小尺度涡旋则通过亚网格模型

来处理。这种方法能够捕捉到湍流的主要特征，同时减少计算成本。

1.2.1.1LES的数学模型

LES基于Navier-Stokes方程，但对这些方程进行滤波处理，以消除小尺度

涡旋的影响。滤波后的方程称为滤波Navier-Stokes方程。例如，对于速度场u，

其滤波值U定义为：

U(x)=\intu(x)G(x-x)dx

其中，−是滤波函数，它决定了滤波的尺度。

1.2.1.2亚网格模型

在LES中，小尺度涡旋的影响通过亚网格模型来描述。常见的亚网格模型

包括：

Smagorinsky模型：基于湍流的耗散率来估计亚网格应力。

WALE模型：考虑了涡旋的拉伸和压缩效应。

动态模型：通过计算过程中的信息来动态调整模型参数。

1.2.2示例：使用Python进行LES的简单模拟

虽然LES的完整模拟通常需要专业的软件和高性能计算资源，下面是一个

使用Python进行LES概念性模拟的简化示例，仅用于说明LES的基本思想：

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格和时间步长

nx,ny=128,128

dx,dy=2.0/nx,2.0/ny

dt=0.01

nt=100

#初始化速度场

u=np.zeros((nx,ny))

v=np.zeros((nx,ny))

#定义滤波函数

deffilter_field(field):

returnnp.convolve(field,np.ones((3,3))/9,mode=same)

2

#LES主循环

forninrange(nt):

u_filtered=filter_field(u)

v_filtered=filter_field(v)

#更新速度场（简化