Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Scalability analysis of heavy-duty gas turbines using data-driven machine learning

With the increasing integration of variable renewable energy sources into power systems, the role of flexible power generation technologies like gas turbines (GT) in rapid grid balancing remains crucial. This sustained importance underscores the need for scaled and precise modeling of GT to ensure effective integration within evolving energy frameworks. While physics-driven GT models integrate thermodynamics, fluid dynamics, and combustion principles, they often rely on approximate mathematical representations to accommodate scaling that may not capture the actual complex dynamics for GTs and inertial effects associated to GTs with different ratings. In this study, a data-driven model is proposed using machine learning (ML) techniques to conduct GT scalability analysis and performance evaluation with high accuracy. The ML model, trained on data from various operating conditions and performance parameters, aims to uncover intricate relationships and patterns, resembling GT characteristics at different scales (ratings). The model is developed to capture complex system interaction and to adapt to changing operational scenarios at different capacities, providing valuable insights of power system dynamics. In this study, the real-time digital simulator platform was employed to generate training data for the ML model and assess its dynamic characteristics. The ultimate objective was to develop a detailed modeling framework based on governing equations and data-driven ML capable of predicting key performance indicators, in thermal systems such as GTs, including power output, speed, fuel consumption, and exhaust temperature under diverse operating conditions at different scales. The developed ML framework demonstrated high accuracy, with mean relative errors for GT power prediction, reference speed, exhaust temperature, and compressor pressure ratio (CPR) parameters consistently below 0.1% across typical load fluctuation scenarios. Maximum deviations were limited to approximately 0.5 K for exhaust temperature and 0.009 for CPR, underscoring the model’s ability to replicating dynamic GT behavior with high precision. The adaptability of the ML model enables its application across diverse operational conditions and its extension to other thermal systems. By leveraging advanced ML techniques, this study presents a robust and scalable modeling framework that enhances GT simulation precision, facilitating improved integration into evolving power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Structurally Driven Selective Adsorption of Hydrocarbons by Metal Substitution in Isostructural Rare-Earth Metal–Organic Frameworks

The design and realization of highly selective nanoporous materials are necessary to target critical separations across industries. By leveraging pore size, pore shape, and linker functionalization, the design of nanoporous solid adsorbents will enable the rapid production of energy efficient separation materials for high-value gas mixtures. This study uses a combination of modeling, synthesis, and gas adsorption testing to investigate a new class of small-pore isostructural rare-earth (RE) 2,5-dihydroxyterephthalic acid (DOBDC) metal–organic frameworks (MOFs) (RE: Pr-, Gd-, Er-, Yb; DOBDC = 2,5-dihydroxyterephthalic acid) and their adsorption selectivity for acetylene/ethylene mixtures. Density functional theory simulations identified that selective binding of acetylene over ethylene in the Gd-, Er-, and Yb-DOBDC MOFs was due to hydrogen-bonding between acetylene and the linker hydroxyl. Adsorption experiments validated the computational results by identifying mechanisms that control the acetylene/ethylene adsorption selectivity and high acetylene adsorption. Furthermore, dynamic column breakthrough experiments with the Gd-DOBDC MOF validated the simulations and indicated that ethylene can be separated from acetylene in a mixture containing 1 vol % acetylene and 39 vol % ethylene (balance argon). In conclusion, the results highlight the complexity of gas binding in functional porous materials and how combining modeling and experiment enables a fundamental understanding of gas–framework interactions that can be leveraged for the design of future separation materials.

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

Optimal Linear Water Density for Proton Transport in Tunnel Oxides

Proton intercalation and transport is key to pseudocapacitive energy storage of oxide electrodes in acidic electrolytes. Although proton transport in water confined in layered or 2D materials is attracting great interest, much less is known about how the proton is transported in the 1D channels of tunnel oxides such as hexagonal WO 3 (h-WO 3 ). Here, we use first-principles molecular dynamics to reveal an optimal linear density of four water molecules per nanometer that yields the highest proton diffusivity. The volcano shape of proton diffusivity versus linear water density is a result of balancing the linear hydrogen-bond chain and the rotation of the water molecules to enable the Grotthuss mechanism. Furthermore, this insight provides a unifying view of proton transport along a single file of water molecules confined in hydrophilic 1D channels.

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

Long-range exciton diffusion in molecular non-fullerene acceptors

Abstract The short exciton diffusion length associated with most classical organic semiconductors used in organic photovoltaics (5-20 nm) imposes severe limits on the maximum size of the donor and acceptor domains within the photoactive layer of the cell. Identifying materials that are able to transport excitons over longer distances can help advancing our understanding and lead to solar cells with higher efficiency. Here, we measure the exciton diffusion length in a wide range of nonfullerene acceptor molecules using two different experimental techniques based on photocurrent and ultrafast spectroscopy measurements. The acceptors exhibit balanced ambipolar charge transport and surprisingly long exciton diffusion lengths in the range of 20 to 47 nm. With the aid of quantum-chemical calculations, we are able to rationalize the exciton dynamics and draw basic chemical design rules, particularly on the importance of the end-group substituent on the crystal packing of nonfullerene acceptors.

14 SOLAR ENERGY↗
DOE OSTI2024

Synthesis of the land carbon fluxes of the Amazon region between 2010 and 2020

The Amazon is the largest continuous tropical forest in the world and plays a key role in the global carbon cycle. Human-induced disturbances and climate change have impacted the Amazon carbon balance. Here we conduct a comprehensive synthesis of existing state-of-the-art estimates of the contemporary land carbon fluxes in the Amazon using a set of bottom-up methods (i.e., dynamic vegetation models and bookkeeping models) and a top-down inversion (atmospheric inversion model) over the Brazilian Amazon and the whole Biogeographical Amazon domain. Over the whole biogeographical Amazon region bottom-up methodologies suggest a small average carbon sink over 2010-2020, in contrast to a small carbon source simulated by top-down inversion (2010-2018). However, these estimates are not significantly different from one another when accounting for their large individual uncertainties, highlighting remaining knowledge gaps, and the urgent need to reduce such uncertainties. Nevertheless, both methodologies agreed that the Brazilian Amazon has been a net carbon source during recent climate extremes and that the south-eastern Amazon was a net land carbon source over the whole study period (2010-2020). Overall, our results point to increasing human-induced disturbances (deforestation and forest degradation by wildfires) and reduction in the old-growth forest sink during drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

No evidence of canopy-scale leaf thermoregulation to cool leaves below air temperature across a range of forest ecosystems

Understanding and predicting the relationship between leaf temperature (T leaf ) and air temperature (T air ) is essential for projecting responses to a warming climate, as studies suggest that many forests are near thermal thresholds for carbon uptake. Based on leaf measurements, the limited leaf homeothermy hypothesis argues that daytime T leaf is maintained near photosynthetic temperature optima and below damaging temperature thresholds. Specifically, leaves should cool below T air at higher temperatures (i.e., > ~25–30°C) leading to slopes <1 in T leaf /T air relationships and substantial carbon uptake when leaves are cooler than air. This hypothesis implies that climate warming will be mitigated by a compensatory leaf cooling response. A key uncertainty is understanding whether such thermoregulatory behavior occurs in natural forest canopies. We present an unprecedented set of growing season canopy-level leaf temperature (T can ) data measured with thermal imaging at multiple well-instrumented forest sites in North and Central America. Our data do not support the limited homeothermy hypothesis: canopy leaves are warmer than air during most of the day and only cool below air in mid to late afternoon, leading to T can /T air slopes >1 and hysteretic behavior. We find that the majority of ecosystem photosynthesis occurs when canopy leaves are warmer than air. Using energy balance and physiological modeling, we show that key leaf traits influence leaf-air coupling and ultimately the T can /T air relationship. Canopy structure also plays an important role in T can dynamics. Future climate warming is likely to lead to even greater T can , with attendant impacts on forest carbon cycling and mortality risk.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

CoREC: Scalable and Resilient In-memory Data Staging for In-situ Workflows

The dramatic increase in the scale of current and planned high-end HPC systems is leading new challenges, such as the growing costs of data movement and IO, and the reduced mean time between failures (MTBF) of system components. In-situ workflows, i.e., executing the entire application workflows on the HPC system, have emerged as an attractive approach to address data-related challenges by moving computations closer to the data, and staging-based frameworks have been effectively used to support in-situ workflows at scale. However, the resilience of these staging-based solutions has not been addressed, and they remain susceptible to expensive data failures. Furthermore, naive use of data resilience techniques such as n-way replication and erasure codes can impact latency and/or result in significant storage overheads. In this article, we present CoREC, a scalable and resilient in-memory data staging runtime for large-scale in-situ workflows. CoREC uses a novel hybrid approach that combines dynamic replication with erasure coding based on data access patterns. It also leverages multiple levels of replications and erasure coding to support diverse data resiliency requirements. Furthermore, the article presents optimizations for load balancing and conflict-avoiding encoding, and a low overhead, lazy data recovery scheme. We have implemented the CoREC runtime and have deployed with the DataSpaces staging service on leadership class computing machines and present an experimental evaluation in the article. Here, the experiments demonstrate that CoREC can tolerate in-memory data failures while maintaining low latency and sustaining high overall storage efficiency at large scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Insights into Methodologies and Stochastic Optimization of Thermal Energy Storage-Coupled Advanced Reactor Systems: A Comparison of Methods for Accessing Long-Term Sub-System Sizing Adequacy

This paper investigates the potential of coupling Thermal Energy Storage (TES) with Advanced Reactors (ARs) to address uncertainties posed by climate change in deep decarbonized power systems. The TES Use-case Team at Idaho National Laboratory (INL) has examined the potential of storing thermal energy from ARs during low demand periods and optimizing discharge during peak-priced hours, in both steady-state and transient conditions. Building on this groundwork, this study bridges the gaps in optimal sizing of the sub-system of TES-coupled AR systems using Risk Analysis Virtual Environment (RAVEN) and Holistic Energy Resource Optimization Network (HERON), INL?s framework for grid optimization. By applying this framework, we present statistically-robust optimal charge, discharge including balance of plant (BOP), and storage sizing for the High-Temperature Gas-Cooled Reactor (HTGR) with 203 MWth output. To this end, we generated synthetic price samples for 30 years using 2018 ? 2021 real-time market data from ERCOT, PJM and MISO. Our results reveals that the TES-coupled HTGR system is highly effective in maximizing revenue from electricity sales. We observed a substantial increase of 40 % in ERCOT and a noteworthy 15 % increase in PJM and MISO when compared to the conventional BOP without TES. This improvement is achieved through regionally-tailored sub-system sizing, which ranges from 398 to 416 MWth for the discharge system and 610 to 1029 MWth for the TES. We find that the average electricity price directly impacts the overall economics, while price volatility influences storage size. Additional sensitivity analyses were performed to access the impact of key assumptions on system economics and sizing, differentiating the optimization window (i.e., 24 ? 219 hours of chronological observations) and by imposing storage continuity condition in tracking TES cycles. We observed that at the 120-hour of the optimization window, a reasonable balance between computation time and accuracy was achieved. Our analysis also highlights the significance of conducting multi-day cycle analysis (> 120-hour) for TES to capture interaction between electricity prices and storage dynamics, providing a comprehensive understanding of TES behavior that AR developers should integrate into their plant designs.

25 ENERGY STORAGE↗
DOE OSTI2026

Adaptive_Evolution_2026

Thermoanaerobacterium saccharolyticum, an anaerobic and thermophilic bacterium capable of metabolizing sugar monomers and soluble oligomers into ethanol, has been proposed for use in consolidated bioprocessing in coculture with compatible cellulolytic bacteria such asClostridium thermocellum. Although the mixed acid fermentation of both of these strains has been engineered to produce ethanol as the major fermentation product, the maximum titer produced thus far byC. thermocellum, about 3% (w/v) ethanol, is half that produced byT. saccharolyticum. There is thus motivation to understand the mechanistic basis of the robust ethanol pathway inT. saccharolyticumso that key features can be recapitulated inC. thermocellumand other organisms.Previously, we characterized theindividualrole of the main genes responsible for electron transfer in the ethanol production ofT. saccharolyticum. However, the consequences of thecombinedloss of function of all these genes have not been investigated, nor has the way in which fermentative metabolism adapts to such constraints. In this work, we combined knockouts of ferredoxin nicotinamide oxidoreductase (fnor) genes (nfnAandnfnB) and hydrogenase genes (hydAandhfsD) and studied their effects on fermentation and grow. We showed that these genetic modifications together impair growth and decrease electron transfer from reduced ferredoxin, thereby redirecting flux from the pyruvate ferredoxin oxidoreductase enzyme to the pyruvate formate lyase enzyme. We also performed adaptive evolution of these mutants to rescue their growth, and most notably, we observed a single nucleotide variation in the alcohol dehydrogenaseadhAgene. Through molecular dynamics simulations and enzymatic assays, we determined that this point mutation causes a structural change that impairs the AdhA specificity for the NADPH cofactor and increases NADH-linked activity to restore redox balance. These findings consolidate our understanding of the functioning of electron transfer pathways in this organism.

alcohol dehydrogenase↗
DOE OSTI2024

HarDWR - Harmonized Water Rights Records

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - Recalculated based on data sourced from WestDAAT - Changed using a Site ID column to identify unique records to using aa combination of Site ID and Allocation ID - Removed the Water Management Area (WMA) column from the harmonized records. The replacement is a separate file which stores the relationship between allocations and WMAs. This allows for allocations to contribute to water right amounts to multiple WMAs during the subsequent cumulative process. - Added a column describing a water rights legal status - Added "Unspecified" was a water source category - Added an acre-foot (AF) column - Added a column for the classification of the right's owner v1.02 - Added a .RData file to the dataset as a convenience for anyone exploring our code. This is an internal file, and the one referenced in analysis scripts as the data objects are already in R data objects. v1.01 - Updated the names of each file with an ID number less than 3 digits to include leading 0s v1.0 - Initial public release Description Here we present an updated database of Western U.S. water right records. This database provides consistent unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of seven broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of inter-sectoral water allocation changes though water markets (Grogan et al., *in review*). Specifically, the data were formatted for use as input to a process-based hydrologic model, Water Balance Model (WBM), with a water rights module (Grogan et al., *in review*). While this specific study motivated the development of the database presented here, water management in the U.S. West is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al., 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. West. We produced the water rights database presented here in four main steps: (1) data collection, (2) data quality control, (3) data harmonization, and (4) generation of cumulative water rights curves. Each of steps (1)-(3) had to be completed in order to produce (4), the final product that was used in the modeling exercise in Grogan et al. (*in review*). All data in each step is associated with a spatial unit called a Water Management Area (WMA), which is the unit of water right administration utilized by the state in which the right came from. Steps (2) and (3) required use to make assumptions and interpretation, and to remove records from the raw data collection. We describe each of these assumptions and interpretations below so that other researchers can choose to implement alternative assumptions an interpretation as fits their research aims. Motivation for Changing Data Sources The most significant change has been a switch from collecting the raw water rights directly from each state to using the water rights records presented in WestDAAT, a product of the Water Data Exchange (WaDE) Program under the Western States Water Council (WSWC). One of the main reasons for this is that each state of interest is a member of the WSWC, meaning that WaDE is partially funded by these states, as well as many universities. As WestDAAT is also a database with consistent categorization, it has allowed us to spend less time on data collection and quality control and more time on answering research questions. This has included records from water right sources we had previously not known about when creating v1.0 of this database. The only major downside to utilizing the WestDAAT records as our raw data is that further updates are tied to when WestDAAT is updated, as some states update their public water right records daily. However, as our focus is on cumulative water amounts at the regional scale, it is unlikely most records updates would have a significant effect on our results. The structure of WestDAAT led to several important changes to how HarWR is formatted. The most significant change is that WaDE has calculated a field known as `SiteUUID`, which is a unique identifier for the Point of Diversion (POD), or where the water is drawn from. This separate from `AllocationNativeID`, which is the identifier for the allocation of water, or the amount of water associated with the water right. It should be noted that it is possible for a single site to have multiple allocations associated with it and for an allocation to be able to be extracted from multiple sites. The site-allocation structure has allowed us to adapt a more consistent, and hopefully more realistic, approach in organizing the water right records than we had with HarDWR v1.0. This was incredibly helpful as the raw data from many states had multiple water uses within a single field within a single row of their raw data, and it was not always clear if the first water use was the most important, or simply first alphabetically. WestDAAT has already addressed this data quality issue. Furthermore, with v1.0, when there were multiple records with the same water right ID, we selected the largest volume or flow amount and disregarded the rest. As WestDAAT was already a common structure for disparate data formats, we were better able to identify sites with multiple allocations and, perhaps more importantly, allocations with multiple sites. This is particularly helpful when an allocation has sites which cross WMA boundaries, instead of just assigning the full water amount to a single WMA we are now able to divide the amount of water between the number of relevant WMAs. As it is now possible to identify allocations with water used in multiple WMAs, it is no longer practical to store this information within a single column. Instead the stAllocationToWMATab.csv file was created, which is an allocation by WMA matrix containing the percent Place of Use area overlap with each WMA. We then use this percentage to divide the allocation's flow amount between the given WMAs during the cumulation process to hopefully provide more realistic totals of water use in each area. However, not every state provides areas of water use, so like HarDWR v1.0, a hierarchical decision tree was used to assign each allocation to a WMA. First, if a WMA could be identified based on the allocation ID, then that WMA was used; typically, when available, this applied to the entire state and no further steps were needed. Second was the spatial analysis of Place of Use to WMAs. Third was a spatial analysis of the POD locations to WMAs, with the assumption that allocation's POD is within the WMA it should belong to; if an allocation still had multiple WMAs based on its POD locations, then the allocation's flow amount would be divided equally between all WMAs. The fourth, and final, process was to include water allocations which spatially fell outside of the state WMA boundaries. This could be due to several reasons, such as coordinate errors / imprecision in the POD location, imprecision in the WMA boundaries, or rights attached with features, such as a reservoir, which crosses state boundaries. To include these records, we decided for any POD which was within one kilometer of the state's edge would be assigned to the nearest WMA. Other Changes WestDAAT has Allowed In addition to a more nuanced and consistent method of assigning water right's data to WMAs, there are other benefits gained from using the WestDAAT dataset. Among those is a consistent categorization of a water right's legal status. In HarDWR v1.0, legal status was effectively ignored, which led to many valid concerns about the quality of the database related to the amounts of water the rights allowed to be claimed. The main issue was that rights with legal status' such as "application withdrawn", "non-active", or "cancelled" were included within HarDWR v1.0. These, and other water rights status' which were deemed to not be in use have been removed from this version of the database. Another major change has been the addition of the "unspecified water source category. This is water that can come from either surface water or groundwater, or the source of which is unknown. The addition of this source category brings the total number of categories to three. Due to reviewer feedback, we decided to add the acre-foot (AF) column so that the data may be more applicable to a wider audience. We added the ownerClassification column so that the data may be more applicable to a wider audience. File Descriptions The dataset is a series of various files organized by state sub-directories. In addition, each file begins with the state's name, in case the file is separate from its sub-directory for some reason. After the state name is the text which describes the contents of the file. Here is each file described in detail. Note that st is a placeholder for the state's name. stFullRecords_HarmonizedRights.csv: A file of the complete water records for each state. The column headers for each of this type of file are: state - The name of the state to which the allocations belong to. FIPS - The two digit numeric state ID code. siteID - The site location ID for POD locations. A site may have multiple allocations, which are the actual amount of water which can be drawn. In a simplified hypothetical, a farm stead may have an allocation for "irrigation" and an allocation for "domestic" water use, but the water is drawn from the same pumping equipment. It should be noted that many of the site ID appear to have been added by WaDE, and therefore may not be recognized by a given state's water rights database. allocationID - The allocation ID for the water right. For most states this is the water right ID, and what is recommended to use should a right be looked up on a given state's water rights database. The water amounts associated with these IDs tend to be finer scaled than those associated with siteID. It should be noted that some allocations may be extracted from multiple sites, particularly for larger Places of Use. ownerClassification - A classification of the types of owners for water rights. The most common is `Private` which incorporates a wide range of entities. Several classifications would be grouped into a government category, most of which are for the U.S. Federal Government. These allocations could be listed as "Federal", "United States of America", or as the names of any number of federal agencies. The last major grouping of entities is for "Native American"s. priorityDate - The date we use as the water right priority date for our modeling analysis. This is the legal priority date when it is available. However, for some rights, specifically from California and New Mexico, we used a pseudo priority date (e.g. well completion date or start of well drilling date) when a legal priority date was not available. The most questionable dates come from New Mexico, where the only date associated with certain water right records was the date the allocation was recorded in the database. As the allocation record creation tended to be within a few months of the filing of the application of the water right, from manually double checking the water rights, and our analysis focuses on aggregating water rights on the timescale of years, we determined it was acceptable to use such dates to include as many records as possible. primaryBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories WestDAAT. This column is the original WaDE category for the primary water use at the PoD site. allocationBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories for WestDAAT. This column is the original WaDE category

Economics↗
DOE OSTI2023

Nonlinear dynamics, bifurcations, and multi-stability in a vibro-impact system with geometric and multi-segmented freeplay nonlinearities

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. Here, the effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

42 ENGINEERING↗
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Optimizing the fuel efficiency of an opposed piston engine for electric power generation

This paper investigates the optimal crankshaft motion for an opposed piston (OP) engine in a novel hybrid architecture to maximize fuel efficiency. The OP engine was selected for this work due to its inherent thermodynamic benefits and the balanced nature of the engine which can achieve downsizing through reducing the number of cylinders rather than the individual cylinder volume. The typical geartrain required on an OP engine was exchanged for two electric motors, reducing friction loss and decoupling the crankshafts. Using the motors to control the crankshaft motion profiles, this architecture introduces capabilities to dynamically vary compression ratio, combustion volume, and scavenging dynamics. To leverage these opportunities, an optimization scheme was developed utilizing nonlinear optimization of a 0-D model to compute the crankshaft motion profile that maximizes the work generated by the system. This optimization was then iteratively coupled with a high fidelity model which supplies the cylinder flow boundary conditions. This iterative approach reduces the model complexity used in the optimal control problem (OCP) while capturing the gas exchange dynamics critical to the 2-stroke cycle of the OP engine. By using the rate of change of motor torque as the input to the OCP, the torque fluctuation in a single cycle can be limited to ensure tracking feasibility. The results show crankshaft velocity slows during the compression stroke and conversely accelerates during the expansion stroke, reducing the peak motor torque required for control and thus reducing the motor losses. The extended residence time at top dead center, however, leads to an increase in heat transfer, illustrating the trade-off between the work extraction efficiency and the indicated engine efficiency.

Engineering↗
DOE OSTI2023

Modeling simultaneous particle shrinkage, dissolution and breakage using the modified moving grid technique

Simultaneous shrinkage, dissolution and breakage are important particle size reduction phenomena that characterize processes like the reactive degradation of solid chemicals. The dynamics of the particle size distribution (PSD) for such processes are non-trivial to model due to the number expending processes brought about by the eventual dissolution of particles. To this end, Population Balance Model (PBM) resolved through the sectional techniques is the natural approach. Here, we introduce a modified Moving Grid technique (m-MGT) to accurately resolve the particle size reduction phenomena. Our technique mimics the perpetual particle shrinkage through a continuously left-moving size grid and incorporates a strategic grid removal routine to capture the disappearance of particles. Coupled with the Fixed Pivot (FP) discretization for breakage, our m-MGT not only preserves the moment-related properties, but also benchmarked very well against the analytical number densities and exhibited a minimum of first-order convergence in all assessed case studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Neural Network for Principle of Least Action

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager–Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the “ground-truth” molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ 38 from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

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

Laisk measurements in the nonsteady state: Tests in plants exposed to warming and variable CO2 concentrations

Abstract Light respiration (RL) is an important component of plant carbon balance and a key parameter in photosynthesis models. RL is often measured using the Laisk method, a gas exchange technique that is traditionally employed under steady-state conditions. However, a nonsteady-state dynamic assimilation technique (DAT) may allow for more rapid Laisk measurements. In 2 studies, we examined the efficacy of DAT for estimating RL and the parameter Ci* (the intercellular CO2 concentration where Rubisco's oxygenation velocity is twice its carboxylation velocity), which is also derived from the Laisk technique. In the first study, we compared DAT and steady-state RL and Ci* estimates in paper birch (Betula papyrifera) growing under control and elevated temperature and CO2 concentrations. In the second, we compared DAT-estimated RL and Ci* in hybrid poplar (Populus nigra L. × P. maximowiczii A. Henry “NM6”) exposed to high or low CO2 concentration pre-treatments. The DAT and steady-state methods provided similar RL estimates in B. papyrifera, and we found little acclimation of RL to temperature or CO2; however, Ci* was higher when measured with DAT compared to steady-state methods. These Ci* differences were amplified by the high or low CO2 pre-treatments. We propose that changes in the export of glycine from photorespiration may explain these apparent differences in Ci*.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Granal thylakoid structure and function: explaining an enduring mystery of higher plants

Summary In higher plants, photosystems II and I are found in grana stacks and unstacked stroma lamellae, respectively. To connect them, electron carriers negotiate tortuous multi‐media paths and are subject to macromolecular blocking. Why does evolution select an apparently unnecessary, inefficient bipartition? Here we systematically explain this perplexing phenomenon. We propose that grana stacks, acting like bellows in accordions, increase the degree of ultrastructural control on photosynthesis through thylakoid swelling/shrinking induced by osmotic water fluxes. This control coordinates with variations in stomatal conductance and the turgor of guard cells, which act like an accordion's air button. Thylakoid ultrastructural dynamics regulate macromolecular blocking/collision probability, direct diffusional pathlengths, division of function of Cytochrome b 6 f complex between linear and cyclic electron transport, luminal pH via osmotic water fluxes, and the separation of pH dynamics between granal and lamellar lumens in response to environmental variations. With the two functionally asymmetrical photosystems located distantly from each other, the ultrastructural control, nonphotochemical quenching, and carbon‐reaction feedbacks maximally cooperate to balance electron transport with gas exchange, provide homeostasis in fluctuating light environments, and protect photosystems in drought. Grana stacks represent a dry/high irradiance adaptation of photosynthetic machinery to improve fitness in challenging land environments. Our theory unifies many well‐known but seemingly unconnected phenomena of thylakoid structure and function in higher plants.

59 BASIC BIOLOGICAL SCIENCES↗