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 1,441 records · Page 80

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

Direct Numerical Simulations↗

Experiments with phase change thermal energy storage canisters for Space Station Freedom

The solar dynamic power module proposed for the Space Station Freedom (SSF) uses the heat of fusion of a phase change material (PCM) to efficiently store thermal energy for use during eclipse periods. The PCM, a LiF-20CaF2 salt, is contained in annular, metal canisters located in a heat receiver at the focus of a solar concentrator. PCM canister ground-based experiments and analytical heat transfer studies are discussed. The hardware, test procedures, and test results from these experiments are discussed. After more than 900 simulated SSF orbital cycles, no canister cracks or leaks were observed and all data were successfully collected. The effect of 1-g test orientation on canister wall temperatures was generally small while void position was strongly dependent on test orientation and canister cooling. In one test orientation, alternating wall temperature data were measured that supports an earlier theory of oscillating vortex flow in the PCM melt. Analytical canister wall temperatures compared very favorably with experimental temperature data. This illustrates that ground-based canister thermal performance can be predicted well by analyses that employ straight-forward, engineering models of void behavior and liquid PCM free convection.

Kerslake, Thomas W.↗

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

Seasonal Evolution of Saturn's Polar Temperatures and Composition

The seasonal evolution of Saturn's polar atmospheric temperatures and hydrocarbon composition is derived from a decade of Cassini Composite Infrared Spectrometer (CIRS) 7-16 micrometers thermal infrared spectroscopy. We construct a near-continuous record of atmospheric variability poleward of 60 deg from northern winter/southern summer (2004, Ls = 293 deg) through the equinox (2009, Ls= 0 deg) to northern spring/southern autumn (2014, Ls = 56 deg). The hot tropospheric polar cyclones that are entrained by pro-grade jets within 2-3 deg of each pole, and the hexagonal shape of the north polar belt, are both persistent features throughout the decade of observations. The hexagon vertices rotated westward by approx. equal to 30 deg longitude between March 2007 and April 2013, confirming that they are not stationary in the Voyager-defined System III longitude system as previously thought. Tropospheric temperature contrasts between the cool polar zones (near 80-85 deg) and warm polar belts (near 75-80 deg) have varied in both hemispheres, resulting in changes to the vertical wind shear on the zonal jets in the upper troposphere and lower stratosphere. The extended region of south polar stratospheric emission has cooled dramatically poleward of the sharp temperature gradient near 75 deg S (by approximately -5 K/yr), coinciding with a depletion in the abundances of acetylene (0030 +/- 0.005 ppm/yr) and ethane (0.35 +/- 0.1 ppm/yr), and suggestive of stratospheric upwelling with vertical wind speeds of w approx. equal to +0.1 mm/s. The upwelling appears most intense within 5 deg latitude of the south pole. This is mirrored by a general warming of the northern polar stratosphere (+5 K/yr) and an enhancement in acetylene (0.030 +/- 0.003 ppm/yr) and ethane (0.45 +/- 0.1 ppm/yr) abundances that appears to be most intense poleward of 75 deg N, suggesting subsidence at w approx. equal to -0.15 mm/ s. However, the sharp gradient in stratospheric emission expected to form near 75 deg N by northern summer solstice (2017, Ls = 90 deg) has not yet been observed, so we continue to await the development of a northern summer stratospheric vortex. The peak stratospheric warming in the north occurs at lower pressure levels (p less than 1 mbar) than the peak stratospheric cooling in the south (p greater than 1 mbar). Vertical motions are derived from both the temperature field (using the measured rates of temperature change and the deviations from the expectations of radiative equilibrium models) and hydrocarbon distributions (solving the continuity equation). Vertical velocities tend towards zero in the upper troposphere where seasonal temperature contrasts are smaller, except within the tropospheric polar cyclones where w approx. equal to +0.02 mm/s. North polar minima in tropospheric and stratospheric temperatures were detected in 2008-2010 (lagging one season, or 6-8 years, behind winter solstice); south polar maxima appear to have occurred before the start of the Cassini observations (1-2 years after summer solstice), consistent with the expectations of radiative climate models. The influence of dynamics implies that the coldest winter temperatures occur in the 75-80 deg region in the stratosphere, and in the cool polar zones in the troposphere, rather than at the poles themselves. In addition to vertical motions, we propose that the UV-absorbent polar stratospheric aerosols entrained within Saturn's vortices contribute significantly to the radiative budget at the poles, adding to the localized enhancement in the south polar cooling and north polar warming poleward of +/-75 deg.

polar temperatures↗

Final Report: Nanoscale Dynamical Heterogeneity in Complex Magnetic Materials

A magnetic object can be demagnetized by dropping it on a hard surface, but what does ‘demagnetized’ actually mean? In 1919 Heinrich Barkhausen proved the existence of magnetic domains, which are regions of uniform magnetization that are much larger than atoms but much smaller than a macroscopic object. A material is fully magnetized when domain magnetizations are aligned, while it is demagnetized when the domain magnetizations are randomly oriented and the net magnetization is zero. The heterogeneity of a demagnetized object leads to interesting questions. Magnets are unstable when their poles align, and stable when their poles anti-align, so why is the magnetized state ever stable? What do domains look like? What is the structure of a domain wall? How does the magnetized state transform to the demagnetized state? How do domains appear and disappear? What are the statistical properties of domains and how do these vary as the domain pattern evolves? Some of these questions remain the focus of intense study nearly a century after Barkhausen’s discovery. For example, just a few years ago a new kind of magnetic texture called a skyrmion was discovered. A skyrmion is a magnetic domain that is a nanometer-scale, topologically protected vortex. ‘Topologically protected’ means that skyrmions are hard to destroy and so are stable for extended periods. Skyrmions are characterized by integral quantum numbers and are observed to move with little dissipation and so could store and process information with very low power input. Our research project uses soft x-rays, which offer very high magnetic contrast, to probe magnetic heterogeneity and to measure how it evolves in time under external influences. We will condition a soft x-ray beam so that the wave fronts will be coherent, that is, they will be smooth and well-defined. When coherent soft x-ray beam interacts with a magnetic material, the magnetic heterogeneity is imprinted onto the wave fronts and projected into a diffraction pattern. These patterns will be analyzed to understand the structure, motion, and statistical properties of magnetic textures and their boundaries. Over the period covered by this grant we will study a) the structure, phase behaviors, and motion of skyrmions in various thin film materials, and 2) the statistical properties of Barkhausen cascades, which are a key factor in how magnetization varies.

36 MATERIALS SCIENCE↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

HART-II Acoustic Predictions using a Coupled CFD/CSD Method

This paper documents results to date from the Rotorcraft Acoustic Characterization and Mitigation activity under the NASA Subsonic Rotary Wing Project. The primary goal of this activity is to develop a NASA rotorcraft impulsive noise prediction capability which uses first principles fluid dynamics and structural dynamics. During this effort, elastic blade motion and co-processing capabilities have been included in a recent version of the computational fluid dynamics code (CFD). The CFD code is loosely coupled to computational structural dynamics (CSD) code using new interface codes. The CFD/CSD coupled solution is then used to compute impulsive noise on a plane under the rotor using the Ffowcs Williams-Hawkings solver. This code system is then applied to a range of cases from the Higher Harmonic Aeroacoustic Rotor Test II (HART-II) experiment. For all cases presented, the full experimental configuration (i.e., rotor and wind tunnel sting mount) are used in the coupled CFD/CSD solutions. Results show good correlation between measured and predicted loading and loading time derivative at the only measured radial station. A contributing factor for a typically seen loading mean-value offset between measured data and predictions data is examined. Impulsive noise predictions on the measured microphone plane under the rotor compare favorably with measured mid-frequency noise for all cases. Flow visualization of the BL and MN cases shows that vortex structures generated in the prediction method are consist with measurements. Future application of the prediction method is discussed.

Boyd, D. Douglas, Jr.↗

Arcjet thruster research and technology. Phase 1: Executive summary

The principal objective of this two phase program is to conduct the development research required to make the low power arcjet a flight ready technology. Many important results were obtained during Phase 1 to move closer to this objective. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. These aided the conceptual understanding of the arcjet and guided design work. A hydrazine (N2H4) arcjet was designed that combined a flight qualified catalyst bed with a modular arcjet. Extensive testing was performed which demonstrated the feasibility of using this propellant in an arcjet for the first time. Startup techniques were developed, stability maintained, material compatibility tests conducted, and performance mapping tests performed. Specific impulse values from 400 to 730 seconds were produced with a non-optimized design. These levels are higher than were originally thought possible and proved that extremely high enthalpy values can be obtained with constricted arc technology. Erosion rate data are promising for lifetime extensions to meet flight application requirements. Power control unit (PCU) development was started with the design and fabrication of a laboratory high switching frequency supply. Valuable data were obtained on PCU operation and on the interaction with the dynamic arc. Phase 2 efforts presently underway are resolving key issues for multi-hundred hour lifetimes, are continuing to investigate arcjet/PCU interactions, and will demonstrate duty cycle N2H4 arcjet/PCU operation in a simulated flight mode for lifetimes consistent with initial applications.

Source record↗

Features of steady magnetospheric convection

The large-scale patterns of ionospheric convection and particle precipitation are described during two intervals of steady magnetospheric convection (SMC) on November 24, 1981. The unique data set used in the analysis includes recordings from the worldwide network of magnetometers and all-sky cameras, global auroral images from the Dynamics Explorer (DE) 1 spacecraft, and particle precipitation data from low-altitude National Oceanic and Atmospheric Administration (NOAA) 6 and NOAA 7 spacecraft. The data show that intense magnetospheric convection continued during more than 10 hours under the steady southward interplanetary magnetic field without any distinct substorm signatures. All data sets available confirmed the stable character of the large-scale magnetospheric configuration during this period. In particular, the magnetic flux threading the polar cap was stable (within 10%) during 3.5 hours of continued DE 1 observations. The dayside cusp was located at an unusually low latitude (70 deg CGL). The nightside auroral pattern consisted of two distinct regions. The diffuse aurora in the equatorward half of the expanded (10 deg wide) auroral oval was well-separated from the bright, active auroral forms found in the vicinity of the poleward boundary of the oval. The twin-vortex convection pattern had no signature of the Harang discontinuity; its nightside 'convection throat' was spatially coincident with the poleward active auroras. This region of the auroral oval was identified as the primary site of the short-lived transient activations during the SMC intervals. The energetic particle observations show that the auroral precipitation up to its high-latitude limit is on closed field lines and that particle acceleration up to greater than 30-keV energy starts close to this limit. The isotropic boundaries of the greater than 30-keV protons and electrons were found close to each other, separating regions of discrete and diffuse precipitation. This suggests that these precipitation types originate on the very taillike and very dipolelike field lines, respectively.

Yahnin, A.↗

The Evolution of Friction Stir Welding Theory at Marshall Space Flight Center

From 1995 to the present the friction stir welding (FSW) process has been under study at Marshall Space Flight Center (MSFC). This is an account of the progressive emergence of a set of conceptual tools beginning with the discovery of the shear surface, wiping metal transfer, and the invention of a kinematic model and making possible a treatment of both metallurgical structure formation and process dynamics in friction stir welding from a unified point of view. It is generally observed that the bulk of the deformation of weld metal around the FSW pin takes place in a very narrow, almost discontinuous zone with high deformation rates characteristic of metal cutting. By 1999 it was realized that this zone could be treated as a shear surface like that in simple metal cutting models. At the shear surface the seam is drawn out and compressed and pressure and flow conditions determine whether or not a sound weld is produced. The discovery of the shear surface was followed by the synthesis of a simple 3- flow kinematic model of the FSW process. Relative to the tool the flow components are: (1) an approaching translational flow at weld speed V, (2) a rotating cylindrical plug flow with the angular velocity of the tool , and (3) a relatively slow ring vortex flow (like a smoke ring) encircling the tool and driven by shoulder scrolls and pin threads. The rotating plug flow picks up an element of weld metal, rotates it around with the tool, and deposits it behind the tool ( wiping metal transfer ); it forms plan section loops in tracers cut through by the tool. Radially inward flow from the ring vortex component retains metal longer in the rotating plug and outward flow expels metal earlier; this interaction forms the looping weld seam trace and the tongue and groove bimetallic weld contour. The radial components of the translational and ring vortex flows introduce parent metal intrusions into the small grained nugget material close to the tool shoulder; if this feature is pronounced, nugget collapse may result. Certain weld features, in particular internal banding seen in transverse section as onion rings and associated surface ridges called tool marks , have long implied an oscillation flow component, but have only recently been attributed in the literature to tool eccentricity. Rotating plug shape, typically a hollow cylinder flared at the end where it sticks to the shoulder, varies as pressure distribution on the tool determines where sticking occurs. Simplified power input estimates balanced against heat loss estimates give reasonable temperature estimates, explain why the power requirement changes hardly at all over a wide range of RPM s, and yield isotherms that seem to fall along boundaries of parameter windows of operation.

Nunes, Arthur C.↗

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↗

Numerical design of advanced multi-element airfoils

The current study extends the application of computational fluid dynamics to three-dimensional high-lift systems. Structured, overset grids are used in conjunction with an incompressible Navier-Stokes flow solver to investigate flow over a two-element high-lift configuration. The computations were run in a fully turbulent mode using the one-equation Baldwin-Barth turbulence model. The geometry consisted of an unswept wing which spanned a wind tunnel test section. Flows over full and half-span Fowler flap configurations were computed. Grid resolution issues were investigated in two dimensional studies of the flapped airfoil. Results of the full-span flap wing agreed well with experimental data and verified the method. Flow over the wing with the half-span was computed to investigate the details of the flow at the free edge of the flap. The results illustrated changes in flow streamlines, separation locations, and surface pressures due to the vortex shed from the flap edge.

Mathias, Donovan L.↗

Computational Aeroacoustic Analysis System Development

Many industrial and commercial products operate in a dynamic flow environment and the aerodynamically generated noise has become a very important factor in the design of these products. In light of the importance in characterizing this dynamic environment, Rocketdyne has initiated a multiyear effort to develop an advanced general-purpose Computational Aeroacoustic Analysis System (CAAS) to address these issues. This system will provide a high fidelity predictive capability for aeroacoustic design and analysis. The numerical platform is able to provide high temporal and spatial accuracy that is required for aeroacoustic calculations through the development of a high order spectral element numerical algorithm. The analysis system is integrated with well-established CAE tools, such as a graphical user interface (GUI) through PATRAN, to provide cost-effective access to all of the necessary tools. These include preprocessing (geometry import, grid generation and boundary condition specification), code set up (problem specification, user parameter definition, etc.), and postprocessing. The purpose of the present paper is to assess the feasibility of such a system and to demonstrate the efficiency and accuracy of the numerical algorithm through numerical examples. Computations of vortex shedding noise were carried out in the context of a two-dimensional low Mach number turbulent flow past a square cylinder. The computational aeroacoustic approach that is used in CAAS relies on coupling a base flow solver to the acoustic solver throughout a computational cycle. The unsteady fluid motion, which is responsible for both the generation and propagation of acoustic waves, is calculated using a high order flow solver. The results of the flow field are then passed to the acoustic solver through an interpolator to map the field values into the acoustic grid. The acoustic field, which is governed by the linearized Euler equations, is then calculated using the flow results computed from the flow solver.

Hadid, A.↗

Observed increase in the peak rain rates of monsoon depressions

Most extreme precipitation in the densely populated region of central India is produced by atmospheric vortices called monsoon lows and monsoon depressions. Here we use satellite and gauge-based precipitation estimates with atmospheric reanalyses to assess 40-year trends in the rain rates of these storms, which have remained unknown. We show that rain rates increased in the rainiest quadrant of monsoon depressions, southwest of the vortex center; precipitation decreased in eastern quadrants, yielding no clear trend in precipitation averaged over the entire storm diameter. In an atmospheric reanalysis, ascent increased in the region of amplifying precipitation, but we could not detect trends in the intensity of rotational winds around the storm center. These storm changes occurred in a background environment where humidity increased rapidly over land while warming was more muted. Monsoon lows, which we show produce less precipitation than depressions, exhibit weaker trends that are less statistically robust.

54 ENVIRONMENTAL SCIENCES↗

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↗