US 12,347,532 B2
Boltzmann-based method for simulating CVI densification process of composite material
Aijun Li, Zhejiang (CN); Dan Zhang, Zhejiang (CN); Jingchao Yuan, Zhejiang (CN); and Meihua Shi, Zhejiang (CN)
Assigned to SHAOXING RESEARCH INSTITUTE OF SHANGHAI UNIVERSITY, Zhejiang (CN); and SHANGHAI UNIVERSITY, Shanghai (CN)
Appl. No. 17/291,742
Filed by SHAOXING RESEARCH INSTITUTE OF SHANGHAI UNIVERSITY, Zhejiang (CN); and SHANGHAI UNIVERSITY, Shanghai (CN)
PCT Filed Jun. 4, 2020, PCT No. PCT/CN2020/094447
§ 371(c)(1), (2) Date May 6, 2021,
PCT Pub. No. WO2020/244597, PCT Pub. Date Dec. 10, 2020.
Claims priority of application No. 201910487883.6 (CN), filed on Jun. 5, 2019.
Prior Publication US 2022/0130497 A1, Apr. 28, 2022
Int. Cl. G16C 60/00 (2019.01); G06F 30/10 (2020.01); G06F 30/25 (2020.01); G06F 113/26 (2020.01); G16C 20/10 (2019.01)
CPC G16C 60/00 (2019.02) [G06F 30/10 (2020.01); G06F 30/25 (2020.01); G06F 2113/26 (2020.01)] 3 Claims
 
1. A Boltzmann-based method for simulating a Chemical Vapor Infiltration (CVI) densification process of a composite material, comprising the steps of:
a) modelling on a computer, a three-dimensional model of an actual preform, producing a three-dimensional matrix by scanning pixels one by one, recording a component number of each pixel at the same time, and storing component information of the pixel is in a matrix form in one-to-one correspondence with a spatial position, which is referred to as a component matrix;
b) assigning a phase component: comprising extracting target attributes according to a spatial distribution relationship and attribute relationships recorded by the component matrix, establishing a matrix based on the target attributes, performing normalization processing to obtain a phase matrix, and counting volume occupancy of a certain phase in the spatial position from the phase matrix, wherein the volume occupancy is a ratio of an amount of substances in the phase component to an amount of all accommodated in space;
c) dividing grid comprising re-dividing the phase matrix according to an actual precision requirement, and establishing an independent phase matrix for each of different phase components, wherein all the independent phase matrices are superimposed into a re-divided grid three-dimensional model of the actual preform;
d) assigning material attributes comprising assigning material attributes according to a re-divided grid phase matrix by way of one-to-one multiplication of the independent phase matrix and the established matrix based on the target attributes;
e) setting boundary conditions according to different actual conditions;
f) calculating a gas-phase flow field by using a Lattice Boltzmann Method (LBM), comprising two sub-steps:
1) Taking a virtual time step so that all particles move without restriction on a set of virtual grids; and
2) Then, taking a modified time step, releasing the particles entering the boundary according to the boundary conditions described by each independent phase matrix, and modifying the gas-phase flow field;
g) calculating a chemical reaction by a phase transformation algorithm, comprising scaling the chemical reaction according to an actual reaction relationship after calculation under specified conditions, and adding and subtracting on a basis of each independent phase matrix;
h) circulating the steps f) and g) until the reaction ends;
i) processing results-comprising outputting and counting calculated independent phase matrixes, and comparing each calculated independent phase matrix with an original matrix to obtain relevant information after CVI densification; and j) performing a CVI densification process on the actual preform according to the relevant information after CVI densification, and forming the composite material comprising the actual preform and a matrix.