Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Linear Model”

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 505 records · Page 28

Time-lapse approach to correct deficiencies of 2D soil zymography

Membrane zymography is commonly used in rhizosphere ecology to map enzyme activities in intact soil samples and plant roots. The method consists of incubating a membrane saturated with an enzyme-specific fluorogenic substrate on the soil/root surface followed by measurements of fluorescence intensity of the product in the membrane. The traditional zymography is based on assumptions that fluorescence on membrane images is linearly increasing with time during zymography and an increase in product content is numerically equal to enzyme activity in soil below the membrane. These assumptions are unlikely to hold in experimental settings. In this work we introduce a new zymography technique, time-lapse zymography (TLZ); the approach that eliminates the need for assumptions of the traditional zymography and provides more realistic estimates of enzymatic activities. We assessed the new technique in a series of laboratory and modeling experiments, including quantification of the fluorescent product diffusion (e.g. MUF: 4-methylumbelliferone) from plant roots, enzyme activity measurements on the roots, and HYDRUS-2D & HP2 software calibration with obtained data. The calibrated model was used to analyze the processes governing spatial and temporal dynamics of MUF contents in the membrane. The results indicated that the enzyme diffusion within the membrane-soil system was negligible, and measured zymograms were adequately reproduced solely by accounting for substrate and product diffusions and for catalytic enzyme reaction described by the Michaelis-Menten equation. TLZ enabled identifying and using linear parts on MUF time series and accounting for MUF losses in each zymogram pixel, considerably improving the accuracy as compared to the traditional method. Results demonstrated that enzymatic activity from only a thin soil layer (~0.2 mm) is adequately represented in zymograms.

54 ENVIRONMENTAL SCIENCES↗

Many-Body Dynamical Localization in a Kicked Lieb-Liniger Gas

The kicked rotor system is a textbook example of how classical and quantum dynamics can drastically differ. The energy of a classical particle confined to a ring and kicked periodically will increase linearly in time whereas in the quantum version the energy saturates after a finite number of kicks. The quantum system undergoes Anderson localization in angular-momentum space. Conventional wisdom says that in a many-particle system with short-range interactions the localization will be destroyed due to the coupling of widely separated momentum states. Here we provide evidence that for an interacting one-dimensional Bose gas, the Lieb-Liniger model, the dynamical localization can persist at least for an unexpectedly long time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DeepPhenoMem V1.0: deep learning modelling of canopy greenness dynamics accounting for multi-variate meteorological memory effects on vegetation phenology

Abstract. Vegetation phenology plays a key role in controlling the seasonality of ecosystem processes that modulate carbon, water and energy fluxes between the biosphere and atmosphere. Accurate modelling of vegetation phenology in the interplay of Earth's surface and the atmosphere is thus crucial to understand how the coupled system will respond to and shape climatic changes. Phenology is controlled by meteorological conditions at different timescales: on the one hand, changes in key meteorological variables (temperature, water, radiation) can have immediate effects on the vegetation development; on the other hand, phenological changes can be driven by past environmental conditions, known as memory effects. However, the processes governing meteorological memory effects on phenology are not completely understood, resulting in their limited performance of vegetation phenology represented in land surface models. A deep learning model, specifically a long short-term memory network (LSTM), has the potential to capture and model the meteorological memory effects on vegetation phenology. Here, we apply the LSTM to model the vegetation phenology using meteorological drivers and high-temporal-resolution canopy greenness observations through digital repeat photography by the PhenoCam network. We compare a multiple linear regression model, a no-memory-effect LSTM model and a full-memory-effect LSTM model to predict the whole seasonal greenness trajectory and the corresponding phenological transition dates across 50 sites and 317 site years during 2009–2018, covering deciduous broadleaf forests, evergreen needleleaf forests and grasslands. Results show that the deep learning model outperforms the multiple linear regression model, and the full-memory-effect LSTM model performs better than the no-memory-effect model for all three plant function types (median R2 of 0.878, 0.957 and 0.955 for broadleaf forests, evergreen needleleaf forests and grasslands). We also find that the full-memory-effect LSTM model is capable of predicting the seasonal dynamic variations of canopy greenness and reproducing trends in shifting phenological transition dates. We also performed a sensitivity analysis of the full-memory-effect LSTM model to assess its plausibility, revealing its coherence with established knowledge of vegetation phenology sensitivity to meteorological conditions, particularly changes in temperature. Our study highlights that (1) multi-variate meteorological memory effects play a crucial role in vegetation phenology, and (2) deep learning opens up new avenues for improving the representation of vegetation phenological processes in land surface models via a hybrid modelling approach.

Geology↗

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗

Advanced Models and Fast Simulation Algorithms for MARS Plants

With the successful commercialization of modular multilevel converters (MMC) and rapid increase in photovoltaic (PV) penetration and energy storage systems (ESS), one of the research areas being explored dynamically is integration of PV systems and ESS systems to MMCs. On this front, an integrated development of PV plants, ESS, and alternating- and direct-current (ac/dc) systems through Multi-Port Autonomous Reconfigurable Solar (MARS) plant is proposed and being researched upon. However, simulation of MARS systems on traditional existing software takes very long time (order of days) due to the presence of a large number of states and non-linear devices. Therefore, development of fast high fidelity models is of paramount importance to test the hardware design, control design and planning of power system expansions. In this context, two research objectives are addressed in this paper: firstly, an ultra fast single- or multi- CPU simulation algorithm to simulate the MARS system based on state-space models, hybrid discretization algorithm with a relaxation technique that reduces the imposed computational burden and a multi-rate method is presented. Secondly, a reduced order model of MARS is also presented. Both the developed techniques are validated with respect to reference PSCAD/EMTDC model.

Marthi, Phani Ratna Vanamali↗

Simulation of high-frequency dissolved oxygen dynamics in a shallow estuary, the Corsica River, Chesapeake Bay

Understanding shallow water biogeochemical dynamics is a challenge in coastal regions, due to the presence of highly variable land-water interface fluxes, tight coupling with sediment processes, tidal dynamics, and diurnal variability in biogeochemical processes. While the deployment of continuous monitoring devices has improved our understanding of high-frequency (12 - 24 hours) variability and spatial heterogeneity in shallow regions, mechanistic modeling of these dynamics has lagged behind conceptual and empirical models. The inherent complexity of shallow water systems is represented in the Corsica River estuary, a small basin within the Chesapeake Bay ecosystem, where abundant monitoring data have been collected from long-term monitoring stations, continuous monitoring sensors, synoptic sensor surveys, and measurements of sediment-water fluxes. A state-of-the-art modeling system, the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), was applied to the Corsica domain with a high-resolution grid and nutrient loads from the most recent version of the Chesapeake Bay watershed model. The Corsica SCHISM model reproduced observed high-frequency variability in dissolved oxygen, as well as seasonal variability in chlorophyll-a and sediment-water fluxes. Time-series signal analyses using Empirical Model Decomposition and spectral analysis revealed that the diurnal and M2 tide frequencies are the dominant high-frequency modes and physical transport contributes a larger share to dissolved oxygen budgets than biogeochemical processes on an hourly time scale. Heterogeneity and patchiness in dissolved oxygen resulting from phytoplankton distributions and geometry-driven eddies amplify the physical transport effect, and on longer time scales oxygen is controlled more by photosynthesis and respiration. Our simulation demonstrates that interactions among physical and biological dynamics generate complex high-frequency variability in water quality and non-linear reposes to nutrient loading and environmental forcing in shallow water systems.

54 ENVIRONMENTAL SCIENCES↗

Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach

Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Outcomes and Insights From Simplified Analytic Trajectory Optimization for a Tethered Underwater Kite

This letter formulates and solves a periodic trajectory optimization problem for a tethered underwater kite. The goal is to maximize the average mechanical power harvested by the kite. The type of kite considered in this letter extracts electricity from ocean currents by moving cross-current as its reel away from its base station, and consumes electricity to reel back. The problem of optimizing this kite’s trajectory is challenging due to the high dimensionality and nonlinearity of its dynamics. To tackle this challenge, the literature often separates the problem into two subproblems focusing on optimizing the cross-current and the reel-in/reel-out components of the trajectory, respectively, which may be sub-optimal. In contrast, this letter solves for the combined cross-current and reel-in/reel-out trajectory by linearizing the dynamics of the kite around a zero-power reference equilibrium trajectory in spherical coordinates. This allows the trajectory optimization problem to be solved analytically for simple sinusoidal input perturbations from equilibrium. Here we use linear quadratic regulation to enable the nonlinear kite model to track the optimized trajectory. The result is a computationally efficient approach that achieves an attractive Loyd factor of 19.9%, while providing important insights into the nature of the optimal trajectory.

42 ENGINEERING↗

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Investigation of the thermal conductivity of molten LiF-NaF-KF with experiments, theory, and equilibrium molecular dynamics

Molten salts are being proposed for numerous advanced energy applications, including advanced nuclear reactors, concentrating solar power plants, thermal energy storage, and fusion reactors. Accurate knowledge of the thermophysical properties of molten salts directly impact the performance of these energy systems and are essential for design and safety analyses. Thermal conductivity data for fluoride molten salts and mixtures are especially lacking. In this work, experimental measurements of thermal conductivity using the steady-state variable gap technique were performed on eutectic LiF-NaF-KF from 834 to 1195 K. The experiment accounts for radiative, convective, and conductive heat losses. In addition, theoretical and molecular dynamics models are used, from 750 K up to 1300 K, to estimate the thermal conductivity for comparison with the experimental results. The results of experiments show a weak negative deviation of thermal conductivity with temperature, unlike previous experimental results in the literature. The measured thermal conductivity magnitudes agree with the theoretical and molecular dynamics predictions, aside from the data above 1100 K, where heat losses and radiative errors are the most significant, having a 16% maximum deviation from theory. These experimental results provide new thermal conductivity data for the LiF-NaF-KF system and further validation of the predictive models. The theoretical model was used to map the composition and temperature dependent thermal conductivity of LiF-NaF-KF and the mapping’s deviation from a linear additivity estimation of thermal conductivity. Additionally, this mapping showed the highest deviations from linearity for KF-LiF rich mixtures and increasing deviation with temperature. Notably, the deviation from linearity near the LiF-NaF-KF eutectic composition was around 25%.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extinction Dynamics of the Foot-and-Mouth Disease Virus Carrier State Under Natural Conditions

Foot-and-mouth disease (FMD) is one of the most economically important livestock diseases worldwide. Following the clinical phase of FMD, a large proportion of ruminants remain persistently infected for extended periods. Although extinction of this carrier state occurs continuously at the animal and population levels, studies vary widely in their estimates of the duration of persistent infection. There is a need for robust statistical models to capture the dynamics of persistent infection for the sake of guiding FMD control and trade policies. The goal of the current study was to develop and assess statistical models to describe the extinction of FMD virus (FMDV) persistent infection using data from primary longitudinal studies of naturally infected cattle and Asian buffalo in Vietnam and India. Specifically, accelerated failure time (AFT) models and generalized linear mixed models (GLMM) were developed to predict the probability of persistent infection in seropositive animals and identified carriers at the individual animal level at sequential time points after outbreaks. The primary studies were analyzed by country and combined using an individual-participant data meta-analysis approach. The models estimated similar trends in the duration of persistent infection for the study/species groups included in the analyses, however the significance of the trends differed between the models. The overall probabilities of persistent infection were similar as predicted by the AFT and GLMM models: 6 months: 99% (AFT) /80% (GLMM), 12 months: 51% (AFT) /32% (GLMM), 18 months: 6% (AFT) /5% (GLMM), 24 months: 0.8% (AFT) /0.6% (GLMM). These models utilizing diverse and robust data sets predict higher probabilities of persistence than previously published, suggesting greater endurance of carriers subsequent to an outbreak. This study demonstrates the utility of statistical models to investigate the dynamics of persistent infection and the importance of large datasets, which can be achieved by combining data from several smaller studies in meta-analyses. Results of this study enhance current knowledge of the FMDV carrier state and may inform policy decisions regarding FMDV persistent infection.

59 BASIC BIOLOGICAL SCIENCES↗

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees Celsius or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics. The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees C or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

Capacitive response of biological membranes

We present a minimal model to analyze the capacitive response of a biological membrane subjected to a step voltage via blocking electrodes. Through a perturbative analysis of the underlying electrolyte transport equations, we show that the leading-order relaxation of the transmembrane potential is governed by a capacitive timescale, τ_{C}=λ_{D}L/D(2+Γδ^{M}/L/4+Γδ^{M}/λ_{D}), where λ_{D} is the Debye screening length, L is the electrolyte width, Γ is the ratio of the permittivity of the electrolyte to the membrane, δ^{M} is the membrane thickness, and D is the ionic diffusivity. This timescale is considerably shorter than the traditional RC timescale λ_{D}L/D for a bare electrolyte due to the membrane's low permittivity and finite thickness. Beyond the linear regime, however, salt diffusion in the bulk electrolyte drives a secondary, nonlinear relaxation process of the transmembrane potential over a longer timescale τ_{L}=L^{2}/4π^{2}D. A simple equivalent-circuit model accurately captures the linear behavior, and the perturbation expansion remains applicable across the entire range of observed physiological transmembrane potentials. Together, these findings underscore the importance of the faster capacitive timescale and nonlinear effects on the bulk diffusion timescale in determining transmembrane potential dynamics for a range of biological systems.

Farhadi, Jafar↗