Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Time”

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,711 records · Page 95

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗

Second-order Kinetics of DTPA and Plutonium in Rat Plasma

We report that in 2008, Serandour et al. reported on their in vitro experiment involving rat plasma samples obtained after an intravenous intake of plutonium citrate. Different amounts of DTPA were added to the plasma samples and the percentage of low-molecular-weight plutonium measured. Only when the DTPA dosage was three orders of magnitude greater than the recommended 30 μmol/kg was 100% of the plutonium apparently in the form of chelate. These data were modeled assuming three competing chemical reactions with other molecules that bind with plutonium. Here, time-dependent second-order kinetics of these reactions are calculated, intended eventually to become part of a complete biokinetic model of DTPA action on actinides in laboratory animals or humans. The probability distribution of the ratio of stability constants for the reactants was calculated using Markov Chain Monte Carlo. In conclusion, these calculations substantiate that the inclusion of more reactions is needed in order to be in agreement with known stability constants.

61 RADIATION PROTECTION AND DOSIMETRY↗

Simplified Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two time step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting rates of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx are obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

Conditions and mechanisms for the formation of nano-sized Delafossite (CuFeO{sub 2}) at temperatures ≤90 °C in aqueous solution

In this study, we present the mechanism of CuFeO{sub 2} formation in aqueous solution at low temperatures ≤90 °C, using sulfate salts as reactants. Furthermore, we demonstrate the influence of experimental conditions (alkalization, reaction and ageing temperature and time) on the synthesized nanoparticles. In all cases, GR–SO{sub 4}, a Fe(II–III) layered double hydroxysulphate (Fe{sup 2+}{sub 4}Fe{sup 3+}{sub 2}(OH){sub 12}·SO{sub 4}) and Cu{sub 2}O precipitate first. During further OH{sup −} supply GR–SO{sub 4} oxidizes and forms Fe{sub 10}O{sub 14}(OH){sub 2}, Cu{sub 2}O and CuFeO{sub 2} crystals. Due to the high pH further CuFeO{sub 2} crystals grow at the cost of the unstable intermediate products. The reaction rate increases with increasing ageing temperature, reaction pH and, in particular, NaOH concentrations in the solution. As a result, highly crystalline CuFeO{sub 2} (3R and 2H polytypes) nanoparticles showing hexagonal morphology can be synthesized at 70 °C within 10 h or at 50 °C within 1 week. The formation of 2H polytype is favored by additional OH{sup −} supply during the pH-stat time and rather low temperatures. - Highlights: • We solve the formation mechanism of pure CuFeO{sub 2} using sulfates as reactants. • CuFeO{sub 2} nanoparticles crystallize on cost of green rust, Fe{sub 10}O{sub 14}(OH){sub 2} and Cu{sub 2}O. • The reaction rate increases with increasing temperature and OH- concentration. • CuFeO{sub 2} nanoparticles form at 50 °C within one week and at 70 °C within 10 h. • 2H-polytype of CuFeO{sub 2} is favored by additional NaOH supply during pH-stat-time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deactivation of Sn-Beta zeolites caused by structural transformation of hydrophobic to hydrophilic micropores during aqueous-phase glucose isomerization

The structural changes underlying the deactivation of Sn-Beta zeolites under aqueous-phase reaction conditions at elevated temperatures (373 K) are investigated using spectroscopic characterization and site titration techniques together with turnover rates for glucose isomerization, a well-understood probe reaction for which changes in measured rates can be ascribed to specific changes in catalyst structure. In the case of hydrophobic, low-defect Sn-Beta zeolites (Sn-Beta-F), treatment in hot liquid water (373 K) for short times (<1 h) prior to reaction causes glucose–fructose isomerization turnover rates (per open Sn site, 373 K) to increase, while longer-term exposure (>3 h) to hot liquid water causes turnover rates to decrease and approach values characteristic of hydrophilic, defect-rich Sn-Beta zeolites (Sn-Beta-OH). However, turnover rates on hydrophilic Sn-Beta-OH zeolites are insensitive to the duration of hot liquid water exposure prior to reaction. Activation and deactivation phenomena on Sn-Beta-F zeolites occur concomitantly with the formation of silanol defects (by ~2–10×) with increasing durations (0–24 h) of hot water treatment, despite negligible differences in open and closed Sn site speciation as quantified ex situ by CD 3 CN IR spectra. Mechanistic interpretations of these phenomena suggest that silanol groups present at low densities serve as binding sites for water molecules and clusters, which confer enthalpic stability to kinetically-relevant hydride-shift transition states and increase turnover rates, while silanol groups present in higher densities stabilize extended hydrogen-bonded water networks, which entropically destabilize kinetically-relevant transition states and decrease turnover rates. Intraporous voids within hydrophobic Sn-Beta-F zeolites become increasingly hydrophilic as silanol groups are formed by hydrolysis of framework siloxane bridges with increasing durations of water treatment, thereby decreasing aqueous-phase glucose isomerization turnover rates (per open Sn site). These findings imply design strategies that suppress framework hydrolysis would attenuate the deactivation of Lewis acid zeolites in aqueous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Efficient Electrochemical Tanks-in-Series Model for Lithium Sulfur Batteries

This article applies and efficiently implements the Tanks-in-Series methodology ( J. Electrochem. Soc. , 167 , 013534 (2020)) to generate a computationally efficient electrochemical model for Lithium-Sulfur batteries. The original Tank model approach for Lithium-ion batteries is modified to account for porosity changes with time. In addition, an exponential scaling method is introduced that enables efficient simulation of the model equations to address the wide range of time constants present for different reactions in the Lithium-Sulfur system. The Tank Model achieves acceptable voltage error even for transport-limited discharged conditions. Predictions of internal electrochemical variables are examined, and electrochemical implications of the approximations discussed. This suggests significant potential for real-time applications such as optimal charging, cell-balancing, and estimation, and represents a step forward in efforts to incorporate detailed electrochemical models in advanced Battery Management Systems for Lithium-Sulfur batteries.

Parke, Caitlin D. (ORCID:0000000177324044)↗

Computational study on the removal of photolabile protecting groups by photochemical reactions

The application and versatility of photolabile protecting groups (PPGs) have continued to expand since their inception almost 50 years ago. Here, in the present work we employ density functional theory (DFT) to characterize ground and excited state properties of several compounds with PPGs and to explore possible, controllable photodissociation reactions that remove these PPGs. The reaction mechanism is examined by analyzing the DFT results as well as dynamical simulations of nonadiabatic electron transfer reactions using multilayer multiconfiguration time-dependent Hartree theory. The photoexcited compounds appear to undergo ultrafast back electron transfer to the ground electronic state via internal conversion. For some compounds, the electronic-nuclear couplings facilitate sufficient photo energy transfer to the nuclear degrees of freedom that dominate the reaction coordinate, leading to the dissociation reaction that removes the PPGs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transitioning NWChem to the Next Generation of Manycore Machines

The NorthWest Chemistry (NWChem) modeling software is a popular molecular chemistry simulation software that was designed from the start to work on massively parallel processing supercomputers[6, 28, 49]. It contains an umbrella of modules that today includes Self Consistent Field (SCF), second order Mller-Plesset perturbation theory (MP2), Coupled Cluster, multi-conguration selfconsistent eld (MCSCF), selected conguration interaction (CI), tensor contraction engine (TCE) many body methods, density functional theory (DFT), time-dependent density functional theory (TDDFT), real time time-dependent density functional theory, pseudopotential plane-wave density functional theory (PSPW), band structure (BAND), ab initio molecular dynamics, Car-Parrinello molecular dynamics, classical molecular dynamics (MD), QM/MM, AIMD/MM, GIAO NMR, COSMO, COSMO-SMD, and RISM solvation models, free energy simulations, reaction path optimization, parallel in time, among other capabilities[ 22]. Moreover new capabilities continue to be added with each new release.

Environmental Molecular Sciences Laboratory↗

Pathways of anaerobic organic matter decomposition in tundra soils from Barrow, Alaska

Arctic tundra soils store a large quantity of organic carbon that is susceptible to decomposition and release to the atmosphere as methane (CH 4 ) and carbon dioxide (CO 2 ) under a warming climate. Anaerobic processes that generate CH 4 and CO 2 remain unclear because previous studies have focused on aerobic decomposition pathways. To predict releases of CO 2 and CH 4 from tundra soils, it is necessary to identify pathways of soil organic matter decomposition under the anoxic conditions that are prevalent in Arctic ecosystems. Here molecular and spectroscopic techniques were used to monitor biological degradation of water-extractable organic carbon (WEOC) during anoxic incubation of tundra soils from a region of continuous permafrost in northern Alaska. Organic and mineral soils from the tundra active layer were incubated at –2, +4, or +8°C for up to 60 days to mimic the short-term thaw season. Results suggest that, under anoxic conditions, fermentation converted complex organic molecules into simple organic acids that were used in concomitant Fe-reduction and acetoclastic methanogenesis reactions. Nonaromatic compounds increased over time as WEOC increased. Organic acid metabolites initially accumulated in soils but were mostly depleted by day 60 because organic acids were consumed to produce Fe(II), CO 2 , and CH 4 . We conclude that fermentation of nonprotected organic matter facilitates methanogenesis and Fe reduction reactions, and that the proportion of organic acids consumed by methanogenesis increases relative to Fe reduction with increasing temperature. As a result, the decomposition pathways observed in this study are important to consider in numerical modeling of greenhouse gas production in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Advances in NASA research on nuclear-pumped lasers

NASA has been primarily interested in nuclear-pumped lasers using the He-3 or U-235F6 reaction for lasant excitation. With He-3 excitation, a large volume, multiple-path He-3-Ar nuclear laser has produced an output of 1 kilowatt. Power deposition was shown to be homogeneous over this volume. The CO laser has been pumped for the first time using the He-3 reaction, producing approximately 200 Watts. Using a boron-10 coating to excite N2, nuclear lasing has been achieved in CO2 in a transfer laser configuration. Nuclear lasing of Ar-Xe has been demonstrated using fission fragment excitation from U-235F6. Research on the gas core reactor has resulted in a steady state operational power of 30 kilowatts with flowing U-235F6 in an argon vortex.

De Young, R. J.↗

Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions

Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗

The mechanisms of flame holding in the wake of a bluff body

The flame holding mechanism for lean methane- and lean propane-air flames is examined under conditions where the recirculation zone is absent. The main objective of this work is to study the holding process in detail in an attempt to determine the mechanism of flame holding and also the conditions where this mechanism is viable and when it fails and blow-off occurs. Inverted flames held in the wake of a flat strip were studied. Experiments with different sizes of flame holders were performed. The velocity flow field was determined using a laser Doppler velocimetry technique. Equation of continuity was used to calculate the flame temperature from the change in area of flow streamlines before and after the flame. Observations of the inverted flame itself were obtained using schlieren and direct photography. Results show that there are different mechanisms operative at the time of blow-off for lean propane and methane flames. Blow-off or extinction occurs for lean propane-air flame in spite of the reaction going to completion and the disparity between the heat loss and the gain in mass diffusion in the reaction zone i.e., Le 1.0 causes the flame to blow-off. For methane-air flame the controlling factor or blow-off is incomplete reaction due to higher blowing rate leading to reduced residence time in the reaction zone.

Strehlow, R. A.↗

Research on the photoluminescence and up-conversion luminescence properties of Y2Mo4O15: Yb, Ho under 454 and 980 nm excitation

Highlights: • In this paper, Y{sub 2}Mo{sub 4}O{sub 15}:Yb, Ho phosphors were synthesised using the conventional solid-state reaction method for the first time, and its photoluminescence and up-conversion emission properties were characterised for the first time under 454 and 980 nm excitation. The special red UCL mechanism was also analysed for the first time. • We observe an abnormal red UCL in Y{sub 2}Mo{sub 4}O{sub 15}:Yb, Ho under 980 nm excitation. • The reducing of decay time of {sup 5}I{sub 6} states is responsible for the red UCL. • Very weak {sup 5}I{sub 6} → {sup 5}I{sub 8} IR emission further confirm the short decay time of {sup 5}I{sub 6} in Y{sub 2}Mo{sub 4}O{sub 15}:Yb, Ho. - Abstract: Y{sub 1.8-x}Yb{sub 0.2}Ho{sub x}Mo{sub 4}O{sub 15} (x = 0.005, 0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 0.64) phosphors were synthesised at 700 °C by using a solid-state reaction method. The samples were characterised by X-ray diffraction, diffuse reflectance spectra, photoluminescence and up-conversion emission spectra. The photoluminescence spectra of Y{sub 2}Mo{sub 4}O{sub 15}:Yb, Ho shows extremely strong green down-shift luminescence corresponding to the {sup 5}F{sub 4}, {sup 5}S{sub 2} → {sup 5}I{sub 8} transition of Ho{sup 3+} ions with I{sub g}/I{sub r} > 20 under 454 nm excitation. While at 980 nm excitation, Y{sub 2}Mo{sub 4}O{sub 15}:Yb, Ho sample exhibits unusual red up-conversion luminescence along with a very weak green emission (I{sub r}/I{sub g} = 6.4), which is resulted from the severe non-radiative relaxation ({sup 5}I{sub 6} → {sup 5}I{sub 7}) process due to the very short decay time of {sup 5}I{sub 6} states of Ho{sup 3+} ions. The mechanism of down-shift luminescence is the electrical dipole – dipole transitions. The quenching concentration for down-shift luminescence is higher than that for up-conversion luminescence, which results from the effects of different luminescence processes and large distance between doped rare earth ions in Y{sub 2}Mo{sub 4}O{sub 15}. The different luminescence mechanism and color performance are also studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Minimum-Time and Vibration Avoidance Attitude Maneuver for Spacecraft with Torque and Momentum Limit Constraints in Redundant Reaction Wheel Configuration

In this paper, we present an optimal open-loop slew trajectory algorithm developed at GSFC for the so-called "Yardstick design" of the James Webb Space Telescope (JWST). JWST is an orbiting infrared observatory featuring a lightweight, segmented primary mirror approximately 6 meters in diameter and a sunshield approximately the size of a tennis court. This large, flexible structure will have significant number of lightly damped, dominant flexible modes. With very stringent requirements on pointing accuracy and image quality, it is important that slewing be done within the required time constraint and with minimal induced vibration in order to maximize observing efficiency. With reaction wheels as control actuators, initial wheel speeds as well as individual wheel torque and momentum limits become dominant constraints in slew performance. These constraints must be taken into account when performing slews to ensure that unexpected reaction wheel saturation does not occur, since such saturation leads to control failure in accurately tracking commanded motion and produces high frequency torque components capable of exciting structural modes. A minimum-time constraint is also included and coupled with reaction wheel limit constraints in the optimization to minimize both the effect of the control torque on the flexible body motion and the maneuver time. The optimization is on slew command parameters, such as maximum slew velocity and acceleration, for a given redundant reaction wheel configuration and is based on the dynamic interaction between the spacecraft and reaction wheel motion. Analytical development of the slew algorithm to generate desired slew position, rate, and acceleration profiles to command a feedback/feed forward control system is described. High-fidelity simulation and experimental results are presented to show that the developed slew law achieves the objectives.

Ha, Kong Q.↗

Constraining the Neutron Star Compactness: Extraction of the Al 23 ( p , γ ) Reaction Rate for the r p Process

The 23 Al(p,γ) 24 Si reaction is among the most important reactions driving the energy generation in type-I x-ray bursts. However, the present reaction-rate uncertainty limits constraints on neutron star properties that can be achieved with burst model-observation comparisons. In this work, we present a novel technique for constraining this important reaction by combining the GRETINA array with the neutron detector LENDA coupled to the S800 spectrograph at the National Superconducting Cyclotron Laboratory. The 23 Al(d,n) reaction was used to populate the astrophysically important states in 24 Si. This enables a measurement in complete kinematics for extracting all relevant inputs necessary to calculate the reaction rate. For the first time, a predicted close-lying doublet of a $2^{+}_{2}$ and ($4^{+}_{1}, 0^{+}_{2}$) state in 24 Si was disentangled, finally resolving conflicting results from two previous measurements. Furthermore, it was possible to extract spectroscopic factors using GRETINA and LENDA simultaneously. This new technique may be used to constrain other important reaction rates for various astrophysical scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Time-dependent 2.2 MeV and 0.5 MeV lines from solar flares

The time dependences of the 2.2 MeV and 0.51 MeV gamma ray lines from solar flares are calculated and the results are compared with observations of the 1972, August 4 and 7 flares. Time lag between the nuclear reactions and the formation of these two lines are caused, respectively, by capture of the neutrons, and by deceleration of the positrons and decay of the radioactive nuclei. Results show that the calculation is consistent with the observed rise of the 2.2 MeV line on August 4, and it does not require different time dependences for the accelerated protons and electrons in the flare region. The above lags can explain the delayed gamma ray emission observed on August 7. Positrons of energies greater than about 10 MeV could be detected in interplanetary space following large solar flares.

Wang, H. T.↗

Time-dependent 2.2-MeV and 0.5-MeV lines from solar flares

The time dependences of the 2.2- and 0.51-MeV gamma-ray lines from solar flares are calculated, and the results are compared with observations of the 1972 August 4 and 7 flares. The time lag between the nuclear reactions and the formation of these two lines is caused by capture of the neutrons and subsequent deceleration of the positrons and decay of the radioactive nuclei. Our main results are that the calculation is consistent with the observed rise of the 2.2-MeV line on August 4, and it does not require different time dependences for the accelerated protons and high-energy electrons in the flare region. The above lags can explain the delayed gamma-ray emission observed on August 7. Positrons of energies greater than about 10 MeV could be detected in interplanetary space following large solar flares.

Wang, H. T.↗