Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “yield point”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Three-dimensional phenotyping of peach tree-crown architecture utilizing terrestrial laser scanning

Tree training systems for temperate fruit have been developed throughout history by pomologists to improve light interception, fruit yield, and fruit quality. These training systems direct crown and branch growth to specific configurations. Quantifying crown architecture could aid the selection of trees that require less pruning or that naturally excel in specific growing/training system conditions. Regarding peaches [Prunus persica (L.) Batsch], access tools such as branching indices have been developed to characterize tree-crown architecture. However, the required branching data (BD) to develop these indices are difficult to collect. Traditionally, BD have been collected manually, but this process is tedious, time-consuming, and prone to human error. These barriers can be circumnavigated by utilizing terrestrial laser scanning (TLS) to obtain a digital twin of the real tree. TLS generates three-dimensional (3D) point clouds of the tree crown, wherein every point contains 3D coordinates (x, y, z). To facilitate the use of these tools for peach, we selected 16 young peach trees scanned in 2021 and 2022. These 16 trees were then modeled and quantified using the open-source software TreeQSM. As a result, “in silico” branching and biometric data for the young peach trees were calculated to demonstrate the capabilities of TLS phenotyping of peach tree-crown architecture. The comparison and analysis of field measurements (in situ) and in silico BD, biometric data, and quantitative structural model branch uncertainty data were utilized to determine the reconstructive model’s reliability as a source substitute for field measurements. Mean average deviation when comparing young tree (YT) height was approx. 5.93%, with crown volume was approx. 13.26% across both 2021 and 2022. All point clouds of the YTs in 2022 showed residuals lower than 12 mm to cylinders fitted to all branches, and mean surface coverage greater than 40% for both the trunk and primary branching orders.

09 BIOMASS FUELS↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The DESI DR1 Peculiar Velocity Survey: Global Zero-point and H 0 Constraints

The Dark Energy Spectroscopic Instrument (DESI) in its first Data Release (DR1) already provides more than 100,000 galaxies with relative distance measurements. The primary purpose of this paper is to perform the calibration of the zero-point for the DESI Fundamental Plane and Tully–Fisher relations, which allows us to measure the Hubble constant, H 0 . This sample has a lower statistical uncertainty than any previously used to measure H 0 , and we investigate the systematic uncertainties in absolute calibration that could limit the accuracy of that measurement. We improve upon the DESI Early Data Release Fundamental Plane H 0 measurement by (a) using a group catalog to increase the number of calibrator galaxies and (b) investigating alternative calibrators in the nearby Universe. Our baseline measurement calibrates to the SH0ES/Pantheon+ type Ia supernovae, and finds H 0 = 73.7 ± 0.06 (stat.) ± 1.1 (syst.) km s −1 Mpc −1 . Calibrating to surface brightness fluctuation distances yields a similar H 0 . We explore measurements using other calibrators, but these are currently less precise since the overlap with DESI peculiar velocity tracers is much smaller. In future data releases with an even larger peculiar velocity sample, we plan to calibrate directly to Cepheids and the tip of the red giant branch, which will enable the uncertainty to decrease towards a percent-level measurement of H 0 . This will provide an alternative to supernovae as the Hubble flow sample for H 0 measurements.

Carr, Anthony [Korea Astronomy and Space Science I↗
DOE OSTI2022

Plastic Waste Product Captures Carbon Dioxide in Nanometer Pores

We discuss, plastic waste (PW) and increasing atmospheric carbon dioxide (CO 2 ) levels are among the top environmental concerns presently facing humankind. With an ambitious 2050 zero-CO 2 emissions goal, there is a demand for economical CO 2 capture routes. Here we show that the thermal treatment of PW in the presence of potassium acetate yields an effective carbon sorbent with pores width of 0.7–1.4 nm for CO 2 capture. The PW to carbon sorbent process works with single or mixed streams of polyolefin plastics. The CO 2 capacity of the sorbent at 25 °C is 17.0 ± 1.1 wt % (3.80 ± 0.25 mmol g –1 ) at 1 bar and 5.0 ± 0.6 wt % (1.13 ± 0.13 mmol g –1 ) at 0.15 bar, and it regenerates upon reaching 75 ± 5 °C. The CO 2 capture cost from flue gas via this technology is estimated to be <$21 ton –1 CO 2 , much lower than competing CO 2 capture technologies. Hence, this PW-derived carbon material should find utility in the capture of CO 2 from point sources of high CO 2 emissions while providing a use for otherwise deleterious PW.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Rutherford-in-Copper-Channel Conductor for the MARCO Solenoidal Detector Magnet

MARCO will be a superconducting solenoid for a new particle physics detector of the upcoming Electron Ion Collider (EIC) at Brookhaven National Laboratory (NY, USA). The design field at the interaction point is 2.0 T with a nominal current of about 4 kA at 4.5 K. For this magnet, an aluminumstabilized conductor was considered at the very preliminary design phase. After that, since the lack of manufacturers able to deal with aluminum extrusion, copper has been chosen as stabilizer. Thanks to a dedicated design effort, the copper stabilizer and accompanying material choices provided acceptable hadronic interaction length. The conductor is based on a NbTi Rutherford cable that is soldered in a U-shaped copper profile. Here, this paper goes through the definition of the conductor crosssection according to the project requirements and the supplier capabilities. Its characteristics will be detailed and discussed, in particular the RRR and the yield strength of the copper channel needed for protection and mechanics respectively. Finally, a proposal on how to make the joint between two conductors is presented.

Stacchi, Francesco [Commissariat a l'Energie Atomi↗
DOE OSTI2023

Synthesis of Hf 0.75 Ta 0.25 B 2 for self-coating TPS

Ultra-high temperature ceramic materials (UHTCs) are important for designing high-performance aerospace vehicles that can withstand repeated exposures to high temperatures. UHTCs mixed with silicides often have well-controlled oxidation due to the formation of protective silicates, which create a regenerative outer protective layer. Borides are known to have high melting points, high hardness and reasonably good oxidation resistance. In this study, Hf 0.75 Ta 0.25 B 2 (HTB) is proposed as a potential alternative to YSZ protective coatings and ZrB 2- SiC composites via the formation of Hf 6 Ta 2 O 17 (HTO) passivation layer. The objective of this paper is to explore the parameter space of HTB synthesis via borocarbothermal (BCTR). Effects on particle size, phase purity, and residual oxygen content were analyzed with parameters of atmosphere composition, reactant grain sizes, and differing reaction pathways. The BCTR of HfO 2 and Ta 2 O 5 were analyzed to predict HTB behavior. It was shown that both one-step and two-step reaction routes can yield HTB, but two-step yields a purer product. Nano B 4 C produced finer HTB and facilitated reaction completion. Using a reducing atmosphere also enhanced reaction completion

36 MATERIALS SCIENCE↗
DOE OSTI2022

Climate-Driven Legacies in Simulated Microbial Communities Alter Litter Decomposition Rates

The mechanisms underlying diversity-functioning relationships have been a consistent area of inquiry in biogeochemistry since the 1950s. Though these mechanisms remain unresolved in soil microbiomes, many approaches at varying scales have pointed to the same notion—composition matters. Confronting the methodological challenge arising from the complexity of microbiomes, this study used the model DEMENTpy, a trait-based modeling framework, to explore trait-based drivers of microbiome-dependent litter decomposition. We parameterized DEMENTpy for five sites along a climate gradient in Southern California, United States, and conducted reciprocal transplant simulations analogous to a prior empirical study. The simulations demonstrated climate-dependent legacy effects of microbial communities on plant litter decomposition across the gradient. This result is consistent with the previous empirical study across the same gradient. An analysis of community-level traits further suggests that a 3-way tradeoff among resource acquisition, stress tolerance, and yield strategies influences community assembly. Simulated litter decomposition was predictable with two community traits (indicative of two of the three strategies) plus local environment, regardless of the system state (transient vs. equilibrium). Although more empirical confirmation is still needed, community traits plus local environmental factors (e.g., environment and litter chemistry) may robustly predict litter decomposition across spatial-temporal scales. In conclusion, this study offers a potential trait-based explanation for climate-dependent community effects on litter decomposition with implications for improved understanding of whole-ecosystem functioning across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Projecting the number of new HIV infections to formulate the "Getting to Zero" strategy in Illinois, USA

Objectives: Getting to Zero (GTZ) initiatives focus on expanding use of antiretroviral treatment (ART) and pre-exposure prophylaxis (PrEP) to eliminate new HIV infections. Computational models help inform policies for implementation of ART and PrEP continuums. Such models, however, vary in their design, and may yield inconsistent predictions. Using multiple approaches can help assess the consistency in results obtained from varied modeling frameworks, and can inform optimal implementation strategies. Methods: A study using three different modeling approaches is conducted. Two approaches use statistical time series analysis techniques that incorporate temporal HIV incidence data. A third approach uses stochastic stimulation, conducted using an agent-based network model (ABNM). All three approaches are used to project HIV incidence among a key population, young Black MSM (YBMSM), over the course of the GTZ implementation period (2016-2030). Results: All three approaches suggest that simultaneously increasing PrEP and ART uptake is likely to be more effective than increasing only one, but increasing ART and PrEP by 20% points may not eliminate new HIV infections among YBMSM. The results further suggest that a 20% increase in ART is likely to be more effective than a 20% increase in PrEP. All three methods consistently project that increasing ART and PrEP by 30% simultaneously can help reach GTZ goals. Conclusions: Increasing PrEP and ART uptake by about 30% might be necessary to accomplish GTZ goals. Such scale-up may require addressing psychosocial and structural barriers to engagement in HIV and PrEP care continuums. ABNMs and other flexible modeling approaches can be extended to examine specific interventions that address these barriers and may provide important data to guide the successful intervention implementation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Attaining the promise of plant gene editing at scale

Crop improvement relies heavily on genetic variation that arises spontaneously through mutation. Modern breeding methods are very adept at combining this genetic variation in ways that achieve remarkable improvements in plant performance. Novel traits have also been created through mutation breeding and transgenesis. The advent of gene editing, however, marks a turning point: With gene editing, synthetic variation will increasingly supplement and, in some cases, supplant the genetic variation that occurs naturally. We are still in the very early stages of realizing the opportunity provided by plant gene editing. At present, typically only one or a few genes are targeted for mutation at a time, and most mutations result in loss of gene function. New technological developments, however, promise to make it possible to perform gene editing at scale. RNA virus vectors, for example, can deliver gene-editing reagents to the germ line through infection and create hundreds to thousands of diverse mutations in the progeny of infected plants. With developmental regulators, edited somatic cells can be induced to form meristems that yield seed-producing shoots, thereby increasing throughput and shrinking timescales for creating edited plants. As these approaches are refined and others developed, they will allow for accelerated breeding, the domestication of orphan crops and the reengineering of metabolism in a more directed manner than has ever previously been possible.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Data compression and covariance matrix inspection: Cosmic shear

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and “removing” 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the massively optimized parameter estimation and data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like Ω m , the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as S 8 . For the case considered—cosmic shear from the first year of data from the Dark Energy Survey—only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter A IA is the most susceptible to inaccuracies in the covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2023

Every-other-layer dipolar excitons in a spin-valley locked superlattice

Monolayer semiconducting transition metal dichalcogenides possess broken inversion symmetry and strong spin-orbit coupling, leading to a unique spin-valley locking effect. In 2H stacked pristine multilayers, spin-valley locking yields an electronic superlattice structure, where alternating layers correspond to barriers and quantum wells depending on the spin-valley indices. Here, in this work, we show that the spin-valley locked superlattice hosts a kind of dipolar exciton with the electron and hole constituents separated in an every-other-layer configuration: that is, either in two even or two odd layers. Such excitons become optically bright via hybridization with intralayer excitons. This effect is also manifested by the presence of multiple anti-crossing patterns in the reflectance spectra, as the dipolar exciton is tuned through the intralayer resonance by an electric field. The reflectance spectra further reveal an excited state orbital of the every-other-layer exciton, pointing to a sizable binding energy in the same order of magnitude as the intralayer exciton. As layer thickness increases, the dipolar exciton can form a one-dimensional Bose–Hubbard chain displaying layer number-dependent fine spectroscopy structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of divertor material on neutral recycling and discharge fueling in DIII-D

We report that experiments with the lower divertor of DIII-D during the Metal Rings Campaign (MRC) show that the fraction F of atomic D in the total recycling flux is material-dependent and varies through the ELM cycle, which may affect divertor fueling. Between ELMs, F C ~ 10% and FW ~ 40%, consistent with expectations if all atomic recycling is due to reflections. During ELMs, F C increases to 50% and F W to 60%. In contrast, the total D recycling coefficient including atoms and molecules R stays close to unity near the strike point where the surface is saturated with D. During ELMs, R can deviate from unity, increasing during high energy ELM-ion deposition (net D release) and decreasing at the end of the ELM which leads to ability of the target to trap the ELM-deposited D. The increase of R > 1 in response to an increase in ion impact energy Ei has been studied with small divertor target samples using Divertor Materials Evaluation System (DiMES). An electrostatic bias was applied to DiMES to change E i by 90 eV. On all studied materials including C, Mo, uncoated and W-coated TZM (>99% Mo, Ti, and Zr alloy), W, and W fuzz, an increase of E i transiently increased the D yield (and R) by ~10%. On C there was also an increase in the molecular D 2 yield, probably due to ion-induced D 2 desorption. Despite the measured increase in F on W compared to C, attached H-mode shots with OSP on W during MRC did not demonstrate a higher pedestal density. About 8% increase in the edge density could be seen only in attached L-mode scenarios. The difference can be explained by higher D trapping in the divertor and lower divertor fueling efficiency in H- versus L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coupling Between Electrons and Charge Density Wave Fluctuation and its Possible Role in Superconductivity

In most charge density wave (CDW) systems of different material classes, ranging from traditional correlated systems in low-dimension to recent topological systems with Kagome lattice, superconductivity emerges when the system is driven toward the quantum critical point (QCP) of CDW via external parameters of doping and pressure. Despite this rather universal trend, the essential hinge between CDW and superconductivity has not been established yet. Here, the evidence of coupling between electron and CDW fluctuation is reported, based on a temperature- and intercalation-dependent kink in the angle-resolved photoemission spectra of 2H-Pd x TaSe 2 . Kinks are observed only when the system is in the CDW phase, regardless of whether a long- or short-range order is established. Notably, the coupling strength is enhanced upon long-range CDW suppression, albeit the coupling energy scale is reduced. Interestingly, the estimation of the superconducting critical temperature by incorporating the observed coupling characteristics into McMillan's equation yields results closely resembling the known values of the superconducting dome. The results thus highlight a compelling possibility that this new coupling mediates Cooper pairs, which provides new insights into the competing relationship not only for CDW but also for other competing orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Microstructure effects on high velocity microparticle impacts of copper

Constitutive models can fail to predict high-rate deformation behavior due to their inability to account for microstructural effects. In part, this is because of a dearth of experimental benchmarking data in the high strain-rate, low pressure regime, since many high-rate experiments also probe a region of strong shockwaves, at which point microstructure effects no longer play a primary role. This work uses laser-induced particle impact testing to quantitatively study high velocity impacts of small, rigid alumina microspheres on flat copper substrates with varying amounts of initial cold work in the weak shock regime, but at very high strain rates up to ~10 7 s –1 . Through paired experiments and numerical simulations, this work shows that the initial microstructure condition can have significant influence on dynamical mechanical properties in this range. Specifically, prior work hardening of the copper substrate leads to increased rebounding of the microparticles (i.e., less plastic dissipation in the impact) as well as smaller craters. Each of these experimental measurables can be converted into a strength measure, i.e., the dynamic yield strength or dynamic hardness, respectively, neither of which is well predicted consistently by existing constitutive laws. The general trend of hardening can be captured by such models by incorporating an existing “pre-strain,” suggesting that future calibration of the materials parameters may yield a good fit over a broader range of conditions. Our results emphasize the importance of reporting the microstructural condition in dynamic studies, as well as the necessity of accounting for these factors when formulating and optimizing constitutive models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of a Genome-Scale Metabolic Model of Clostridium thermocellum and Its Applications for Integration of Multi-Omics Datasets and Computational Strain Design

Solving environmental and social challenges such as climate change requires a shift from our current non-renewable manufacturing model to a sustainable bioeconomy. To lower carbon emissions in the production of fuels and chemicals, plant biomass feedstocks can replace petroleum using microorganisms as biocatalysts. The anaerobic thermophile Clostridium thermocellum is a promising bacterium for bioconversion due to its capability to efficiently degrade lignocellulosic biomass. However, the complex metabolism of C. thermocellum is not fully understood, hindering metabolic engineering to achieve high titers, rates, and yields of targeted molecules. In this study, we developed an updated genome-scale metabolic model of C. thermocellum that accounts for recent metabolic findings, has improved prediction accuracy, and is standard-conformant to ensure easy reproducibility. We illustrated two applications of the developed model. We first formulated a multi-omics integration protocol and used it to understand redox metabolism and potential bottlenecks in biofuel (e.g., ethanol) production in C. thermocellum. Second, we used the metabolic model to design modular cells for efficient production of alcohols and esters with broad applications as flavors, fragrances, solvents, and fuels. The proposed designs not only feature intuitive push-and-pull metabolic engineering strategies, but also present novel manipulations around important central metabolic branch-points. We anticipate the developed genome-scale metabolic model will provide a useful tool for system analysis of C. thermocellum metabolism to fundamentally understand its physiology and guide metabolic engineering strategies to rapidly generate modular production strains for effective biosynthesis of biofuels and biochemicals from lignocellulosic biomass.

59 BASIC BIOLOGICAL SCIENCES↗