Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗
NASA NTRSDate not supplied

The INSTEP Monitoring Network: Merging High-and-Low Cost Measurements to Characterize California Wildfires

Despite challenges with data quality and scope, low-cost sensor networks have skyrocketed in popularity over the last 15 years, making air quality data available on refined spatial scales. More recently, studies have leveraged both high and low-quality instruments to create stronger “hybrid” models, with most studies focusing on particulate matter. Low-cost measurements typically represent ground-level emissions only, providing context for human health issues from climate change-driven events such as wildfires. Since low-cost sensors’ capabilities are localized, daily events and microclimates tend to dominate the data rather than larger regional or atmospheric trends. Likewise, their low cost explains their high uncertainty. In contrast, some regulatory-grade instruments produce column measurements as well, providing reliable information on a broader scope. To bridge this gap while expanding into gas-phase measurements, we deployed 12 air quality sensor packages in California, USA during the 2022 wildfire season. These INSTEP (Inexpensive Network Sensor Technology Exploring Pollution) monitors measure carbon monoxide (CO), carbon dioxide (CO2), ozone (O3), nitrogen dioxide (NO2), and several hydrocarbons including methane (CH4) and formaldehyde (HCHO). Half of the monitors were co-located with remote sensing spectrometers: NASA Pandora and Total Column Carbon Observing Network (TCCON). The overlap in pollutants includes NO2, O3, and HCHO between the INSTEP monitors and the Pandora column measurements. TCCON covers column CO, CO2, and CH4, rounding out our comparison. Most of the monitors were distributed throughout the San Francisco Bay area, and an additional three were located within 100 km of Los Angeles. The sites ranged in geographic and population characteristics, including desert, mountainous, coastal, and urban locations. Since varying environmental conditions such as temperature and pressure are known to challenge sensor performance, we will apply newer sensor “calibration” techniques meant to combat this. We will normalize our sensor signals by z-scoring them prior to applying a single calibration model in the form of multivariate linear regression or an artificial neural network. While this technique has been validated for the hydrocarbon and ozone sensor types (metal oxide), it has not yet been tested on electrochemical and non-dispersive infrared sensors, which are also used in the INSTEP monitors. This will serve as a test to see if this normalization technique – or another – is most effective in accounting for environmental differences among sensors. Related data analysis efforts have found success with a variety of geospatial analysis techniques, including weighted network models in which high-quality instruments are given higher weights than their low-cost counterparts. Our preliminary analysis will focus on kriging, which uses a Gaussian algorithm to assign weights, providing estimated pollution levels at locations between monitors. Smoke trajectory and evolution will also be considered using both measurement types. We also aim to baseline subtract our emission estimates from each region to determine which portion of emissions are regional and local, further characterizing burn differences in northern and southern California fires. Future directions include using INSTEP jointly with TEMPO satellite data, and mobile deployments on aircraft and uncrewed aerial vehicles (UAV).

Low-cost sensors↗
NASA NTRS2000

Power System Options Evaluated for the Radiation and Technology Demonstration Mission

The Radiation and Technology Demonstration (RTD) Mission is under joint study by three NASA Centers: the NASA Johnson Space Center, the NASA Goddard Space Flight Center, and the NASA Glenn Research Center at Lewis Field. This Earth-orbiting mission, which may launch on a space shuttle in the first half of the next decade, has the primary objective of demonstrating high-power electric thruster technologies. Secondary objectives include better characterization of Earth's Van Allen trapped-radiation belts, measurement of the effectiveness of the radiation shielding for human protection, measurement of radiation effects on advanced solar cells, and demonstration of radiation-tolerant microelectronics. During the mission, which may continue up to 1 year, the 2000-kg RTD spacecraft will first spiral outward from the shuttle-deployed, medium-inclination, low Earth orbit. By the phased operation of a 10-kW Hall thruster and a 10-kW Variable Specific Impulse Magneto-Plasma Rocket, the RTD spacecraft will reach a low-inclination Earth orbit with a radius greater than five Earth radii. This will be followed by an inward spiraling orbit phase when the spacecraft deploys 8 to 12 microsatellites to map the Van Allen belts. The mission will conclude in low Earth orbit with the possible retrieval of the spacecraft by the space shuttle. A conceptual RTD spacecraft design showing two photovoltaic (PV) array wings, the Hall thruster with propellant tanks, and stowed microsatellites is presented. Early power system studies assessed five different PV array design options coupled with a 120-Vdc power management and distribution system (PMAD) and secondary lithium battery energy storage. Array options include (1) state-of-the-art 10-percent efficient three-junction amorphous SiGe thin-film cells on thin polymer panels deployed with an inflatable (or articulated) truss, (2) SCARLET array panels, (3) commercial state-of-the-art, planar PV array rigid panels with 25-percent efficient, three-junction GaInP2/GaAs/Ge solar cells, (4) rigid panels with 25-percent efficient, three-junction GaInP2/GaAs/Ge solar cells, in a 2 -concentrator trough configuration, and (5) thin polymer panels with 25-percent efficient, three-junction GaInP2/GaAs/Ge solar cells deployed with an inflatable (or articulated) truss. To assess the relative merits of these PV array design options, the study group developed a dedicated Fortran code to predict power system performance and estimate system mass. This code also modeled Earth orbital environments important for accurately predicting PV array performance. The most important environmental effect, solar cell radiation degradation, was calculated from electron-proton fluence input from the industry standard AE8/AP8 trapped radiation models and the concept of damage equivalence. Power systems were sized to provide 10 kW of thruster power and approximately 1 kW of spacecraft power at end of life. Of the five PV array design options, the option 1 (thin-film cells) power system was the most massive 590 kg, whereas the option 4 (trough concentrator) power system was the lightest 260 kg. Arguably, the lowest cost would come from the option 3 (commercial array panels) power system with an acceptable, albeit greater, system mass of 320 kg. Predicted power system performance during the spiral-out mission phase is shown the preceding graph for the option 5 (flexible-panel) array. From the results, the radiation-induced power loss over time is evident as the spacecraft slowly spirals outward through the trapped proton belt. The importance of the spiral trip time is also evident in the two curves representing 74-day and 182-day spiral-out periods. The longer spiral time introduces a beginning-of-life power oversizing penalty greater than 1 kW. Future studies will analyze power system performance and mass with a 50-Vdc power management and distribution architecture favorable to the VASIMR thruster and longer missions.

Kerslake, Thomas W.↗
NASA NTRS1972

Information from Jupiter's decametric radiation

Radiation associated with decametric bursts are observed from about 39.5 MHz downwards, and its occurrence probability depends partly on the central meridian longitude of Jupiter. Three central meridian longitude regions of emission are designated as Sources A, B, and C. Emissions from B and C regions are also dependent on Io's phase. The frequency band of noise storms drifts slowly up or down the spectrum. Pulses with durations of 20 microseconds display drifts at the rate of 1000 MHz/sec. Decametric data are used to estimate Jupiter's surface magnetic field strength to be 10 gauss; however, if the maximum emission, about 10 MHz, is the electron gyrofrequency at Io's orbit, the surface field strength would be about 800 gauss. A model of the emission beam control exercised by Io is summarized. Data are presented which indicate that the width of source A emissions results from the changes in the angle of Jovicentric declination of the earth. One of two recent theories places the plasmapause at least 50 Jupiter radii from the surface on the sunward side, and the other places it within Io's orbit.

Carr, T. D.↗
NASA NTRS1980

Fragmentation of rotating protostellar clouds

With a three-dimensional hydrodynamic computer code, the behavior of rotating, isothermal gas clouds as they collapse from Jeans unstable configurations is examined in order to determine whether they are susceptible to fragmentation during the initial dynamic collapse phase of evolution. It is found that a gas cloud will not fragment unless (1) it begins collapsing from a radius much smaller than the Jeans radius (i.e., the cloud initially encloses many Jeans masses) and (2) irregularities in the cloud's initial structure (specifically, density inhomogeneities) enclose more than one Jeans mass of material. Instead of fragmenting, most of the models collapse to a ring configuration. The rings appear to be less susceptible to fragmentation from arbitrary perturbations in their structure than has previously been indicated in other work. Because the models, which include the effects of gas pressure, do not readily fragment during a phase of dynamic collapse, it is suggested that gas clouds in the galactic disk undergo fragmentation only during quasi-equilibrium phases of their evolution.

Tohline, J. E.↗
NASA NTRS1993

Low-frequency source parameters of twelve large earthquakes

A global survey of the low-frequency (1-21 mHz) source characteristics of large events are studied. We are particularly interested in events unusually enriched in low-frequency and in events with a short-term precursor. We model the source time function of 12 large earthquakes using teleseismic data at low frequency. For each event we retrieve the source amplitude spectrum in the frequency range between 1 and 21 mHz with the Silver and Jordan method and the phase-shift spectrum in the frequency range between 1 and 11 mHz with the Riedesel and Jordan method. We then model the source time function by fitting the two spectra. Two of these events, the 1980 Irpinia, Italy, and the 1983 Akita-Oki, Japan, are shallow-depth complex events that took place on multiple faults. In both cases the source time function has a length of about 100 seconds. By comparison Westaway and Jackson find 45 seconds for the Irpinia event and Houston and Kanamori about 50 seconds for the Akita-Oki earthquake. The three deep events and four of the seven intermediate-depth events are fast rupturing earthquakes. A single pulse is sufficient to model the source spectra in the frequency range of our interest. Two other intermediate-depth events have slower rupturing processes, characterized by a continuous energy release lasting for about 40 seconds. The last event is the intermediate-depth 1983 Peru-Ecuador earthquake. It was first recognized as a precursive event by Jordan. We model it with a smooth rupturing process starting about 2 minutes before the high frequency origin time superimposed to an impulsive source.

Harabaglia, Paolo↗
NASA NTRS1993

A minimum cost tolerance allocation method for rocket engines and robust rocket engine design

Rocket engine design follows three phases: systems design, parameter design, and tolerance design. Systems design and parameter design are most effectively conducted in a concurrent engineering (CE) environment that utilize methods such as Quality Function Deployment and Taguchi methods. However, tolerance allocation remains an art driven by experience, handbooks, and rules of thumb. It was desirable to develop and optimization approach to tolerancing. The case study engine was the STME gas generator cycle. The design of the major components had been completed and the functional relationship between the component tolerances and system performance had been computed using the Generic Power Balance model. The system performance nominals (thrust, MR, and Isp) and tolerances were already specified, as were an initial set of component tolerances. However, the question was whether there existed an optimal combination of tolerances that would result in the minimum cost without any degradation in system performance.

Gerth, Richard J.↗
NASA NTRS1960

Noise Considerations for Manned Reentry Vehicles

Noise measurements pertaining mainly to the static firing, launch, 0 and exit flight phases are presented for three rocket-powered vehicles 4 in the Project Mercury test program. Both internal and external data 4 from onboard recordings are presented for a range of Mach numbers and dynamic pressures and for different external vehicle shapes. The main sources of noise are noted to be the rocket engines during static firing and launch and the aerodynamic boundary layer during the high-dynamic-pressure portions of the flight. Rocket-engine noise measurements along the surface of the Mercury Big Joe vehicle were noted to correlate well with data from small models and available data for other large rockets. Measurements have indicated that the aerodynamic noise pressures increase approximately as the dynamic pressure increases and may vary according to the external shape of the vehicle, the highest noise levels being associated with conditions of flow separation. There is also a trend for the aerodynamic noise spectra to peak at higher frequencies as the flight Mach number increases.

Hilton, David A.↗
NASA NTRS1998

The Sound Parameter Effect in Metal-Rich Globular Clusters

Recent Hubble Space Telescope observations have found that the horizontal branches (HBs) in the metal-rich globular clusters NGC 6388 and NGC 6441 slope upward with decreasing B - V. Such a slope is not predicted by canonical HB models and cannot be produced by either a greater cluster age or enhanced mass loss along the red giant branch (RGB). The peculiar HB morphology in these clusters may provide an important clue for understanding the second-parameter effect. We have carried out extensive evolutionary calculations and numerical simulations in order to explore three noncanonical scenarios for explaining the sloped HBs in NGC 6388 and NGC 6441: (1) a high cluster helium abundance scenario, in which the HB evolution is characterized by long blue loops; (2) a rotation scenario, in which internal rotation during the RGB phase increases the HB core mass; and (3) a helium-mixing scenario, in which deep mixing on the RGB enhances the envelope helium abundance. All of these scenarios predict sloped HBs with anomalously bright RR Lyrae variables. We compare this prediction with the properties of the two known RR Lyrae variables in NGC 6388. Additional observational tests of these scenarios are suggested.

Hall, D. K↗
NASA NTRS2015

North Atlantic Simulations in Coordinated Ocean-Ice Reference Experiments Phase II (CORE-II) : Inter-Annual to Decadal Variability - Part II

Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958−2007 period from twenty global ocean - sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958−2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.

Ocean model comparisons↗
NASA NTRS1981

Closed coronal structures. II - Generalized hydrostatic model

Numerical computations of stationary solar coronal loop atmospheres are used to extend earlier analytical work. Two classes of loops are examined, namely symmetric loops with a temperature maximum at the top but now having a length greater than the pressure scale height and loops which have a local temperature minimum at the top. For the first class, new scaling laws are found which relate the base pressure and loop length to the base heating, the heating deposition scale height, and the pressure scale height. It is found that loops for which the length is greater than about two to three times the pressure scale height do not have stable solutions unless they have a temperature minimum at the top. Computed models with a temperature inversion at the top are permitted in a wider range of heating deposition scale height values than are loops with a temperature maximum at the top. These results are discussed in relation to observations showing a dependence of prominence formation and stability on the state of evolution of magnetic structures, and a general scenario is suggested for the understanding of loop evolution from emergence in active regions through the large-scale structure phase to opening in coronal holes.

Serio, S.↗
NASA NTRS1982

Solar flare shocks in interplanetary space and solar flare particle events

Data on particle events resulting from three solar flares are examined to determine the effect of shocks on the population of protons at energies not less than 30 MeV. The arrival of the shock is found to mark the start of the decay phase of the particle event, and changes in particle anisotropy are found to be related to features in the shock structure. Two distinct modes of particle propagation are found to exist: a direct component, and a second component which has a time scale of days and is seen in two of the flares which also produced interplanetary shock. In addition, a model is described in which the protons are accelerated by the shocks, and it is shown to be in agreement with observations.

Evenson, P.↗
NASA NTRS1992

A multilayer perceptron solution to the match phase problem in rule-based artificial intelligence systems

In rule-based AI planning, expert, and learning systems, it is often the case that the left-hand-sides of the rules must be repeatedly compared to the contents of some 'working memory'. The traditional approach to solve such a 'match phase problem' for production systems is to use the Rete Match Algorithm. Here, a new technique using a multilayer perceptron, a particular artificial neural network model, is presented to solve the match phase problem for rule-based AI systems. A syntax for premise formulas (i.e., the left-hand-sides of the rules) is defined, and working memory is specified. From this, it is shown how to construct a multilayer perceptron that finds all of the rules which can be executed for the current situation in working memory. The complexity of the constructed multilayer perceptron is derived in terms of the maximum number of nodes and the required number of layers. A method for reducing the number of layers to at most three is also presented.

Sartori, Michael A.↗
NASA NTRS1997

RXTE observations of AGN

The Rossi X-ray Timing Explorer (RXTE) observed three active galaxies during its in-orbit verification phase: NGC 4151; NGC 4945, and MCG 8-11-11. All three were detected from 2 keV to more than 100 keV by a combination of the proportional counter array (PCA) and the high energy X-ray timing experiment (HEXTE). The PCA contains five, xenon/methane, multilayer, multiwire, gas proportional counters covering the 2 to 60 keV range, while HEXTE is an array of eight NaI/CsI phoswich scintillation counters covering the 15 to 250 keV range. The three active galaxies represent the classes of Seyfert 1, Seyfert 2 and intermediate Seyfert galaxies. The results of the fitting of various models containing partial covering fractions, Compton reflection components and high energy spectral breaks are discussed.

Rothschild, R. E.↗
NASA NTRS2003

Three-Gorge Reservoir: A 'Controlled Experiment' for Calibration/Validation of Time-Variable Gravity Signals Detected from Space

With the advances of measurements, modern space geodesy has become a new type of remote sensing for the Earth dynamics, especially for mass transports in the geophysical fluids on large spatial scales. A case in point is the space gravity mission GRACE (Gravity Recovery And Climate Experiment) which has been in orbit collecting gravity data since early 2002. The data promise to be able to detect changes of water mass equivalent to sub-cm thickness on spatial scale of several hundred km every month or so. China s Three-Gorge Reservoir has already started the process of water impoundment in phases. By 2009,40 km3 of water will be stored behind one of the world s highest dams and spanning a section of middle Yangtze River about 600 km in length. For the GRACE observations, the Three-Gorge Reservoir would represent a geophysical controlled experiment , one that offers a unique opportunity to do detailed geophysical studies. -- Assuming a complete documentation of the water level and history of the water impoundment process and aided with a continual monitoring of the lithospheric loading response (such as in area gravity and deformation), one has at hand basically a classical forwardinverse modeling problem of surface loading, where the input and certain output are known. The invisible portion of the impounded water, i.e. underground storage, poses either added values as an observable or a complication as an unknown to be modeled. Wang (2000) has studied the possible loading effects on a local scale; we here aim for larger spatial scales upwards from several hundred km, with emphasis on the time-variable gravity signals that can be detected by GRACE and follow-on missions. Results using the Green s function approach on the PREM elastic Earth model indicate the geoid height variations reaching several millimeters on wavelengths of about a thousand kilometers. The corresponding vertical deformations have amplitude of a few centimeters. In terms of long-wavelength spherical harmonics, the induced geoid height variations are very close to the accuracy of GRACE- recoverable gravity field, while the low-degree (2 to 5) harmonics should be detectable. With a large regional time-variable gravity signal, the Three-Gorge experiment can serve as a useful calibration/verification for GRACE (including the elastic loading effects), and future gravity missions (especially for visco-elastic yielding as well as underground water variations).

Chao, Benjamin F.↗
NASA NTRS2008

NASA's PEM Fuel Cell Power Plant Development Program for Space Applications

A three-center NASA team led by the Glenn Research Center in Cleveland, Ohio is completing a five-year PEM fuel cell power plant development program for future space applications. The focus of the program has been to adapt commercial PEM fuel cell technology for space applications by addressing the key mission requirements of using pure oxygen as an oxidant and operating in a multi-gravity environment. Competing vendors developed breadboard units in the 1 to 5 kW power range during the first phase of the program, and a single vendor developed a nominal 10-kW engineering model power pant during the second phase of the program. Successful performance and environmental tests conducted by NASA established confidence that PEM fuel cell technology will be ready to meet the electrical power needs of future space missions.

Hoberecht, Mark A.↗
NASA NTRS2017

Variation in Coronal Activity from Solar Cycle 24 Minimum to Maximum Using Three-Dimensional Reconstructions of the Coronal Electron Density from STEREO/COR1

Three-dimensional electron density distributions in the solar corona are reconstructed for 100 Carrington rotations (CR 2054 - 2153) during 2007/03 - 2014/08 using the spherically symmetric method from polarized white-light observations with the inner coronagraph (COR1) onboard the twin Solar Terrestrial Relations Observatory (STEREO). These three-dimensional electron density distributions are validated by comparison with similar density models derived using other methods such as tomography and a magnetohydrodynamics (MHD) model as well as using data from the Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO)-C2. Uncertainties in the estimated total mass of the global corona are analyzed based on differences between the density distributions for COR1-A and -B. Long-term variations of coronal activity in terms of the global and hemispheric average electron densities (equivalent to the total coronal mass) reveal a hemispheric asymmetry during the rising phase of Solar Cycle 24, with the northern hemisphere leading the southern hemisphere by a phase shift of 7 - 9 months. Using 14 CR (approx. equal to 13-month) running averages, the amplitudes of the variation in average electron density between Cycle 24 maximum and Cycle 23/24 minimum (called the modulation factors) are found to be in the range of 1.6 - 4.3. These modulation factors are latitudinally dependent, being largest in polar regions and smallest in the equatorial region. These modulation factors also show a hemispheric asymmetry: they are somewhat larger in the southern hemisphere. The wavelet analysis shows that the short-term quasi-periodic oscillations during the rising and maximum phases of Cycle 24 have a dominant period of 7 - 8 months. In addition, it is found that the radial distribution of the mean electron density for streamers at Cycle 24 maximum is only slightly larger (by approx. equal to 30%) than at cycle minimum.

Wang, Tongjiang↗
NASA NTRS1985

A study on radiative damping of planetary waves utilizing stratospheric observations

Satellite ozone observations made by the Stratospheric Aerosol and Gas Experiment (SAGE) and corresponding meteorological temperature data are used to study the radiative damping processes associated with planetary waves during stratospheric warmings. Ramanathan's model has been adapted for the radiative heating and cooling calculations. The derived infrared damping coefficients, based on observed stratospheric ozone and temperature, are compared with the Newtonian cooling coefficients of Dickinson and Fels. It is also shown that the negative correlation between temperature and ozone solar heating in the upper stratosphere accelerates the damping rate due to infrared cooling alone, in agreement with the theoretical analysis and an earlier report based on observations. In addition, it is also found in this analysis that the phase relationship between ozone and solar heating waves is characterized by its behavior in three distinct layers. In the regions above about 2 mb and also below about 10 mb, the waves are closely in-phase. Between approximately 2 and 10 mb, they show a departure from the in-phase relationship which can be attributed to the so-called 'opacity effect'. This effect significantly determines the magnitude of radiative damping.

Ghazi, A.↗