Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient analysis”

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

CNRS Molten Salt Reactor Benchmark Analysis Using Griffin-Pronghorn Coupled Multi-physics Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time, since the fuel salt is circulating through the whole primary loop. This movement of fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and the corresponding redistribution of the DNPs in the active region of the core. To capture this effect, neutronics and thermal-hydraulics (T/H) codes need to be able to handle the movement of DNPs and their decay [1, 2]. Idaho National Laboratory (INL) is actively working on developing the neutronics and T/H codes to model the MSRs with flowing fuel. As one of the efforts, the neutronics code Griffin [3] and the T/H code Pronghorn [4] that are built upon the MOOSE framework [5] have been extended to handle the flowing fuel involving the drift of DNPs, and this capability is being tested and verified with MSR benchmark problems including the CNRS benchmark problem [6]. Additionally, a multi-physics analysis of the molten salt reactor experiment (MSRE) is being conducted to verify and validate the codes against available experimental data [7]. This work presents the analysis of the CNRS benchmark problem using the Griffin-Pronghorn coupled multi-physics code system of INL. Both steady-state and transient problems are analyzed, and the results are compared with those of other participants presented in [6]. The following section provides a description of the benchmark, along with input parameters and observables for each step. Then, the results for each step are presented in comparison with other participants’ results, followed by summary and conclusion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Kinetic and X-ray Absorption Spectroscopic Analysis of Catalytic Redox Cycles over Highly Uniform Polymetal Oxo Clusters

Metal–organic framework materials (MOFs) offer an opportunity for investigating catalytic properties of polymetal oxo clusters that are highly well defined and uniform in nature, in contrast to other classes of catalysts that may exhibit a propensity toward active site heterogeneity. We report herein a kinetic and X-ray absorption spectroscopy (XAS) analysis of the two-electron oxidation of CO over divalent metal sites in MIL-100(M = Fe, Cr) (MIL = Materials of Institut Lavoisier) materials carrying μ 3 -oxo bridged trimers, and connect observations about the kinetic relevance of redox steps to density functional theory (DFT) predictions published previously. The high degree of uniformity evident from in situ titration measurements leads to a congruence in mechanistic inferences made from steady-state catalytic, transient stoichiometric, isotopic exchange, and isotopic tracer data that all point to a sequential mechanism comprised of separate oxidation and reduction half-cycle steps conjoined by an active oxygen intermediate. In situ XAS data reinforce mechanistic conclusions derived from kinetic analysis, and suggest that the active oxygen intermediate may be more appropriately characterized as an iron-oxyl (Fe 3+ –O – ) rather than an iron-oxo (Fe 4+ =O 2– ) species. The Cr analogue of MIL-100 exhibits contrasting rate features that can be rationalized using an identical sequence of steps as MIL-100(Fe), but with a highly dissimilar set of kinetic parameters that can also be validated using transient stoichiometric experiments. The larger coverages of active oxygen intermediates on MIL-100(Cr) are consistent with predictions from prior DFT studies that suggest more stable and less reactive active oxygen species for metals with lower d-electron counts, and point to metal identity as a lever for precisely controlling the kinetic relevance of oxidation and reduction half-cycles in catalytic redox sequences over polymetal oxo clusters. Importantly, the results presented point to the utility of developing broadly applicable structure–catalytic property relationships over MOF nodes specifically, and highly uniform catalysts more generally.

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

Integrated Steady-State System Package for Nuclear Thermal Propulsion Analysis Using Multi-Dimensional Thermal Hydraulics and Dimensionless Turbopump Treatment

Nuclear thermal propulsion is an evolving technology that can be utilized for long-distance space travel. This technology yields the advantage of a high thrust and specific impulse, but requires an examination of the potential design adjustments necessary to enhance its feasibility. The development of nuclear thermal propulsion requires a comprehensive understanding of the system-level behavior during transient and steady-state operation. This paper extends our previous research by including the proper handling of turbomachinery with multi-channel thermal hydraulic simulations only for steady-state solutions. The system-level approach presented here enables the treatment of the turbopump components through non-dimensional analysis that eliminates the assumption of constant efficiencies. All the other components within the system (e.g., reflector and core) can be discretized to multiple channels and layers, in which the full thermal hydraulic solution is established. The approach chosen here enables the realistic modeling of the propellant flow within the expander cycle by capturing the pressure losses, mass flow rate splits, and enthalpy gain for various operational conditions. The verification of the package is completed through point comparisons of previous investigations into similar system designs. Furthermore, sensitivity studies are used to benchmark the capabilities of the package and investigate solution variations due to the perturbation of operational conditions and regimes. The sensitivity studies performed here are important to capture variation in flow characteristics (e.g., temperature, pressure, mass flow rates) for different design objectives such as the thrust and specific impulse. This work demonstrates that system-level simulations lacking multi-channel capability and proper turbomachinery treatment may yield higher uncertainties in understanding the engine’s response and characteristics to changing various requirements. This is extremely important when screening the design space of such propulsion systems and when transient simulations are required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Cause identification of electromagnetic transient events using spatiotemporal feature learning

This paper presents a spatiotemporal feature learning method for cause identification of electromagnetic transient events in power grids. The proposed method is formulated based on the availability of time-synchronized high-frequency measurements and using the convolutional neural network as the spatiotemporal feature representation along with softmax function for the classification. Despite the existing threshold-based, or energy-based events analysis methods, such as support vector machine autoencoder, and tapered multi-layer perceptron neural network, the proposed feature learning is carried out with respect to both time and space. The effectiveness of the proposed feature learning and the subsequent cause identification is validated through the Electromagnetic Transients Program (EMTP) simulation of different events such as line energization, capacitor bank energization, lightning, fault, and high-impedance fault in the IEEE 30-bus, and the real-time digital simulation of the Western System Coordinating Council (WSCC) 9-bus system.

42 ENGINEERING↗
DOE OSTI2021

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

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

ECAR-6585 MARVEL Core Support Structure Analysis: Project #33526

Idaho National Laboratory (INL) is developing a small microreactor to produce electrical power utilizing a small nuclear core and stirling engines for power output under the Microreactor Application Research Validation and Evaluation (MARVEL) program. The MARVEL Reactor (Figure 1) will be installed in the Transient Reactor Test (TREAT) Facility. Walsh Engineering Services, PC (WES) has contracted with the INL to detail the final design. There are six main systems in the MARVEL, one of which is the Fuel and Core System (FCS). The FCS contains a subsystem–the Core Support Structures Subsystem (CSS)–that supports and orients each of the 36 nuclear fuel pins in the MARVEL Reactor System. The FCS (apart from the fixed Beryllium Oxide radial reflectors) is shown in light green and is depicted in more detail in Figure 1 and Figure 2. The main subject of this analysis is the lower grid plate, which is the principal item in the CSS which supports the fuel, see Figure 3. Additionally, the upper grid plate thermal expansion of the upper grid plate is presented in Appendix A.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Critical Evaluation of a Novel Analysis Technique for Assessment of Printed Circuit Heat Exchangers in High-Temperature Nuclear Service

Application of printed circuit heat exchangers (PCHEs) to very high-temperature reactors (VHTRs) requires mechanical performance assessment methodologies. The PCHE morphology consists of thousands of millimeter-scale channels, for enhanced thermal efficiency, enclosed in a meter-scale PCHE core. PCHE geometry under thermomechanical creep-fatigue transients results in multi-axial interactions between its different segments, such as channeled core, walls, and headers. These global-level interactions influence the local channel-level responses. Hence, developing a PCHE performance assessment methodology, following the ASME Code, Section III, Division 5 provisions, is a critical gap to be filled. There is no analysis or design methodology available in ASME Code to assess a PCHE for its global and local level performances under high temperature and pressure loadings. This paper critically evaluates a recently proposed two-step analysis technique to estimate global interactions and local channel-level responses of PCHEs. In this novel analysis technique, the channeled PCHE core is replaced with orthotropic solid blocks of representative stiffness properties for the global thermomechanical analysis. Here, subsequent channel scale submodel analysis with detailed channel geometry, loading, and elastic-perfectly plastic (EPP) material model estimates the local responses for PCHE performance assessment. This paper critically evaluates this novel technique for its effectiveness in PCHE performance assessment. Finite element (FE) models imitating various analysis issues are developed, and FE analysis results are scrutinized. An important outcome of this study is the validation of the novel two-step PCHE analysis technique for application to the performance assessment of PCHEs in VHTRs.

42 ENGINEERING↗
DOE OSTI2024

Selenium interaction with iron minerals: Quantitative comparison of sorption and coprecipitation impacts on mobility

Given the significance of selenium (Se) as a micronutrient, the radioactive nature of some of its isotopes, and its affinity to iron (Fe) minerals, extensive research has been conducted on the sorption mechanisms between Se and these minerals. Here, in this study, we employ sorption data sourced from the L-SCIE database and coprecipitation data from available literature to achieve the following objectives: i) establish coherence between adsorption and coprecipitation processes, ii) quantitatively evaluate the importance of these processes in nuclear waste repository science, and iii) propose a forward-looking approach for integrating coprecipitation into reactive transport models. Our findings indicate that a correlation between Se adsorption and coprecipitation can be established using the λ formalism. The comparable log(λ Se(IV) /λ Se(VI) ) ratios derived from adsorption and coprecipitation experiments suggest that these processes can be quantitatively compared and evaluated using our numerical approach. Across all iron oxide phases examined, coprecipitation leads to significantly greater immobilization of Se compared to adsorption. Specifically, for hydrous ferric oxide, hematite, and goethite, coprecipitation is predicted to result in 100–1000 times more Se immobilization compared to adsorption, irrespective of the Se oxidation state (Se(IV) or Se(VI)); notably stronger immobilization potential via coprecipitation was observed for magnetite. The modeling approach and quantitative analysis presented herein clearly highlight the importance of including coprecipitation processes when simulating Se (and other elements) transport, particularly under conditions where mineral compositions are transient or evolving with time. Neglecting coprecipitation in models is likely to lead to significant overestimates of migration.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

The Impact of Carotenoid Energy Levels on the Exciton Dynamics and Singlet–Triplet Annihilation in Isolated Bacterial Light-Harvesting 2 Complexes

The light-harvesting 2 (LH2) complex of purple phototrophic bacteria plays a critical role in absorbing solar energy and distributing the excitation energy. Exciton dynamics within LH2 complexes are controlled by the structural arrangement and energy levels of the bacteriochlorophyll (BChl) and carotenoid (Car) pigments. However, there is still debate over the competing light-harvesting versus energy-dissipation pathways. In this work, we compared five variants of the LH2 complex from genetically modified strains of Rhodobacter sphaeroides, all containing the same BChls but different Cars with increasing conjugation: zeta-carotene (N = 7; LH2 Zeta ), neurosporene (N = 9; LH2 Neu ), spheroidene (N = 10; LH2 Spher ), lycopene (N = 11; LH2 Lyco ), and spirilloxanthin (N = 13; LH2 Spir ). Absorption measurements confirmed that the Car excited-state energy decreased with increasing conjugation. Similarly, fluorescence spectra showed that the B850 BChl emission peak had an increasing red shift from LH2 Zeta →(LH2 Neu /LH2 Spher )→LH2 Lyco →LH2 Spir . In contrast, time-resolved fluorescence and ultrafast transient absorption (fs-TA) revealed similar excited-state lifetimes (∼1 ns) for all complexes except LH2 Spir (∼0.7 ns). From fs-TA analysis, an additional ∼7 ps nonradiative dissipation step from B850 BChl was observed for LH2 Zeta . Further, singlet− singlet and singlet−triplet annihilation studies showed a ∼50% average fluorescence lifetime reduction in LH2Zeta at high laser power and high repetition rate, compared to ∼10−15% reductions in LH2 Neu /LH2 Spher /LH2 Lyco and minimal lifetime change in LH2 Spir . In LH2 Zeta , the fastest decay component (<50 ps) became prominent at high repetition rates, consistent with strong singlet−triplet annihilation. Nanosecond TA measurements revealed long-lived (>40 μs) BChl triplet states in LH2 Zeta and signs of damage caused by singlet oxygen, whereas other LH2s showed faster triplet quenching (∼18 ns) by Cars. These findings highlight a key design principle of LH2 complexes: the Car triplet energy must be significantly lower than the BChl triplet energy to efficiently quench BChl triplets that otherwise act as potent “trap states,” causing exciton annihilation in laser-based experiments or photodamage in native membranes.

Absorption↗
DOE OSTI2022

A New Type of Weak Grid IBR Oscillations

This letter presents the discovery of a new type oscillations associated with inverter-based resources (IBRs) in weak grids. While typical weak grid oscillations become worse when the grid strength decreases and the power exporting level increases, this new type of oscillations appear when the grid becomes weaker and the power exporting level decreases. Electromagnetic transient (EMT) simulation results are presented, along with eigenvalue loci for a varying power. The notable signatures of such oscillations are also remarked. In further analysis, a feedback system consisting of two blocks is derived with one block associated with the phase-locked-loop (PLL) and the other associated with the rest of the system. Here, this feedback system quantitatively demonstrates why insufficient damping in the PLL may be an issue for stability. Open-loop analysis shows that a larger phase lag in angle tracking in the low-frequency range due to high power exporting level at the weak grid condition is beneficial to the open-loop poles introduced by the PLL. And this is the reason why such oscillations may appear when power is ramping down.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Impact of Inverter Based Resource Negative Sequence Current Injection on Transmission System Protection

This report documents the results of analysis performed to investigate the impact of inverter-based resource (IBR) response to unbalanced faults on transmission system protective relay dependability and security. Electromagnetic transient (EMT) simulations were performed to simulate IBR response to these faults using existing manufacturer-developed EMT models for four separate IBRs. The study team was composed of IBR manufacturers, relay manufacturers, transmission providers, reliability coordinators and industry consultants with experience in EMT simulation and system protection. The results indicate that under certain conditions, IBR response can result in protective relay misoperations if current protection practices, which were developed based on conventional power sources, are not adapted to the characteristics of IBRs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Impact of Wide-Area Oscillation Damping Control using Measurement-Driven Approach on System Separation - Saudi Grid Case Study

In any interconnected power grid, low-frequency oscillations is a major problem that can limit the power transfer capability and deteriorate power system security due to potential low-damped or even undamped oscillations. Synchronized measurements provided by PMUs enable the design and development of wide-area oscillation damping controllers (WADC) based on measurement-driven model to overcome the limitations of traditional controllers. This work focuses on the impact of a W ADC based on a measurement-driven approach on system separation for the Saudi Electricity Company (SEC) grid. A grid model is provided by SEC for this study. Modal analysis is performed to identify the oscillation mode of interest for which the controller is designed. Damping the dominant oscillation mode helps slightly to improve transient stability of the system. Considering bus frequency at different locations in SEC's system as the candidate observation signals, the optimal observation signal is chosen for the W ADC by using the FFT method. The system transfer function model is constructed by utilizing probing measurements. Then, the W ADC parameters are calculated based on the identified model. The performance of the designed W ADC is tested under different contingencies.

Altarjami, Ibrahim↗
DOE OSTI2021

Excited-State Bond Contraction and Charge Migration in a Platinum Dimer Complex Characterized by X-ray and Optical Transient Absorption Spectroscopy

Interactions between metal centers in dimeric transition metal complexes (TMCs) play important roles in their excited-state energetics and pathways, and, thus, affect their photophysical properties relevant to their applications, e.g., photoluminescent materials and photocatalysis. In this study, we report electronic and nuclear structural dynamics studies of two photoexcited pyrazolate-bridged [Pt(ppy)(µ-R 2 pz)] 2 -type Pt(II) dimers (ppy = 2-phenylpyridine, µ-R 2 pz = 3,5-substituted pyrazolate): [Pt(ppy)(µ-H 2 pz)] 2 (1) and [Pt(NDI-ppy)(µ-Ph 2 pz)] 2 (2, NDI = 1,4,5,8-naphthalenediimide), both of which have distinct ground state Pt-Pt distances. X-ray transient absorption spectroscopy at the Pt L III -edge revealed a new d-orbital vacancy due to the one-electron oxidation of the Pt centers in 1 and 2. However, while a transient Pt-Pt contraction was observed in 2, such an effect was completely absent in 1, demonstrating how the excited states of these complexes are determined by the overlap of the Pt (d z 2 ) orbitals, which is tuned by steric bulk of the pyrazolate R-groups in the 3- and 5-positions. In tandem with analysis of the Pt-Pt distance structural parameter, we observed photoinduced charge separation in 2 featuring a covalently linked NDI acceptor on the ppy ligand. The formation and subsequent decay of the NDI radical anion absorption signals were detected upon photoexcitation using optical transient absorption spectroscopy. The NDI radical anion decayed on the same timescale, hundreds of picoseconds, as that of the d-orbital vacancy signal of the oxidized Pt-Pt core observed in the XTA measurements. The data indicated ultrafast formation of the charge-separated state and subsequent charge recombination to the original Pt(II-II) species.

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

Polarization-resolved core exciton dynamics in LiF using attosecond transient absorption spectroscopy

The ability to control absorption by modifying the polarization of light presents an exciting opportunity to experimentally determine the orbital alignment of absorption features. Here, attosecond extreme ultraviolet (XUV) transient absorption spectroscopy is used to investigate the polarization dependence of core exciton dynamics in LiF thin films at the Li + K edge. XUV pulses excite electrons from the Li 1s core level into the conduction band, allowing for the formation of a p-orbital-like core exciton, aligned along the XUV light polarization axis. A sub-5 fs near-infrared (NIR) probe pulse then arrives at variable time delays, perturbing the XUV-excited states and allowing the coherence decay of the core exciton to be mapped. The coherence lifetimes are found to be ≈ 2.4 ± 0.4 fs, which is attributed to a phonon-mediated dephasing mechanism as in previous core exciton studies. The differential absorption features are also shown to be sensitive to the relative polarization of the XUV and NIR fields. The parallel NIR probe induces couplings between the initial XUV-excited p-like bright exciton and s-like dark excitons. When crossed pump and probe polarizations are used, the coupling between the bright and dark states is no longer dipole-allowed, and the transient absorption signal associated with the coupling is suppressed by approximately 90%. This interpretation is supported by simulations of a few-level model system, as well as analysis of the calculated band structure. Here, the results indicate that laser polarization can serve as a powerful experimental tool for exploring the orbital alignment of core excitonic states in solid-state materials.

Attosecond laser spectroscopy↗
DOE OSTI2022

Searching for ultra-light bosons and constraining black hole spin distributions with stellar tidal disruption events

Stars that pass close to the supermassive black holes located in the center of galaxies can be disrupted by tidal forces, leading to flares that are observed as bright transient events in sky surveys. The rate for these events to occur depends on the black hole spins, which in turn can be affected by ultra-light bosons due to superradiance. We perform a detailed analysis of these effects and show that searches for stellar tidal disruptions have the potential to uncover the existence of ultra-light bosons. In particular, we find that upcoming stellar tidal disruption rate measurements by the Vera Rubin Observatory’s Legacy Survey of Space and Time can be used to either discover or rule out bosons with masses ranging from 10 –20 to 10 –18 eV. Our analysis also indicates that these measurements may be used to constrain a variety of supermassive black hole spin distributions and determine if close-to maximal spins are preferred.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Assessment and Performance of Flexible Quench Antenna Array Diagnostics for Superconducting Magnets

FNAL has been developing multiple versions of flexible quench antennas (flex-QA), including some specially optimized for high sensitivity and/or high resolution, to characterize quench events and transients during current ramping in superconducting magnets. A fundamental feature in our use of these is the creation of grid-like structures of sensitive elements to cover coil surfaces, with the aim of getting precise localization of magnetic flux-change events. The flex-QA are coupled with fast data-acquisition, allowing comprehensive analysis of signals at the relevant fine time scales. In addition to arrays of various flex-QA types being used during cryogenic testing of superconducting magnets, we also are utilizing a newly developed room temperature test stand to better understand QA response characteristics. The data from actual superconducting magnet tests, warm test stand measurements, and simulation data on the same QA designs allows us to draw conclusions on operational feasibility and plan better for improvements of our sensors. In this paper we present data from the multiple tests performed and analysis results. Flex-QA designs are compared, and their features, options, and optimization discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗