Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pathway analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Projected income data under different shared socioeconomic pathways for Washington state

Abstract High-resolution income projections under different Shared Socioeconomic Pathways (SSPs) are essential for the climate change research communities to devise climate change adaptation and mitigation strategies. To generate income projections for Washington state, we obtain state-level GDP per capita projections and convert them into projected annual household income. The resulting state-level income projections are subsequently downscaled to the census block-level based on the Longitudinal Origin-Destination Employment Statistics (LODES) dataset. For accuracy assessment, we downscale historical income data from state- level to block- and block group-level and compare the downscaled results against the actual income data from LODES. County-level accuracy assessment is also conducted based on American Community Survey. The results demonstrate a good agreement (Average R 2 of 0.67, 0.8, and 0.99 for block-, block group-, and county-level, respectively) between the downscaled income data and the reference data, thereby validating the methodology employed. Our approach is applicable to other states for income projections, which can be utilized by a broader audience, including those involved in demographic analysis, economic research, and urban planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Electrolyte Reactivity on the MgV 2 O 4 Cathode Surface

Predictive understanding of the solvation-dependent reactivity and molecular interaction of electrolyte ions and solvent molecules on reactive electrodes has been a major challenge but is essential for addressing instabilities and surface passivation that occur at electrode-electrolyte interface (EEI) of multivalent Mg batteries. In this work, the isolated intrinsic reactivities of prominent chemical species present in magnesium bis(trifluoromethanesulfonimide) (Mg(TFSI) 2 ) in diglyme (G2) electrolytes, including ionic (TFSI - , [Mg(TFSI)] + , [Mg(TFSI):G2] + , [Mg(TFSI):2G2] + ) as well as neutral molecules (G2) on magnesium vanadate cathode (MgV 2 O 4 ) surface has been studied using a combination of first-principles calculations and multimodal analysis of well-defined cathode electrolyte interphase (CEI) layers. Here, our calculations show that non-solvated [Mg(TFSI)] + is the strongest adsorbing species on the MgV 2 O 4 surface compared to all other ions while fully solvated [Mg(TFSI):2G2] + are least favorable to decomposition. The cleavage of C-S bonds in TFSI - to form CF 3 - is predicted to be most desired pathway for all ionic species, which is followed by the cleavage of C-O bonds of G2 to yield CH 3 + or OCH 3 - species. The strong stabilization and electron transfer between ionic electrolyte species and MgV 2 O 4 is found to significantly favor these decomposition reactions on the surface compared to intrinsic gas phase dissociation. Experimentally, we used state-of-the-art ion soft landing to selectively deposit mass-selected TFSI - , [Mg(TFSI):G2] + and [Mg(TFSI):2G2] + on MgV 2 O 4 thin film to form well-defined electrolyte-MgV 2 O 4 interface. Analysis of soft-landed interphase using X-ray photoelectron, X-ray absorption near edge structure, electron energy-loss spectroscopies as well as transmission electron microscopy confirmed the presence of decomposition species (e.g., MgFx, carbonates) formed in the interfacial region and the higher amount of MgFx with [Mg(TFSI):G2] + , which corroborates the theoretical observation. Overall, we established the mechanistic pathway for the electrolyte-induced formation of passivating fluorides on MgV 2 O 4 cathode facilitated by the surface adsorption and charge transfer, which provided essential knowledge for rational design of stable electrolytes for multivalent cathodes.

Cathode-electrolyte interphase formation↗
DOE OSTI2024

Covalent labeling of the Arabidopsis plasma membrane H + ‐ ATPase reveals 3D conformational changes involving the C‐terminal regulatory domain

The plasma membrane proton pump is the primary energy transducing, electrogenic ion pump of the plasma membrane in plants and fungi. Compared to its fungal counterpart, the plant plasma membrane proton pump's regulatory C‐terminal domain (CTD) contains an additional regulatory segment that links multiple sensory pathways regulating plant cell length through phosphorylation and recruitment of regulatory 14‐3‐3 proteins. However, a complete structural model of a plant proton pump is lacking. Here, we performed covalent labeling with mass spectrometric analysis (CL‐MS) on the Arabidopsis pump AHA2 to identify potential interactions between the CTD and the catalytic domains. Our results suggest that autoinhibition in the plant enzyme is much more structurally complex than in the fungal enzyme.

Blackburn, Matthew R. [Department of Biochemistry ↗
DOE OSTI2020

Future projections and uncertainty assessment of precipitation extremes in the Korean peninsula from the CMIP5 ensemble

Abstract Projections of changes in extreme climate are sometimes predicted by multimodel ensemble methods that combine forecasts from individual simulation models using weighted averaging. One method to assign weight to each model is the Bayesian model averaging (BMA) in which posterior probability is used. For the cases of extreme climate, the generalized extreme value distribution (GEVD) is typically used. We applied the approach of GEV‐embedded BMA to a series of 35 years of the annual maximum daily precipitation data (both historical data and data gathered from simulation experiments for future periods) over the Korean peninsula as simulated by the models in the Coupled Model Intercomparison Project Phase 5 (CMIP5). Simulation data under two Representative Concentration Pathway (RCP) scenarios, namely RCP4.5 and RCP8.5, were used. Observed data and 17 CMIP5 models for 12 gird cells in Korea have been examined to predict future changes in precipitation extremes. A simple regional frequency analysis of pooling observations from three stations in each cell was employed to reduce the estimation variance and local fluctuations. A bias correction technique using the regression‐type transfer function was applied to these simulation data. Return levels spanning over 20 and 50 years, as well as the return periods relative to the reference years (1971–2005), were estimated for two future periods, namely Period 1 (2021–2050) and Period 2 (2066–2095). From these analyses, relative increase observed in the spatially averaged 20‐year (50‐year) return level was approximately 23% (16%) and 45% (36%) in the RCP4.5 and RCP8.5 experiments, respectively, by the end of the 21st century. We concluded that extreme rainfalls will likely occur two times and four times more frequently in the RCP4.5 and RCP8.5 scenarios, respectively, as compared to in the reference years by the end of the 21st century.

Lee, Youngsaeng↗
DOE OSTI2026

Upcycling mixed cathode materials to high-energy-density LiFe 0.75 Mn 0.25 PO 4

To address the demand for next-generation cathode materials with high energy density, upcycling LiFePO 4 into LiFe 0.75 Mn 0.25 PO 4 has attracted considerable attention. Nevertheless, existing strategies have yet to achieve both morphology and full elemental recovery under mild ambient conditions. Here, we report an upcycling route that can address this issue by combining leaching and a high-temperature treatment process. The upcycled LiFe 0.75 Mn 0.25 PO 4 exhibits enlarged lattice spacing and a high discharge plateau, and it delivers an energy density of 563.7 Wh/kg, 40.3 Wh/kg higher than recycled LiFePO 4 , which indicates the high value of the proposed upcycling strategy. At 1 C, LiFe 0.75 Mn 0.25 PO 4 also exhibits excellent cycling stability of 91% over 700 cycles. Techno-economic analysis also indicates impressive economic and environmental benefits, including 10.4% less raw materials usage and 12.2% less energy consumption and wastewater generation. This work demonstrates a scalable and economic upcycling strategy and provides a promising pathway for sustainable battery upcycling compatible with industrial conditions.

LMFP↗
DOE OSTI2026

Enhancing SWAT with mechanistic plant hydraulics: development and application in the Hanjiang River Basin

Plant transpiration plays a critical role in global water and energy cycles, requiring better process understanding as climate change intensifies drought stress and alters plant responses. Most hydrological models such as the widely-used SWAT lack representation of plant hydraulics, the mechanistic processes controlling plant water regulation and transpiration. Here, this study developed SWAT-PHS by integrating a plant hydraulics scheme (PHS) into SWAT hydrological model, enabling explicit simulation of root water uptake, sap flow, storage and transpiration at 30-minute timescales for watershed-scale modeling. In the Hanjiang River Basin, SWAT-PHS mitigated overestimation of runoff during the rainy season and underestimation during the dry season, reducing the overall simulation error by 29% across the entire simulation period. The model can simulate reasonable plant water dynamics, including diurnal transpiration patterns and drought responses showing declining transpiration flux, hydraulic buffering through stem water storage, and depth-dependent root water uptake strategies. Sensitivity analysis shows that SWAT-PHS captured mechanistic relationships between plant hydraulic traits and transpiration, with root distribution and stem capacitance positively affecting annual transpiration while vulnerability parameters showed negative effects. This work provides a pathway for improving hydrologic modeling and water resource management by better representing plant water regulation under climate change and expected intensifying water stress conditions.

China↗
DOE OSTI2021

Quantum Sensing of Electron Transfer Pathways in Natural Photosynthesis Using Time-Resolved High-Field Electron Paramagnetic Resonance/Electron–Nuclear Double Resonance Spectroscopy

Photosynthetic integral membrane reaction center (RC) proteins capture and convert sunlight into chemical energy via efficient charge separation achieved through a series of rapid, photoinitiated electron transfer steps. These fast electron transfers create an entangled spin qubit (radical) pair that contains detailed information about the weak magnetic interactions, structure, and dynamics of localized protein environments involved in charge separation events. In this work, we investigate how these entangled electron spin qubits interact with nuclear spins of the protein environment using the high spectral resolution of 130 GHz electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR). Spectroscopic interrogation enabled the observation and probing of protons located in the electron transfer pathway between the membrane-spanning electron pair P + Q A - (where P + is the primary donor, a special pair of bacteriochlorophylls, and Q A is the primary quinone acceptor) in the bacterial RC. Spectroscopic analysis reveals hydrogen-bonding interactions involved in regulating the route that light-induced electrons travel through the RC protein during charge separation.

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

Ethane Dehydrogenation on Single and Dual Centers of Ga-modified γ-Al 2 O 3

Density-functional theory calculations and microkinetic analysis are used to investigate the efficacy of Ga-modified γ-Al 2 O 3 (110) surfaces for the catalytic dehydrogenation of ethane and elucidate the synergy between Ga and Al sites. The model surfaces are modified by either Ga grafting or doping. We consider and analyze numerous active sites and rank them using microkinetic analysis. The kinetic parameters obtained from microkinetic modelling are compared with experimental values for ethane dehydrogenation over Ga 2 O 3 -Al 2 O 3 mixed oxides prepared by co-precipitation. The dominant reaction pathway proceeds via heterolytic C-H bond dissociation to a surface proton and a metal-carbanion intermediate that undergoes β-hydride elimination. We find that grafted Ga sites are catalytically inactive. In contrast, Ga-doped sites exhibit fivefold enhancement in catalytic activity when compared to the sites on pristine Al 2 O 3 , owed to the synergy between neighboring AlIII and GaIV sites. Furthermore, we model and investigate the effect of surface hydroxylation, demonstrate how surface water interferes with the aforementioned synergy between AlIII and GaIV sites and discuss the implications for the catalytic activity of the modified surfaces. Here, an increase in the partial pressure of H 2 O significantly increases the apparent activation energies of dehydrogenation and interestingly changes the most active site.

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

Life Cycle Assessment and Techno-Economic Assessment of Lithium Recovery from Geothermal Brine

Lithium-ion batteries (LIB) play an essential role in the electrification of the transportation sector, and battery demand for lithium compounds will see a significant increase in the coming decades. This has raised concerns on the supply of lithium, and as a result, technologies are being developed to process unconventional lithium sources. One promising technology is to extract lithium from geothermal brine using lithium-aluminum-layered double hydroxide chloride (LDH) sorbent and forward osmosis. A combined life cycle assessment (LCA) and techno-economic assessment (TEA) is conducted to evaluate the environmental and economic performance of this technology. It is assumed that the lithium extraction unit is an add-on to a 50 MW geothermal power plant located in California. The analysis is based on lab-scale experimental data and stoichiometry while considering the economy of scale for an industrial system. LCA results suggest that, compared with conventional LiOH and Li 2 CO 3 production pathways, the new technology achieves 1–95% reduction in environmental impacts. Even higher reduction can be achieved for LiOH produced via electrolysis. This add-on unit for lithium extraction could achieve a payback period of less than 1 year and reach net present values of $454M and $315M and internal rates of return of 792 and 1130% for LiOH and Li 2 CO 3 production pathways, respectively. The favorable environmental and economic performance suggests that LDH sorption coupled with forward osmosis has great potential to enable the domestic production of battery lithium compounds and that further development should be carried out.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Assessment of Accelerated Stress Testing Data for Silicon Photovoltaics Using Tensor Decomposition Methods

In this work, we examine the use of high-order tensor decompositions to analyze degradation pathways emerging from accelerated stress testing of silicon photovoltaic (PV) modules. Matrix-based decompositions are powerful tools for studying two-dimensional data arrays and form the foundation of a host of classical data analysis techniques. Tensors are high-order extrapolations of matrices that are able to account for more parameter dimensions, and a variety of tensor decomposition methods have been developed that similarly seek to extend insights from matrix decompositions to higher dimensions. Applying and interpreting tensor decomposition methods to sequences of PV module image data, we seek to uncover and isolate different degradation modes occurring from accelerated stress testing procedures. Further, we consider the contributions of different modes to PV module performance degradations.

data analysis↗
NASA NTRS1991

Studies of proton irradiated H2O + CO2 and H2O + CO ices and analysis of synthesized molecules

Infrared spectra of H2O + CO2 and H2O + CO ices before and after proton irradiation showed that a major reaction in both mixtures was the interconversion of CO2 yields CO. Radiation synthesized organic compounds such as carbonic acid were identified in the H2O + CO2 ice. Different chemical pathways dominate in the H2O + CO ice in which formaldehyde, methanol, ethanol, and methane were identified. Sublimed material was also analyzed using a mass spectrometer. Implications of these results are discussed in reference to comets.

Moore, M. H.↗
DOE OSTI2022

Pathways to a Sustainable Aviation Ecosystem: Accelerated Sustainable Aviation Fuel Properties Modeling

With world-class models and supercomputing capabilities combined with accurate property measurement, the National Renewable Energy Laboratory (NREL) can turn even the most complex sustainable aviation fuel (SAF) property characterization challenges into achievable successes. NREL provides holistic analysis and modeling to guide investments and accelerate the broad adoption of SAF.

aviation↗
DOE OSTI2022

Achieving 50% Energy Savings in Chicago Homes: A Case Study for Advancing Equity and Climate Goals

Since 2019, the National Renewable Energy Laboratory (NREL) and Elevate have collaborated to identify pathways to deep energy retrofits in Chicago's housing stock, document equity implications and co-benefits of this transition, and validate the findings by implementing retrofits in real Chicago homes. This document summarizes our analysis process to model advanced retrofit packages that lead to greater than 50% energy savings in Chicago homes. Based on these findings, we have also developed a roadmap with the City to guide implementation, and are deploying the recommended retrofit packages in real Chicago homes to realize these energy savings. This work was developed in collaboration with two key stakeholders - the City of Chicago and Commonwealth Edison (ComEd) - and funded by the U.S. Department of Energy (DOE). NREL's Residential Buildings team maintains the best-in-class ResStock™ energy model of the U.S. residential building stock. For this work, we calibrated ResStock to Chicago's unique local housing stock to accurately simulate energy use in Chicago homes both for current conditions and with various retrofit scenarios. We simulated a wide range of potential building retrofits covering all aspects of residential energy use and then grouped these into packages based on energy and utility bill savings. ResStock can model diverse building types and housing characteristics, so we're able to observe the range of outcomes that might occur when these upgrades are deployed across the entire housing stock. We can then estimate potential energy savings from an advanced retrofit program on Chicago's housing stock by comparing the modeled energy use before versus after a retrofit. This novel version of ResStock, calibrated to Chicago with data from Elevate, can help City officials, ComEd, and other partners plan for community-scale decarbonization via residential retrofits. Specifically, this work contributes the following project goals: Develop a building retrofit prioritization strategy for Chicago single-family and 2- to 4-unit buildings; Identify neighborhoods and home types that have the highest potential for savings from electrification; and Assess the impact of advanced building retrofits on energy use, utility bills, and CO2 emissions at the city and building level. Although this study is specific to Chicago, its methods and learnings are applicable across the United States. These findings are especially notable for heat pumps and electrification retrofits in cold climates.

building energy modeling↗
DOE OSTI2026

Comparative techno-economic and environmental analysis of using nuclear energy and fossil energy with carbon sequestration in U.S. Gulf Coast petroleum refineries

U.S. Gulf Coast refineries have processing capacity of approximately 9.4 million barrels of crude oil per day—about 50% of U.S. national capacity—while consuming natural gas, electricity, and hydrogen, resulting in approximately 100 million metric tonnes (MMT) of onsite CO₂ emissions in 2022. This study evaluates two alternative refinery energy supply pathways: nuclear energy (NE) and fossil energy-derived hydrogen with carbon capture and sequestration (FE-CCS; autothermal reforming), intermittent renewable sources are not considered as they cannot reliably meet continuous industrial heat and hydrogen demands. Using publicly available data for refineries, we conducted a bottom-up, facility-level assessment to estimate energy use by type, well-to-gate refinery emissions, and associated costs of integration. Compared to current refinery operations using natural gas energy supply and conventional hydrogen production (steam methane reforming without sequestration), NE and FE-CCS could produce, respectively, average emissions reductions of 37% and 42%, total abatement costs of $\$$52–$\$$221 and $\$$211–$\$$612/MT CO₂, and additional costs of about $\$$0.2–$\$$6 and $\$$3–$\$$11 per barrel of crude, respectively. Our analysis indicates that 25 out of 27 refineries have lower total additional cost ($\$$/bbl.) for the NE scenario than the FE-CCS scenario, making a strong case for NE integration with petroleum refineries. This work's main contribution is a detailed bottom-up refinery-level analysis method that can be utilized by the worldwide refining industry and stakeholders as they assess different technological options, along with their costs and environmental impacts for a specific refinery operation.

Carbon capture and sequestration↗
DOE OSTI2020

Techno-Economic Analysis and Global Warming Potential of a Novel Offshore Macroalgae Biorefinery

The success of a large scale macroalgae-based biorefinery is dependent on the demonstration of favorable system economics and environmental sustainability. This study uses detailed process modeling to quantify the mass and energy flows through the various unit operations required for a novel free-floating macroalgae biorefinery concept. The modular process model served as the foundation for the techno-economic and global warming potential analyses used to quantify the sustainability of the proposed concept. This work includes detailed techno-economic results for a complete macroalgae cultivation and conversion system with multiple hatchery configurations and several emerging technologies. System optimization was achieved through the evaluation of various technology options for each unit operation. Technologies considered include traditional twine and textile substrate hatchery configurations, drone assisted seeding and biomass transport, mechanized line seeding and harvesting, adhesive spore mixtures that simplify seeding operations and improve hatchery energetics, and hydrothermal liquefaction to produce upgradable biocrude. Outputs from the system include renewable diesel (R100), naphtha, biochar, nitrogen and phosphorus fertilizers, and aqueous/solid waste streams. Three different system pathways were explored, yielding a biomass production cost ranging from $210 to $565 per dry metric ton and a minimum fuel selling price from $1.35 to $2.91 per liter of gasoline equivalent. Stochastic manipulation of the process model and sensitivity analyses support these results. The global warming potential analysis shows net greenhouse gas emissions ranging from 14 to 29 gCO2-eq MJ-1, supported by stochastic and sensitivity analyses. The recommendations from this work highlight critical areas for research and development investment such that a sustainable macroalgae cultivation and conversion system can be realized.

NOMAD seaweed cultivation, techno-economic analysi↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2021

Maize specialized metabolome networks reveal organ-preferential mixed glycosides

Despite the scientific and economic importance of maize, little is known about its specialized metabolism. Here, five maize organs were profiled using different reversed-phase liquid chromatography-mass spectrometry methods. The resulting spectral metadata, combined with candidate substrate-product pair (CSPP) networks, allowed the structural characterization of 427 of the 5,420 profiled compounds, including phenylpropanoids, flavonoids, benzoxazinoids, and auxin-related compounds, among others. Only 75 of the 427 compounds were already described in maize. Analysis of the CSPP networks showed that phenylpropanoids are present in all organs, whereas other metabolic classes are rather organ-enriched. Frequently occurring CSPP mass differences often corresponded with glycosyl- and acyltransferase reactions. The interplay of glycosylations and acylations yields a wide variety of mixed glycosides, bearing substructures corresponding to the different biochemical classes. For example, in the tassel, many phenylpropanoid and flavonoid-bearing glycosides also contain auxin-derived moieties. The characterized compounds and mass differences are an important step forward in metabolic pathway discovery and systems biology research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗