Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cavity flow”

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,081 records · Page 60

Gridded electron reversal ionizer

A gridded electron reversal ionizer forms a three dimensional cloud of zero or near-zero energy electrons in a cavity within a filament structure surrounding a central electrode having holes through which the sample gas, at reduced pressure, enters an elongated reversal volume. The resultant negative ion stream is applied to a mass analyzer. The reduced electron and ion space-charge limitations of this configuration enhances detection sensitivity for material to be detected by electron attachment, such as narcotic and explosive vapors. Positive ions may be generated by generating electrons having a higher energy, sufficient to ionize the target gas and pulsing the grid negative to stop the electron flow and pulsing the extraction aperture positive to draw out the positive ions.

Chutjian, Ara↗

Intracavity Rayleigh/Mie Scattering for Multipoint, Two-Component Velocity Measurement

A simultaneous multi-point two-component Doppler velocimeter is described. The system uses two optical cavities: a Fabry-Perot etalon and an optical cavity for collecting and re-circulating the Rayleigh/Mie scattered light that is collected from the measurement volume in two parallel, but opposite directions. Single-pulse measurements of two orthogonal components of the velocity vector in a supersonic free jet were performed to demonstrate the technique. The re-circulation of the light rejected by the interferometer input mirror also increased the signal intensity by a factor of 3.5. 2005 Optical Society of America Interferometric Rayleigh scattering has previously been used for single-point velocity measurements in unseeded gas flow. However, this past work has generally been limited to probing with continuous-wave lasers resulting in time-averaged measurements of velocity. Multiple velocity components have been measured simultaneously by separate instruments.1,2 It has also been demonstrated that two orthogonal velocity components can be measured simultaneously at one point using one interferometer by reflecting back the probing laser beam, although this approach results in directional ambiguity of the flow velocity vector.3 This measurement ambiguity was removed by prior knowledge of the approximate magnitude and sign of the velocity components. Furthermore, it was shown that multiple points could be measured simultaneously with a Rayleigh scattering interferometric approach, but only one component of velocity was measured.4 Another method of performing multiple component velocity measurements with Rayleigh scattering uses a pair of cameras to image the flow, one of which views the flow through an iodine gas filter. This iodine-filter technique has the advantage of allowing high-resolution velocity imaging, but it generally has a lower dynamic range.

Bivolaru, Daniel↗

Dynamics and X-ray emission of a galactic superwind interacting with disk and halo gas

There is a general agreement that the conspicuous extranuclear X-ray, optical-line, and radio-contiuum emission of starbursts is associated with powerful galactic superwinds blowing from their centers. However, despite the significant advances in observational studies of superwinds, there is no consensus on the nature of the emitting material and even on the emission mechanisms themselves. This is to a great extent a consequence of a poor understanding of dynamical processes in the starburst superwind regions. To address this issue, we have conducted two-dimensional hydrodynamical simulations of galactic superwinds. While previous similar studies have used a single (disk) component to represent the ISM of the starburst galaxy, we analyze the interaction of the wind with a two-component disk-halo ambient interstellar medium and argue that this two-component representation is crucial for adequate modeling of starbursts. The emphasis of this study is on the geometry and structure of the wind region and the X-ray emission arising in the wind material and the shocked gas in the disk and the halo of the galaxy. The simulation results have shown that a clear-cut bipolar wind can easily develop under a range of very different conditions. On the other hand, a complex 'filamentary' structure associated with the entrained dense disk material is found to arise within the hot bubble blown out by the wind. The flow pattern within the bubble is dominated equally by the central biconic outflow and a system of whirling motions r elated to the origin and development of the 'filaments'. The filament parameters make them a good candidate for optical-emission-line filamentary gas observed in starburst halos. We find that the history of mass and energy deposition in the starburst region of the galaxy is crucial for wind dynamics. A 'mild' early wind, which arises as a result of the cumulative effect of stellar winds from massive stars, produces a bipolar vertical cavity in the disk and halo gas without strongly affecting the gaseous disk, thus creating conditions for virtually free vertical escape of the hot gas at the later, much more violent supernova-dominated phases of the starburst. We calculate the luminosity, mass, and effective temperature of the X-ray emitting gas in the 'soft' (0.1 to 0.7 keV, 0.7 to 2.2 keV, and 0.1 to 2.2 keV) and 'hard' (1.6 to 8.3 keV) energy bands and estimate the contribution of different gaseous components to the X-ray flux in these bands. Analysis of these parameters enables us to make conclusions regarding the nature of the X-ray-emitting material. We have inferred that the bulk of the soft thermal X-ray emission from starbursts arises in the wind-shocked material of the disk and halo gas rather than in the wind material itself. This enables us to predict that the integrated soft X-ray spectra of starbursts need not show an overabundance of heavy elements which are believed to be produced copiously in the centers of starbursts. Unlike soft X-ray emission, the hard component of thermal X-ray emission is found to originate in the wind material ejected from the starburst region. However, the derived ratio of hard-to-soft X-ray luminosities is too small compared to that observed in starbursts. We conclude therefore that the observed hard X-ray emission of starbursts is probably not associated with the thermal emission of hot wind or ambient shocked gas. Typical temperatures of the bulk of the soft X-ray-emitting material in our very different models have been found to agree well with the ones estimated on the basis of the ROSAT data for the soft component of X-ray emission of nearby starbursts. We predict that temperatures of the extranuclear soft X-ray-emitting gas in starburst galaxies with heavy element abundances near solar should be close to T(sub Xs = 2 to 5 x 10(exp 6)K.

Suchkov, Anatoly A.↗

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗

Design and Operating Characteristics of High-Speed, Small-Bore, Angular-Contact Ball Bearings

The computer program SHABERTH was used to analyze 35-mm-bore, angular-contact ball bearings designed and manufactured for high-speed turbomachinery applications. Parametric tests of the bearings were conducted on a high-speed, high-temperature bearing tester and were compared with the computer predictions. Four bearing and cage designs were studied. The bearings were lubricated either by jet lubrication or through the split inner ring with and without outer-ring cooling. The predicted bearing life decreased with increasing speed because of increased operating contact stresses caused by changes in contact angle and centrifugal load. For thrust loads only, the difference in calculated life for the 24 deg. and 30 deg. contact-angle bearings was insignificant. However, for combined loading, the 24 deg. contact-angle bearing gave longer life. For split-inner-ring bearings, optimal operating conditions were obtained with a 24 deg. contact angle and an inner-ring, land-guided cage, using outer-ring cooling in conjunction with low lubricant flow rates. Lower temperature and power losses were obtained with a single-outer-ring, land-guided cage for the 24 deg. contact-angle bearing having a relieved inner ring and partially relieved outer ring. Inner-ring temperatures were independent of lubrication mode and cage design. In comparison with measured values, reasonably good engineering correlation was obtained using the computer program SHABERTH for predicted bearing power loss and for inner- and outer-ring temperatures. The Parker formula for XCAV (used in SHABERTH, a measure of oil volume in the bearing cavity) may need to be refined to reflect bearing lubrication mode, cage design, and location of cage-controlling land.

Pinel, Stanley I.↗

Stable water isotopes and tritium data from porewater at Elkhorn Slough

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we collected geochemical tracers across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. The experimental transect is located at the Elkhorn Slough National Estuarine Research Reserve.Here we present a data table (.csv) with the results of electrical conductivity, stable water isotopes, and tritium in porewater across several elevations in the hillslope-to-salt marsh continuum (information in the data set). We analyzed the stable water isotope samples by cavity ring-down spectroscopy in a Picarro L2130-i at the University of California Santa Cruz (accuracy of 0.025‰ and 0.1‰ for δ18O and δ2H, respectively). We measured the electric conductivity of the samples with an Orion Star™ A329 multiparameter meter (0.5% accuracy, Thermo Fisher Scientific, Massachusetts, USA). We analyzed 3H samples at Lawrence Livermore National Laboratory by helium-3 accumulation (Clarke et al., 1976; Surano et al., 1992).

54 ENVIRONMENTAL SCIENCES↗

Daily water stable isotopes, transpiration, and matrix potential data for an aspen and engelmann stand in the East River Watershed (version 2)

We provide daily stable isotope (2H & 18O) ratios in soil water and xylem (plant stem) water, as well as the sap flow (transpiration) and the soil's matric potential at a forested site near Gothic, Colorado, in the East River catchment. We measured the stable isotopic composition of the transpiration and the daily transpiration flux sum of three aspen and three engelmann spruce. In both forest stands, we installed a soil profile and measured the soil matric potential at 15, 30, and 60 cm depth as well as the stable isotopes of soil pore water at 5, 10, 30, 60, and 90 cm depths. All isotope measurements were done in situ via vapor probes connected to a cavity ring down spectrometer (Picarro L1240i).We further report the daily meteorological data observed at billy barr near our study site. We also provide for each tree the relative share of root water uptake derived from the isotope measurements via a Bayesian mixing model (MixSIAR).The daily data is provided as a time series in "Iso_MP_Sap_DataDaily_ESSDiveUpload.csv" and the units are provided in "dd.csv"; the location of the instrumented trees and soil profiles are given as latitude and longitude coordinates saved as CSV and KMZ files; and a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata.The data was gathered to investigate the short-term changes of the water sources (i.e., variation of root water uptake from different soil depths) of the studied subalpine trees.Update 07/16/2025: The relative and absolute plant water uptake depths were grouped to ensure that the MixSIAR model was applied with endmembers that differed in their d2H value by at least 3 permill and at least 1 permill for d18O. Whenever the difference between observed d2H values for two or more probes at neighboring depths was less than 3 permill, we used the average value for the source water endmember. For days at which probes that were not next to each other measurements did not differ at least 3 permill, the average of all probes between these two depths was used as the water source endmember representing the depth range between these two probes.

54 ENVIRONMENTAL SCIENCES↗

Circumferential pressure probe

A probe for measuring circumferential pressure inside a body cavity is disclosed. In the preferred embodiment, a urodynamic pressure measurement probe for evaluating human urinary sphincter function is disclosed. Along the length of the probe are disposed a multiplicity of deformable wall sensors which typically comprise support tube sections with flexible side wall areas. These are arranged along the length of the probe in two areas, one just proximal to the tip for the sensing of fluid pressure inside the bladder, and five in the sensing section which is positioned within the urethra at the point at which the urinary sphincter constricts to control the flow of urine. The remainder of the length of the probe comprises multiple rigid support tube sections interspersed with flexible support tube sections in the form of bellows to provide flexibility.

Holmes, Harlan K.↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An improved offset generator developed for Allan deviation measurement of ultra stable frequency standards

Measuring the performance of ultra stable frequency standards such as the Superconducting Cavity Maser Oscillator (SCMO) necessitates improvement of some test instrumentation. The frequency stability test equipment used at JPL includes a 1 Hz Offset Generator to generate a beat frequency between a pair of 100 MHz signals that are being compared. The noise floor of the measurement system using the current Offset Generator is adequate to characterize stability of hydrogen masers, but it is not adequate for the SCMO. A new Offset Generator with improved stability was designed and tested at JPL. With this Offset Generator and a new Zero Crossing Detector, recently developed at JPL, the measurement flow was reduced by a factor of 5.5 at 1 second tau, 3.0 at 1000 seconds, and 9.4 at 10,000 seconds, compared against the previous design. In addition to the new circuit designs of the Offset Generator and Zero Crossing Detector, tighter control of the measurement equipment environment was required to achieve this improvement. The design of this new Offset Generator are described, along with details of the environment control methods used.

R. L. Hamell↗

The interaction of radio frequency electromagnetic fields with atmospheric water droplets and applications to aircraft ice prevention

The feasibility of computerized simulation of the physics of advanced microwave anti-icing systems, which preheat impinging supercooled water droplets prior to impact, was investigated. Theoretical and experimental work performed to create a physically realistic simulation is described. The behavior of the absorption cross section for melting ice particles was measured by a resonant cavity technique and found to agree with theoretical predictions. Values of the dielectric parameters of supercooled water were measured by a similar technique at lambda = 2.82 cm down to -17 C. The hydrodynamic behavior of accelerated water droplets was studied photograhically in a wind tunnel. Droplets were found to initially deform as oblate spheroids and to eventually become unstable and break up in Bessel function modes for large values of acceleration or droplet size. This confirms the theory as to the maximum stable droplet size in the atmosphere. A computer code which predicts droplet trajectories in an arbitrary flow field was written and confirmed experimentally. The results were consolidated into a simulation to study the heating by electromagnetic fields of droplets impinging onto an object such as an airfoil. It was determined that there is sufficient time to heat droplets prior to impact for typical parameter values. Design curves for such a system are presented.

Hansman, R. J., Jr.↗

Vented compression molding

A new process, vented compression molding, is developed for molding low density particulates, such as ablative materials used in the Space Shuttle. This process is a modification of standard compression molding that uses a perforated surface as the molding surface. An excess of material is used to fill the mold cavity and this material is forced through the openings in the surface of the mold as the mold closes. This process creates frictional forces, which can be determined approximately by Mohr-Coulomb yield criterion, generating the pressure that compacts the material. The pressure in the mold is shown to be governed by the geometry of the openings and the mechanical properties of the particulate, and not by the displacement or the volume of the mold. It is determined that this uncoupling of the pressure-displacement relationship allows uniform material distribution and compaction without requiring large-scale flow of the particulates.

Mccree, J. O.↗

Rotordynamic forces in labyrinth seals: Theory and experiment

A theoretical and experimental investigation of the aerodynamic forces generated by a single gland labyrinth seal executing a simultaneous spinning/whirling motion has been conducted. A lumped parameter model for a single gland seal with coupling to an upstream cavity with leakage is developed along with an appropriate solution technique. From this theory, it is shown that the presence of the upstream cavity can, in some cases, augment the cross-stiffness and direct damping by a factor of four. The parameters that govern the coupling are presented along with predictions on their influence. A simple uncoupled model is used to identify the mechanisms responsible for cross force generation. This reduced system is nondimensionalized and the physical significance of the reduced parameters is discussed. Closed form algebraic formulas are given for some simple limiting cases. It is also shown that the total cross-force predicted by the uncoupled model can be represented as the sum of an ideal component due to an inviscid flow with entry swirl and a viscous part due to the change in swirl created by friction inside the gland. The frequency dependent ideal part is solely responsible for the rotordynamic direct damping. The facility designed and built to measure these frequency dependent forces is described. Experimental data confirm the validity and usefulness of this ideal/viscous decomposition. A method for calculating the damping coefficients based on the force decomposition using only the static measurements is presented. Experimental results supporting the predicted cross force augmentation due to the effect of upstream coupling are presented.

Millsaps, Knox T.↗

On some structure-turbulence interaction problems

The interactions between a turbulent flow structure; responding to its excitation were studied. The turbulence was typical of those associated with a boundary layer, having a cross-spectral density indicative of convection and statistical decay. A number of structural models were considered. Among the one-dimensional models were an unsupported infinite beam and a periodically supported infinite beam. The fuselage construction of an aircraft was then considered. For the two-dimensional case a simple membrane was used to illustrate the type of formulation applicable to most two-dimensional structures. Both the one-dimensional and two-dimensional structures studied were backed by a cavity filled with an initially quiescent fluid to simulate the acoustic environment when the structure forms one side of a cabin of a sea vessel or aircraft.

Maekawa, S.↗

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental modeling of jet-ring turbine disk cooling

The experimental facility and some early results are described from a current research program modeling turbine disk cooling with multiple impinging jets, such as employed on the Space Shuttle Main Engine oxygen turbopump. The study is designed to obtain detailed local convection heat transfer rates on specially constructed turbine disk models that employ either multiple cooling jet impingement near the disk outer radius from a jet ring, or alternatively, single entrance coolant supply into the center of the disk cavity. Jet impingement is an effective scheme for cooling of turbine disks at or near the blade attachment region, but the heat transfer mechanisms and merits relative to other schemes are not well understood. The present study employs two specially constructed full scale plastic model disks, contoured and plane, together with the corresponding stators. Local heat transfer rates are determined by a computer vision system from the response of thin liquid crystal coatings applied to the disk test faces. The present results indicate that multiple jet impingement directed at the blade attachment region results in higher cooling rates in that region than does the same flow supplied to the disk center, but this conclusion must be regarded as tentative.

Metzger, D. E.↗

The performance of smooth-wall drag reducing outer-layer devices in rough-wall boundary layers

The boundary layer drag on a 0.28 m x 0.91 m flat-plate was measured to evaluate the performance of smooth-wall drag reducing outer-layer devices in rough-wall boundary layers up to a chord Reynolds number, Re(c) of 17,000. The devices were a pair of thin flat ribbons placed in tandem. The wall conditions examined were smooth, transverse-groove roughnesses having cavity width/depth ratios of 0.7 (d-type) and 3 (k-type), and sandpaper roughness (k-type) of grit sizes 180, 50 and 36. The area-averaged wall-drag was reduced compared to the respective normal levels for all wall conditions. All k-type rough walls exhibited a comparable level of drag reduction (about 7 percent) which was lower than that in the smooth-wall. The d-type rough-wall exhibited a transitional behavior in that the drag reduction dropped from the smooth-wall level to that of the k-type roughnesses with increase in the flow speed. In further experiments on a 6 meter long sandpaper roughened wall and with a pair of symmetric airfoil devices in tandem at an Re(c) of 77,500, it was possible to recover the device drag penalty fully, but no net drag reduction was achieved. The apparent loss of performance of the devices in a rough-wall is explained in terms of their inability to alter the pressure component of the drag in a significant way.

Bandyopadhyay, P. R.↗

A new model of resonance in the winter stratosphere

It is generally accepted that the planetary waves observed in the winter stratosphere are primarily a response to dynamical forcing from the troposphere. Nevertheless, the mechanism by which wave amplitudes sometimes become larger remains uncertain. It is possible that anomalously large waves in the stratosphere might simply be the result of anomalously large tropospheric forcing. However, it has also been suggested that they are a response to a stratospheric-tropospheric cavity being in a near-resonant configuration. It has been suggested that nonlinear self-tuning effects could play an important role in the behavior of such a cavity. Self-tuning may occur when a system starts to one side of resonance, such that the mean-state change induced by growing waves brings the system closer to resonance. A new model of the stratospheric cavity is introduced and is then used to re-examine the possibility of wave growth in the real atmosphere and in atmospheric models due to self-tuning effects. The new model is based on the picture of the winter-time stratosphere which has been revealed by the observations of Ertel's potential vorticity, Q. Isentropic maps of Q show two rather distinct regions, the first containing the circumpolar vortex, where gradients of Q are large and Rossby waves may propagate easily. Surrounding this is a second, low-latitude region where the gradients are generally weak and where, because the Eulerian-mean flow is comparable with their phase speed, Rossby waves must be continually breaking. As the waves are observed to grow the relative sizes of these two regions change in time. This leads to the interesting possibility that self-tuning, mainly due to irreversible changes in the size of the polar vortex, is taking place.

Haynes, P. H.↗