Parameter determination and validation for a mechanistic model of the enzymatic saccharification of cellulose-Iß.

Cost-effective production of fuels and chemicals from lignocellulosic biomass often involves enzymatic saccharification, which has been the subject of intense research and development. Recently, a mechanistic model for the enzymatic saccharification of cellulose has been developed that accounts for distribution of cellulose chain lengths, the accessibility of insoluble cellulose to enzymes, and the distinct modes of action of the component cellulases [Griggs et al. (2012) Biotechnol. Bioeng., 109(3):665-675; Griggs et al. (2012) Biotechnol. Bioeng., 109(3):676-685]. However, determining appropriate values for the adsorption, inhibition, and rate parameters required further experimental investigation. In this work, we performed several sets of experiments to aid in parameter estimation and to quantitatively validate the model. Cellulosic materials differing in degrees of polymerization and crystallinity (α-cellulose-Iß and highly crystalline cellulose-Iß ) were digested by component enzymes (EGI/CBHI/ßG) and by mixtures of these enzymes. Based on information from the literature and the results from these experiments, a single set of model parameters was determined, and the model simulation results using this set of parameters were compared with the experimental data of total glucan conversion, chain-length distribution, and crystallinity. Model simulations show significant agreement with the experimentally derived glucan conversion and chain-length distribution curves and provide interesting insights into multiple complex and interacting physico-chemical phenomena involved in enzymatic hydrolysis, including enzyme synergism, substrate accessibility, cellulose chain length distribution and crystallinity, and inhibition of cellulases by soluble sugars.

A residual-potential boundary for time-dependent, infinite-domain problems in computational acoustics.

A theoretically exact computational boundary is introduced, that is based on modal residual potentials for the spherical geometry. The boundary produces a set of first-order, uncoupled ordinary differential equations for nodal boundary responses, and a set of uncoupled time-stepping equations for modal boundary responses. The two sets are coupled through nodal-modal transformation based on the orthogonal surface functions for the spherical boundary. Numerical results generated with the boundary are presented for a step-wave-excited, elastic, spherical shell submerged in an infinite acoustic medium. Extension of the method to other separable geometries for partial differential equations defined in unbounded domains is mentioned.