# American Institute of Mathematical Sciences

## Modeling, approximation, and time optimal temperature control for binder removal from ceramics

 1 Department of Mathematics, University of Missouri, Columbia, MO 65211, USA 2 Department of Mechanical and Aerospace Engineering, and Department of Chemical Engineering, Columbia, MO 65211, USA 3 Department of Chemical Engineering, University of Missouri, Columbia, MO 65211, USA

* Corresponding author: Stephen J. Lombardo

Received  September 2020 Revised  November 2020 Published  February 2021

The process of binder removal from green ceramic components-a reaction-gas transport problem in porous media-has been analyzed with a number of mathematical techniques: 1) non-dimensionalization of the governing decomposition-reaction ordinary differential equation (ODE) and of the reaction gas-permeability partial differential equation (PDE); 2) development of a pseudo steady state approximation (PSSA) for the PDE, including error analysis via $L^2$ norm and singular perturbation methods; 3) derivation and analysis of a discrete model approximation; and 4) development of a time optimal control strategy to minimize processing time with temperature and pressure constraints. Theoretical analyses indicate the conditions under which the PSSA and discrete models are viable approximations. Numerical results indicate that under a range of conditions corresponding to practical binder burnout conditions, utilization of the optimal temperature protocol leads to shorter cycle times as compared to typical industrial practice.

Citation: Carmen Chicone, Stephen J. Lombardo, David G. Retzloff. Modeling, approximation, and time optimal temperature control for binder removal from ceramics. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021034
##### References:

show all references

##### References:
Top panel: Centerline pressure versus temperature for linear temperature rise test case $\beta = 5.6 \times 10^{-4}$ K/s (that is $T(t) = \beta t+300$) and green body length $L = 0.01$ m using numerical experiments with PDE model (5), PSSA (46), and three-station model (64). Bottom three panels: Relative centerline PSSA error for green body lengths (from top to bottom) $0.001$ m, $0.01$ m and $0.1$ m
Numerically approximated graphs of $L^2$ norms of difference between $\epsilon \rho$ and corresponding PSSA versus temperature in physical units ([$\rho \sqrt{L/2}$]) = mol/m$^{5/2}$ for the linear temperature rise (57) test case where $\beta = 5.6\times 10^{-4}$ and green body lengths $L = 0.001, 0.01, 0.1$ m (depicted in order from bottom to top). The $L^2$ error is obtained using finite-difference approximation (MacCormack's method) of the solution of the model equations (5) and numerical quadrature for the $L^2$ norm, and the theoretical bound is computed likewise from the solution of the ODE (55)
Numerically approximated graphs of the absolute value of the quality function and the absolute difference between $u: = \epsilon \rho$ and the corresponding PSSA versus temperature (both computed at the centerline) for the linear temperature rise (57) test case where $\beta = 5.6\times 10^{-4}$ and three green body lengths $L = 0.001, 0.01, 0.1$ m (depicted in order from bottom to top). The solution of the model equations (5) is approximated using MacCormack's method, and the quality function is computed using (44). Note: To benchmark the values, $\rho_0 = 40$ mol/m$^3$
Graphs of numerical approximations of pressure versus temperature for the linear temperature rise test case $\beta = 5.6 \times 10^{-4}$ K/s (that is $T(t) = \beta t+300$) are depicted for a solid cylinder of height $0.03$ cm and radius $0.005$ cm using the PDE model (5) and corresponding PSSA at the geometric center computed using one (PSSA$_1$), two (PSSA$_2$), and ten (PSSA$_{10}$) of its (transport scaled) power series representation (29) with post process conversion to corresponding physical units
Top panel depicts optimal temperature protocol versus time for green body length $0.01$ m, maximum temperature rise rate of 10 K per minute and maximum gas pressure $200\, 000$ Pa for PDE model, PSSA, and three-station model. Bottom panel shows pressure versus temperature during the optimal temperature protocol for the three models
Top panel depicts for the three-station model an approximation of the time-optimal temperature protocol for a green body of length $0.01$ m for maximum temperature rise rates $0.015$, $0.2$, $0.5$, $10.0$ K per minute and maximum allowed gas pressure $200\,000$ Pa. Bottom panel depicts corresponding pressure versus temperature
Top panel depicts a comparison of the time optimal temperature protocols for the solid cylinder of height 3 cm and radius 0.5 cm as computed using the PDE model and ten terms of the series expansion representation (29) of its corresponding PSSA for the case of maximum heating rate 10 K per minute and maximum allowed gas pressure $150\,000$ Pa, where the maximum is computed at the geometric center of the cylinder. Bottom panel depicts corresponding geometric center gas pressure versus temperature
Model Parameters
 Name Value Description $T_0$ 300.0 K initial temperature $P_0$ $1.0\times 10^5$ Pa ambient gas pressure $R$ 8.314 J/mol-K ideal gas constant $M$ $4.4\times 10^{-2}$ kg/mol molecular weight of gas $E$ $2.22\times 10^5$ J/mol activation energy $A$ $1.67 \times 10^{16}$ /s pre-exponential factor $\mu$ $2.5\times 10^{-5}$ Pa-s viscosity of gas $S$ $6\times 10^6$/m specific surface $k$ $5$ tortuosity $\epsilon_{20}$ 0.4 initial binder volume fraction $\epsilon_3$ 0.5 ceramic volume fraction $\rho_2$ $1.0\times 10^3$ kg/m$^3$ polymer density $L$ 0.001–0.1 m green body length
 Name Value Description $T_0$ 300.0 K initial temperature $P_0$ $1.0\times 10^5$ Pa ambient gas pressure $R$ 8.314 J/mol-K ideal gas constant $M$ $4.4\times 10^{-2}$ kg/mol molecular weight of gas $E$ $2.22\times 10^5$ J/mol activation energy $A$ $1.67 \times 10^{16}$ /s pre-exponential factor $\mu$ $2.5\times 10^{-5}$ Pa-s viscosity of gas $S$ $6\times 10^6$/m specific surface $k$ $5$ tortuosity $\epsilon_{20}$ 0.4 initial binder volume fraction $\epsilon_3$ 0.5 ceramic volume fraction $\rho_2$ $1.0\times 10^3$ kg/m$^3$ polymer density $L$ 0.001–0.1 m green body length
 [1] Luke Finlay, Vladimir Gaitsgory, Ivan Lebedev. Linear programming solutions of periodic optimization problems: approximation of the optimal control. Journal of Industrial & Management Optimization, 2007, 3 (2) : 399-413. doi: 10.3934/jimo.2007.3.399 [2] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [3] Guirong Jiang, Qishao Lu. The dynamics of a Prey-Predator model with impulsive state feedback control. Discrete & Continuous Dynamical Systems - B, 2006, 6 (6) : 1301-1320. doi: 10.3934/dcdsb.2006.6.1301 [4] A. K. Misra, Anupama Sharma, Jia Li. A mathematical model for control of vector borne diseases through media campaigns. Discrete & Continuous Dynamical Systems - B, 2013, 18 (7) : 1909-1927. doi: 10.3934/dcdsb.2013.18.1909 [5] Diana Keller. Optimal control of a linear stochastic Schrödinger equation. Conference Publications, 2013, 2013 (special) : 437-446. doi: 10.3934/proc.2013.2013.437 [6] Tomáš Roubíček. An energy-conserving time-discretisation scheme for poroelastic media with phase-field fracture emitting waves and heat. Discrete & Continuous Dynamical Systems - S, 2017, 10 (4) : 867-893. doi: 10.3934/dcdss.2017044 [7] Shanjian Tang, Fu Zhang. Path-dependent optimal stochastic control and viscosity solution of associated Bellman equations. Discrete & Continuous Dynamical Systems - A, 2015, 35 (11) : 5521-5553. doi: 10.3934/dcds.2015.35.5521 [8] Matthias Erbar, Jan Maas. Gradient flow structures for discrete porous medium equations. Discrete & Continuous Dynamical Systems - A, 2014, 34 (4) : 1355-1374. doi: 10.3934/dcds.2014.34.1355 [9] Guillaume Bal, Wenjia Jing. Homogenization and corrector theory for linear transport in random media. Discrete & Continuous Dynamical Systems - A, 2010, 28 (4) : 1311-1343. doi: 10.3934/dcds.2010.28.1311 [10] Zhihua Zhang, Naoki Saito. PHLST with adaptive tiling and its application to antarctic remote sensing image approximation. Inverse Problems & Imaging, 2014, 8 (1) : 321-337. doi: 10.3934/ipi.2014.8.321 [11] Hyeong-Ohk Bae, Hyoungsuk So, Yeonghun Youn. Interior regularity to the steady incompressible shear thinning fluids with non-Standard growth. Networks & Heterogeneous Media, 2018, 13 (3) : 479-491. doi: 10.3934/nhm.2018021 [12] Yves Dumont, Frederic Chiroleu. Vector control for the Chikungunya disease. Mathematical Biosciences & Engineering, 2010, 7 (2) : 313-345. doi: 10.3934/mbe.2010.7.313 [13] Yunfei Lv, Rong Yuan, Yuan He. Wavefronts of a stage structured model with state--dependent delay. Discrete & Continuous Dynamical Systems - A, 2015, 35 (10) : 4931-4954. doi: 10.3934/dcds.2015.35.4931 [14] Daoyuan Fang, Ting Zhang. Compressible Navier-Stokes equations with vacuum state in one dimension. Communications on Pure & Applied Analysis, 2004, 3 (4) : 675-694. doi: 10.3934/cpaa.2004.3.675 [15] Y. Latushkin, B. Layton. The optimal gap condition for invariant manifolds. Discrete & Continuous Dynamical Systems - A, 1999, 5 (2) : 233-268. doi: 10.3934/dcds.1999.5.233 [16] Martin Bohner, Sabrina Streipert. Optimal harvesting policy for the Beverton--Holt model. Mathematical Biosciences & Engineering, 2016, 13 (4) : 673-695. doi: 10.3934/mbe.2016014 [17] Xingchun Wang, Yongjin Wang. Variance-optimal hedging for target volatility options. Journal of Industrial & Management Optimization, 2014, 10 (1) : 207-218. doi: 10.3934/jimo.2014.10.207 [18] J. Frédéric Bonnans, Justina Gianatti, Francisco J. Silva. On the convergence of the Sakawa-Shindo algorithm in stochastic control. Mathematical Control & Related Fields, 2016, 6 (3) : 391-406. doi: 10.3934/mcrf.2016008 [19] Alberto Bressan, Ke Han, Franco Rampazzo. On the control of non holonomic systems by active constraints. Discrete & Continuous Dynamical Systems - A, 2013, 33 (8) : 3329-3353. doi: 10.3934/dcds.2013.33.3329 [20] Cécile Carrère, Grégoire Nadin. Influence of mutations in phenotypically-structured populations in time periodic environment. Discrete & Continuous Dynamical Systems - B, 2020, 25 (9) : 3609-3630. doi: 10.3934/dcdsb.2020075

2019 Impact Factor: 1.27

## Tools

Article outline

Figures and Tables