Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Library of Advanced Materials for Engineering (LAMÉ) 5.20

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Library of Advanced Materials for Engineering (LAMÉ) 5.18

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Library of Advanced Materials for Engineering (LAMÉ) 5.16

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Library of Advanced Materials for Engineering (LAMÉ) 5.22

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Library of Advanced Materials for Engineering (LAMÉ) (V.5.24)

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options an

36 MATERIALS SCIENCE↗
DOE OSTI2026

Library of Advanced Materials for Engineering (LAMÉ) 5.28

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Library of Advanced Materials for Engineering (LAMÉ) 5.30

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimal operation of solid-oxide electrolysis cells considering long-term chemical degradation

Optimizing the performance of solid oxide electrolysis cells (SOECs) for long-term hydrogen (H 2 ) production at high temperatures is crucial, as prolonged operation leads to efficiency losses and shorter cell lifespans due to chemical degradation. Here, in this work, we adopt a quasi-steady state approach for dynamic optimization over extended operational periods to address the disparity in timescales between cell operation and degradation. Integrating a 2-D non-isothermal SOEC model with balance-of-plant (BOP) equipment, we explore three optimization objectives: minimizing terminal degradation, maximizing integral efficiency, and minimizing the levelized cost of H 2 (LCOH). Our dynamic optimization algorithm reduces LCOH by 9.5% and 16% compared to strategies focusing solely on terminal degradation and integral efficiency, respectively. For electricity prices of 0.03 $\$$/mWh and 0.3 $\$$ mWh optimal replacement schedules range from 5 to 2 years, depending on the operational mode. Furthermore, a flexible operational mode yields additional improvements in LCOH over traditional galvanostatic and potentiostatic modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The Effect of Cholesterol on Membrane Binding and Self-Assembly of Collagen Fibrils

Collagen is a skeleton of native extracellular matrix (ECM) that is known to provide mechanical and structural stability. In an attempt to develop a new connective cellular model with the surrounding ECM without further experimental complications, such as the reconstitution of ECM receptors, we designed the experiments and discovered that the fibrillogenesis of membrane-bound collagen is not spontaneous as it is in the form of free collagen in bulk solution. Furthermore, the confocal microscopic results suggest that cholesterol is a crucial component that facilitates the fibril formation on the membrane surface. In situ X-ray and neutron reflectivity on Langmuir monolayer and solid-supported lipid bilayer models, respectively, reveal two features of cholesterol effects on the collagen fibril formation. Mainly, cholesterol increases the lateral lipid headgroup separation on the membrane surface, which promotes the association degree of collagen monomers. It also enhances the elastic modulus of the membrane to impede membrane filtration by the collagen assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Comments on the SURFplus reactive burn model

The SURFplus reactive burn model is intended for heterogeneous solid HE that produce an excess amount of carbon. It utilizes a fast hotspot reaction (SURF model) plus a second slow reaction for the energy release from carbon clustering. The fast reaction dominates shock initiation, though the fast rate parameters are slightly affected by the late energy release from the slow reaction. The slow reaction has a significant affect on propagating detonation waves. Propagating detonation waves are characterized by the curvature effect; i.e., detonation speed $D_n$ as function of the local front curvature κ. The $D_n(κ)$ curves have a qualitative different shape. This results from the sensitivity of the slope of the $D_n(κ)$ curve to the reaction-zone width up to the sonic point, and the sonic point moving as κ increases from the end of the slow reaction to the neighborhood of the end of the fast reaction. For PBX 9502, the steep slope for small curvature is determined by the slow reaction, while the lower slope for larger curvature is determined by the fast reaction rate in the high detonation pressure regime. The remainder of this paper is organized as follows. In section 2 we review the assumptions of the Shaw and Johnson carbon cluster model used by the SURFplus model for the slow reaction. Section 3 describes in detail the mixture model for a partly burned HE. The form of the slow rate used by the SURFplus model is described in section 4. We conclude with a brief discussion of applying a burn model with 2 rates to either a metalized or a composite PBX.

42 ENGINEERING↗
DOE OSTI2021

Predicting phase transformation kinetics during metal additive manufacturing using non-isothermal Johnson-Mehl-Avrami models: Application to Inconel 718 and Ti-6Al-4V

Here, a computational model was developed to predict solid-state phase transformation kinetics within mechanical parts during metal additive manufacturing processes. This model is a modified version of the Johnson-Mehl-Avrami model for non-isothermal phase transformations that can be applied to various material systems undergoing solid-state phase transformations. Using the thermal history of an additive manufacturing fabricated mechanical part, along with the necessary thermodynamic data and kinetic information as inputs, the model outputs the history of phase fraction evolution during the build process. The model was applied to an Inconel 718 part built by powder bed fusion and a Ti-6Al-4V part built by directed energy deposition. Microstructure characterization and mechanical testing were performed for the validation of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Deriving the Infrared Complex Refractive Indices of Organic Powders for Optical Modeling: Comparison of Methods

Identification of solids via infrared reflection spectroscopy requires a spectral library of all solids likely to be encountered. A confounding factor in populating such a spectral library is that the reflectance spectra of solids vary with their form, including particle size, film thickness, and substrate. To reduce the efforts of experimentally constructing such a library, an alternate strategy is to use the wavelength-dependent optical constants, n and k, of a solid to calculate a series of reflectance spectra corresponding to each scenario or morphology. Because most n/k measurements are best performed on mm-sized crystals, however, the challenge of determining the optical constants increases when a solid is only readily available as a powder, as is often the case. Some organic solids, such as caffeine, are both unavailable in large crystals and difficult to press into pellets. In this study, the infrared optical constants, or complex refractive indices, of caffeine were determined using three different methods: single-angle reflectance, infrared spectroscopic ellipsometry, and quantitative absorbance measurements of KBr pellets. The n and k values derived through each method were used to model the hyperspectral imaging reflectance spectrum of a caffeine film on a steel planchet. Over 1,110 – 870 cm-1, the single-angle reflectance-derived n and k had the best correlation with the experimental spectrum. These results suggest different organic solids may require different methods to determine the most accurate infrared complex refractive indices for synthetic spectral libraries.

FTIR, hyperspectral imaging, infrared, infrared (I↗
DOE OSTI2020

High Temperature Erosion Modeling in Particle Based CSP Systems

1. Introduction. Wear and erosion damage of materials from solid particle and surface interactions is a major issue in various industries. Although more common in mining and oil and gas production, erosion is becoming a critical issue in renewable technologies as well such as particle based concentrated solar thermal power (CSP) systems. In particle based CSP systems, solid particles are used to absorb solar energy and as thermal storage. However, these particles may cause significant amount of wear to system components while through the system at high temperatures. This damage can be costly, and therefore, requires a greater understanding of solid particle erosion in CSP systems. Throughout the years, models and tools have been developed to predict and control erosion in industries such as oil and gas production. However, these erosion models and erosion prediction tools have been mainly developed based on erosion data for much higher velocities and lower temperatures, i.e. the operating conditions for which the existing models have been built are not pertinent to those expected in CSP system. It is known that erosion depends on many factors such as material properties, erodent particle properties, and particle impact speed and angle. However, the effect of temperature on erosion is not vastly investigated. In this work, a temperature-based correlation is introduced that will be used to modify the existing erosion models to predict erosion rates at velocities and temperatures relevant to Gen3 CSP systems. The existing models are first validated against erosion experiments run at low temperatures and low velocities. Subsequently, a temperature correction term is developed that can extend the existing models to high temperatures, based on the available experimental data showing the effect of temperature on impact erosion. 2. Erosion Models. Continuing sequence of impacts from solid particles on surfaces would result in loss of material due to mechanical interaction between solid surface and particles. Erosion equations have been developed to predict erosion under different conditions and erosion mechanism, including cutting and deformation erosion. These models are mainly either mechanistic, empirical, and semi-mechanistic models, when the latter combine the theory of the erosion mechanism and particle motion with the available experimental data. One of the first empirical correlation was introduced by American Petroleum Institute (API) Recommended Practice (RP) 14E [1]. This correlation was very conservative, and several improved empirical and semi-mechanistic equations were developed in the years following. More recently, the models introduced by Erosion/Corrosion Research Center (E/CRC) at the University of Tulsa and Arabnejad. et al. [2] are frequently used in the literature and industry, as they account for many parameters affecting erosion including particle impact speed and angle, material density and hardness, and particle size and shape. However, both these models were developed based on data at relatively high velocities and also do not take into account any mechanical changes in material associated with thermal cycling at high temperatures. 3. Erosion Prediction at High Temperatures. In this work, the accuracy of Arabnejad et al. model and E/CRC model is investigated to predict impact erosion at conditions relevant to Gen3 CSP systems. Two significant deviations expected in Gen3 systems compared to operating conditions used to build these models are particle impact velocities and system temperature. As a first step, the performance of these models was validated at low velocities and temperatures. Impact erosion experiments were conducted on SS316 coupons using HSP 40/70 ceramic particles. Assuming particle rate of 1 kg/s/m, 0.0254 m of particle curtain thickness, and particle velocity of 1.5 m/s in the system, an overall erosion of 0.59 mm/year was calculated. The experimental results were subsequently compared to those from computational simulations and erosion of 0.47 mm/year was obtained using the Arabnejad et al. model. The results indicated that the models work well at low velocities and low temperature conditions. To account for temperature effect, a mathematical correlation was developed using data published by DUCOM [3] for Inconel 600 eroded by alumina particles at high velocities. The correlation was applied to both Arabnejad et. al model and the E/CRC model. The prediction results from these modified models were within 20% of this experimental data. Predictions of of erosion by the modified E/CRC model at three different temperatures are made. Similar to the calculation at low temperature, a typical CSP system with 1 kg/s/m of particle flow rate per unit length of the particle curtain and a curtain thickness of 0.0254 m is assumed. The annual thickness loss calculations were run for different particle impact velocities, assuming a uniform particle impact area equal to the cross-sectional area of the curtain (particle-particle interactions and dispersion of particles are not considered in obtaining the results). It is observed that erosion increases exponentially as temperature increases. Furthermore, it is also noted that, erosion changes non-linearly with impact velocity. Currently, experiments are also being conducted to measure erosion of SS316 at 800 ºC with HSP 40/70 particles at low impact velocities. We expect to use the results from high temperature testing to further improve the temperature correlation function. Similar models are also being developed for abrasion erosion resulting from particle sliding along the surfaces as well as attrition from particle to particle and particle to surface interactions. References. [1] Institute, A. P. (1991). API Recommended Practice for Design and Installation of Offshore Production Platform Piping System, API RP 14E. [2] Arabnejad, H., Mansouri, A., Shirazi, S. A., and McLaury, B. S. (2015a). Development of mechanistic erosion equation for solid particles. Wear, 332–333, 1044–1050. http://doi.org/10.1016/j.wear.2015.01.031. [3] https://ducom.com/high-temperature-erosion-evaluating-sample-wear/

14 SOLAR ENERGY↗
DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗