Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Gaussian process based surrogate approach for the optimization of cylindrical targets

Simulating direct-drive inertial confinement experiments presents significant computational challenges, both due to the complexity of the codes required for such simulations and the substantial computational expense associated with target design studies. Machine learning models, and in particular, surrogate models, offer a solution by replacing simulation results with a simplified approximation. In this study, we apply surrogate modeling and optimization techniques that are well established in the existing literature to one-dimensional simulation data of a new cylindrical target design containing deuterium–tritium fuel. These models predict yields without the need for expensive simulations. We find that Bayesian optimization with Gaussian process surrogates enhances sampling efficiency in low-dimensional design spaces but becomes less efficient as dimensionality increases. Nonetheless, optimization routines within two-dimensional and five-dimensional design spaces can identify designs that maximize yield, while also aligning with established physical intuition. Optimization routines, which ignore constraints on hydrodynamic instability growth, are shown to lead to unstable designs in 2D, resulting in yield loss. However, routines that utilize 1D simulations and impose constraints on the in-flight aspect ratio converge on novel cylindrical target designs that are stable against hydrodynamic instability growth in 2D and achieve high yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Star Formation in Splash Bridges

Splash bridges are created from the direct collision of two gas-rich disk galaxies. These direct collisions can eject gas masses on the order of 10 10 M ⊙ stripped from the stellar disks of each galaxy. The Taffy Galaxy system (UGC 1294/5) is a prototypical example of a splash bridge system. CO observations of the Taffy revealed that its splash bridge contains a mass of H 2 equal to the Milky Way's H 2 mass. However, the little visible star formation occurring within the bridge highlights the need for models of direct gas-rich disk collisions. The Arp 194 system displays what may be another splash bridge resulting from the collision between two disk galaxies. The region between the two stellar disks contains bright clumps of active star formation. We aim to better understand the conditions for star formation in splash bridges by employing a Jeans criterion to determine where gravitational instabilities occur in the shocked and cooling gas of gas-rich disk collisions. The splash bridge results are obtained from our previous work using a sticky particle code and post-processed. We find that the inclination between the gas disks and the collision velocity with which the gas collides strongly affects the fraction of gas that will become gravitationally unstable. Low inclinations between gas disks produce starbursts whereas high inclinations result in steady-rate star formation. The offset of the gas disks at impact will determine how many gas elements directly collide but does not strongly affect the resulting star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hydro-scaling of direct-drive cylindrical implosions at the OMEGA and the National Ignition Facility

Deceleration-phase Rayleigh–Taylor instability (RTI) growth during inertial confinement fusion capsule implosions significantly affects the performance as it mixes cold ablator material into the fuel. Precise measurements of such instability growth are essential for both validating the existing simulation codes and improving our predictive capability. RTI measurements on the inner surface of a spherical shell are limited and are often inferred indirectly at limited convergence. In contrast, cylindrical implosions allow for direct diagnostic access to the converging interface by imaging down the cylinder axis while retaining the effects of convergence. We have performed direct-drive cylindrical implosion experiments at both the OMEGA and the NIF laser facilities using scaled targets. RTI growth is demonstrated to be scale-invariant between the cylindrical targets at OMEGA and similar targets at the NIF, which are scaled up by a factor of three in the radial dimension. Single-mode (m = 0) instability growth factors of ~14 are measured at a convergence ratio (CR) ~ 2.5 with nearly identical mode growth at both scales. The measurements are in agreement with xRAGE radiation-hydrodynamics simulations. In addition, we have developed the Bayesian-inference-engine method to account for the variations in the target alignment, magnification, and the parallax effect in the measurement, allowing a more precise comparison between the experimental data and the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

A Study of Cavitation-Ignition Bubble Combustion

We present the results of an experimental and computational study of the physics and chemistry of cavitation-ignition bubble combustion (CIBC), a process that occurs when combustible gaseous mixtures are ignited by the high temperatures found inside a rapidly collapsing bubble. The CIBC process was modeled using a time-dependent compressible fluid-dynamics code that includes finite-rate chemistry. The model predicts that gas-phase reactions within the bubble produce CO and other gaseous by-products of combustion. In addition, heat and mechanical energy release through a bubble volume-expansion phase are also predicted by the model. We experimentally demonstrate the CIBC process using an ultrasonically excited cavitation flow reactor with various hydrocarbon-air mixtures in liquid water. Low concentrations (< 160 ppm) of carbon monoxide (CO) emissions from the ultrasonic reactor were measured, and found to be proportional to the acoustic excitation power. The results of the model were consistent with the measured experimental results. Based on the experimental findings, the computational model, and previous reports of the "micro-diesel effect" in industrial hydraulic systems, we conclude that CIBC is indeed possible and exists in ultrasonically- and hydrodynamically-induced cavitation. Finally, estimates of the utility of CIBC process as a means of powering an idealized heat engine are also presented.

Nguyen, Quang-Viet↗
DOE OSTI2021

The AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo

In this work, we present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 10 12 M⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, and Gizmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Benchmarking numerical methods for impact and cratering applications

Large scale computational models are important for studying impact cratering events that are prevalent both on Earth and, more broadly, in this solar system. To address these problems, models must reliably account for both large length scales (e.g., kilometers) and relatively long time scales (hundreds of seconds). This work benchmarks two such approaches, a more traditional hydrodynamics approach and a finite-discrete element method (FDEM), for impact cratering applications. Both 2D and 3D results are discussed for two different impact velocities, 5 km/s and 20 km/s, striking normal to the target and, for 3D simulations, 45° from vertical. In addition, comparisons to previously published data are presented. Finally, differences in how these methods model damage are discussed. Ultimately, both approaches show successful modeling of several different impact scenarios.

58 GEOSCIENCES↗
DOE OSTI2024

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Multi-variable, Multi-value Hydrodynamic Modeling Approach to Support Marine Energy and Coastal Resilience Applications

Funded by the Marine Energy Seedling Laboratory Program, the study reported on here aimed to develop a multi-variable, multi-value hydrodynamic modeling approach based on the Finite-Volume Community Ocean Model (FVCOM) to support a wide range of marine renewable energy and coastal resilience applications. We migrated the marine hydrokinetic energy module to simulate tidal energy extraction by tidal turbines under wave-current interactions. We explored the particle tracking module in FVCOM and applied it to Sequim Bay, Washington, to illustrate the surface current field and mimic the pollutant transport near the surface. A tracer-water age module was added to FVCOM to examine water exchange through the entrance. Lastly, a user-specified vertical coordinate option was added to FVCOM to improve simulation of the effect of floating turbines on surface currents and scalar transport. The above-described model development has been tested using a simplified model domain and demonstrated in real-world coastal systems. The outcome of this project will provide a useful modeling code and approach for marine renewable energy and coastal resilience research and applications.

16 TIDAL AND WAVE POWER↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Thermocatalytic Heat Pipes for Geothermal Resource Recovery

Heat pipes are an important technology that allow orders of magnitude faster heat transfer than simple conduction. However, operating principles in heat pipes place fundamental bounds on their performance (critical heat flux and efficiency). Conventional heat pipe functionality is inherently tied to vaporization and condensation of the working fluid charged in the heat pipe. These fluids each have different operating temperature ranges based on the capillary, entrainment, sonic, and boiling limits of the heat pipe design. These limits, typically the capillary limit, dictate the maximum heat flux a heat pipe can carry, and most importantly for geothermal systems, the distance over which the pipes can operate (100 to 200 m maximum under optimum conditions). A thermocatalytic heat pipe breaks the inherent limitations of phase change thermo- and hydrodynamics and can transform heat pipe technology as a potentially more efficient means of extracting heat from a geothermal resource. The thermocatalytic heat pipe uses a working fluid to transport both sensible and chemical heat. An endothermic chemical reaction at depth removes heat from the reservoir and produces reactive intermediates, which are transported to the surface and used to run a reverse exothermic reaction that releases heat for use in power generation or other useful purposes. This technology offers two distinct advantages over conventional geothermal heat recovery technologies: (1) lower heat loss to the rock outside of the geothermal reservoir, and (2) higher heat transfer rates to the well field within the geothermal reservoir. Both advantages offer opportunity to reduce risks and lower costs of geothermal energy recovery. In this report, we discuss an initial effort to assess the efficacy and limitations of this technology for extracting heat from both porous/permeable and nominally impermeable geothermal reservoirs. Numerical simulation capabilities of the STOMP-GT code were enhanced to enable simulations of thermochemical heat pipes traversing geothermal reservoirs. An array of potential thermochemical reaction systems was evaluated and screened. Of these, an ethanol dehydration reaction was most promising in the vapor-liquid reaction set. A solid-phase dehydration reaction (CuSO4·5H2O) showed the highest reaction enthalpy per unit volume but would require development of a nonaqueous carrier fluid to implement it in a heat pipe. Subsurface reservoir simulations predicted long-term performance of the heat pipes for each geothermal reservoir type. The performance of U-shaped wells and coaxial wells was evaluated for a suite of reactions for both hydrothermal and hot dry rock reservoirs and was compared with a baseline case of simply pumping water through the wells. The heat pipe technology was additionally evaluated for an enhanced geothermal system (EGS) with an injection borehole, production borehole, and intervening hydraulically conductive fracture. All reservoir types showed significant improvement in heat recovered over a 20-year operating period ranging from a 1.8X increase for the hot dry rock case to more than 2.5X more energy recovered for the EGS case.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
NASA NTRS1986

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗
NASA NTRS2011

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗