Search NASA⌕ Search

SEARCH · Search NASA

Results for “System Code”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Pseudouridinylation of mRNA coding sequences alters translation

Chemical modifications of RNAs have long been established as key modulators of nonprotein-coding RNA structure and function in cells. There is a growing appreciation that messenger RNA (mRNA) sequences responsible for directing protein synthesis can also be posttranscriptionally modified. The enzymatic incorporation of mRNA modifications has many potential outcomes, including changing mRNA stability, protein recruitment, and translation. In this work, we tested how one of the most common modifications present in mRNA coding regions, pseudouridine (Ψ), impacts protein synthesis using a fully reconstituted bacterial translation system and human cells. Our work reveals that replacing a single uridine nucleotide with Ψ in an mRNA codon impedes amino acid addition and EF-Tu GTPase activation. A crystal structure of the Thermus thermophilus 70S ribosome with a tRNA Phe bound to a ΨUU codon in the A site supports these findings. We also find that the presence of Ψ can promote the low-level synthesis of multiple peptide products from a single mRNA sequence in the reconstituted translation system as well as human cells, and increases the rate of near-cognate Val-tRNA Val reacting on a ΨUU codon. The vast majority of Ψ moieties in mRNAs are found in coding regions, and our study suggests that one consequence of the ribosome encountering Ψ can be to modestly alter both translation speed and mRNA decoding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Modern gyrokinetic particle-in-cell simulation of fusion plasmas on top supercomputers

The Gyrokinetic Toroidal Code at Princeton (GTC-P) is a highly scalable and portable particle-in-cell (PIC) code. It solves the 5D Vlasov-Poisson equation featuring efficient utilization of modern parallel computer architectures at the petascale and beyond. Motivated by the goal of developing a modern code capable of dealing with the physics challenge of increasing problem size with sufficient resolution, new thread-level optimizations have been introduced as well as a key additional domain decomposition. GTC-P's multiple levels of parallelism, including inter-node 2D domain decomposition and particle decomposition, as well as intra-node shared memory partition and vectorization have enabled pushing the scalability of the PIC method to extreme computational scales. In this paper, we describe the methods developed to build a highly parallelized PIC code across a broad range of supercomputer designs. This particularly includes implementations on heterogeneous systems using NVIDIA GPU accelerators and Intel Xeon Phi (MIC) co-processors and performance comparisons with state-of-the-art homogeneous HPC systems such as Blue Gene/Q. New discovery science capabilities in the magnetic fusion energy application domain are enabled, including investigations of Ion-Temperature-Gradient (ITG) driven turbulence simulations with unprecedented spatial resolution and long temporal duration. Performance studies with realistic fusion experimental parameters are carried out on multiple supercomputing systems spanning a wide range of cache capacities, cache-sharing configurations, memory bandwidth, interconnects and network topologies. Furthermore, these performance comparisons using a realistic discovery-science-capable domain application code provide valuable insights on optimization techniques across one of the broadest sets of current high-end computing platforms worldwide.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Optimizing the orbital occupation in the multiple minima problem of magnetic materials from the metaheuristic firefly algorithm

Here, we present the use and implementation of the firefly algorithm to help in scanning the multiple metastable minima of orbital occupations in density functional theory (DFT) plus Hubbard U correction and to identify the ground state occupations in strongly correlated materials. We show the application of this implementation with the Abinit code on KCoF 3 and UO 2 crystals, which are typical d and f electron systems with numerous occupation minima. We demonstrate the validity and performance of the method by comparing with previous methodologies. The method is general and can be applied to any code using constrained occupation matrices.

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

Coupling laser physics to radiation-hydrodynamics

In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Symbiotic bacteriophages exhibit multiple adaptive strategies in activated sludge flocs and contribute to floc stability

Despite the importance of phages for the dynamics and functions of microbial communities, it remains largely unexplored how symbiotic phages adapt in activated sludge systems and influence microbial aggregate stability, which is critical for clarification and reliable performance. Here, based on 12,127 phage contigs recovered from 12 activated sludge microbiomes, the symbiotic phages exhibited broad host ranges and infected dominant prokaryotes in activated sludge. Moreover, these phage communities exhibited high lysogenicity (49% to 66% lysogenic phages) and harbored diverse auxiliary metabolic genes (AMGs) that could enhance microbial aggregate stability, such as genes coding for lipopolysaccharide, peptidoglycan, and extracellular polysaccharide biosynthesis. Additionally, prokaryotic antiviral systems were widely distributed in dominant prokaryotes, particularly the Restriction-Modification (RM) and CRISPR-Cas systems, which could mitigate potential harmful phage infections. Overall, this study reveals the contemporary beneficial relationship between phages and prokaryotic hosts in activated sludge microbiomes and the potential benefits of symbiotic phages on floc stability.

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

Extension of OpenMC for Fixed Source Transmutation Calculations

This report documents work performed under a Strategic Partnership Project between Argonne National Laboratory (ANL) and the United Kingdom Atomic Energy Authority (UKAEA). The overall goal of this project is to extend OpenMC [1], a community-developed Monte Carlo particle transport code, to be able to perform fixed-source transmutation calculations. In fission and fusion reactors, the high flux of energetic neutrons causes materials within the reactor to “transmute,” or undergo a nuclear reaction that results in the addition/removal of neutrons and protons from the nucleus of an atom. If subjected to these reactions for long enough, the overall composition and physical properties of the material itself begin to change as a result of transmutation. Such a feature is vital for predicting the decrease in tritium production rate within a breeder blanket during the lifetime of a fusion reactor. OpenMC is capable of simulating neutron transport in fission/fusion systems, thereby allowing it to estimate the flux that causes transmutation. It is also capable of solving the transmutation equations, which determine how the composition of a material changes over time due to neutron irradiation and radioactive decay. However, solving the transmutation equations was previously only possible when the source of neutrons came from a fission system. In a fusion system, the source of neutrons is typically determined by a separate code and then given as an input to the particle transport simulation. This is known as a fixed source calculation. Through this project, we have extended OpenMC to solve the transmutation equations for a fixed source calculation. Evaluating the change in material compositions due to transmutation and its effect on physical properties is of key importance to a range of engineering analyses for fission and fusion systems. For example, in a fusion reactor, estimating the dose rate at different physical locations resulting from irradiated materials in the reactor allows designers to ensure that workers are not exposed to doses beyond applicable regulations. In order to properly dispose of irradiated materials, designers also need to estimate the radiotoxicity, which again relies on knowledge of the material composition at some future time. The specific tasks for this project that were agreed to between ANL and UKAEA were as follows: 1. Make changes and additions in the openmc.deplete and related modules in OpenMC to support transmutation calculations following a fixed source transport simulation. 2. Make necessary changes to OpenMC to model transmutation due to an arbitrary set of reactions needed for fusion applications. Use this new capability to generate a depletion chain file based on the TENDL nuclear data library. 3. Improve the openmc.deplete module in OpenMC to keep track of gases produced as a by-product of nuclear reactions during transmutation calculations. 4. Validate the new capabilities by carrying out fixed-source transmutation calculations on a suitable benchmark problem using OpenMC and a comparable Monte Carlo neutron transport code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Expanding the genetic toolkit: adenine and cytosine base editors for gene disruption in Aspergillus niger

Despite revolutionizing fungal genetic engineering, conventional CRISPR/Cas9-mediated knockouts rely on DNA double-strand breaks (DSBs), which can cause unwanted insertions and deletions, chromosomal abnormalities, and cytotoxicity. Base editors such as adenine base editors (ABEs), which convert A‧T to G‧C, and cytosine base editors (CBEs), which convert C‧G to T‧A, offer a safer alternative by enabling predictable, target-specific single-nucleotide changes without introducing DSBs. To overcome the limitations of traditional genome editing in filamentous fungi, we developed efficient base-editing systems in Aspergillus niger . For the first time, we constructed an ABE in A. niger , achieving up to 80% editing efficiency and inducing predictable A-to-G mutations at the intended intron sites, disrupting gene function through mRNA mis-splicing. We also developed a highly efficient CBE system, capable of introducing premature stop codons with 50–100% efficiency. To broaden the editing scope, we implemented a Cas9-NG variant recognizing a relaxed PAM sequence requiring only a single guanine (G), enabling editing at start codons and splice sites. Leveraging this expanded scope, we established gene disruption approaches by targeting start codons via ABE-mediated A-to-G conversions (ATG-to-GTG and ATG-to-ACG) and CBE-mediated C-to-T conversion (ATG-to-ATA). Additionally, our base-editing systems enable multiplex gRNA delivery and marker-free editing of multiple genes. Collectively, the scope-expanding strategies increase the number of genes targetable for disruption by base-editing in A. niger by 26.3% and enable near-complete coverage of 96% of the coding genes. Overall, this work demonstrates the potential of ABE and CBE systems as versatile, efficient, and safer alternatives to DSBs-based gene disruption in filamentous fungi.

Aspergillus↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Decayheatml

This code is designed to predict and analyze the decay heat generated in molten salt reactors (MSRs) using a hybrid approach that combines machine learning and segmented polynomial fitting. The accurate prediction of decay heat is essential for reactor safety and the optimization of spent fuel storage. The code operates through several key components: 1) Data Architecture: It incorporates a modular data architecture that handles various MSR-specific operational parameters such as power density, humidity content, and air ingress. These parameters are sampled using Sobol sequences to ensure comprehensive coverage of operational uncertainties. 2) Machine Learning Framework: The code employs a diverse set of machine learning models, including polynomial regression, decision trees, random forests, gradient boosting, support vector regression, k-nearest neighbors, multi-layer perceptrons, and symbolic regression. These models are trained to predict decay heat over a wide temporal range, from immediate shutdown up to 10,000 years. 3) Region-Optimized Training: The temporal domain is divided into multiple regions, each modeled separately to capture distinct decay heat characteristics across different time scales. This approach significantly improves the accuracy and interpretability of predictions. 4) Segmented Polynomial Interpretation (SPI): The SPI method translates machine learning predictions into piecewise polynomial equations. These equations are physically interpretable and can be directly integrated into existing engineering workflows and safety analyses. 5) Front-End Interfaces: The code includes both a Jupyter notebook interface for research development and a Streamlit web application for operational deployment. These interfaces allow users to interactively explore decay heat predictions, adjust operational parameters, and visualize results in real-time. 6) Applications: The framework supports various applications, including safety system validation and spent fuel container optimization. It enables real-time evaluation of worst-case decay heat scenarios, informing the design of passive safety systems and optimizing container designs for long-term storage. Overall, this code provides a robust, accurate, and user-friendly tool for predicting decay heat in MSRs, enhancing reactor safety, and optimizing spent fuel management.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2022

Implementation of real‐time TDDFT for periodic systems in the open‐source PySCF software package

Abstract We present a new implementation of real‐time time‐dependent density functional theory (RT‐TDDFT) for calculating excited‐state dynamics of periodic systems in the open‐source Python‐based PySCF software package. Our implementation uses Gaussian basis functions in a velocity gauge formalism and can be applied to periodic surfaces, condensed‐phase, and molecular systems. As representative benchmark applications, we present optical absorption calculations of various molecular and bulk systems and a real‐time simulation of field‐induced dynamics of a (ZnO) 4 molecular cluster on a periodic graphene sheet. We present representative calculations on optical response of solids to infinitesimal external fields as well as real‐time charge‐transfer dynamics induced by strong pulsed laser fields. Due to the widespread use of the Python language, our RT‐TDDFT implementation can be easily modified and provides a new capability in the PySCF code for real‐time excited‐state calculations of chemical and material systems.

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

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Stochastic GW -GPU: Rapid Quasi-Particle Energies for Molecules beyond 10,000 Atoms

StochasticGW is a code for computing accurate quasi-particle (QP) energies of molecules and material systems in the GW approximation. StochasticGW utilizes the stochastic Resolution of the Identity (sROI) technique to enable a massively parallel implementation with computational costs that scale semilinearly with system size, allowing the method to access systems with tens of thousands of electrons. Here, we introduce a new implementation, StochasticGW-GPU, for which the main bottleneck steps have been ported to GPUs and give substantial performance improvements over previous versions of the code. We showcase the new code by computing band gaps of hydrogenated silicon clusters (Si x H y ) containing up to 10,001 atoms and 35,144 electrons, and we obtain individual QP energies with a statistical precision of better than ±0.03 eV with times-to-solution of less than 1 h.

Thomas, Phillip S. [Lawrence Berkeley National Lab↗
DOE OSTI2018

A novel risk classification system for 30-day mortality in children undergoing surgery

A simple, objective and accurate way of grouping children undergoing surgery into clinically relevant risk groups is needed. The purpose of this study, is to develop and validate a preoperative risk classification system for postsurgical 30-day mortality for children undergoing a wide variety of operations. The National Surgical Quality Improvement Project-Pediatric participant use file data for calendar years 2012–2014 was analyzed to determine preoperative variables most associated with death within 30 days of operation (D30). Risk groups were created using classification tree analysis based on these preoperative variables. The resulting risk groups were validated using 2015 data, and applied to neonates and higher risk CPT codes to determine validity in high-risk subpopulations. A five-level risk classification was found to be most accurate. The preoperative need for ventilation, oxygen support, inotropic support, sepsis, the need for emergent surgery and a do not resuscitate order defined non-overlapping groups with observed rates of D30 that vary from 0.075% (Very Low Risk) to 38.6% (Very High Risk). When CPT codes where death was never observed are eliminated or when the system is applied to neonates, the groupings remained predictive of death in an ordinal manner.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

The Mycobacterium tuberculosis protein O -phosphorylation landscape

Bacterial phosphosignaling has long been synonymous with the histidine kinases of the two component systems (TCSs), but many bacteria, including Mycobacterium tuberculosis (Mtb), also code for Ser/Thr protein kinases (STPKs). STPKs are the main phosphosignaling enzymes in eukaryotes, but the full extent of phosphorylation on protein Ser/Thr and Tyr (O-phosphorylation) in bacteria is unclear. Here, we explored the global signaling capacity of the STPKs in Mtb. We generated a panel of STPK loss- and gain-of-function strains and measured the resulting O-phosphorylation and transcriptional changes. This deep phosphoproteome shows that O-phosphorylation in Mtb is a vastly underexplored protein modification that affects >70% of the proteome. The substrate-kinase interactions of STPKs show an extensive interface with the transcriptional machinery, resulting in STPK regulation of gene expression of over 30% of Mtb genes. Finally, Mtb O-phosphorylation gives rise to an expansive, distributed, and cooperative network of a complexity that has previously only been associated with eukaryotic phosphosignaling networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Sparsified Time-dependent PDEs FNO (STFNO) v1.0.0

STFNO (Sparsified Time-dependent PDEs FNO code) is an extension of the popular Fourier Neural Operator (FNO) architecture to the solution of coupled systems of time-dependent partial differential equations. STFNO leverages the sparsified dependencies on the field quantities based on the semi-discretiezed form of the PDEs, enabling significant reduction in the number of model parameters. STFNO has been extensively tested on two fusion simulation codes, NIMROD and GTC, and can be easily tailored to other systems of PDEs.

Rahman, Mustafa [Lawrence Berkeley National Labora↗
DOE OSTI2024

Dialogue Between Lighting and HVAC Systems: Improving building system integration

Lighting systems have long been capable of sensing when someone enters or exits a room and using that knowledge to turn lights on or off. More recently, connected lighting systems with sensors integrated into every luminaire have become broadly available, facilitating highly granular occupancy detection. Similarly, HVAC systems have long been able to use an understanding of building occupancy to adjust temperature setpoints and reduce energy use without significant impacts to occupant comfort. Energy codes (e.g., ASHRAE/IES Standard 90.1, IECC, Title 24) now require “occupied standby HVAC control,” whereby systems adjust both temperature and ventilation setpoints in zones that are determined to be unoccupied during normal occupancy hours. Here, this article discusses current issues that stymie the integration of Lighting and HVAC systems, and DOE activities focused on addressing them.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗