Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Towards Robust Calibration of the AWSD Reactive Burn Model

Calibration of a reactive burn model for detonation of high explosive is an important step towards predictive hy drodynamic simulations of detonation. A typical calibration consists of varying model parameters (e.g., rate constants, activation energies) until results of hydrodynamic simulations match the experimental data for a certain set of ex periments. Hydrodynamic simulations of the dependence of steady detonation velocity on the radius of a cylindrical high-explosive charge - often used in such calibrations - can be computationally expensive. In this work, we propose a method where such expensive simulations are performed infrequently, and only to parameterize and refine a surrogate model for the dependence of the detonation velocity on calibrated parameters. The method is developed, implemented and applied to an example problem - calibration of the AWSD reactive burn model for important high explosive PBX 9502. Two different flavors of the surrogate model are investigated, and the calibration is performed successfully.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Adjoint DSMC for nonlinear spatially-homogeneous Boltzmann equation with a general collision model

We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in Caflisch et al., which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2022

A correction procedure for self-induced velocity of a finite-sized particle in two-way coupled Euler–Lagrange simulations

The importance of incorporating a correction to undo the self-induced perturbation velocity of a particle, when its size becomes comparable to the Eulerian grid, in a two-way coupled Euler–Lagrange (EL) simulation is now well appreciated. The present work improves upon the prior correction procedures in a few important ways. First, the past correction procedures have been scalar-based with the assumption that the quasi-steady force is the source of self-induced velocity perturbation. Here we generalize to a vector correction procedure and thereby the directions of feedback force and relative velocity can be different. This allows the correction procedure to be used even in the presence of added-mass, history, and lift forces. Second, the effect of a nearby wall has been systematically included in the correction procedure. The correction procedure depends on fundamental Oseen solutions of streamwise and transverse regularized feedback forces. We present a Fourier transform-based analytical approach to obtaining these regularized Oseen solutions. We also present a step-by-step numerical procedure for obtaining the Oseen solutions in any EL code. With the analytical or numerical Oseen functions, the correction procedure can be easily implemented in any EL code. Iterations are required in solving the implicit correction equations and it is demonstrated that the correction procedure converges rapidly within three or four iterations. In conclusion, a simple empirical approach is also presented to account for unsteady effects in the correction procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SUGAR: An improved empirical model of Type Ia supernovae based on spectral features

Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Influence of Cs Loading on Pt/m-ZrO 2 Water–Gas Shift Catalysts

Certain alkali metals (Na, K) at targeted loadings have been shown in recent decades to significantly promote the LT-WGS reaction. This occurs at alkali doping levels where a redshift in the C-H band of formate occurs, indicating electronic weakening of the bond. The C-H bond breaking of formate is the proposed rate-limiting step of the formate associative mechanism, lending support to the occurrence of this mechanism in H2-rich environments of the LT-WGS stage of fuel processors. Continuing in this vein of research, 2%Pt/m-ZrO 2 was promoted with various levels of Cs in order to explore its influence on the rate of formate intermediate decomposition, as well as that of LT-WGS in a fixed bed reactor. In situ DRIFTS experiments revealed that Cs promoter loadings of 3.87% to 7.22% resulted in significant acceleration of the forward formate decomposition in steam at 130 °C. Of all of the alkali metals tested to date, the redshift in the formate ν(CH) band with the incorporation of Cs was the greatest. XANES difference experiments at the Pt L 2 and L 3 edges indicated that the electronic effect was not likely due to an enrichment of electronic density on Pt. CO 2 TPD experiments revealed that, unlike Na and K promoters, Cs behaves more like Rb in that the decomposition of the second intermediate in LT-WGS, carbonate species, is hindered due to (1) increased basicity of Cs, (2) the tendency of Cs to cover Pt sites that facilitate CO 2 decomposition, and (3) the tendency of Cs to increase Pt particle size as shown by EXAFS results, resulting in fewer Pt sites that facilitate CO2 decomposition. As such, the LT-WGS rate was hindered overall and the rate-limiting step shifted to carbonate decomposition (CO 2 removal). Like its Rb counterpart, low levels of added Cs (e.g., 0.72%Cs) were found to improve the stability of the catalyst relative to the unpromoted catalyst; the stability comparison was made at similar CO conversion level as well as similar space velocity.

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

Plasma surface ionization wave interactions with single channels

Abstract The study of plasma surface ionization waves (SIWs) in recent years has primarily focused on planar surfaces and periodic two dimensional structures. In application, substrates are likely to have non-planar morphology such as cracks, pores, and steps. Additionally, targets for the applications of medicine or catalysis may have targets with heterogeneous composition. This classification of targets are brought under the umbrella of complex interfaces. In this work, plasma SIWs were incident on a complex target consisting of a single channel cut into glass slides. The SIW velocities for the in-channel portion of the wave and radially propagating portion of the wave were tracked. It was found that surface wave velocities are not significantly affected by channel geometry, but primarily increase with pulse energy. A third propagation direction for the SIW is characterized in the azimuthal direction relative to the radial portion of the wave. Channel geometry is found to greatly effect the area treated by the plasma but not the propagation velocity of the surface wave. Surface wave morphology and the impact on application is also discussed. A simple model was introduced to understand the mechanisms behind SIW escape from a channel. It was found that the ratio of pulse energy to a geometry dependent minimum energy can predict the escape angle of a SIW from single channels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗
DOE OSTI2025

Catalytic water splitting kinetics and optimization

This project was motivated primarily by the opportunity of turning photocatalysis into a cost efficient hydrogen generation method for domestic energy and chemical production, and by the need to provide detailed mechanisms of photocatalytic reactions of broad scientific and technological interest. Specifically, we focused on the development of a new technique for analysis of (photo)catalytic reaction mechanisms to guide photocatalyst design. Our approach was based on applying velocity map imaging (VMI) to surface catalyzed reactions in order to identify photocatalytic reaction products, intermediates, and relate these to photocatalyst composition, dimensions, size and other characteristics. By coupling VMI with other spectroscopic and microstructural characterization, we aimed to reveal the features that limit performance of photocatalysis. As a first step to achieve this goal, we constructed a custom VMI microscope specifically designed for investigation of surface photocatalytic reactions and validated it’s design using a photodissociatinon of MnBr(CO) 5 adsorbed on a metallic surface. Using both 230 nm and 266 nm pulsed laser excitations, we observed that there are likely three main pathways for photodissociation of Mn(CO)5Br: from the gas phase, the surface, and ejected Mn complex from the surface increasing the gas phase concentration.

08 HYDROGEN↗
DOE OSTI2021

Overview of the TCV tokamak program: Scientific progress and facility upgrades

The TCV tokamak is augmenting its unique historical capabilities (strong shaping, strong electron heating) with ion heating, additional electron heating compatible with high densities, and variable divertor geometry, in a multifaceted upgrade program designed to broaden its operational range without sacrificing its fundamental flexibility. The TCV program is rooted in a three-pronged approach aimed at ITER support, explorations towards DEMO, and fundamental research. A 1 MW, tangential neutral beam injector (NBI) was recently installed and promptly extended the TCV parameter range, with record ion temperatures and toroidal rotation velocities and measurable neutral-beam current drive. ITER-relevant scenario development has received particular attention, with strategies aimed at maximizing performance through optimized discharge trajectories to avoid MHD instabilities, such as peeling-ballooning and neoclassical tearing modes. Experiments on exhaust physics have focused particularly on detachment, a necessary step to a DEMO reactor, in a comprehensive set of conventional and advanced divertor concepts. The specific theoretical prediction of an enhanced radiation region between the two X-points in the low-field-side snowflake-minus configuration was experimentally confirmed. Fundamental investigations of the power decay length in the scrape-off layer (SOL) are progressing rapidly, again in widely varying configurations and in both D and He plasmas; in particular, the double decay length in L-mode limited plasmas was found to be replaced by a single length at high SOL resistivity. Experiments on disruption mitigation by massive gas injection and electron-cyclotron resonance heating (ECRH) have begun in earnest, in parallel with studies of runaway electron generation and control, in both stable and disruptive conditions; a quiescent runaway beam carrying the entire electrical current appears to develop in some cases. Developments in plasma control have benefited from progress in individual controller design and have evolved steadily towards controller integration, mostly within an environment supervised by a tokamak profile control simulator. TCV has demonstrated effective wall conditioning with ECRH in He in support of the preparations for JT-60SA operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Sap Velocity Data for Urban Trees in Chicago, Illinois (2024-2025)

This dataset contains uncorrected sap velocity measurements using the heat ratio method (HRM) collected using ICT International SFM1x sensors at five urban sites in Chicago, Illinois, as part of the DOE CROCUS project. The data includes continuous monitoring of sap velocity from various tree species, including Maples (Acer spp.): Sugar Maple (Acer saccharum), Silver Maple (Acer saccharinum), and Red Maple (Acer rubrum); Oaks (Quercus spp.): Swamp White Oak (Quercus bicolor); American Elm (Ulmus americana); Honey Locust (Gleditsia triacanthos); Cottonwood (Populus deltoides); and Tree of Heaven (Ailanthus altissima) across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), and West Woodlawn "Blacks in Green" (BIG). These include both street trees and those in urban park locations. Measurements were collected at 15-20 minute intervals, depending on the sensor, and transmitted via Long Range Wide Area Network (LoRaWAN) protocols. The wireless data was collected by Sage Network (https://sagecontinuum.org/) nodes. The dataset includes sensor ID, Global Positioning System (GPS) coordinates, tree species (common and scientific names), tree identification number, diameter at breast height (DBH in cm), uncorrected sap velocity measurements (cm/hr) from both inner and outer probes, and Sage Node identifiers so the data can be mapped to related variables such as air quality and wind speed that were collected on the Sage nodes. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 3-bit binary system indicating physical range violations (< -10 or > 60 cm/hr), step spikes (absolute difference > 36 cm/hr), and stuck sensor conditions (> 10 consecutive identical values). These are raw data, not corrected for wood anatomy or species-specific characteristics. Data is provided in comma separated (CSV) format. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Multi-Function Research LoRaWAN (MFR) Nodes. DOIs for the supporting data are provided as part of this data package.

Chicago↗
DOE OSTI2025

Regularizing the linearly extrapolated BDF2 scheme for incompressible flows with time relaxation

This paper presents a highly-efficient finite element scheme for the time relaxation model (TRM). The efficiency is achieved through the second-order BDF2 time-stepping scheme with linear extrapolation (BDF2LE). The accuracy of the scheme is also greatly enhanced through the use of the divergence-free Scott-Vogeulis finite elements, and van Cittert approximate deconvolution. A complete finite element analysis is provided, which includes rigorous proofs for the stability, well-possessedness, and convergence of both velocity and pressure solutions. Furthermore, we also demonstrate that the inclusion of the linear time relaxation term preserves the long-time stability of the unregularized BDF2LE scheme. Finally, numerical experiments are presented that demonstrate the added stability and accuracy that time relaxation can provide.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

TEAMER – Enhanced Flow Measurement for Aquantis Tidal Turbine Test

The AQ10 is a floating, two-bladed, passive yawing tidal turbine developed by Aquantis that has a 10-meterrotor diameter, 160 kW rating, and employs reliable off-the-shelf powertrain and power conversion hardware. Aquantis is planning on-water turbine power performance and loads (blade loading and thrust)testing, where the turbine will be pushed through still water up to 4 knots and placed in a ‘station keeping’ tow in a tidal race up to its rated speed. On-water testing will be conducted using vessels and floating platforms on the sea surface to improve ease of testing and reduce disturbance to the environment. In this TEAMER project, Pacific Northwest National Laboratory (PNNL) will conduct water velocity and turbulence measurements in front of the turbine during on-water testing using acoustic Doppler instrumentation. By measuring both the tidal current flowing past the turbine and the resulting electrical power output, test results will provide a power curve (power vs flow speed) for the turbine up to rated power. Measurements of turbulence and velocity shear in front of the rotor will also provide information to assess the structural response of the rotor blades. With this analysis, Aquantis can use the performance and loads data to validate Tidal Bladed and OpenFAST simulations of the measured operating conditions. Measuring the power performance of a prototype turbine is a valuable step to improving device development and conducting a complete power performance assessment to IEC/TS 62600-200 standards in the future.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Integrated phononic waveguide on thin-film lithium niobate on diamond

Here, we demonstrate wavelength-scale phononic waveguides formed by transfer-printed thin-film lithium niobate (LN) on bulk diamond, a material stack that combines the strong piezoelectricity of LN with the high acoustic velocity and color-center compatibility of diamond. We characterize a delay line based on a 100 μm long phononic waveguide at room and cryogenic temperatures. The total insertion loss through the device at 4 K is −5.8 dB, corresponding to a $>$ 50 % transducer efficiency, at a frequency of 2.8 GHz. Our work represents a step toward phonon-mediated hybrid quantum systems consisting of strain-sensitive color centers in diamond.

Malik, Sultan [Stanford Univ., CA (United States)]↗