Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference models”

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,063 records · Page 59

Laser Disdrometer (LDIS) Instrument Handbook

The second-generation particle size velocity (Parsivel2, Figure 1) disdrometer, or laser disdrometer, referred to as “LDIS” within the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) User Facility, is a modern, laser-based disdrometer designed for comprehensive and reliable measurement of all types of precipitation. LDIS captures detailed information on the size and velocity of individual hydrometeors that fall to the ground, classifying them into a range of 32 categories. These raw measurements are processed by a fast signal processor (provided by the vendor) to calculate various parameters, including precipitation type, amount, intensity, kinetic energy, visibility in the precipitation, and equivalent radar reflectivity. These measurements have served as the primary data set and a valuable reference for a variety of studies, ranging from observational analyses of precipitation processes and climate model evaluations to monitoring or validating radar/satellite retrievals (e.g., Wang et al. 2018, Giangrande et al. 2019, Jackson et al. 2020).

47 OTHER INSTRUMENTATION↗

SCALE Lattice Physics Code Assessments of Accident Tolerant Fuel

This report highlights accident-tolerant fuel (ATF) code assessment activities performed under Project NRCHQ- 60-17-T-0017, Lattice Physics Enhancements and Assessment. ATF covers a broad range of advanced fuel and clad designs for light water reactors (LWRs) to enhance performance under several accident conditions. Several ATF concepts are anticipated to be deployed as lead test rods (LTRs) or lead test assemblies (LTAs) within the next five years. The purpose of this work is to assess the predictive capabilities of NRC neutronics codes that underpin various licensing calculations. ATF designs use different fuel and clad materials compared to standard UO 2 fuel and zirconium-alloy claddings (hereafter UO 2 -Zry). These new materials and geometrical designs need to be assessed to quantify the impact of nuclear data uncertainties on quantities of interest (QOIs) in licensing calculations and the impact of modeling approximations which may be valid for UO 2 -Zry but not for ATF. This report outlines a systematic approach for ATF neutronics code assessment which includes sensitivity and uncertainty (S/U) analysis of nuclear data, identification of experimental benchmark and gaps for code validation, investigation of modeling approximations, and code-to-code comparisons of calculated QOIs against high-fidelity reference continuous energy (CE) Monte Carlo (MC) calculations. This report focuses on the assessment of the SCALE/Polaris lattice physics code for reactor safety analysis. Polaris lattice physics calculations generate few group (FG) cross sections for PARCS full-core calculations. (Full-core analysis, spent fuel analysis, and severe accident analysis will be performed in future work.) The selected ATF concepts for this report include Cr 2 O 3 and Al 2 O 3 -Cr 2 O 3 -doped UO 2 fuel, U 3 Si 2 fuel, FeCrAl cladding, SiC cladding, and Cr-coated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Highly Resolved Projections of Passenger Electric Vehicle Charging Loads for the Contiguous United States: Results From and Methods Behind Bottom-Up Simulations of County-Specific Household Electric Vehicle Charging Load (Hourly 8760) Profiles Projected Through 2050 for Differentiated Household and Vehicle Types

This report documents enhancements made to the TEMPO (Transportation Energy & Mobility Pathway Options TM ) model to project spatially, demographically, and temporally resolved national-scale EV charging load profiles and describes three scenarios and corresponding datasets created for the NREL demand-side grid (dsgrid) project in support of bulk power systems modeling. In brief, TEMPO was enhanced to disaggregate national and annual energy demand projections into household and county-level projections of passenger electric vehicle (EV) hourly charging load profiles (8760 profiles), accounting for consumer, travel, and temperature variations that impact EV energy demand. In alignment with NREL's forward-looking grid modeling, three scenarios for EV adoption covering 2020-2050 were created-- Annual Energy Outlook (AEO) Reference Case, Electrification Futures Study (EFS) High Electrification, and All EV Sales by 2035 --and associated datasets have been included in the dsgrid platform for public use.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

VA Community Determinants of Health Data Curation Documentation FY26-Q2

The U.S. Department of Veterans Affairs (VA) places the health and well-being of our nation’s veterans as its top priority. VA is dedicated to offering timely access to high-quality, evidence-based mental health care that meets the needs of veterans and supports their reintegration into society. One of our core missions is to prevent suicide among veterans through innovative approaches and resources. With funding from the VA Office of Mental Health and Suicide Prevention (OMHSP), the Community Determinants of Health (EDH) Data project has developed innovative datasets associated with specific health outcomes, a methodology for transforming spatiotemporal data from one spatial reference (e.g., a 1 km grid) to another (e.g., U.S. Census Tracts), and capabilities for modeling health outcomes. These datasets represent an enhancement of the Agency for Healthcare Research and Quality (AHRQ), addressing key gaps by introducing finer spatial resolution (Census Tract) and additional geographical covariates into existing data. The curation and standardization of these datasets is a complex task since they often originate from various sources and are measured at different spatial and temporal resolutions. For example, U.S. Census data products typically use census blocks, block groups, or counties, while data like weather data are available on 1 km grids. Some economic data may only be available at the ZIP code level. In this context, ‘standardized’ means that all datasets share the same spatial extent (e.g., U.S. Census Tract and/or County), and ‘curated’ implies a repeatable process with data provenance and the use of appropriate methodologies for covariate conversion. The Community Determinants of Health datasets draw from multiple sources, resulting in variables with varying degrees of availability, patterns of missing data, and methodological considerations across different sources, geographies, and years.

99 GENERAL AND MISCELLANEOUS↗

A System for Standardizing and Combining U.S. Environmental Protection Agency Emissions and Waste Inventory Data

The U.S. Environmental Protection Agency (USEPA) provides databases that agglomerate data provided by companies or states reporting emissions, releases, wastes generated, and other activities to meet statutory requirements. These databases, often referred to as inventories, can be used for a wide variety of environmental reporting and modeling purposes to characterize conditions in the United States. Yet, users are often challenged to find, retrieve, and interpret these data due to the unique schemes employed for data management, which could result in erroneous estimations or double-counting of emissions. To address these challenges, a system called Standardized Emission and Waste Inventories (StEWI) has been created. The system consists of four python modules that provide rapid access to USEPA inventory data in standard formats and permit filtering and combination of these inventory data. When accessed through StEWI, reported emissions of carbon dioxide to air and ammonia to water are reduced approximately two- and four-fold, respectively, to avoid duplicate reporting. StEWI will greatly facilitate the use of USEPA inventory data in chemical release and exposure modeling and life cycle assessment tools, among other things. To date, StEWI has been used to build the recent USEEIO model and the baseline electricity life cycle inventory database for the Federal LCA Commons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Porcine Arterivirus, Influenza B, and Their Coinfection on Antiviral Response in the Porcine Lung

Interferon (IFN) cytokines induce an autonomous antiviral state in cells of the infected site to restrict virus spreading and critically regulate overall antiviral response. The antiviral state leads to host protection through expression of hundreds of IFN-stimulated genes that restrict viral infection through multiple mechanisms, for example, directly in viral genome degradation and indirectly through cellular metabolic inhibition. Young pigs were split into four treatment groups: control, porcine reproductive and respiratory syndrome virus (PRRSV, also known as porcine arterivirus) infected, influenza B virus (IBV) infected, and IBV/PRRSV coinfection. Lung tissue was collected at 3, 5, and 7 days post infection (dpi) for control, PRRSV and IBV/PRRSV coinfection, and at 3 and 5 dpi for IBV. Transcriptomic analysis, using usegalaxy.org tools, was performed against the S.scrofa 11.1 reference genome. Differentially expressed gene (DEG) analysis was carried out using DeSeq2 based on the model treatment + dpi + treatment:dpi + E. Downstream analysis examined the interaction of DEG at each dpi for over-enriched gene ontology (G.O.) terms and pathways. Comparisons of the infected groups vs. the controls yielded a total of (n = 1412) DEGs for the PRRSV group and (n = 1578) for the IBV/PRRSV group across all timepoints. The IBV group had (n = 64) total DEGs across 3 and 5 dpi. Expression data were considered statistically significant based on false discovery rate (FDR) ⫹ 0.1. Venn diagram comparisons of the DEGs across dpi showed that groups shared only 16 DEGs at 3 dpi, no DEGs were shared at 5 dpi, and for 7 dpi, only the PRRSV and IBV/PRRSV groups were compared and shared a total of 43 DEGs. Across the comparisons, differential expression was observed in antiviral genes such as IRF1, MX1, and OAS2. The IBV and IBV/PRRSV groups showed higher expression of antiviral genes at earlier dpi than the PRRSV group. Additionally, downregulated genes from the comparisons clustered around Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways effecting lung development and cellular integrity. Early expression of host IFN and antiviral genes may lead to viral RNA degradation, and assembly and transcription inhibition in the IBV infections. In comparison, expression of antiviral genes in the PRRSV group decreased across time. The decrease may explain why PRRSV infections persist, while IBV clears. Moreover, all infected groups showed prolonged upregulation in neutrophil degranulation pathway activity, possibly exacerbating symptomatic lung lesion pathology seen in these respiratory infections.

59 BASIC BIOLOGICAL SCIENCES↗

Stochastic Analysis for Long Term Capital Structures, Systems, and Components Refurbishment and Replacement

As commercial Nuclear Power Plants (NPPs) pursue extended plant operation in the form of Second License Renewal (SLR), opportunities exist for these plants to provide capital investments to ensure long-term safe and economic performance. At the current time, several utilities have announced an intention to pursue extended operation for one or more of their NPPs via SLR . The goal of this research is to develop a risk-informed approach to evaluate and prioritize plant capital investments made in preparation for, and during the period of, extended plant operations to support decisions for NPP operations. Since the capital investments are influenced by various factors, such as markets, safety and regulatory, the decision-making process of NPP operations should take into account relevant factors for balancing risks, costs and profits. The traditional method of capital budgeting is based on the priority list of candidate projects using economic measures such as benefit-investment ratio, net present value (NPV) and internal rate of return. In the literatures, the problem of capital budgeting or the variant can be represented by an appropriate knapsack problem. The knapsack approach to capital budgeting takes as input as investment, along with the cost and profit of each project. The objective of capital budgeting is to find the combination of the binary decisions for every investment such that the overall profit is as large as possible. The output is a collection of projects to be carried out, and we refer this selected collection of projects as a project portfolio. One limitation of traditional optimization models for capital budgeting is that they do not account for risk/uncertainty in profit and cost streams associated with individual projects, they do not account for risk in resource availability in future years [1,2,3]. Projects can incur cost over-runs, especially when projects are large, performed infrequently, and when there is risk regarding technical viability, external contractors, and/or suppliers of requisite parts and materials. Occasionally, projects are performed ahead of schedule and with cost savings. Planned budgets for capital improvements can be cut and key personnel may be lost. Or, there may be surprise windfalls in budgets for maintenance activities due to decreased costs for “unplanned” maintenance. In these cases, how should we resolve capital budgeting when we have risk forecasts for costs, profits and budgets? One approach we proposed in this summary is to re-solve the optimization models based on assumed statistical distributions of given parameters. If these distributions were not available, a two-stage stochastic optimization approach can be used to provide priority lists to decision-makers to support better risk-informed decisions [4, 5]. In this summary, we will only focus on the first approach.

42 ENGINEERING↗

Aspects of Crystal Growth and Atomic-Scale Characterization of U/Th Age-Dating Particles

Trace analysis of particulate material collected on environmental swipe samples continues to be one of the cornerstones of the IAEA’s process of verifying member state compliance as part of the Non-Proliferation Treaty (NPT). While complementary to traditional bulk analysis, the ability to measure individual particles can be far more impactful from a radiochronometry perspective. The 234U/230Th radiochronometer is being evaluated as a means of establishing model ages for discrete particles with sufficient confidence. To support the extension of current capabilities, reference materials with a suitable particle form factor and well-documented purification date are actively sought for development with large geometry secondary ion mass spectroscopy.

Nuclear Safeguards, IAEA, Particle Standards↗

Availability of Critical Benchmark Experiments for the Pebble Tanker Transportation Model for Nuclear Criticality Safety Validation of TRISO Pebbles

This study addresses the need for comprehensive investigations into TRi-structural ISOtropic (TRISO) fuel pebble transportation validation. In this work, an exploratory model, the pebble tanker(PT), was developed with the aim of facilitating the validation of nuclear criticality safety calculations in the context of industrial-scale transportation of TRISO fuel. The PT model was designed to investigate the availability and applicability of critical benchmark experiments crucial for assessing the transportation of these pebbles. This work incorporated sensitivity/uncertainty (S/U) similarity studies to quantify the applicability of critical benchmark experiments and to address nuclear data uncertainties in the context of TRISO transportation. Two container models were investigated: one for the Hermes-type pebble and one for the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering fuel, containment, and either water or air, to enable a focus on the underlying physics of applications involving TRISO fuel pebbles using the PT model. A crucial aspect under consideration was the capacity of the transport package to hold pebbles while ensuring subcriticality in the flooded state. An approach in the criticality validation process involves assessing the similarity between systems through an integral index parameter evaluation. This involves calculating a correlation coefficient (referred to as c k ) based on shared nuclear data–induced uncertainty between a benchmark experiment and the application of the PT model. To facilitate this analysis, the SCALE tools, particularly the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences, were employed for comprehensive studies in neutronics and S/U analysis. Our findings showed that there are sufficient critical experimental benchmarks to perform this validation of the PT model in the most reactive state, i.e. when the tanker is flooded. This paper provides valuable insights into validating a transport package for Generation IV TRISO fuel pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Output Capacitance Loss of GaN HEMTs in Steady-State Switching

The output capacitance (C OSS ) loss, a loss produced when the device’s output capacitor is charged and discharged, has become a concern for GaN high electron mobility transistors (HEMTs) in high-frequency applications. This work presents a new, easy-to-implement method for the C OSS loss characterization based on the unclamped inductive switching (UIS) setup. As compared to prior approaches, this method involves the device’s ON-state conduction and could measure the C OSS loss in a single pulse and the steady-state switching. The C OSS loss of three types of mainstream commercial GaN HEMTs is characterized, which exhibit some common dependencies including a non-monotonic relation with the dv/dt (or resonance frequency), a linear relation with the ON-state current, a power-law relation with the peak blocking voltage, and little temperature dependence. In addition, their C OSS losses all show minimal distinctions in a single pulse and the steady-state switching, despite the increased on-resistance in the steady-state switching. Furthermore, this suggests that the traps accounting for the C OSS loss possess different de-trapping time constants as compared to the traps governing the dynamic on-resistance. Finally, a unified model is established to describe the C OSS loss of all three types of GaN HEMTs. These results provide important references for the high-frequency application of GaN HEMTs and new insights into the physical origin of their C OSS loss.1

42 ENGINEERING↗

A Chromosome-Scale Genome Assembly of the Flax Rust Fungus Reveals the Two Unusually Large Effector Proteins, AvrM3 and AvrN

Rust fungi comprise thousands of species, many of which cause disease on important crop plants. The flax rust fungus Melampsora lini has been a model species for the genetic dissection of plant immunity since the 1940s; however, the highly fragmented and incomplete reference genome has so far hindered progress in effector gene discovery. Here, we generated a fully phased, chromosome-scale assembly of the two nuclear genomes of M. lini strain CH5, resolving an additional 320 Mbp of the sequence. The 482-Mbp dikaryotic genome is at least 79% repetitive, with a large proportion (approximately 40%) of the genome comprising young, highly similar transposable elements. The assembly resolves the known effector gene loci, some of which carry complex duplications that were collapsed in the previous assembly. Using a genetic map followed by manual correction of gene models, we identified the AvrM3 and AvrN genes, which encode unusually large fungal effector proteins and trigger defense responses when co-expressed with the corresponding resistance genes. We located the genes linked to the tetrapolar mating system on chromosomes 4 and 9, but in contrast to the cereal rusts that have one pheromone receptor gene per haplotype, in flax rust, three pheromone receptor genes were found, with two of them closely linked on one haplotype. Taken together, we show that a high-quality assembly is crucial for resolving complex gene loci, and given the increasing number of fungal effectors of large size, the commonly applied criterion for effector candidates of being small proteins needs to be reconsidered.

Melampsora↗

Estimating scalar turbulent fluxes with slow-response sensors in the stable atmospheric boundary layer

Conventional and recently developed approaches for estimating turbulent scalar fluxes under stable atmospheric conditions are evaluated, with a focus on gases for which fast sensors are not readily available. First, the relaxed eddy accumulation (REA) classical approach and a recently proposed mixing length parameterization, labeled A22, are tested against eddy-covariance computations. Using high-frequency measurements collected from two contrasting sites (the frozen tundra near Utqiaġvik, Alaska, and a sparsely vegetated grassland in Wendell, Idaho, during winter), it is shown that the REA and A22 models outperform the conventional Monin–Obukhov similarity theory (MOST) utilized widely to infer fluxes from mean gradients. Second, scenarios where slow trace gas sensors are the only viable option in field measurements are investigated using digital filtering applied to fast-response sensors to simulate their slow-response counterparts. With a filtered scalar signal, the observed filtered eddy-covariance fluxes are referred to here as large-eddy-covariance (LEC) fluxes. A virtual eddy accumulation (VEA) approach, akin to the REA model but not requiring a mechanical apparatus to separate the gas flows, is also formulated and tested. A22 outperforms VEA and LEC in predicting the observed unfiltered (total) eddy-covariance (EC) fluxes; however, VEA can still capture the LEC fluxes well. This finding motivates the introduction of a sensor response time correction into the VEA formulation to offset the effect of sensor filtering on the underestimated net averaged fluxes. The only needed parameter for this correction is the mean velocity at the instrument height, a surrogate of the advective timescale. The VEA approach is very suitable and simple to use with gas sensors of intermediate speed (∼ 0.5 to 1 Hz) and with conventional open- or closed-path setups.

58 GEOSCIENCES↗

Updates to the Hanford Soil Inventory Model (SIM Version 2) for FY 2023

The purpose of this environmental calculation file (ECF) is to document the updates made to the Soil Inventory Model version 2.1 (SIM-v2.1) (ECF-HANFORD-21-0073, Updates to the Hanford Soil Inventory Model (SIM Version 2) for FY 2021) in FY 2023. These updates to SIM-v2.1 are referred to as the SIM-v2.2 version. The following updates have been made: • Include inventory estimate for a new analyte • hexavalent chromium (Cr(VI)), separate from total chromium inventory • Revise the inventory discharged to the 216-U-10 and 216-T-4 Pond systems based on partitioning of waste streams discharges over time and space among influent ditches and ponds • Enhance the preprocessing and postprocessing of the input and output files Note that the methodology for estimating Cr(VI) inventories is based on ECF-200W-23-0040, Recommendations for Updating Liquid Discharged Inventory and Transport Modeling Parameters for Cumulative Impacts Evaluation of Hexavalent Chromium in the 200 West Area while the revision of inventory discharged to the 216-U-10 and 216-T-4 Pond systems is based on ECF-HANFORD-19-0032, Distribution of Infiltration in the 216-U-10, 216-B-3 Pond, and 216-T-4 Pond Systems 1944-1997.

54 ENVIRONMENTAL SCIENCES↗

Multimodal CO2 Transportation Cost Model

This model provides a cost estimate for transporting CO2 via truck or rail in the United States using commercially available equipment and technologies. The model includes an analysis of direct and indirect CO2 emissions to determine costs per net tonne of CO2 transported. Publicly available data and methods published in the peer-reviewed literature are used to the extent possible; references are available at the bottom of the "Calculations" sheet. Upstream (i.e., liquefaction, buffer storage) and downstream (i.e., buffer storage, reconditioning) activities are included in the model of emissions and costs. All capital and operating expenditures are estimated in the "Calculations" sheet. The cost of financing the project is determined in the "FINEX" sheet. All user inputs are done via drop down menus in on the "User Interface" sheet. Summary results are also provided on the "User Interface" sheet.

Myers, CoreyA↗

A Latent-Variable Formulation of the Poisson Canonical Polyadic Tensor Model: Maximum Likelihood Estimation and Fisher Information

We establish parameter inference for the Poisson canonical polyadic (PCP) tensor model through a latent-variable formulation. Our approach exploits the observation that any random PCP tensor can be derived by marginalizing an unobservable random tensor of one dimension larger. The loglikelihood of this larger dimensional tensor, referred to as the “complete” loglikelihood, is comprised of multiple rank one PCP loglikelihoods. Using this methodology, we first derive maximum likelihood estimators for the PCP model and demonstrate that several existing algorithms for fitting non-negative matrix and tensor factorizations are Expectation-Maximization algorithms. Next, we derive the observed and expected Fisher information matrices for the PCP model. The Fisher information provides us crucial insights into the well-posedness of the tensor model, such as the role that tensor rank plays in identifiability and indeterminacy. For the special case of rank one PCP models, we demonstrate that these results are greatly simplified.

97 MATHEMATICS AND COMPUTING↗

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗

DECOVALEX-2023: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX initiative is an international research collaboration (www.decovalex.org), initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of the complex perturbations and coupled processes in geologic repositories and how these impact long-term performance predictions. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2023, started in early 2020 and ended in late 2023. More than fifty research teams associated with 17 international DECOVALEX partner organizations participated in the comparative evaluation of eight modeling tasks covering a wide range of spatial and temporal scales, geological formations, and coupled processes. This Virtual Special Issue on DECOVALEX-2023 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes. While primarily focused on radioactive waste, much of the work included here has wider application to many geoengineering topics.

Coupled processes↗