A biogeography-based optimization algorithm with mutation strategies for model parameter estimation of solar and fuel cells

Mathematical models are useful tools for simulation, evaluation, optimal operation and control of solar cells and proton exchange membrane fuel cells (PEMFCs). To identify the model parameters of these two type of cells efficiently, a biogeography-based optimization algorithm with mutation strategies (BBO-M) is proposed. The BBO-M uses the structure of biogeography-based optimization algorithm (BBO), and both the mutation motivated from the differential evolution (DE) algorithm and the chaos theory are incorporated into the BBO structure for improving the global searching capability of the algorithm. Numerical experiments have been conducted on ten benchmark functions with 50 dimensions, and the results show that BBO-M can produce solutions of high quality and has fast convergence rate. Then, the proposed BBO-M is applied to the model parameter estimation of the two type of cells. The experimental results clearly demonstrate the power of the proposed BBO-M in estimating model parameters of both solar and fuel cells.

An improved TLBO with elite strategy for parameters identification of PEM fuel cell and solar cell models

Clean and renewable energy generation and supply has drawn much attention worldwide in recent years, the proton exchange membrane (PEM) fuel cells and solar cells are among the most popular technologies. Accurately modeling the PEM fuel cells as well as solar cells is critical in their applications, and this involves the identification and optimization of model parameters. This is however challenging due to the highly nonlinear and complex nature of the models. In particular for PEM fuel cells, the model has to be optimized under different operation conditions, thus making the solution space extremely complex. In this paper, an improved and simplified teaching-learning based optimization algorithm (STLBO) is proposed to identify and optimize parameters for these two types of cell models. This is achieved by introducing an elite strategy to improve the quality of population and a local search is employed to further enhance the performance of the global best solution. To improve the diversity of the local search a chaotic map is also introduced. Compared with the basic TLBO, the structure of the proposed algorithm is much simplified and the searching ability is significantly enhanced. The performance of the proposed STLBO is firstly tested and verified on two low dimension decomposable problems and twelve large scale benchmark functions, then on the parameter identification of PEM fuel cell as well as solar cell models. Intensive experimental simulations show that the proposed STLBO exhibits excellent performance in terms of the accuracy and speed, in comparison with those reported in the literature.

Has the time for renewables finally arrived? – An event study on EPA’s Clean Power Plan

The recent proposals presented by EPA aimed to reduce the dependency of fossil fuels and to lower current emissions levels, hoping to gradually shift electric generation units to renewable energy sources. Actually, the Final Rule Proposal announcement day exhibited a negative Abnormal Return on Fossil Fuels but the following days had positive Abnormal Returns, mostly due to legislative change perceived by financial markets which eased up implementation periods of the proposed measures in the Final Rule when compared to the Draft Rule. Oppositely, Renewables and Solar Portfolios exhibited negative Cumulative Abnormal Returns over the period surrounding the Final Rule.

Temperature- and Pressure-induced Proton Transfer in the 1:1 Adduct Formed between Squaric Acid and 4,4 '-Bipyridine

We have applied a combination of spectroscopic and diffraction methods to study the adduct formed between squaric acid and bypridine, which has been postulated to exhibit proton transfer associated with a single-crystal to single-crystal phase transition at ca. 450 K. A combination of X-ray single-crystal and very-high flux powder neutron diffraction data confirmed that a proton does transfer from the acid to the base in the high-temperature form. Powder X-ray diffraction measurements demonstrated that the transition was reversible but that a significant kinetic energy barrier must be overcome to revert to the original structure. Computational modeling is consistent with these results. Modeling also revealed that, while the proton transfer event would be strongly discouraged in the gas phase, it occurs in the solid state due to the increase in charge state of the molecular ions and their arrangement inside the lattice. The color change is attributed to a narrowing of the squaric acid to bipyridine charge-transfer energy gap. Finally, evidence for the possible existence of two further phases at high pressure is also presented.

Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.

On the evolution of the solar wind between 1 and 5 AU at the time of the Cassini Jupiter flyby: Multispacecraft observations of interplanetary coronal mass ejections including the formation of a merged interaction region

The Cassini flyby of Jupiter occurred at a time near solar maximum. Consequently, the pre-Jupiter data set reveals clear and numerous transient perturbations to the Parker Spiral solar wind structure. Limited plasma data are available at Cassini for this period due to pointing restrictions imposed on the instrument. This renders the identification of the nature of such structures ambiguous, as determinations based on the magnetic field data alone are unreliable. However, a fortuitous alignment of the planets during this encounter allowed us to trace these structures back to those observed previously by the Wind spacecraft near the Earth. Of the phenomena that we are satisfactorily able to trace back to their manifestation at 1 AU, two are identified as being due to interplanetary coronal mass ejections. One event at Cassini is shown to be a merged interaction region, which is formed from the compression of a magnetic cloud by two anomalously fast solar wind streams. The flux-rope structure associated with this magnetic cloud is not as apparent at Cassini and has most likely been compressed and deformed. Confirmation of the validity of the ballistic projections used here is provided by results obtained from a one-dimensional magnetohydrodynamic projection of solar wind parameters measured upstream near the Earth. It is found that when the Earth and Cassini are within a few tens of degrees in heliospheric longitude, the results of this one-dimensional model predict the actual conditions measured at 5 AU to an impressive degree. Finally, the validity of the use of such one-dimensional projections in obtaining quasi-solar wind parameters at the outer planets is discussed.

A solar storm observed from the Sun to Venus using the STEREO, Venus Express, and MESSENGER spacecraft

The suite of SECCHI optical imaging instruments on the STEREO-A spacecraft is used to track a solar storm, consisting of several coronal mass ejections (CMEs) and other coronal loops, as it propagates from the Sun into the heliosphere during May 2007. The 3-D propagation path of the largest interplanetary CME (ICME) is determined from the observations made by the SECCHI Heliospheric Imager (HI) on STEREO-A (HI-1/2A). Two parts of the CME are tracked through the SECCHI images, a bright loop and a V-shaped feature located at the rear of the event. We show that these two structures could be the result of line-of-sight integration of the light scattered by electrons located on a single flux rope. In addition to being imaged by HI, the CME is observed simultaneously by the plasma and magnetic field experiments on the Venus Express and MESSENGER spacecraft. The imaged loop and V-shaped structure bound, as expected, the flux rope observed in situ. The SECCHI images reveal that the leading loop-like structure propagated faster than the V-shaped structure, and a decrease in in situ CME speed occurred during the passage of the flux rope.We interpret this as the result of the continuous radial expansion of the flux rope as it progressed outward through the interplanetary medium. An expansion speed in the radial direction of ~30 km s-1 is obtained directly from the SECCHI-HI images and is in agreement with the difference in speed of the two structures observed in situ. This paper shows that the flux rope location can be determined from white light images, which could have important space weather applications.

Numerical simulation of the 12 May 1997 CME Event: The role of magnetic reconnection

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the 12 May 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance (≈1R ) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

Magnetic discontinuities in the near-earth solar wind: evidence of in-transit turbulence or remnants of coronal structure?

Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.

Solar cycle 24: implications for energetic particles and long-term space climate change

The recent solar minimum was the longest and deepest of the space age, with the lowest average sunspot numbers for nearly a century. The Sun appears to be exiting a grand solar maximum (GSM) of activity which has persisted throughout the space age, and is headed into a significantly quieter period. Indeed, initial observations of solar cycle 24 (SC24) continue to show a relatively low heliospheric magnetic field strength and sunspot number (R), despite the average latitude of sunspots and the inclination of the heliospheric current sheet showing the rise to solar maximum is well underway. We extrapolate the available SC24 observations forward in time by assuming R will continue to follow a similar form to previous cycles, despite the end of the GSM, and predict a very weak cycle 24, with R peaking at ∼65–75 around the middle/end of 2012. Similarly, we estimate the heliospheric magnetic field strength will peak around 6nT. We estimate that average galactic cosmic ray fluxes above 1GV rigidity will be ∼10% higher in SC24 than SC23 and that the probability of a large SEP event during this cycle is 0.8, compared to 0.5 for SC23. Comparison of the SC24 R estimates with previous ends of GSMs inferred from 9300 years of cosmogenic isotope data places the current evolution of the Sun and heliosphere in the lowest 5% of cases, suggesting Maunder Minimum conditions are likely within the next 40 years.

A 27 day persistence model of near-Earth solar wind conditions: a long lead-time forecast and a benchmark for dynamical models

Geomagnetic activity has long been known to exhibit approximately 27 day periodicity, resulting from solar wind structures repeating each solar rotation. Thus a very simple near-Earth solar wind forecast is 27 day persistence, wherein the near-Earth solar wind conditions today are assumed to be identical to those 27 days previously. Effective use of such a persistence model as a forecast tool, however, requires the performance and uncertainty to be fully characterized. The first half of this study determines which solar wind parameters can be reliably forecast by persistence and how the forecast skill varies with the solar cycle. The second half of the study shows how persistence can provide a useful benchmark for more sophisticated forecast schemes, namely physics-based numerical models. Point-by-point assessment methods, such as correlation and mean-square error, find persistence skill comparable to numerical models during solar minimum, despite the 27 day lead time of persistence forecasts, versus 2–5 days for numerical schemes. At solar maximum, however, the dynamic nature of the corona means 27 day persistence is no longer a good approximation and skill scores suggest persistence is out-performed by numerical models for almost all solar wind parameters. But point-by-point assessment techniques are not always a reliable indicator of usefulness as a forecast tool. An event-based assessment method, which focusses key solar wind structures, finds persistence to be the most valuable forecast throughout the solar cycle. This reiterates the fact that the means of assessing the “best” forecast model must be specifically tailored to its intended use.

Evidence for solar wind modulation of lightning

The response of lightning rates over Europe to arrival of high speed solar wind streams at Earth is investigated using a superposed epoch analysis. Fast solar wind stream arrival is determined from modulation of the solar wind V y component, measured by the Advanced Composition Explorer spacecraft. Lightning rate changes around these event times are determined from the very low frequency arrival time difference (ATD) system of the UK Met Office. Arrival of high speed streams at Earth is found to be preceded by a decrease in total solar irradiance and an increase in sunspot number and Mg II emissions. These are consistent with the high speed stream's source being co-located with an active region appearing on the Eastern solar limb and rotating at the 27 d period of the Sun. Arrival of the high speed stream at Earth also coincides with a small (~1%) but rapid decrease in galactic cosmic ray flux, a moderate (~6%) increase in lower energy solar energetic protons (SEPs), and a substantial, statistically significant increase in lightning rates. These changes persist for around 40 d in all three quantities. The lightning rate increase is corroborated by an increase in the total number of thunder days observed by UK Met stations, again persisting for around 40 d after the arrival of a high speed solar wind stream. This result appears to contradict earlier studies that found an anti-correlation between sunspot number and thunder days over solar cycle timescales. The increase in lightning rates and thunder days that we observe coincides with an increased flux of SEPs which, while not being detected at ground level, nevertheless penetrate the atmosphere to tropospheric altitudes. This effect could be further amplified by an increase in mean lightning stroke intensity that brings more strokes above the detection threshold of the ATD system. In order to remove any potential seasonal bias the analysis was repeated for daily solar wind triggers occurring during the summer months (June to August). Though this reduced the number of solar wind triggers to 32, the response in both lightning and thunder day data remained statistically significant. This modulation of lightning by regular and predictable solar wind events may be beneficial to medium range forecasting of hazardous weather.

Detection of lower tropospheric responses to solar energetic particles at midlatitudes

Solar energetic particles (SEPs) occasionally contribute additional atmospheric ionization beyond that arising from the usual galactic cosmic ray background. During an SEP event associated with a solar flare on April 11, 2013, the vertical ionization rate profile obtained using a balloon-borne detector showed enhanced ionization with a 26% increase at 20 km, over Reading, United Kingdom. Fluctuations in atmospheric electrical parameters were also detected at the surface, beneath the balloon’s trajectory. As no coincident changes in geomagnetism occurred, the electrical fluctuations are very likely to be associated with increased ionization, as observed by the balloon measurements. The lack of response of surface neutron monitors during this event indicates that energetic particles that are not detected at the surface by neutron monitors can nevertheless enter and influence the atmosphere’s weather-generating regions.

The Solar Stormwatch CME catalogue: results from the first space weather citizen science project

Solar Stormwatch was the first space weather citizen science project, the aim of which was to identify and track coronal mass ejections (CMEs) observed by the Heliospheric Imagers aboard the STEREO satellites. The project has now been running for approximately 4 years, with input from >16000 citizen scientists, resulting in a dataset of >38000 time-elongation profiles of CME trajectories, observed over 18 pre-selected position angles. We present our method for reducing this data set into aCME catalogue. The resulting catalogue consists of 144 CMEs over the period January-2007 to February-2010, of which 110 were observed by STEREO-A and 77 were observed by STEREO-B. For each CME, the time-elongation profiles generated by the citizen scientists are averaged into a consensus profile along each position angle that the event was tracked. We consider this catalogue to be unique, being at present the only citizen science generated CME catalogue, tracking CMEs over an elongation range of 4 degrees out to a maximum of approximately 70 degrees. Using single spacecraft fitting techniques, we estimate the speed, direction, solar source region and latitudinal width of each CME. This shows that, at present, the Solar Stormwatch catalogue (which covers only solar minimum years) contains almost exclusively slow CMEs, with a mean speed of approximately 350 kms−1. The full catalogue is available for public access at www.met.reading.ac.uk/spate/stormwatch. This includes, for each event, the unprocessed time-elongation profiles generated by Solar Stormwatch, the consensus time-elongation profiles and a set of summary plots, as well as the estimated CME properties.

Solar wind-magnetosphere coupling functions on timescales of 1 day to 1 year

There are no direct observational methods for determining the total rate at which energy is extracted from the solar wind by the magnetosphere. In the absence of such a direct measurement, alternative means of estimating the energy available to drive the magnetospheric system have been developed using different ionospheric and magnetospheric indices as proxies for energy consumption and dissipation and thus the input. The so-called coupling functions are constructed from the parameters of the interplanetary medium, as either theoretical or empirical estimates of energy transfer, and the effectiveness of these coupling functions has been evaluated in terms of their correlation with the chosen index. A number of coupling functions have been studied in the past with various criteria governing event selection and timescale. The present paper contains an exhaustive survey of the correlation between geomagnetic activity and the near-Earth solar wind and two of the planetary indices at a wide variety of timescales. Various combinations of interplanetary parameters are evaluated with careful allowance for the effects of data gaps in the interplanetary data. We show that the theoretical coupling, P�, function first proposed by Vasyliunas et al. is superior at all timescales from 1-day to 1-year.