Results for “goodness of fit”

DOE OSTI2024

Is Clean Hydrogen Production a Good Fit for Questa? (Final Economic Impact Results) [Slides]

The Village of Questa, New Mexico is aiming to become a regional clean energy hub with robust and diverse employment opportunities for the local community supported by the energy sector and by other businesses inspired or attracted by abundant clean energy, outdoor recreation, and cultural opportunities. A coalition of stakeholders in the Village of Questa, comprising the Village, Kit Carson Electric Cooperative (KCEC), Questa Economic Development Fund, and Chevron, is exploring options to develop hydrogen production facilities as an opportunity to create jobs, provide reliable clean energy, and utilize former mine resources. Questa is home to a molybdenum mine owned by Chevron that closed in 2014. Several residents in Questa and surrounding communities lost their jobs when the mine closed and transitioned from active operations into environmental remediation. Although remediation efforts have been ongoing since 2014 and are expected to continue for at least 16 more years, the number of jobs with Chevron is much smaller now than it was before the closure. Between available workforce, brownfield land, and water rights formerly supporting mine operations but now in a transition period, there are considerable local resources that could be directed toward clean energy generation. Questa's electricity supply is already 100% solar during daylight hours thanks to Kit Carson Electric Cooperative's (KCEC's) strategic decision-making and partnering over the last decade. Now, Questa, KCEC, and Chevron are exploring the potential costs and benefits of siting an electrolytic hydrogen production facility and additional solar photovoltaic (PV) capacity in Questa to further advance the region's clean energy economy. In this report, we estimated the potential economic impacts (i.e., jobs, value added, gross output, tax revenue) of constructing and operating a combined hydrogen (32 MW polymer electrolyte membrane electrolizer + 7.5 MW fuel cell) and solar facility (22.5 MW) in the Village of Questa, as well as the resulting economic spillovers to Taos County and the state of New Mexico. We employ an input-output model that leverages IMPLAN's economic data for the region complemented by construction and operating expenses estimated by NREL and feedback from the local coalition to evaluate the direct, indirect and induced effects of the project construction (transient impacts) and operation (more permanent impacts). Based on the area's average trade profile, feedback from the coalition and current market conditions, these projects are expected to support 487 full-time equivalent jobs during construction, generating $\$24$ million in income for those workers and $\$82$ million in local economic activity in the state. Of those jobs, 106 are expected to be construction sector jobs. These investments are also estimated to add $\$36.5$ million to New Mexico's gross state product (GSP). In the Village of Questa, we estimate 16 jobs will be supported in construction and transportation industries, generating $\$0.9$ million in earnings. In Taos County, the construction phase is expected to support 285 jobs primarily in construction and professional services, while manufacturing jobs dominate the results for the Rest of New Mexico. The Village is also estimated to receive $\$0.9$ million in tax revenue from the construction phase alone. Once in operation, the project continues to impact the state and Questa. Around 20 jobs (full-time equivalent for each year of operation) are supported across New Mexico, with approximately 11 directly employed in Questa by both facilities. The total annual local economic activity supported by ongoing operations is just over $\$1.3$ million/yr, generating $\$1.6$ million/yr in additional income in the state. Annual operations are estimated to add $\$2.1$ million to the state's GSP. The Village is expected to receive around $\$43,000$/yr in tax revenue. Impacts vary significantly depending on which businesses are supplying materials, equipment and services, and where construction workers reside. Choosing local suppliers will most benefit Questa and the New Mexico economy, adding up to 500 jobs during construction and 13 long-term jobs. Local and state governments may consider ways to incentivize local businesses in order to maximize economic benefits.

08 HYDROGEN↗

DOE OSTI2024

Is Clean Hydrogen Production a Good Fit for Questa? Final Economic Impact Results

The Village of Questa, New Mexico is aiming to become a regional clean energy hub with robust and diverse employment opportunities for the local community supported by the energy sector and by other businesses inspired or attracted by abundant clean energy, outdoor recreation, and cultural opportunities. A coalition of stakeholders in the Village of Questa, comprising the Village, Kit Carson Electric Cooperative (KCEC), Questa Economic Development Fund, and Chevron, is exploring options to develop hydrogen production facilities as an opportunity to create jobs, provide reliable clean energy, and utilize former mine resources. Questa is home to a molybdenum mine owned by Chevron that closed in 2014. Several residents in Questa and surrounding communities lost their jobs when the mine closed and transitioned from active operations into environmental remediation. Although remediation efforts have been ongoing since 2014 and are expected to continue for at least 16 more years, the number of jobs with Chevron is much smaller now than it was before the closure. Between available workforce, brownfield land, and water rights formerly supporting mine operations but now in a transition period, there are considerable local resources that could be directed toward clean energy generation. Questa's electricity supply is already 100% solar during daylight hours thanks to Kit Carson Electric Cooperative's (KCEC's) strategic decision-making and partnering over the last decade. Now, Questa, KCEC, and Chevron are exploring the potential costs and benefits of siting an electrolytic hydrogen production facility and additional solar photovoltaic (PV) capacity in Questa to further advance the region's clean energy economy. In this report, we estimated the potential economic impacts (i.e., jobs, value added, gross output, tax revenue) of constructing and operating a combined hydrogen (32 MW polymer electrolyte membrane electrolizer + 7.5 MW fuel cell) and solar facility (22.5 MW) in the Village of Questa, as well as the resulting economic spillovers to Taos County and the state of New Mexico. We employ an input-output model that leverages IMPLAN's economic data for the region complemented by construction and operating expenses estimated by NREL and feedback from the local coalition to evaluate the direct, indirect and induced effects of the project construction (transient impacts) and operation (more permanent impacts). Based on the area's average trade profile, feedback from the coalition and current market conditions, these projects are expected to support 487 full-time equivalent jobs during construction, generating $\$24$ million in income for those workers and $\$82$ million in local economic activity in the state. Of those jobs, 106 are expected to be construction sector jobs. These investments are also estimated to add $\$36.5$ million to New Mexico's gross state product (GSP). In the Village of Questa, we estimate 16 jobs will be supported in construction and transportation industries, generating $\$0.9$ million in earnings. In Taos County, the construction phase is expected to support 285 jobs primarily in construction and professional services, while manufacturing jobs dominate the results for the Rest of New Mexico. The Village is also estimated to receive $\$0.9$ million in tax revenue from the construction phase alone. Once in operation, the project continues to impact the state and Questa. Around 20 jobs (full-time equivalent for each year of operation) are supported across New Mexico, with approximately 11 directly employed in Questa by both facilities. The total annual local economic activity supported by ongoing operations is just over $\$1.3$ million/yr, generating $\$1.6$ million/yr in additional income in the state. Annual operations are estimated to add $\$2.1$ million to the state's GSP. The Village is expected to receive around $\$43,000$/yr in tax revenue. Impacts vary significantly depending on which businesses are supplying materials, equipment and services, and where construction workers reside. Choosing local suppliers will most benefit Questa and the New Mexico economy, adding up to 500 jobs during construction and 13 long-term jobs. Local and state governments may consider ways to incentivize local businesses in order to maximize economic benefits.

08 HYDROGEN↗

DOE OSTI2022

Is Clean Hydrogen Production a Good Fit for Questa? Project Introduction

This presentation introduces the Questa Communities Local Energy Action Program (C-LEAP) project to local stakeholders, provides general information about green hydrogen and presents some preliminary information about how a green hydrogen production facility could be developed in Questa.

Communities LEAP↗

DOE OSTI2023

Is Clean Hydrogen Production a Good Fit for Questa? Intermediate Feasibility Study Results

NREL conducted a feasibility study for a clean hydrogen production facility that would be located in Questa, New Mexico. This slide deck presents the history of the coalition of stakeholders that received NREL technical assistance through the Department of Energy's Communities Leading Energy Action Program. It describes how hydrogen is generated and its end uses. It describes two conceptual applications explored in the feasibility analysis, a grid-tied power generation and long-term storage application and a heavy-duty vehicle fueling application. The presentation details the configuration of these to applications in the Questa pilot facility, including system components, required inputs (water, electricity, and land for pv generation), and outputs in kilograms of hydrogen, megawatt hours of power delivered, and hours of vehicle hours of operation. The presentation addresses key concerns raised by groups engaged through a public outreach process, including safety, health, and environmental issues; water use; visual impacts; and costs and alternatives.

heavy-duty vehicle fueling↗

DOE OSTI2023

Monte Carlo goodness-of-fit tests for degree corrected and related stochastic blockmodels

Abstract We construct Bayesian and frequentist finite-sample goodness-of-fit tests for three different variants of the stochastic blockmodel for network data. Since all of the stochastic blockmodel variants are log-linear in form when block assignments are known, the tests for the latent block model versions combine a block membership estimator with the algebraic statistics machinery for testing goodness-of-fit in log-linear models. We describe Markov bases and marginal polytopes of the variants of the stochastic blockmodel and discuss how both facilitate the development of goodness-of-fit tests and understanding of model behaviour. The general testing methodology developed here extends to any finite mixture of log-linear models on discrete data, and as such is the first application of the algebraic statistics machinery for latent-variable models.

Karwa, Vishesh↗

DOE OSTI2024

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Practical guide to understanding goodness-of-fit metrics used in chemical state modeling of x-ray photoelectron spectroscopy data by synthetic line shapes using nylon as an example

Chemical state analysis of a sample surface through fitting bell-shaped curves to x-ray photoelectron spectroscopic polymer data is reviewed using nylon to introduce and discuss aspects of data analysis. Different strategies for modeling chemistry in nylon spectra are presented and in so doing, a case is made to include in published science the design logic and implementation in terms of line shapes and optimization parameter constraints between components in a peak model. Imperfections in line shape relative to the true shape for photoemission lines, when compensated for using constraints to optimization parameters, are shown to provide chemical state information about a sample that justify, for peak models constructed with these limitations, metrics for goodness-of-fit different from those expected for pulse-counted data.

Materials Science↗

DOE OSTI2022

An improved representation of the relationship between photosynthesis and stomatal conductance leads to more stable estimation of conductance parameters and improves the goodness-of-fit across diverse data sets

Stomata play a central role in surface-atmosphere exchange by controlling the flux of water and CO 2 between the leaf and the atmosphere. Representation of stomatal conductance (g sw ) is therefore an essential component of models that seek to simulate water and CO 2 exchange in plants and ecosystems. For given environmental conditions at the leaf surface (CO 2 concentration and vapor pressure deficit or relative humidity), models typically assume a linear relationship between g sw and photosynthetic CO 2 assimilation (A). However, measurement of leaf-level g sw response curves to changes in A are rare, particularly in the tropics, resulting in only limited data to evaluate this key assumption. Here, we measured the response of g sw and A to irradiance in six tropical species at different leaf phenological stages. We showed that the relationship between g sw and A was not linear, challenging the key assumption upon which optimality theory is based-that the marginal cost of water gain is constant. Our data showed that increasing A resulted in a small increase in g sw at low irradiance, but a much larger increase at high irradiance. We reformulated the popular Unified Stomatal Optimization (USO) model to account for this phenomenon and to enable consistent estimation of the key conductance parameters g 0 and g 1 . Our modification of the USO model improved the goodness-of-fit and reduced bias, enabling robust estimation of conductance parameters at any irradiance. In addition, our modification revealed previously undetectable relationships between the stomatal slope parameter g 1 and other leaf traits. We also observed nonlinear behavior between A and g sw in independent datasets that included data collected from attached and detached leaves, and from plants grown at elevated CO 2 concentration. We propose that this empirical modification of the USO model can improve the measurement of g sw parameters and the estimation of plant and ecosystem-scale water and CO 2 fluxes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Hidden conformal symmetry from the lattice

We analyze newly expanded and refined data from lattice studies of an SU(3) gauge theory with eight Dirac fermions in the fundamental representation. We focus on the light composite states emerging from these studies, consisting of a set of pseudoscalars and a single light scalar. We first consider the view that this theory is just outside the conformal window. In this case, the pseudoscalars arise from spontaneous breaking of chiral symmetry. Identifying the scalar in this case as an approximate dilaton, we fit the lattice data to a dilaton effective field theory, finding that it yields a good fit even at lowest order. For comparison, we then consider the possibility that the theory is inside the conformal window. The fermion mass provides a deformation, triggering confinement. We employ simple scaling laws to fit the lattice data, and find that it is of lesser quality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Return of the templates: Revisiting the Galactic Center excess with multimessenger observations

The Galactic center excess (GCE) remains one of the most intriguing discoveries from the Fermi Large Area Telescope (LAT) observations. We revisit the characteristics of the GCE by first producing a new set of high-resolution galactic diffuse gamma-ray emission templates. This diffuse emission, which accounts for the bulk of the observed gamma rays, is ultimately due to cosmic-ray interactions with the interstellar medium. Using recent high-precision cosmic-ray observations, in addition to the continuing Fermi-LAT observations and observations from lower energy photons, we constrain the properties of the galactic diffuse emission. We describe a large set of diffuse gamma-ray emission templates which account for a very wide range of initial assumptions on the physical conditions in the inner galaxy. The broad properties of the GCE that we find in this work are qualitatively unchanged despite the introduction of this new set of templates, though its quantitative features appear mildly different than those obtained in previous analyses. In particular, we find a high-energy tail at higher significance than previously reported. This tail is very prominent in the northern hemisphere, and less so in the southern hemisphere. This strongly affects one prominent interpretation of the excess: known millisecond pulsars are incapable of producing this high-energy emission, even in the relatively softer southern hemisphere, and are therefore disfavored as the sole explanation of the GCE. The annihilation of dark matter particles of mass 40 - 7 + 10 GeV (95% CL) to b quarks with a cross-section of ⟨ σ A v ⟩ = 1. 4 - 0.3 + 0.6 × 10 - 26 cm 3 s − 1 provides a good fit to the excess especially in the relatively cleaner southern sky. Dark matter of the same mass range annihilating to b quarks or heavier dark matter particles annihilating to heavier Standard Model bosons can combine with millisecond pulsars to provide a good fit to the southern hemisphere emission as well, as can a broken power-law spectrum which would be related to recent cosmic-ray burst activity. As part of this paper, we make publicly available all of our templates and the data covariance matrix we have generated to account for systematic uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Inelastic dark matter scattering off Thallium cannot save DAMA

We study the compatibility of the observed DAMA modulation signal with inelastic scattering of dark matter (DM) off of the 0.1% Thallium (Tl) dopant in DAMA. In this work we test whether there exist regions of parameter space where the Tl interpretation gives a good fit to the most recent data from DAMA, and whether these regions are compatible with the latest constraints from other direct detection experiments. Previously, Chang et al. in 2010, had proposed the Tl interpretation of the DAMA data, and more recently (in 2019) the DAMA/LIBRA collaboration found regions in parameter space of Tl inelastic scattering that differ by more than 10σ from a no modulation hypothesis. We have expanded upon their work by testing whether the regions of parameter space where inelastic DM-Tl scattering gives a good fit to the most recent DAMA data survive the constraints placed by the lack of a DM signal in XENON1T and CRESST-II. In addition, we have tested how these regions change with the main sources of uncertainty: the Tl quenching factor, which has never been measured directly, and the astrophysical uncertainties in the DM distribution. We conclude that inelastic DM scattering off Tl cannot explain the DAMA data in light of null results from other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A new class of multivariate goodness of fit tests for multivariate normal mixtures

Not provided.

Mathematics↗

DOE OSTI2023

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Empirical Estimation of Shortest Route Length along U.S. Interstate Highways Based on Great Circle Distance

In this study, 98 regression models were specified for easily estimating shortest distances based on great circle distances along the U.S. interstate highways nationwide and for each of the continental 48 states. This allows transportation professionals to quickly generate distance, or even distance matrix, without expending significant efforts on complicated shortest path calculations. For simple usage by all professionals, all models are present in the simple linear regression form. Only one explanatory variable, the great circle distance, is considered to calculate the route distance. For each geographic scope (i.e., the national or one of the states), two different models were considered, with and without the intercept. Based on the adjusted R-squared, it was observed that models without intercepts generally have better fitness. Additionally, all these models generally have good fitness with the linear regression relationship between the great circle distance and route distance. At the state level, significant variations in the slope coefficients between the state-level models were also observed. Furthermore, a preliminary analysis of the effect of highway density on this variation was conducted.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗

DOE OSTI2020

Using SPT polarization, Planck 2015, and non-CMB data to constrain tilted spatially-flat and untilted nonflat ΛCDM, XCDM, and ΦCDM dark energy inflation cosmologies

We use six tilted spatially-flat and untilted non-flat dark energy cosmological models in analyses of South Pole Telescope polarization (SPTpol) cosmic microwave background (CMB) data, alone and in combination with Planck 2015 CMB data and non-CMB data, namely, the Pantheon Type Ia supernovae apparent magnitudes, a collection of baryon acoustic oscillation data points, Hubble parameter measurements, and growth rates. Although the cosmological models that best-fit the Planck CMB and non-CMB data do not provide good fits to the SPTpol data, with the χ2's exceeding the expected value, given the uncertainties, in each model the cosmological parameter constraints from the SPTpol data and from the Planck CMB and non-CMB data are largely mutually consistent. Furthermore, when the smaller angular scale SPTpol data are used jointly with either the Planck data alone or with the Planck CMB and the non-CMB data to constrain untilted non-flat models, spatially-closed models remain favored over their corresponding flat limits. When used in conjunction with Planck data, non-CMB data (baryon acoustic oscillation measurements in particular, from six experiments) have significantly more constraining power than the SPTpol data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3

Chang, Hilary↗

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↗