Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic control”

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 127 records · Page 7

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗

Encoding beyond cosmological horizons in de Sitter JT gravity

Black hole event horizons and cosmological event horizons share many properties, making it natural to ask whether our recent advances in understanding black holes generalize to cosmology. To this end, we discuss a paradox that occurs if observers can access what lies beyond their cosmological horizon in the same way that they can access what lies beyond a black hole horizon. In particular, distinct observers with distinct horizons may encode the same portion of spacetime, violating the no-cloning theorem of quantum mechanics. This paradox is due precisely to the observer-dependence of the cosmological horizon — the sharpest difference from a black hole horizon — although we will argue that the gravity path integral avoids the paradox in controlled examples.

2D gravity↗

Kinetics of Shear-Induced Structural Ordering in Dense Colloids

The macroscopic rheological response of a colloidal solution is highly correlated with the local microscopic structure, as revealed by an in situ Rheo-SAXS experiment with a high temporal resolution. Oscillatory shear can induce a strain-controlled ordering-to-disorder transition, resulting in a shear-thickening process that is different from the normal shear-thickening behavior that is driven by hydrodynamics and particle friction. Here, we reveal that there is a complex time-dependent kinetics toward structural ordering under different applied strains. When the strain amplitude reaches a critical value that starts to induce disordering in the system, the pathway toward the dynamic equilibrium can also become highly non-monotonic. Within the same oscillatory cycle, there is a strong correlation of ordering with different phases of the oscillation, with the system oscillating between two dynamic metastable states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D Simulations Capture the Persistent Low-Mode Asymmetries Evident in Laser-Direct-Drive Implosions on OMEGA

Spherical implosions in Inertial Confinement Fusion (ICF) are inherently sensitive to perturbations that may arise from experimental constraints and errors. Control and mitigation of low-mode (long wavelengths) perturbations is a key milestone to improving implosion performances. Here, we present the first 3-D radiation-hydrodynamic simulations of directly driven ICF implosions with an inline package for polarized Crossed-Beam Energy Transfer (CBET). Simulations match bang times, yields (separately accounting for laser-induced high modes and fuel age), hot spot flow velocities and direction, for which polarized CBET contributes to the systematic flow orientation evident in the OMEGA implosion database. Current levels of beam mispointing, imbalance, target offset and asymmetry from polarized CBET degrade yields by more than 40%. The effectiveness of two mitigation strategies for low-modes is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinematic and dynamical origins of mean-p T fluctuations in heavy-ion collisions

Event-by-event fluctuations of the mean transverse momentum (mean-p T ) provide a sensitive probe of collective dynamics beyond single-particle spectra and anisotropic flow. We present a systematic study of mean-p T fluctuation observables using a Bayesian-calibrated multistage hydrodynamic framework, including quantitative comparisons to RHIC measurements and model-based investigations of beam-energy and kinematic-acceptance effects. The experimental definitions employed by the STAR and ALICE Collaborations are implemented explicitly and found to yield consistent results within controlled limits. We study the centrality and beam-energy dependence of the observable, its sensitivity to key soft-sector ingredients, and the impact of the kinematic p T acceptance. By introducing scaled-p T cuts, we demonstrate that a part of the apparent energy dependence arises from kinematic projection effects, while the remaining trends reflect genuine collective dynamics. Our results establish mean-p T fluctuations as a nontrivial and independent validation of calibrated hydrodynamic descriptions of the quark–gluon plasma.

Event-by-event correlations↗

Coupled Aerodynamic and Hydrodynamic Hybrid Simulation of Floating Offshore Wind Turbines

The development and innovation of floating offshore wind energy in the U.S. requires detailed high-fidelity observations and measurements of turbine and platform loading due to wind, waves, and currents. However, full-scale and quasi-full-scale experiments require significant financial and temporal investments for construction, experimental testing, and long-term field campaigns. To support the commercial advancement of the offshore wind energy industry, specialized wind tunnel and wave basin experimental facilities are critical to be able to test FOWT designs at small scale under controlled conditions prior to full-scale deployment. Oregon State University (OSU) is internationally known as a leader in water and energy research, development, and testing. The O.H. Hinsdale Wave Research Laboratory (HWRL) and the Wallace Energy Systems and Renewables Facility (WESRF) at OSU have extensive experience building, modeling, monitoring, controlling, and actuating scaled systems. Experiments on wave-structure interaction have been performed at the HWRL since its establishment in 1972. Studies have included the interaction of waves with coastal structures (breakwaters, seawalls, buildings, cylinders, bridges, fixed foundations of offshore wind turbines, etc.) and with floating structures (e.g., wave energy converters, maneuvering of vessels, etc.). Hinsdale is actively used by marine energy technology developers, both for private testing and OSU-collaborative research projects. However, despite the availability of several large-scale facilities for hydrodynamic testing (at OSU and elsewhere in the U.S.), existing experimental laboratories are generally limited in their ability to accurately generate combined wind and wave conditions. The simulation of both wind and waves in experimental testing is complicated due to a number of constraints, including: [i] incompatible similitude laws governing the wind and waves for scaled experiments, [ii] producing accurate wind over a large enough control volume via fans, and [iii] generating wind that reasonably represents the atmospheric boundary layer in existing wave basins/flumes. Hence, physical test data providing insight into the simultaneous wave- and wind-structure response of floating offshore wind components can be difficult to generate. Given the aforementioned challenges in classic hydrodynamic experiments, the motivation of this project is to establish a real-time hybrid simulation (RTHS) approach that can apply aero- and hydro-dynamic loading by augmenting wave-only experimental facilities with virtual aerodynamic forces through numerical models representing the remaining dynamic forces. RTHS is a physical-numerical approach that partitions a prototype system into physical and numerical sub-assemblies that interact with each other through actuators and sensors in real time. In coupling physical and numerical models, the hybrid simulation approach applied herein is ideal for problems with: (1) structures subjected to different scaling laws, such as floating offshore wind turbines subjected to combined aero/hydro-dynamic loading, (2) structures that are too large or complex to be tested entirely in a laboratory setting, such as deep-water mooring applications, and (3) component testing, where the behavior of a portion of the assembly is uncertain but still interacts with other portions of the structure, such as testing the fatigue life of turbine blades. Few U.S. experimental facilities are able to test simultaneous aero- and hydro-dynamic loading and none can accurately produce aero/hydro-dynamic response on scaled FOWT models due to conflicting similitude laws between the wind (commonly Reynolds) and the waves (commonly Froude). To aid in accelerating the development of the U.S. floating offshore industry, there is a significant need to develop a flexible, modular framework that can expand the capacities of existing wave-only laboratories. The project goal is to demonstrate a hydrodynamic real-time hybrid simulation (hydro-RTHS) framework that couples numerical wind and physical waves acting on a FOWT, thus representing simultaneous aero/hydro-dynamic loading. The FOWT is partitioned into a full-scale numerical sub-assembly associated with the aerodynamics and a model-scale physical sub-assembly associated with the hydrodynamics. The numerical-physical partition associated with hydro-RTHS mitigates scaling constraints by supplying different scaling laws to the physical and numerical sub-assemblies. Herein, length, force, and time are scaled and exchanged between the sub-assemblies using Froude scaling to represent the open-channel flow in the physical sub-assembly. Other similitude laws could also be utilized depending on the problem definition. It is envisioned that the ability to model FOWTs under waves and wind, with mitigation of similitude distortions, would result in reduced development costs (currently, FOWT concept development is performed with full-size pro- totypes at enormous expense and risk) and increase the reliability of the FOWT industry (since extreme wave and wind conditions and contingency events can be tested safely in a controlled environment).

16 TIDAL AND WAVE POWER↗

Few is different: deciphering many-body dynamics in mesoscopic quantum gases

Emergent macroscopic descriptions of matter, such as hydrodynamics, are central to our description of complex physical systems across a wide spectrum of energy scales. The conventional understanding of these many-body phenomena has recently been shaken by a number of experimental findings. Collective behavior of matter has been observed in mesoscopic systems, such as high-energy hadron–hadron collisions, or ultracold gases with only a few strongly interacting fermions. In such systems, the separation of scales between macroscopic and microscopic dynamics (at the heart of any effective theory) is inapplicable. To address the conceptual challenges that arise from these observations and explore the universality of emergent descriptions of matter, the EMMI Rapid Reaction Task Force was assembled. This document summarizes the RRTF discussions on recent theoretical and experimental advances in this rapidly developing field. Leveraging technological breakthroughs in the control of quantum systems, we can now quantitatively explore what it means for a system to exhibit behavior beyond the sum of its individual parts. In particular, the report highlights how the (in)applicability of hydrodynamics and other effective theories can be probed across three principal frontiers: the size frontier, the equilibrium frontier, and the interaction frontier.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamical correlation between ions in liquid seen in real space and time

Controlling ionic transport in liquids is anticipated to provide new scientific and technological opportunities, but it requires accurate knowledge of atomic-scale dynamics of ions beyond the hydrodynamic description. Atomic dynamics in liquids is characterized by strong and dynamical correlations among atoms, which render a conventional approach to describing the dynamics in reciprocal space challenging. We propose an alternative approach to describe liquid dynamics in real space and time using the time-delayed two-body correlation function, Van Hove correlation function, and discuss recent results on the local dynamics in water and electrolytes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗

Robust chain aggregation of low-entropy rigid ladder polymers in solution

Conjugated polymers have been widely investigated where ladder-type conjugated polymers receive more attention due to their rigid backbones and extraordinary properties. However, the understanding of how the rigid conformation of ladder polymers translates to material properties is still limited. Here, we systematically investigated the solution aggregation properties of a carbazole-derived conjugated ladder polymer (LP) and its analogous non-ladder control polymer (CP) via light scattering, neutron scattering, and UV-vis absorption spectroscopy characterization techniques, revealing a highly robust, temperature-insensitive aggregation behavior of the LP. The experimental findings were further validated by computational molecular dynamics simulations. Here we found that the peak positions and intensities of the UV spectra of the LP remained constant between 20 °C and 120 °C in chlorobenzene solution. The polymer also showed a stable hydrodynamic radius measured by dynamic light scattering from 20 °C to 70 °C in the chlorobenzene solution. Using small-angle neutron scattering, no Guinier region was reached in the measured q range down to 0.008 Å –1 , even at elevated temperature. In contrast, the non-ladder control polymer CP was fully soluble in the chlorobenzene solvent without the observation of any notable aggregates. The Brownian dynamics simulation showed that during polymer aggregation, the entropy change of the LP was significantly less negative than that of the non-ladder control polymer. These findings revealed the low entropy nature of rigid conjugated ladder polymers and the low entropy penalty for their aggregation, which is promising for highly robust intermolecular interactions at high temperatures. Such a unique thermodynamic feature of rigid ladder polymers can be leveraged in the design and application of next-generation electronic and optoelectronic devices that function under unconventional high temperature conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of kinetic mix enhancement in thin-shell OMEGA implosions

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D 3 He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fall line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D 3 He-p mix yield of 1.7 × 10 9 as compared to the experimental value of 9.3⁢ (±2.1) × 10 8 . This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D 3 He-p yields of 5×10 5 . This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗

Tunable electron–flexural phonon interaction in graphene heterostructures

Peculiar electron-phonon interaction characteristics underpin the ultrahigh mobility, electron hydrodynamics, superconductivity, and superfluidity observed in graphene heterostructures. Here, the Lorenz ratio (L) between the electronic thermal conductivity and the product of the electrical conductivity and temperature provides unique insight into electron-phonon interactions that is inaccessible to past graphene measurements. Here we show an unusual L peak in degenerate graphene near 60 Kelvin and decreased peak magnitude with increased mobility. When combined with ab initio calculations of the many-body electron-phonon self-energy and analytical models, this experimental observation reveals that broken reflection symmetry in graphene heterostructures can relax a restrictive selection rule to allow quasielastic electron coupling with an odd number of flexural phonons, contributing to the increase of L toward the Sommerfeld limit at an intermediate temperature sandwiched between the low-temperature hydrodynamic regime and the inelastic electron-phonon scattering regime above 120 Kelvin. In contrast to past practices of neglecting flexural phonon contributions to transport in two-dimensional materials, this work suggests that tunable electron-flexural phonon coupling can provide a handle to control quantum matter at the atomic scale, such as magic angle twisted bilayer graphene where low-energy excitations may mediate Cooper pairing of flat-band electrons.

36 MATERIALS SCIENCE↗

Ferromagnetic liquid droplets with adjustable magnetic properties

Significance Structured functional liquids combine mechanical versatility of fluids with solid-state properties, such as ferromagnetism, and offer a route to synthesize and control magnetic liquids for adaptive liquid robotics. Studies on these intriguing materials are only in their nascent state, and a profound understanding of the physical state is still lacking. We use hydrodynamics experiments to probe how the magnetization of ferromagnetic liquid droplets, governed by the assembly and jamming of magnetic nanoparticles at liquid–liquid interfaces, and their response to external stimuli can be tuned by chemical, structural, and magnetic means. Our results highlight the leading role of structural short-range order on magnetic properties, which provide a path toward nano-patterning structured liquids.

36 MATERIALS SCIENCE↗

Simulation of a High-Confinement Tidal Energy Array (Abstract)

When the projected area of a current turbine array occupies a substantial fraction of a tidal or river channel cross-sectional area, the confinement associated with the channel boundaries substantially increases turbine efficiency. This effect can be exploited to substantially reduce the levelized cost of energy for current turbine arrays. However, the presence of such an array will affect inflow velocity, water depth, and circulation in the surrounding area. This has implications for array performance and environmental impacts. This project combines existing numerical modeling capability from PNNL with an understanding of confined array performance from UW. The work will 1) Simulate tidal currents and water surface elevation in a regional circulation model that includes the Knik Arm site at appropriate resolution; 2) Incorporate a varying thrust coefficient into this model that reflects the outcomes of tidal hydrodynamic simulations of confined turbine arrays; 3) Use the model to determine the effects of array operation on water levels and inflow conditions as a function of different blockage ratios and control strategies; and 4) Identify the “far field” environmental implications of operating an array of this type. UW will use the outputs from this model to evaluate the implications of “feedback effects” and variable inflow conditions on levelized cost of energy estimates for a confinement exploiting array of cross-flow turbines.

16 TIDAL AND WAVE POWER↗

Interacting galaxies in the IllustrisTNG simulations - I: Triggered star formation in a cosmological context

ABSTRACT We use the IllustrisTNG cosmological hydrodynamical simulations to investigate how the specific star formation rates (sSFRs) of massive galaxies (M* > 1010 M⊙) depend on the distance to their closest companions. We estimate sSFR enhancements by comparing with control samples that are matched in redshift, stellar mass, local density, and isolation, and we restrict our analysis to pairs with stellar mass ratios of 0.1 to 10. At small separations (∼15 kpc), the mean sSFR is enhanced by a factor of 2.0 ± 0.1 in the flagship (110.7 Mpc)3 simulation (TNG100-1). Statistically significant enhancements extend out to 3D separations of 280 kpc in the (302.6 Mpc)3 simulation (TNG300-1). We find similar trends in the EAGLE and Illustris simulations, although their sSFR enhancements are lower than those in TNG100-1 by about a factor of two. Enhancements in IllustrisTNG galaxies are seen throughout the redshift range explored (0 ≤ $z$ < 1), with the strength of the enhancements decreasing with increasing redshift for galaxies with close companions. In order to more closely compare with observational results, we separately consider 2D projected distances between galaxies in IllustrisTNG. We detect significant sSFR enhancements out to projected separations of 260 kpc in TNG300-1, with projection effects diluting the size of the enhancements by about 20 per cent below 50 kpc. We find similar sSFR enhancements in TNG100-1 and Sloan Digital Sky Survey galaxies, with enhancements extending out to projected separations of about 150 kpc for star-forming galaxies at $z$ < 0.2. Finally, by summing over all separations, we estimate that the presence of closest companions boosts the average sSFR of massive galaxies in TNG100-1 by 14.5 per cent.

Patton, David R.↗

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗