Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High Gradient Acceleration”

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 127 records · Page 7

Thermal Performance of a Conduction-Cooled CCT Dipole ReBCO Magnet: Several Cycles of Cool-Down and Thermal Gradient Measurements

Here, this paper presents experimental results from conduction-cooled thermal testing of a ReBCO canted cosine theta (CCT) magnet (C2), originally designed and fabricated at LBNL using CORC cables. While the performance of the coil under liquid helium and nitrogen environments has been previously established, this study explores its behavior under conduction cooling using a large test cryostat at The Ohio State University. The magnet, measuring 613 mm in length and weighing ∼75 kg, consists of four helical layers wound with ReBCO-based CORC wire and was thermally anchored to a copper cold ring supported by a G-10 strongback. Cooling was provided by two Sumitomo RDK-415D cryocoolers, offering a combined 3 W at 4.2 K and 150 W at 77 K. Multiple thermal cycles were performed, with cooldown durations of up to 45 hours. Final base temperatures of approximately 10.8 K (at the coil edge) and 12.0 K (at the coil center) were achieved, with an axial temperature difference of approximately 1.2 K. The warm-up period extended over approximately 24.6 hours. Voltage measurements from all four layers were recorded during cooldown and warmup. The system demonstrated stable cooldown performance, repeatable gradients, and good thermal anchoring. These results support the feasibility of conduction cooling in large-scale HTS magnets, aligning with broader goals for “green” cryogen-free accelerator technologies and paving the way for more sustainable, scalable, and energy-efficient high-field magnet systems in next-generation particle accelerators.

Canted cosine theta magnets↗

Extending SEER for Extreme Heterogeneity

Heterogeneous and multi-device nodes are increasingly common in high-performance computing and data centers, yet existing programming models often lack simple, transparent, and portable support for these diverse architectures. The main contribution of this work is the development of novel SEER capabilities to address this challenge by providing a descriptive programming model that allows applications to seamlessly leverage heterogeneous nodes across various device types. SEER uses efficient memory management and can select the proper device[s] depending on the computational cost of the applications. This is completely transparent to the programmer, thereby providing a highly productive programming environment. Integrating extreme heterogeneity into the SEER library as shown with the use of NVIDIA and AMD GPUs simultaneously allows it to expand and exploit the performance possibilities. Our analysis based on the well-known Conjugate Gradient algorithm reports accelerations above 1.5 × on computationally demanding steps of such an algorithm by using both architectures simultaneously.

Teranishi, Keita [ORNL] (ORCID:0000000166472690)↗

Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators

Hamiltonian Monte Carlo (HMC) is a powerful and accurate method to sample from the posterior distribution in Bayesian inference. However, HMC techniques are computationally demanding for Bayesian neural networks due to the high dimensionality of the network’s parameter space and the non-convexity of their posterior distributions. Therefore, various approximation techniques, such as variational inference (VI) or stochastic gradient MCMC, are often employed to infer the posterior distribution of the network parameters. Such approximations introduce inaccuracies in the inferred distributions, resulting in unreliable uncertainty estimates. In this work, we propose a hybrid approach that combines inexpensive VI and accurate HMC methods to efficiently and accurately quantify uncertainties in neural networks and neural operators. The proposed approach leverages an initial VI training on the full network. We examine the influence of individual parameters on the prediction uncertainty, which shows that a large proportion of the parameters do not contribute substantially to uncertainty in the network predictions. This information is then used to significantly reduce the dimension of the parameter space, and HMC is performed only for the subset of network parameters that strongly influence prediction uncertainties. This yields a framework for accelerating the full batch HMC for posterior inference in neural networks. We demonstrate the efficiency and accuracy of the proposed framework on deep neural networks and operator networks, showing that inference can be performed for large networks with tens to hundreds of thousands of parameters. Finally, we show that this method can effectively learn surrogates for complex physical systems by modeling the operator that maps from upstream conditions to wall-pressure data on a cone in hypersonic flow.

Bayesian inference↗

Solar flares, the lunar surface, and gas rich meteorites

Investigations on the Fe-group nuclei track density vs depth in lunar rocks and Surveyor 3 TV camera filter glass were critically examined considering more factors than previously. The analysis gives a firmer basis to the observation of the preferential leakage of low energy Fe nuclei from the accelerating region of the sun. The track density gradients in lunar rock 12022 and filter glass are used to determine the lunar erosion rate of 3 angstroms/yr. Track gradients are less steep than predicted from energy spectrum observed in the Surveyor glass, perhaps due to sputtering. High densities of etchable tracks were found at all depths down to 60 cm in fines from Apollo cores and also in thin sections of the Pesjanoe, Pantar, and Fayetteville gas-rich meteorites. It is felt unlikely that suprathermal heavy ions were responsible for the high track densities.

Barber, D. J.↗

Motion of the Heliospheric Termination Shock at High Heliographic Latitude

We expect the mean distance of the heliospheric termination shock to be greater (smaller) at polar latitudes than at equatorial latitudes, depending on whether the mean dynamic pressure of the solar wind is greater or smaller at high latitudes. The heliospheric termination shock is expected to move in response to variation in upstream solar wind conditions, so that at any particular instant the termination shock will resemble a distorted asymmetric balloon with some parts moving inward and others moving outward. If the shock is a gasdynamic or magnetohydrodynamic shock the results of the analysis depend only very weakly on the nature of the upstream disturbance; typical speeds of the disturbed shock are approximately 100 to 200 km/s. In the absence of a significant latitude gradient of the typical magnitude of solar wind disturbances typical motions of the disturbed shock at polar latitudes would be about twice as fast, due to the higher speed of the high-latitude wind. If the dynamics of the termination shock are dominated by acceleration of the anomalous component of the cosmic rays, the motion of the shock in response to a given disturbance is substantially slower than in the gasdynamic case. Conceivably, particle acceleration might be a less important effect at higher latitudes, and we envision the possibility of a termination shock that is dominated by particle acceleration at lower latitudes and is an MHD shock at high latitudes. In this event high latitude solar wind disturbances would produce substantially larger inward and outward motions of the shock in the polar regions.

Barnes, Aaron↗

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗

Surface Cracking and Interface Reaction Associated Delamination Failure of Thermal and Environmental Barrier

In this paper, surface cracking and interface reactions of a ZrO2-8wt%Y2O3 and mullite/BSAS/Si thermal and environmental barrier coating system on SiC/SiC ceramic matrix composites were characterized after long-term combined laser thermal gradient and furnace cyclic tests in a water vapor containing environment. The surface cracking was analyzed based on the coating thermal gradient sintering behavior and thermal expansion mismatch stress characteristics under the thermal cyclic conditions. The interface reactions that were largely enhanced by the coating surface cracking in the water vapor environment were investigated in detail, and the reaction phases were identified for the coating system after the long-term exposure. The accelerated coating delamination failure was attributed to the increased delamination driving force under the thermal gradient cyclic loading and the reduced interface adhesion due to the detrimental interface reactions. The coating design issues will also be discussed based on the observed failure mechanisms under the high-heat-flux test conditions.

Zhu, Dongming↗

Conventional and Bimodal Nuclear Thermal Rocket (NTR) Artificial Gravity Mars Transfer Vehicle Concepts

A variety of countermeasures have been developed to address the debilitating physiological effects of zero-gravity (0-g) experienced by cosmonauts and astronauts during their approximately 0.5 to 1.2 year long stays in low Earth orbit (LEO). Longer interplanetary flights, combined with possible prolonged stays in Mars orbit, could subject crewmembers to up to approximately 2.5 years of weightlessness. In view of known and recently diagnosed problems associated with 0-g, an artificial gravity (AG) spacecraft offers many advantages and may indeed be an enabling technology for human flights to Mars. A number of important human factors must be taken into account in selecting the rotation radius, rotation rate, and orientation of the habitation module or modules. These factors include the gravity gradient effect, radial and tangential Coriolis forces, along with cross-coupled acceleration effects. Artificial gravity Mars transfer vehicle (MTV) concepts are presented that utilize both conventional NTR, as well as, enhanced bimodal nuclear thermal rocket (BNTR) propulsion. The NTR is a proven technology that generates high thrust and has a specific impulse (Isp) capability of approximately 900 s-twice that of today's best chemical rockets. The AG/MTV concepts using conventional Nuclear Thermal Propulsion (NTP) carry twin cylindrical International Space Station (ISS)- type habitation modules with their long axes oriented either perpendicular or parallel to the longitudinal spin axis of the MTV and utilize photovoltaic arrays (PVAs) for spacecraft power. The twin habitat modules are connected to a central operations hub located at the front of the MTV via two pressurized tunnels that provide the rotation radius for the habitat modules. For the BNTR AG/MTV option, each engine has its own closed secondary helium(He)-xenon (Xe) gas loop and Brayton Rotating Unit (BRU) that can generate 10s of kilowatts (kWe) of spacecraft electrical power during the mission coast phase eliminating the need for large PVAs. A single inflatable TransHab-type habitation module is also used with multiple vertical floors oriented radial to the MTV spin axis. The BNTR MTV's geometry-long and linear-is naturally compatible with AG operation. By rotating the vehicle about its center-of-mass (CM) and perpendicular to its flight vector at approximately 3.0 to 5.2 rpm, a centrifugal force and AG environment corresponding to approximately 0.38 to 1.0 g can be established to help maintain crew fitness out to Mars and back. Vehicles using NTP/ Bimodal Nuclear Thermal Propulsion (BNTP) can more readily accommodate the heavier payload mass and increased RCS propellant loading associated with AG operation, and can travel faster to and from Mars thereby reducing the crew's exposure to galactic cosmic radiation and solar flares. Mission scenario descriptions, key vehicle features and operational characteristics for each propulsion option are presented using the lift capability and payload volumes estimated for the Space Launch System (SLS)-1B and followon Heavy Lift Vehicle (HLV).

Nuclear Propulsion↗

Advancements in HF-free bipolar pulsed electropolishing for next-generation superconducting cavities

Hydrofluoric acid (HF)-free bipolar pulsed electropolishing (BPEP) offers an environmentally sustainable alternative to conventional Buffered Chemical Polishing (BCP) and Electropolishing (EP) techniques for superconducting radiofrequency (SRF) cavities. Recent studies at Jefferson Lab have demonstrated that a single-cell 1.3 GHz niobium cavity processed using HF-free BPEP achieved an accelerating gradient Eacc of 35 MV/m with a quality factor Q0 of 1E10 at 2 K, following extensive research and optimization. This talk will highlight the challenges encountered in developing this technique, key insights gained from experimental studies, and ongoing efforts to enhance its capabilities. In particular, we will explore its potential for refining Nb3Sn-coated niobium cavities via vapor diffusion techniques and for electroplating Nb3Sn films onto various cavity substrates, contributing to the advancement of high-performance SRF systems.

Tian, Hui [Thomas Jefferson National Accelerator F↗

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗

Field, frequency, and temperature dependencies of the surface resistance of nitrogen diffused niobium superconducting radio frequency cavities

We investigate the rf performance of several single-cell superconducting radio-frequency cavities subjected to low temperature heat treatment in nitrogen environment. The cavities were treated at temperature 120 – 165 ° C for an extended period of time (24–48 h) either in high vacuum or in a low partial pressure of ultrapure nitrogen. The improvement in Q 0 with a Q rise was observed when nitrogen gas was injected at ∼ 300 ° C during the cavity cooldown from 800 ° C and held at 165 ° C , without any degradation in accelerating gradient over the baseline performance. The treatment was applied to several elliptical cavities with frequency ranging from 0.75 to 3.0 GHz, showing an improved quality factor as a result of low temperature nitrogen treatments. The Q rise feature is similar to that achieved by nitrogen alloying Nb cavities at higher temperature, followed by material removal by electropolishing. The surface modification was confirmed by the change in electronic mean free path and tuned with the temperature and duration of heat treatment. The decrease of the temperature-dependent surface resistance with increasing rf field, resulting in a Q rise, becomes stronger with increasing frequency and decreasing temperature. The data suggest a crossover frequency of ∼ 0.95 GHz above that the Q rise phenomenon occurs at 2 K. Some of these results can be explained qualitatively with an existing model of intrinsic field-dependence of the surface resistance with both equilibrium and nonequilibrium quasiparticle distribution functions. The change in the Q slope below 0.95 GHz may result from masking contribution of trapped magnetic flux to the residual surface resistance. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Conventional and Bimodal Nuclear Thermal Rocket (NTR) Artificial Gravity Mars Transfer Vehicle Concepts

A variety of countermeasures have been developed to address the debilitating physiological effects of "zero-gravity" (0-g) experienced by cosmonauts and astronauts during their approximately 0.5-1.2 year long stays in LEO (Low Earth Orbit). Longer interplanetary flights, combined with possible prolonged stays in Mars orbit, could subject crewmembers to up to approximately 2.5 years of weightlessness. In view of known and recently diagnosed problems associated with 0-g, an artificial gravity spacecraft offers many advantages and may indeed be an enabling technology for human flights to Mars. A number of important human factors must be taken into account in selecting the rotation radius, rotation rate, and orientation of the habitation module or modules. These factors include the gravity gradient effect, radial and tangential Coriolis forces, along with cross-coupled acceleration effects. Artificial gravity (AG) Mars transfer vehicle (MTV) concepts are presented that utilize both conventional NTR, as well as, enhanced "bimodal" nuclear thermal rocket (BNTR) propulsion. The NTR is a proven technology that generates high thrust and has a specific impulse (I (sub sp)) capability of approximately 900 s - twice that of today's best chemical rockets. The AG/MTV concepts using conventional NTP carry twin cylindrical "ISS-type" habitation modules with their long axes oriented either perpendicular or parallel to the longitudinal spin axis of the MTV and utilize photovoltaic arrays (PVAs) for spacecraft power. The twin habitat modules are connected to a central operations hub located at the front of the MTV via two pressurized tunnels that provide the rotation radius for the habitat modules. For the BNTR AG/MTV option, each engine has its own "closed" secondary helium-xenon gas loop and Brayton rotating unit that can generate tens of kilowatts (kW (sub e)) of spacecraft electrical power during the mission coast phase eliminating the need for large PVAs. A single inflatable "TransHab-type" habitation module is also used with multiple vertical floors oriented radial to the MTV spin axis. The BNTR MTV's geometry - long and linear - is naturally compatible with AG operation. By rotating the vehicle about its center-of-mass and perpendicular to its flight vector at approximately 3.0 - 5.2 rpm, a centrifugal force and AG environment corresponding to approximately 0.38 - 1.0 g can be established to help maintain crew fitness out to Mars and back. Vehicles using NTP/BNTP can more readily accommodate the heavier payload mass and increased RCS propellant loading associated with AG operation, and can travel faster to and from Mars thereby reducing the crew's exposure to galactic cosmic radiation and solar flares. Mission scenario descriptions, key vehicle features and operational characteristics for each propulsion options are presented using the lift capability and payload volumes estimated for the SLS-1A and HLV.

Nuclear Propulsion↗

Reducing field emission in superconducting cavities with ozone

Contaminations on the active surface of high-power microwave resonators, such as radio-frequency accelerator cavities, can greatly enlarge electron field emission and reduce their performance. In this work, we demonstrate the use of ozone to clean the surface of niobium superconducting radio-frequency cavities resulting in reduced field emission and enhanced accelerating gradient before onset of field emission. We report results from processing four cavities.

43 PARTICLE ACCELERATORS↗

STS-40 orbital acceleration research experiment flight results during a typical sleep period

The Orbital Acceleration Research Experiment (OARE), an electrostatic accelerometer package with complete on-orbit calibration capabilities, was flown for the first time aboard the Space Shuttle on STS-40. This is also the first time an accelerometer package with nano-g sensitivity and a calibration facility has flown aboard the Space Shuttle. The instrument is designed to measure and record the Space Shuttle aerodynamic acceleration environment from the free molecule flow regime through the rarified flow transition into the hypersonic continuum regime. Because of its sensitivity, the OARE instrument defects aerodynamic behavior of the Space Shuttle while in low-earth orbit. A 2-hour orbital time period on day seven of the mission, when the crew was asleep and other spacecraft activities were at a minimum, was examined. During the flight, a 'trimmed-mean' filter was used to produce high quality, low frequency data which was successfully stored aboard the Space Shuttle in the OARE data storage system. Initial review of the data indicated that, although the expected precision was achieved, some equipment problems occurred resulting in uncertain accuracy. An acceleration model which includes aerodynamic, gravity-gradient, and rotational effects was constructed and compared with flight data. Examination of the model with the flight data shows the instrument to be sensitive to all major expected low frequency acceleration phenomena; however, some erratic instrument bias behavior persists in two axes. In these axes, the OARE data can be made to match a comprehensive atmospheric-aerodynamic model by making bias adjustments and slight linear corrections for drift. The other axis does not exhibit these difficulties and gives good agreement with the acceleration model.

Blanchard, R. C.↗

Progress and challenges in traveling-wave (TW) SRF cavity

Traveling-wave (TW) technology can push the accelerator field gradient of niobium SRF cavity to 70MV/m or higher beyond the limit of 50~60MV/m in Standing-wave (SW) technology. The early stages of TW SRF cavity developments had been funded by several SBIR grants to Euclid Techlabs and completed in collaboration with Fermilab through a 1-cell prototype and a proof-of-principle 3-cell TW cavity. The TW resonance excitation in the 3-cell TW cavity at 2K was demonstrated through the low power RF test in early 2024. A high-power test of the 3-cell in TW mode being prepared. To advance a design and technology to fabricate a novel high gradient TW SRF cavity, FNAL proposed a half-meter TW RF design and R&Ds to realize that are in progress. Here we will report the recent progress in the 3-cell TW cavity and the challenges towards a half-meter scale TW cavity.

43 PARTICLE ACCELERATORS↗

A Parallel Particle Swarm Optimization Algorithm Accelerated by Asynchronous Evaluations

A parallel Particle Swarm Optimization (PSO) algorithm is presented. Particle swarm optimization is a fairly recent addition to the family of non-gradient based, probabilistic search algorithms that is based on a simplified social model and is closely tied to swarming theory. Although PSO algorithms present several attractive properties to the designer, they are plagued by high computational cost as measured by elapsed time. One approach to reduce the elapsed time is to make use of coarse-grained parallelization to evaluate the design points. Previous parallel PSO algorithms were mostly implemented in a synchronous manner, where all design points within a design iteration are evaluated before the next iteration is started. This approach leads to poor parallel speedup in cases where a heterogeneous parallel environment is used and/or where the analysis time depends on the design point being analyzed. This paper introduces an asynchronous parallel PSO algorithm that greatly improves the parallel e ciency. The asynchronous algorithm is benchmarked on a cluster assembled of Apple Macintosh G5 desktop computers, using the multi-disciplinary optimization of a typical transport aircraft wing as an example.

Venter, Gerhard↗

Effects of gravity on sheared and nonsheared turbulent nonpremixed flames

The present numerical study is concerned with the fundamental physics of the multiway interaction between turbulence, chemical reaction, and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous, three-dimensional governing equations. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration. The effects of buoyancy include the generation of local shear, baroclinic production or destruction of vorticity, and countergradient heat and mass transport. Increased vorticity and small-scale turbulence promote further mixing and reaction. However, if the strain-rates become too high, local flame extinction can occur. Our objective is to gain an understanding of the complex interactions between the physical phenomena involved, with particular attention to the effects of buoyancy on the turbulence structure, flame behavior, and factors influencing flame extinction.

Elghobashi, Said↗

An Estimation of the Flying Qualities of the Kaiser Fleetwing All-Wing Airplane from Tests of a 1/7-Scale Model, TED No. NACA 2340

An investigation of a 1/7-scale powered model of the Kaiser Fleetwing all-wing airplane was made in the Langley full-scale tunnel to provide data for an estimation of the flying qualities of the airplane. The analysis of the stability and control characteristics of the airplane has been made as closely as possible in accordance with the requirements of the Bureau of Aeronautics, Navy Department's specifications, and a summary of the more significant conclusions is presented as follows. With the normal center of gravity located at 20 percent of the mean aerodynamic chord, the airplane will have adequate static longitudinal stability, elevator fixed, for all flight conditions except for low-power operation at low speeds where the stability will be about neutral. There will not be sufficient down-elevator deflection available for trim above speeds of about 130 miles per hour. It is probable that the reduction in the up-elevator deflections required for trim will be accompanied by reduced elevator hinge moments for low-power operation at low flight speeds. The static directional stability for this airplane will be low for all rudder-fixed or rudder-free flight conditions. The maximum rudder deflection of 30 deg will trim only about 15 deg yaw for most flight conditions and only 10 deg yaw for the condition with low power at low speeds. Also, at low powers and low speeds, it is estimated that the rudders will not trim the total adverse yaw resulting from an abrupt aileron roll using maximum aileron deflection. The airplane will meet the requirements for stability and control for asymmetric power operation with one outboard engine inoperative. The airplane would have no tendency for directional divergence but would probably be spirally unstable, with rudders fixed. The static lateral stability of the airplane will probably be about neutral for the high-speed flight conditions and will be only slightly increased for the low-power operation in low-speed flight. The airplane will not roll against the ailerons in a side-slip maneuver. Although the airplane would probably meet the minimum requirements of pb/2V of 0.07 at all speeds, there will be a loss in rolling ability of the airplane at high aileron deflections and at low flight speeds. It is estimated that the wing stall will be a gradual movement forward from the trailing edge and will be accompanied by no sudden pitching or rolling accelerations. Some stall warning may be indicated by reduction in the elevator and aileron force gradients and by the shaking of the controls caused by unsteady flow over the surfaces near the stall.

Brewer, Gerald W.↗