Real-time feasibility of nonlinear model predictive control for semi-batch reactors subject to uncertainty and disturbances

This paper presents two nonlinear model predictive control based methods for solving closed-loop stochastic dynamic optimisation problems, ensuring both robustness and feasibility with respect to state output constraints. The first one is a new deterministic approach, using the wait-and-see strategy. The key idea is to specifically anticipate violation of output hard-constraints, which are strongly affected by instantaneous disturbances, by backing off of their bounds along the moving horizon. The second method is a stochastic approach to solve nonlinear chance-constrained dynamic optimisation problems under uncertainties. The key aspect is the explicit consideration of the stochastic properties of both exogenous and endogenous uncertainties in the problem formulation (here-and-now strategy). The approach considers a nonlinear relation between uncertain inputs and the constrained state outputs. The performance of the proposed methodologies is assessed via an application to a semi-batch reactor under safety constraints, involving strongly exothermic reactions

Construction of global optimization constrained NLP test cases from unconstrained problems

This paper presents a novel construction technique for constrained nonconvex Nonlinear Programming Problem (NLP) test cases, derived from the evaluation tree structure of standardized bound constrained problems for which the global solution is known. It is demonstrated in a step-by-step procedure how first an equality constrained problem can be derived from an unconstrained one, with bounds imposed on all variables, using the Directed Acyclic Graph (DAG) of the unconstrained objective function and the use of interval arithmetic to derive bounds for the new variables introduced. An advantage of the proposed methodology is that several standard unconstrained global optimization test cases can be constructed for varying number of optimization variables, thus leading to adjustable size derived NLP’s. Further to this in a second step it is demonstrated how any subset of the equalities derived can be relaxed into inequalities giving an equivalent optimization problem. Finally, in a third step it is demonstrated how, by reducing the number of equality constraints derived, it is possible to obtain more complex expressions in the constraints and objective function. The methodology is highlighted throughout by motivating examples and a sample code in Mathematica TM is provided in the Appendix.This is the author accepted manuscript. The final version is available from Elsevier via https://doi.org/https://doi.org/10.1016/j.cherd.2016.03.01

ICRS-Filter: A randomized direct search algorithm for constrained nonconvex optimization problems

This work presents a novel algorithm and its implementation for the stochastic optimization of generally constrained Nonlinear Programming Problems (NLP). The basic algorithm adopted is the Iterated Control Random Search (ICRS) method of Casares and Banga (1987) with modifications such that random points are generated strictly within a bounding box defined by bounds on all variables. The ICRS algorithm serves as an initial point determination method for launching gradient-based methods that converge to the nearest local minimum. The issue of constraint handling is addressed in our work via the use of a filter based methodology, thus obviating the need for use of the penalty functions as in the basic ICRS method presented in Banga and Seider (1996),which handles only bound constrained problems. The proposed algorithm, termed ICRS-Filter, is shown to be very robust and reliable in producing very good or global solutions for most of the several case studies examined in this contribution.This is the author accepted manuscript. The final version is available from Elsevier via http://dx.doi.org/10.1016/j.cherd.2015.12.00

Optimal Operation Strategy for Biohydrogen Production

Hydrogen produced by microalgae is intensively researched as a potential alternative to conventional energy sources. Scaling-up of the process is still an open issue, and to this end, accurate dynamic modeling is very important. A challenge in the development of these highly nonlinear dynamic models is the estimation of the associated kinetic parameters. This work presents the estimation of the parameters of a revised Droop model for biohydrogen production by Cyanothece sp. ATCC 51142 in batch and fed-batch reactors. The latter reactor type results in an optimal control problem in which the influent concentration of nitrate is optimized which has never been considered previously. The kinetic model developed is demonstrated to predict experimental data to a high degree of accuracy. A key contribution of this work is the prediction that hydrogen productivity can achieve 3365 mL/L through an optimally controlled fed-batch process, corresponding to an increase of 116% over other recently published strategies.Author E. A. del Rio-Chanona would like to acknowledge CONACyT scholarship No. 522530 and the Secretariat of Public Education and the Mexican government for funding this project. Author P. Dechatiwongse is supported by a scholarship from the Royal Thai Government, Thailand. Solar Hydrogen Project was funded by the UK Engineering and Physical Sciences Research Council (EPSRC), project reference EP/F00270X/1.This is the author accepted manuscript. The final version is available from ACS via http://dx.doi.org/10.1021/acs.iecr.5b0061

Echelle long-slit optical spectroscopy of evolved stars

We present echelle long-slit optical spectra of a sample of objects evolving off the AGB, most of them in the pre-planetary nebula (pPN) phase, obtained with the ESI and MIKE spectrographs at Keck-II and Magellan-I, respectively. The total wavelength range covered with ESI (MIKE) is ~3900 to 10900 A (~3600 to 7200A). In this paper, we focus our analysis mainly on the Halpha profiles. Prominent Halpha emission is detected in half of the objects, most of which show broad Halpha wings (up to ~4000 km/s). In the majority of the Halpha-emission sources, fast, post-AGB winds are revealed by P-Cygni profiles. In ~37% of the objects Halpha is observed in absorption. In almost all cases, the absorption profile is partially filled with emission, leading to complex, structured profiles that are interpreted as an indication of incipient post-AGB mass-loss. All sources in which Halpha is seen mainly in absorption have F-G type central stars, whereas sources with intense Halpha emission span a larger range of spectral types from O to G. Shocks may be an important excitation agent of the close stellar surroundings for objects with late type central stars. Sources with pure emission or P Cygni Halpha profiles have larger J-K color excess than objects with Halpha mainly in absorption, which suggests the presence of warm dust near the star in the former. The two classes of profile sources also segregate in the IRAS color-color diagram in a way that intense Halpha-emitters have dust grains with a larger range of temperatures. (abridged)Comment: 68 pages, 14 figures, accepted for publication in ApJS (abstract abridged

Nonlinear electrodynamics and CMB polarization

Recently WMAP and BOOMERanG experiments have set stringent constraints on the polarization angle of photons propagating in an expanding universe: $\Delta \alpha = (-2.4 \pm 1.9)^\circ$. The polarization of the Cosmic Microwave Background radiation (CMB) is reviewed in the context of nonlinear electrodynamics (NLED). We compute the polarization angle of photons propagating in a cosmological background with planar symmetry. For this purpose, we use the Pagels-Tomboulis (PT) Lagrangian density describing NLED, which has the form $L\sim (X/\Lambda^4)^{\delta - 1}\; X$, where $X=1/4 F_{\alpha\beta} F^{\alpha \beta}$, and $\delta$ the parameter featuring the non-Maxwellian character of the PT nonlinear description of the electromagnetic interaction. After looking at the polarization components in the plane orthogonal to the ($x$)-direction of propagation of the CMB photons, the polarization angle is defined in terms of the eccentricity of the universe, a geometrical property whose evolution on cosmic time (from the last scattering surface to the present) is constrained by the strength of magnetic fields over extragalactic distances.Comment: 17 pages, 2 figures, minor changes, references adde

V2324Cyg - an F-type star with fast wind

For the first time high-resolution optical spectroscopy of the variable star V2324Cyg associated with the IR-source IRAS20572+4919 is made. More than 200 absorption features (mostly FeII, TiII, CrII, YII, BaII, and YII) are identified within the wavelength interval 4549-7880AA. The spectral type and rotation velocity of the star are found to be F0III and Vsini=69km/s, respectively. HI and NaID lines have complex PCyg-type profiles with an emission component. Neither systematic trend of radial velocity Vr with line depth Ro nor temporal variability of Vr have been found. We determined the average heliocentric radial velocity Vr=-16.8\pm 0.6km/s. The radial velocities inferred from the cores of the absorption components of the H$\beta$ and NaI wind lines vary from -140 to -225km/s (and the expansion velocities of the corresponding layers, from about 120 to 210km/s). The maximum expansion velocity is found for the blue component of the split H$\alpha$ absorption: 450km/s for December 12, 1995. The model atmospheres method is used to determine the star's parameters: Teff=7500K, log g=2.0, $\xi_t$=6.0km/s, and metallicity, which is equal to the solar value. The main peculiarity of the chemical abundances pattern is the overabundance of lithium and sodium. The results cast some doubt on the classification of V2324Cyg as a post-AGB star.Comment: 28 pages, 6 figure

Red-giant and main-sequence solar-like oscillators in binary systems revealed by ESA Gaia Data Release 3 -- Reconstructing stellar and orbital evolution from binary-star ensemble seismology

Binary systems constitute a valuable astrophysics tool for testing our understanding of stellar structure and evolution. Systems containing a oscillating component are interesting as asteroseismology offers independent parameters for the oscillating component that aid the analysis. About 150 of such systems are known in the literature. To enlarge the sample of these benchmark objects, we crossmatch the Two-Body-Orbit Catalogue (TBO) of Gaia DR3, with catalogs of confirmed solar-like oscillators on the main-sequence and red-giant phase from NASA Kepler and TESS. We obtain 954 new binary system candidates hosting solar-like oscillators, of which 45 and 909 stars are on the main sequence and red-giant, resp., including 2 new red giants in eclipsing systems. 918 oscillators in potentially long-periodic systems are reported. We increase the sample size of known solar-like oscillators in binary systems by an order of magnitude. We present the seismic properties of the full sample and conclude that the grand majority of the orbital elements in the TBO is physically reasonable. 82% of all TBO binary candidates with multiple times with APOGEE are confirmed from radial-velocity measurement. However, we suggest that due to instrumental noise of the TESS satellite the seismically inferred masses and radii of stars with $\nu_\textrm{max}$$\lesssim$30$\mu$Hz could be significantly overestimated. For 146 giants the seismically inferred evolutionary state has been determined and shows clear differences in their distribution in the orbital parameters, which are accounted the accumulative effect of the equilibrium tide acting in these evolved binary systems. For other 146 systems hosting oscillating stars values for the orbital inclination were found in the TBO. From testing the TBO on the SB9 catalogue, we obtain a completeness factor of 1/3.Comment: under review for publication in A&A (22 pages + 4 pages of appendix, 21 figures, 33 pages of tables in the Appendix

Identification of Tetranectin as a Potential Biomarker for Metastatic Oral Cancer

Lymph node involvement is the most important predictor of survival rates in patients with oral squamous cell carcinoma (OSCC). A biomarker that can indicate lymph node metastasis would be valuable to classify patients with OSCC for optimal treatment. In this study, we have performed a serum proteomic analysis of OSCC using 2-D gel electrophoresis and liquid chromatography/tandem mass spectrometry. One of the down-regulated proteins in OSCC was identified as tetranectin, which is a protein encoded by the CLEC3B gene (C-type lectin domain family 3, member B). We further tested the protein level in serum and saliva from patients with lymph-node metastatic and primary OSCC. Tetranectin was found significantly under-expressed in both serum and saliva of metastatic OSCC compared to primary OSCC. Our results suggest that serum or saliva tetranectin may serve as a potential biomarker for metastatic OSCC. Other candidate serum biomarkers for OSCC included superoxide dismutase, ficolin 2, CD-5 antigen-like protein, RalA binding protein 1, plasma retinol-binding protein and transthyretin. Their clinical utility for OSCC detection remains to be further tested in cancer patients