Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortex dynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing ferroelectric domain structures and their switching dynamics in SrBi2Ta2O9 by in-situ electric biasing in transmission electron microscopy

Abstract Lead-free SrBi 2 Ta 2 O 9 (SBT) has been a promising ferroelectric material for various applications such as electronics and data storage due to its outstanding ferroelectric properties including high fatigue endurance and low leakage current. However, the atomic-scale domain structure and switching dynamics of ferroelectric SBT remain elusive. This study reveals that spontaneous polarization arises from canted bismuth-cation displacements, forming 90° and Ising-type 180° domain walls. Interestingly, topological pairs of ferroelectric vortex and antivortex connect ferroelectric boundaries where three domain walls converge. In situ electrical biasing transmission electron microscopy (TEM) reveals the dominance of 180° switching over 90°, where oxygen octahedral connectivity is protected by ferroelastic energy in the 90° domain wall. Consequently, all 180° domain walls and (anti)vortices are removed, leaving only 90° domain walls in the electrically poled states. Chemical deterioration along domain walls highlights vulnerability of SBT to ferroelectric fatigue. This study provides insight into crucial aspects for practical applications of SBT.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Vortical interactions in turbulent thermoacoustic systems

This study examines the dynamics of vortical interactions and their implications for mitigating thermoacoustic instability in a turbulent combustor. The regions of intense vortical interactions are identified as vortical communities in the network space of weighted directed vortical networks constructed from two-dimensional experimental velocity data. One can expect vortical interactions in the combustor to be strongest near the moment of vortex shedding, as the shed vortices gradually weaken due to dissipation while convecting downstream. However, we show that, during the state of thermoacoustic instability, there is a non-trivial consistent phase lag of approximately 52° between the shedding of the coherent structures from the backward-facing step and the time instant when the vortical interactions attain their local maximum value. We explain this phase lag by investigating the correlation between acoustic pressure fluctuations, spatio-temporal dynamics of coherent structures and vortical interactions in the reaction field of the combustor. We also show the aperiodic variation of vortical interactions during the states of combustion noise and aperiodic epochs of intermittency. Furthermore, the spatio-temporal evolution of pairs of vortical communities with the maximum inter-community interactions provides insight into explaining the critical regions detected in the reaction field during the states of intermittency and thermoacoustic instability, also identified in previous studies. As a result, we further show that the most efficient suppression of thermoacoustic instability via air microjet injection is achieved when steady air jets are introduced to disrupt the maximum inter-community interactions present during the state of thermoacoustic instability.

Sahay, Ankit [Indian Institute of Technology Madra↗
DOE OSTI2026

Emergence of vorticity and viscous stress in finite-scale quantum hydrodynamics

The Madelung equations offer a hydrodynamic description of quantum systems, from single particles to quantum fluids. In this formulation, the probability density is mapped onto the fluid density and the phase is treated as a scalar potential generating the velocity field. As examples of potential flows, quantum fluids described in this way are inherently irrotational, but quantum vortices may arise at discrete points where the phase is undefined. In this paper, starting from this irrotational description of a quantum fluid, a coarse-graining procedure is applied to arrive at a macroscopic description of the quantum fluid in terms of a hierarchy of moments in which the role of velocity is played by a Favre average of the microscopic velocity field. This hierarchy is truncated using an explicit closure derived from an expansion in a finite length scale. The resulting coarse-grained fields are shown to allow for finite vorticity at any point in the fluid. Additionally, it is shown that this vorticity obeys a similar equation to the vorticity equation in classical hydrodynamics and includes a vortex-stretching term. The particular closure employed here also gives rise to a novel stress term in the fluid equations, which in the appropriate limit appears analogous to an artificial viscous stress from computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗
DOE OSTI2024

Self-Focused Pulse Propagation Is Mediated by Spatiotemporal Optical Vortices

We show that the dynamics of high-intensity laser pulses undergoing self-focused propagation in a nonlinear medium can be understood in terms of the topological constraints imposed by the formation and evolution of spatiotemporal optical vortices (STOVs). STOVs are born from point phase defects on the sides of the pulse nucleated by spatiotemporal phase shear. These defects grow into closed loops of spatiotemporal vorticity that initially exclude the pulse propagation axis, but then reconnect to form a pair of toroidal vortex rings that wrap around it. STOVs constrain the intrapulse flow of electromagnetic energy, controlling the focusing-defocusing cycles and pulse splitting inherent to nonlinear pulse propagation. Here, we illustrate this in two widely studied but very different regimes, relativistic self-focusing in plasma and nonrelativistic self-focusing in gas, demonstrating that STOVs mediate nonlinear propagation irrespective of the detailed physics.

42 ENGINEERING↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
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↗
DOE OSTI2024

Imaging and structure analysis of ferroelectric domains, domain walls, and vortices by scanning electron diffraction

Direct electron detectors in scanning transmission electron microscopy give unprecedented possibilities for structure analysis at the nanoscale. In electronic and quantum materials, this new capability gives access to, for example, emergent chiral structures and symmetry-breaking distortions that underpin functional properties. Quantifying nanoscale structural features with statistical significance, however, is complicated by the subtleties of dynamic diffraction and coexisting contrast mechanisms, which often results in a low signal-to-noise ratio and the superposition of multiple signals that are challenging to deconvolute. Here we apply scanning electron diffraction to explore local polar distortions in the uniaxial ferroelectric Er(Mn,Ti)O 3 . Using a custom-designed convolutional autoencoder with bespoke regularization, we demonstrate that subtle variations in the scattering signatures of ferroelectric domains, domain walls, and vortex textures can readily be disentangled with statistical significance and separated from extrinsic contributions due to, e.g., variations in specimen thickness or bending. The work demonstrates a pathway to quantitatively measure symmetry-breaking distortions across large areas, mapping structural changes at interfaces and topological structures with nanoscale spatial resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization

While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗
DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Automated Operational Forecasting of Monsoon Low Pressure Systems

Monsoon low pressure systems (LPSs) are the dominant rain-bearing weather system of South Asia, often producing extreme precipitation and hydrological disasters in a region inhabited by nearly two billion people. Despite the importance of these storms, no operational system has automatically identified and tracked LPS in real time in numerical weather prediction model output; many commonly used vortex-tracking algorithms are ill suited for monsoon LPS because of the weak winds and cold cores of these systems. Here, we describe a new system that uses optimized algorithms to identify monsoon LPS in short- to medium-range forecasts from the U.S. Global Ensemble Forecast System (GEFS) and a version of the deterministic Global Forecast System (GFS) adapted and used operationally by the Indian Institute of Tropical Meteorology (IITM). We also assess the historical performance of these models in forecasting South Asian monsoon LPS, comparing this with the performance of the Integrated Forecasting System of the ECMWF. We assess the accuracy of model predictions of LPS genesis, position, intensity, and precipitation rates for forecast lead times of 1–5 days, yielding quantitative information on model biases to guide operational forecasters and disaster managers. The system we introduce here could be extended to other low-latitude regions affected by dynamically weak, heavily precipitating atmospheric vortices that are often not included in tropical cyclone inventories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Manipulating Ferroelectric Topological Polar Structures with Twisted Light

Abstract The dynamic control of non‐equilibrium states represents a central challenge in condensed matter physics. While intense terahertz fields drive metal‐insulator transitions and ferroelectricity via soft phonon modes, recent theory suggests that twisted light with orbital angular momentum (OAM) offers a distinct route to manipulate ferroelectric order and stabilize topological excitations including skyrmions, vortices, and Hopfions. Control of ferroelectric polarization in quasi‐2D CsBiNb 2 O 7 (CBNO) is demonstrated using non‐resonant twisted ultra‐violet (UV) light (375 nm, 800 THz). Combining in situ X‐ray Bragg coherent diffractive imaging (BCDI), twisted optical Raman spectroscopy, and density functional theory (DFT), three‐dimensional (3D) ionic displacements, strain fields, and polarization changes are resolved in single crystals. Operando measurements reveal light‐induced strain hysteresis under twisted light–a hallmark of nonlinear, history‐dependent ferroelastic switching driven by OAM. Discrete, irreversible domain transitions emerge as the topological charge ℓ is cycled, stabilizing non‐trivial domain textures including vortex‐antivortex pairs, Bloch/anti‐Bloch points, and merons. These persist after OAM removal, indicating a memory effect. Competing mechanisms are discussed, including multiphoton absorption, strain‐mediated polarization switching, and defect‐wall interactions. The findings establish structured light as a tool for deterministic, reversible control of ferroic states, enabling optically reconfigurable non‐volatile devices.

Chemistry↗