Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unit dynamics”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Step Towards Machine Learning-based Coherent Generator Grouping for Emergency Control Applications in Modern Power Grid

A dynamic neural network (NN) based multi-class classifier is proposed for improving online prediction of coherent generator groups (CGGs), following the occurrences of various contingencies in the power grid. This is motivated by the increasing availability of the measurements from phasor measurement units (PMUs) and the number of grouping schemes is limited. The proposed method consists of three steps. First, by performing offline simulations, a library of system dynamic responses characterized by post-contingency rotor angles and speeds of individual generators is obtained. To generate sufficient data, up to N-2 contingencies and the uncertain parameters associated with the power grid including type and location of disturbance and fault clearing times are modeled. Secondly, the training data-set is produced by generating labels for individual contingencies using a hierarchical clustering method based on rotor angle and speed data. Lastly, the dynamic NN models are trained for online applications such as emergency controls and controlled islanding. The proposed method is tested on the standard 16-generator 68-bus system to demonstrate its performance. Furthermore, the impact of the sample data lengths on the CGG numbers is evaluated. It is interesting to observe that the time domain stability behaviors can be determined by examining the changes in the CGG numbers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2025

Longitudinal Patterns in Nitrate Revealed Through Drone‐Based Measurements in an Agriculturally Influenced Midwestern U.S. River

Nitrate concentrations in streams and rivers in the Midwestern United States are often elevated, reflecting the predominance of agriculture in the surrounding landscape. Recent advances in technology, including surface water drones and more precise sensors, provide opportunities to investigate nitrate dynamics with high spatial and temporal resolution. We deployed an aquatic drone, the AquaBOT, in a sixth-order, agriculturally influenced river to examine longitudinal patterns in water quality. Our goal was to measure the spatial and temporal heterogeneity in nitrate and nitrate removal processes and determine the influence of tributary inputs on main stem chemistry. We navigated the drone along a 12-km reach of the Des Moines River (Iowa, USA) nine times between June 2021 and August 2022. Across the deployments, mean nitrate concentration was positively related to discharge and was nearly two orders of magnitude higher in spring than summer. We observed contrasting patterns in main stem nitrate, which decreased downstream during some runs (e.g., 3.1–2.7 mg N L −1 in June 2021), demonstrating net nitrate uptake along the reach, and remained constant on other dates. Similarly, tributaries to the Des Moines had a varied influence on riverine nitrate. Tributaries either increased or decreased main stem nitrate concentrations depending on the tributary and the date. Nitrate removal rates were spatially and temporally variable but showed some consistency at the subreach (2 km) scale, with two subreaches often showing elevated rates of nitrate removal across dates. In conclusion, our study reveals nuanced heterogeneity in nitrate dynamics of the Des Moines River despite the homogeneity of agricultural land cover in the watershed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Customized open source renewable energy models validated through PHIL lab experiments

Energy models for power systems require ongoing updates to reflect advancements in equipment technology and the increasing complexity of power electronic devices. This study utilizes a Power Hardware-in-the-Loop (PHIL) experimental setup to validate custom photovoltaic (PV) inverter models, aiming to enhance and expedite the development of advanced renewable energy models. The research compares the performance of a physical inverter with generic Renewable Energy Source (RES) models recommended by the Western Electricity Coordinating Council (WECC). As inverter-based renewable energy sources become more prevalent in modern electrical grids, it is crucial that dynamic models accurately represent their real-world behavior. Accurate models improve our understanding of these energy resources and their interactions with the grid. The proposed model enhancements are designed to better reflect real inverter performance, based on insights from PHIL experiments. These models are developed using the open source Modelica language and the OpenIPSL Modelica Library, allowing integration across various simulation tools without re-implementation. The paper concludes with a thorough assessment, comparing the enhanced models with PHIL experiments on a real PV inverter in a controlled laboratory setting. As a result, the study provides the enhanced WECC RES models and validation data as open source resources, facilitating further research and development.

Modelica↗
NASA NTRS2018

Hardware Assessment in Support of the Dynamic Power Convertor Development Effort

Stirling Radioisotope Power Systems (RPS) are being developed by NASAs RPS Program in collaboration with the U.S. Department of Energy (DOE). Efforts ranging from 2001 to 2015 enabled development of the Technology Demonstration Convertor (TDC) for use in the 110-watt Stirling Radioisotope Generator (SRG-110) and the Advanced Stirling Convertor (ASC) for use in the Advanced Stirling Radioisotope Generator (ASRG). The DOE selected Lockheed Martin Space Systems Company (LMSSC) as the system integration contractor for both flight development efforts. The SRG-110 housed two TDCs fabricated by Infinia and resulted in the production of 16x demonstration units and 2x engineering units. The project was redirected in 2006 to make use of a more efficient and lower mass ASCs under development by Sunpower Inc. The DOE managed the flight contract with LMSSC and subcontractor Sunpower Inc. from 2007 to 2013 to build the ASRG, with support from NASA Glenn Research Center (GRC). Sunpower Inc. held two parallel contracts to produce ASCs, one with Lockheed Martin to produce ASC-F flight units and one with GRC for the production of ASC-E3 engineering unit pathfinders that were used to refine the flight design and production processes. The DOE initiated termination of the ASRG contract in late 2013. After ASRG had ended, GRC completed characterization testing of the ASRG Engineering Unit 2 (EU2) and the GRC contract with Sunpower was also completed. The NASA RPS Program Office has recently initiated a new Dynamic Power Conversion development effort which includes the potential maturation of Stirling, Brayton, and Rankine power convertors for the next generation of RPS. The effort started with the request for proposal and review of submits. Contracts are anticipated for release in 2017 and will initially focus on a design phase prior to fabrication and testing. This new effort will focus on robustness in addition to high efficiency, specific power, and reliability. Also, some requirements introduced during the ASRG contract have also been included in the new effort, such as constant lateral loading. Due to the focus on robustness and new requirements relative to the older TDC design, the Stirling Cycle Development Project has initiated an assessment of government owned hardware to help inform requirements evolution and evaluation of future designs. While lessons learned from the ASRG flight development project have been taken into consideration, the evaluation of the TDC design had not been completed for some existing environments or relatively new requirements. To further assess the TDC design, a series of tasks were initiated to evaluate degradation for units that have operated unattended for over 105,000 hours, demonstrate robustness to a random vibration environment, characterize and evaluate performance for varying lateral load profiles. The status for each task are described.

Wilson, Scott D.↗
DOE OSTI2025

A framework for testing soil carbon dynamics post land-use transition in a multisector dynamics model

Soil carbon plays a crucial role in the global carbon cycle. Changes in land use can determine whether carbon is stored or is emitted into the atmosphere as carbon dioxide, which has broad implications for the human and Earth systems. These feedbacks to the carbon cycle and their socio-economic drivers are modelled by many global multisector dynamics models to project future possibilities for the human-Earth system. One notable model of this class is the Global Change Analysis Model (GCAM), which uses a simplified process to model soil organic carbon (SOC) content after land-use transition across 384 land units. While the current GCAM soil carbon framework is based on scientific principles, it has not been tested against experimental data. This work examines rates of SOC change from GCAM input data. Specifically, first order rate constants derived from model inputs were compared to values from two syntheses to assess GCAM’s accuracy. Welch’s t-tests and linear models were used to determine if rate constants were consistent across all tested geographical areas and land-use transition types. While we found that there was general agreement on the direction and magnitude (i.e., rate) of SOC change, the rate constant derived from GCAM and empirical values differed strongly in a subset of specific instances. These results indicate that GCAM’s current SOC dynamics during land use transition successfully capture broad patterns of change in this critical carbon pool, but should be interpreted with caution at finer spatial scales. One potential cause of these discrepancies is our highly aggregated variable, soil timescale, which could be made more granular to improve accuracy. When using economically rooted multisector dynamics models, such as GCAM, it is critical to understand such model limitations for representing specific Earth system processes.

carbon↗
NASA NTRS2017

Hardware Assessment in Support of the Dynamic Power Convertor Development Effort

Stirling Radioisotope Power Systems (RPS) are being developed by NASA's RPS Program in collaboration with the U.S. Department of Energy (DOE). Efforts ranging from 2001 to 2015 enabled development of the Technology Demonstration Convertor (TDC) for use in the 110-watt Stirling Radioisotope Generator (SRG-110) and the Advanced Stirling Convertor (ASC) for use in the Advanced Stirling Radioisotope Generator (ASRG). The DOE selected Lockheed Martin Space Systems Company (LMSSC) as the system integration contractor for both flight development efforts. The SRG-110 housed two TDCs fabricated by Infinia and resulted in the production of 16x demonstration units and 2x engineering units. The project was redirected in 2006 to make use of a more efficient and lower mass ASCs under development by Sunpower Inc. The DOE managed the flight contract with LMSSC and subcontractor Sunpower Inc. from 2007 to 2013 to build the ASRG, with support from NASA Glenn Research Center (GRC). Sunpower Inc. held two parallel contracts to produce ASCs, one with Lockheed Martin to produce ASC-F flight units and one with GRC for the production of ASC-E3 engineering unit pathfinders that were used to refine the flight design and production processes. The DOE initiated termination of the ASRG contract in late 2013. After ASRG had ended, GRC completed characterization testing of the ASRG Engineering Unit #2 (EU2) and the GRC contract with Sunpower was also completed. The NASA RPS Program Office has recently initiated a new Dynamic Power Conversion development effort which includes the potential maturation of Stirling, Brayton, and Rankine power convertors for the next generation of RPS. The effort started with the request for proposal and review of submits. Contracts are anticipated for release in 2017 and will initially focus on a design phase prior to fabrication and testing. This new effort will focus on robustness in addition to high efficiency, specific power, and reliability. Also, some requirements introduced during the ASRG contract have also been included in the new effort, such as constant lateral loading. Due to the focus on robustness and new requirements relative to the older TDC design, the Stirling Cycle Development Project has initiated an assessment of government owned hardware to help inform requirements evolution and evaluation of future designs. While lessons learned from the ASRG flight development project have been taken into consideration, the evaluation of the TDC design had not been completed for some existing environments or relatively new requirements. To further assess the TDC design, a series of tasks were initiated to evaluate degradation for units that have operated unattended for over 105,000 hours, demonstrate robustness to a random vibration environment, characterize and evaluate performance for varying lateral load profiles. The status for each task are described.

Wilson, Scott D.↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ecological acclimation: A framework to integrate fast and slow responses to climate change

Ecological responses to climate change occur across vastly different time-scales, from minutes for physiological plasticity to decades or centuries for community turnover and evolutionary adaptation. Accurately predicting the range of ecosystem trajectories will require models that incorporate both fast processes that may keep pace with climate change and slower ones likely to lag behind and generate disequilibrium dynamics. However, the knowledge necessary for this integration is currently fragmented across disciplines. We develop ‘ecological acclimation’ as a unifying framework to emphasize the similarity of dynamics driven by processes operating on dramatically different time-scales and levels of biological organization. The framework focuses on ecoclimate sensitivities, measured as the change in an ecological response variable per unit of climate change. Acclimation processes acting at different time-scales cause these sensitivities to shift in magnitude and even direction over time. We highlight shifting ecoclimate sensitivities in case studies from diverse ecosystems, including terrestrial plant communities, coral reefs and soil microbiomes. Models predicting future ecosystem states inevitably make assumptions about acclimation processes; these assumptions must be explicit for users to evaluate whether a model is appropriate for a given forecast horizon. Similarly, decision frameworks that clearly account for multiple acclimation processes and their distinct time-scales will help natural resource managers plan for ecological impacts of climate change from years to many decades into the future. We outline a synthetic research programme focused on the time-scales of ecological acclimation to reduce uncertainty in ecological forecasts.

climate adaptation↗
DOE OSTI2025

Probing Electrode Transformation under Dynamic Operation for Alkaline Water Electrolysis

Alkaline water electrolyzers (AWEs) play a pivotal role in the realm of large-scale hydrogen production. However, AWEs face significant challenges in electrode degradation particularly under dynamic operating conditions, induced by reverse current phenomenon during frequent startup/shutdown. Herein, this study aims to rationalize the degradation mechanisms of AWEs under these conditions. A three-electrode membrane electrode assembly (MEA) setup is first utilized to decouple polarization behaviors of anode and cathode in AWEs. Following a proposed accelerated stress testing protocol, the setup allows for tracking individual electrode performance transformations during frequent reverse current operation. Integrating operando cell studies with in situ and post-mortem characterizations, it is showed that continuous formation of highly active species, nickel (oxy)hydroxides, improves the anode performance for oxygen evolution reaction. On the contrary, irreversible oxidation of nickel to β-nickel hydroxide results in a severe degradation of cathode, leading to material dissolution, poor electrical conductivity and loss of catalytic activity for hydrogen evolution reaction. These results provide insights in nickel-based electrode transformation mechanisms for alkaline water electrolysis and indicate that cathode with higher redox reversibility can potentially improve durability of AWEs under dynamic conditions.

08 HYDROGEN↗
DOE OSTI2026

Alternating salt and freshwater floods of coastal soils impact soil structure, hydraulic properties, and oxygen dynamics

Coastal soils are increasingly impacted by hydrologic intensification in the form of rising sea level and flooding from storm surges and precipitation. Alternating saltwater (SW) and freshwater (FW) exposure has the potential to disperse colloids, which can lead to disintegration of soil structure, clogging of pore spaces, and reduction in ecologically and biogeochemically important functions like infiltration or gas exchange. To investigate how hydrologic intensification affects soil structure and oxygen dynamics, we conducted a series of laboratory-based flood simulations. Intact soil cores of the A and B horizons (22 total) from the toe slope of an upland coastal forest along the western shore of Chesapeake Bay. We subjected the cores to 24 h of saturation with either FW or alternating brackish SW and FW, with a 24-h draining event in between flood events. Oxygen diffusion into soil during draining was reduced by up to 30% for soils flooded alternating SW-FW compared to soils flooded only with FW. We attributed the reduced oxygen diffusion to clogging of smaller pores by colloids, with colloid redistribution observed as an increase in specific surface area down the core profile of up to 59%. After three SW-FW floods (six floods total), there were significant changes in pore size distribution, significant redistribution of colloids, and the A horizon became sodic. We concluded that a small number of SW flooding events can induce a measurable change in soil physical properties that directly impacts the biogeochemical dynamics.

Rod, Kenton A. [Pacific Northwest National Laborat↗
DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗
DOE OSTI2026

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗