Search NASA⌕ Search

SEARCH · Search NASA

Results for “biophysical simulations”

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

ITreeForeCast: An integrated modeling software to simulate tree level growth and forest carbon storage

Healthy trees in forest act as a natural carbon sink, capturing carbon. As they grow, they store carbon in their trunks, leaves and roots. Not all trees store carbon at the same rate, or in the same quantities, as it depends on a variety of biophysical and climatic factors. Furthermore, although carbon estimation in trees can be complex, the precision of estimates is tightly linked to trees growth, both in diameter and height. However, the simulation of carbon uptake by forest and forest growth has each been modeled separately, and independently at differing levels of detail and spatial resolution. In this paper, we introduce ITreeForeCast, a simulation model combining the two types of modeling on a unified platform, enabling the investigation of impacts of management strategies on carbon sequestration and wood products. ITreeForeCast is a user-extendable framework that offers new opportunities to model, simulate, and visualize the dynamics of individual trees in a forest, simulate management strategies over time, and carbon uptake.

09 - BIOMASS FUELS↗
DOE OSTI2022

Deforestation-induced climate change reduces carbon storage in remaining tropical forests

Abstract Biophysical effects from deforestation have the potential to amplify carbon losses but are often neglected in carbon accounting systems. Here we use both Earth system model simulations and satellite–derived estimates of aboveground biomass to assess losses of vegetation carbon caused by the influence of tropical deforestation on regional climate across different continents. In the Amazon, warming and drying arising from deforestation result in an additional 5.1 ± 3.7% loss of aboveground biomass. Biophysical effects also amplify carbon losses in the Congo (3.8 ± 2.5%) but do not lead to significant additional carbon losses in tropical Asia due to its high levels of annual mean precipitation. These findings indicate that tropical forests may be undervalued in carbon accounting systems that neglect climate feedbacks from surface biophysical changes and that the positive carbon–climate feedback from deforestation-driven climate change is higher than the feedback originating from fossil fuel emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulated Land Allocation, Nitrogen Use, and Nitrogen Loss in the Mississippi Atchafalaya River Basin for Various RFS2 (Renewable Fuel Standard) Policy Scenarios

Total nitrogen leaching rates were calculated over the Mississippi Atchafalaya River Basin (MARB) using an integrated economic-biophysical modeling approach. Land allocation for corn production and total nitrogen application rates were calculated for crop reporting districts using the Biofuel and Environmental Policy Analysis Model (BEPAM) for 5 RFS2 policy scenarios. These were used as input in the Integrated BIosphere Simulator-Agricultural Version (Agro-IBIS) and the Terrestrial Hydrologic Model with Biogeochemistry (THMB) to calculate the nitrogen loss. Land allocation and total nitrogen application simulations were simulated for the period 2016-2030 for 303 crop reporting districts (https://www.nass.usda.gov/Data_and_Statistics/County_Data_Files/Frequently_Asked_Questions/county_list.txt). The final 2030 values are reported here. Both are stored in csv files. Units for land allocation are million ha and nitrogen application are million kg. The nitrogen leaching rates were modeled with a spatial resolution of 5' x 5' using the North American Datum of 1983 projection and stored in NetCDF files. The 30-year average is calculated over the last 30 years of the 45 years being simulated. Leaching rates are calculated in kg-N/ha.

bioenergy crops↗
DOE OSTI2024

Asymmetric oligomerization state and sequence patterning can tune multiphase condensate miscibility

Endogenous biomolecular condensates, composed of a multitude of proteins and RNAs, can organize into multiphasic structures with compositionally distinct phases. This multiphasic organization is generally understood to be critical for facilitating their proper biological function. However, the biophysical principles driving multiphase formation are not completely understood. Here we use in vivo condensate reconstitution experiments and coarse-grained molecular simulations to investigate how oligomerization and sequence interactions modulate multiphase organization in biomolecular condensates. We demonstrate that increasing the oligomerization state of an intrinsically disordered protein results in enhanced immiscibility and multiphase formation. Interestingly, we find that oligomerization tunes the miscibility of intrinsically disordered proteins in an asymmetric manner, with the effect being more pronounced when the intrinsically disordered protein, exhibiting stronger homotypic interactions, is oligomerized. Our findings suggest that oligomerization is a flexible biophysical mechanism that cells can exploit to tune the internal organization of biomolecular condensates and their associated biological functions.

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

California annual grass phenology and allometry influence ecosystem dynamics and fire regime in a vegetation demography model

Grass-dominated ecosystems cover wide areas of the land surface yet have received far less attention from the Earth System Model (ESM) community. This limits model projections of ecosystem dynamics in response to global change and coupled vegetation–climate dynamics. We used the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a dynamic vegetation demography model, to determine ecosystem sensitivity to alternate, observed grass allometries and biophysical traits, and evaluated model performance in capturing California C 3 annual grasslands structure and fire regimes. Grass allometry, leaf physiology, plant phenology, and plant mortality all drove the seasonal variation in matter and energy exchange and fire dynamics in California annual grasslands. Allometry influenced grassland structure and function mainly through canopy architecture-mediated space and light competition instead of through carbon partitioning strategy. Regional variation in grassland annual burned area was driven by variation in ecosystem productivity. Our study advances the modeling of grassy ecosystems in ESMs by establishing the importance of grass allometry and plant phenology and mortality in driving C 3 annual grassland seasonal dynamics and fire regime. The calibrated annual grass allometry and biophysical traits presented can be applied in future studies to project climate–vegetation–fire feedbacks in annual grass-dominant ecosystems under global change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Examining metastatic behavior within 3D bioprinted vasculature for the validation of a 3D computational flow model

Understanding the dynamics of circulating tumor cell (CTC) behavior within the vasculature has remained an elusive goal in cancer biology. To elucidate the contribution of hydrodynamics in determining sites of CTC vascular colonization, the physical forces affecting these cells must be evaluated in a highly controlled manner. To this end, we have bioprinted endothelialized vascular beds and perfused these constructs with metastatic mammary gland cells under physiological flow rates. By pairing these in vitro devices with an advanced computational flow model, we found that the bioprinted analog was readily capable of evaluating the accuracy and integrated complexity of a computational flow model, while also highlighting the discrete contribution of hydrodynamics in vascular colonization. This intersection of these two technologies, bioprinting and computational simulation, is a key demonstration in the establishment of an experimentation pipeline for the understanding of complex biophysical events.

42 ENGINEERING↗
DOE OSTI2021

Simulating measurable ecosystem carbon and nitrogen dynamics with the mechanistically defined MEMS 2.0 model

Abstract. For decades, predominant soil biogeochemical models have used conceptual soil organic matter (SOM) pools and only simulated them to a shallow depth in soil. Efforts to overcome these limitations have prompted the development of the new generation SOM models, including MEMS 1.0, which represents measurable biophysical SOM fractions, over the entire root zone, and embodies recent understanding of the processes that govern SOM dynamics. Here we present the result of continued development of the MEMS model, version 2.0. MEMS 2.0 is a full ecosystem model with modules simulating plant growth with above- and belowground inputs, soil water and temperature by layer, decomposition of plant inputs and SOM, and mineralization and immobilization of nitrogen (N). The model simulates two commonly measured SOM pools – particulate and mineral-associated organic matter (POM and MAOM, respectively). We present results of calibration and validation of the model with several grassland sites in the US. MEMS 2.0 generally captured the soil carbon (C) stocks (R2 of 0.89 and 0.6 for calibration and validation, respectively) and their distributions between POM and MAOM throughout the entire soil profile. The simulated soil N matches measurements but with lower accuracy (R2 of 0.73 and 0.31 for calibration and validation of total N in SOM, respectively) than for soil C. Simulated soil water and temperature were compared with measurements, and the accuracy is comparable to the other commonly used models. The seasonal variation in gross primary production (GPP; R2 = 0.83), ecosystem respiration (ER; R2 = 0.89), net ecosystem exchange (NEE; R2 = 0.67), and evapotranspiration (ET; R2 = 0.71) was well captured by the model. We will further develop the model to represent forest and agricultural systems and improve it to incorporate new understanding of SOM decomposition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Long-range structural defects by pathogenic mutations in most severe glucose-6-phosphate dehydrogenase deficiency

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common blood disorder, presenting multiple symptoms, including hemolytic anemia. It affects 400 million people worldwide, with more than 160 single mutations reported in G6PD. The most severe mutations (about 70) are classified as class I, leading to more than 90% loss of activity of the wild-type G6PD. The crystal structure of G6PD reveals these mutations are located away from the active site, concentrating around the noncatalytic NADP + -binding site and the dimer interface. However, the molecular mechanisms of class I mutant dysfunction have remained elusive, hindering the development of efficient therapies. To resolve this, we performed integral structural characterization of five G6PD mutants, including four class I mutants, associated with the noncatalytic NADP + and dimerization, using crystallography, small-angle X-ray scattering (SAXS), cryogenic electron microscopy (cryo-EM), and biophysical analyses. Comparisons with the structure and properties of the wild-type enzyme, together with molecular dynamics simulations, bring forward a universal mechanism for this severe G6PD deficiency due to the class I mutations. We highlight the role of the noncatalytic NADP + -binding site that is crucial for stabilization and ordering two β-strands in the dimer interface, which together communicate these distant structural aberrations to the active site through a network of additional interactions. This understanding elucidates potential paths for drug development targeting G6PD deficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An agent-based modeling tool supporting bioenergy and bio-product community communication regarding cellulosic bioeconomy development

The emergence of the cellulosic bioeconomy requires not only adequate technological, economic, and policy advances, but also effective communication and coordination among the multiple stakeholders in the bioeconomy community. Aiming to facilitate stakeholder communication and collaboration for cellulosic bioeconomy growth, an agent-based model is developed as a community communication tool in this study. The tool simulates the behaviors and interactions of key stakeholders, as well as the responses and interactions between the stakeholders and biophysical, socioeconomic, and environmental conditions. The agent behavior rules are developed and validated based on behavior and economic theories and information attained from a stakeholder focus group study and a land use survey. A pilot study in a Central Illinois watershed shows stakeholder synergies, i.e., positive feedback loops across multiple stakeholders. The modeling tool's role in facilitating community communication is demonstrated via the focus group study. A comparison of policy scenarios suggests that the cellulosic bioeconomy development is fastest under a portfolio of policies combining the views of the multiple stakeholders. Critically, subsidy for small-scale bio-facilities is essential as it rapidly initializes an early market for cellulosic feedstocks and the adoption of Miscanthus. As a result, the communication tool is applied to bringing stakeholders and researchers together to develop shared visions, assess barriers, and explore solution portfolios to initialize and speed up the regional cellulosic bioeconomy.

09 BIOMASS FUELS↗
DOE OSTI2022

cabbi-bio/BECT-ABM

The emergence of the cellulosic bioeconomy requires not only adequate technological, economic, and policy advances, but also effective communication and coordination among the multiple stakeholders in the bioeconomy community. Aiming to facilitate stakeholder communication and collaboration for cellulosic bioeconomy growth, an agent-based model is developed as a community communication tool in this study. The tool simulates the behaviors and interactions of key stakeholders, as well as the responses and interactions between the stakeholders and biophysical, socioeconomic, and environmental conditions. The agent behavior rules are developed and validated based on behavior and economic theories and information attained from a stakeholder focus group study and a land use survey.

Yang, Pan↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Identifying sequence perturbations to an intrinsically disordered protein that determine its phase-separation behavior

Phase separation of intrinsically disordered proteins (IDPs) commonly underlies the formation of membraneless organelles, which compartmentalize molecules intracellularly in the absence of a lipid membrane. Identifying the protein sequence features responsible for IDP phase separation is critical for understanding physiological roles and pathological consequences of biomolecular condensation, as well as for harnessing phase separation for applications in bioinspired materials design. To expand our knowledge of sequence determinants of IDP phase separation, we characterized variants of the intrinsically disordered RGG domain from LAF-1, a model protein involved in phase separation and a key component of P granules. Based on a predictive coarse-grained IDP model, we identified a region of the RGG domain that has high contact probability and is highly conserved between species; deletion of this region significantly disrupts phase separation in vitro and in vivo. We determined the effects of charge patterning on phase behavior through sequence shuffling. We designed sequences with significantly increased phase separation propensity by shuffling the wild-type sequence, which contains well-mixed charged residues, to increase charge segregation. This result indicates the natural sequence is under negative selection to moderate this mode of interaction. We measured the contributions of tyrosine and arginine residues to phase separation experimentally through mutagenesis studies and computationally through direct interrogation of different modes of interaction using all-atom simulations. Finally, we show that despite these sequence perturbations, the RGG-derived condensates remain liquid-like. Together, these studies advance our fundamental understanding of key biophysical principles and sequence features important to phase separation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2020

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗