Gray-box Soft Sensor for Water Content Monitoring in Fluidized Bed Granulation.

Soft sensors are powerful tools for the implementation of process analytical technology (PAT). They are categorized into white-box (first-principle), black-box (statistical), and gray-box models. Gray-box models integrate white-box and black-box models to address each drawback, i.e., prediction accuracy and intuitiveness. Although they have been applied to various industrial processes, their applicability to water content monitoring in fluidized bed granulation has not been reported. In this study, we evaluated three types of gray-box models, i.e., parallel, serial, and combined gray-box models, in terms of prediction accuracy using real operating data on a commercial scale with two formulations. The gray-box models were constructed by integrating the heat and mass balance model (white-box model) and locally weighted partial least squares regression (LW-PLSR) model (black-box model). LW-PLSR was utilized to cope with collinearity and nonlinearity. In the serial gray-box models, LW-PLSR models adjusted the fitting parameters of the white-box model depending on the process parameters for each query. In the parallel gray-box or combined gray-box models, LW-PLSR models compensated for the output error of the white-box or serial gray-box models, respectively. The results demonstrated that all three types of gray-box models improved the prediction accuracy of the white-box models regardless of the formulation. Besides, we proposed the assessment method based on Hotelling's T2 and Q residual for gray-box models using LW-PLSR, which contributes decision support to select gray-box or white-box model. The accurate and descriptive gray-box models are expected to enhance process understanding and precise quality control in fluidized bed granulation.

Screening of sleep apnea based on heart rate variability and long short-term memory.

PURPOSE: Sleep apnea syndrome (SAS) is a prevalent sleep disorder in which apnea and hypopnea occur frequently during sleep and result in increase of the risk of lifestyle-related disease development as well as daytime sleepiness. Although SAS is a common sleep disorder, most patients remain undiagnosed because the gold standard test polysomnography (PSG), is high-cost and unavailable in many hospitals. Thus, an SAS screening system that can be used easily at home is needed. METHODS: Apnea during sleep affects changes in the autonomic nervous function, which causes fluctuation of the heart rate. In this study, we propose a new SAS screening method that combines heart rate measurement and long short-term memory (LSTM) which is a type of recurrent neural network (RNN). We analyzed the data of intervals between adjacent R waves (R-R interval; RRI) on the electrocardiogram (ECG) records, and used an LSTM model whose inputs are the RRI data is trained to discriminate the respiratory condition during sleep. RESULTS: The application of the proposed method to clinical data showed that it distinguished between patients with moderate-to-severe SAS with a sensitivity of 100% and specificity of 100%, results which are superior to any other existing SAS screening methods. CONCLUSION: Since the RRI data can be easily measured by means of wearable heart rate sensors, our method may prove to be useful as an SAS screening system at home.

Prediction Performance and Economic Efficiency of Soft Sensors for in-Line Water Content Monitoring in Fluidized Bed Granulation: PP-Based Model vs. NIRS-Based Model.

Soft sensors play a crucial role as process analytical technology (PAT) tools. They are classified into physical models, statistical models, and their hybrid models. In general, statistical models are better estimators than physical models. In this study, two types of standard statistical models using process parameters (PPs) and near-infrared spectroscopy (NIRS) were investigated in terms of prediction accuracy and development cost. Locally weighted partial least squares regression (LW-PLSR), a type of nonlinear regression method, was utilized. Development cost was defined as the cost of goods required to construct an accurate model of commercial-scale equipment. Eleven granulation lots consisting of three laboratory-scale, two pilot-scale, and six commercial-scale lots were prepared. Three commercial-scale granulation lots were selected as a validation dataset, and the remaining eight granulation lots were utilized as calibration datasets. The results demonstrated that the PP-based and NIRS-based LW-PLSR models achieved high prediction accuracy without using the commercial-scale data in the calibration dataset. This practical case study clarified that the construction of accurate LW-PLSR models requires the calibration samples with the following two features: 1) located near the validation samples on the subspace spanned by principal components (PCs), and 2) having a wide range of variations in PC scores. In addition, it was confirmed that the reduction in cost and mass fraction of active pharmaceutical ingredient (API) made the PP-based models more cost-effective than the NIRS-based models. The present work supports to build accurate models efficiently and save the development cost of PAT.

Approach to Establishment of Control Strategy for Oral Solid Dosage Forms Using Continuous Manufacturing.

As a result of the research activities of the Japan Agency for Medical Research and Development (AMED), this document aims to show an approach to establishing control strategy for continuous manufacturing of oral solid dosage forms. The methods of drug development, technology transfer, process control, and quality control used in the current commercial batch manufacturing would be effective also in continuous manufacturing, while there are differences in the process development using continuous manufacturing and batch manufacturing. This document introduces an example of the way of thinking for establishing a control strategy for continuous manufacturing processes.

Autoencoder-Based Extrasystole Detection and Modification of RRI Data for Precise Heart Rate Variability Analysis.

Heart rate variability, which is the fluctuation of the R-R interval (RRI) in electrocardiograms (ECG), has been widely adopted for autonomous evaluation. Since the HRV features that are extracted from RRI data easily fluctuate when arrhythmia occurs, RRI data with arrhythmia need to be modified appropriately before HRV analysis. In this study, we consider two types of extrasystoles-premature ventricular contraction (PVC) and premature atrial contraction (PAC)-which are types of extrasystoles that occur every day, even in healthy persons who have no cardiovascular diseases. A unified framework for ectopic RRI detection and a modification algorithm that utilizes an autoencoder (AE) type of neural network is proposed. The proposed framework consists of extrasystole occurrence detection from the RRI data and modification, whose targets are PVC and PAC. The RRI data are monitored by means of the AE in real time in the detection phase, and a denoising autoencoder (DAE) modifies the ectopic RRI caused by the detected extrasystole. These are referred to as AE-based extrasystole detection (AED) and DAE-based extrasystole modification (DAEM), respectively. The proposed framework was applied to real RRI data with PVC and PAC. The result showed that AED achieved a sensitivity of 93% and a false positive rate of 0.08 times per hour. The root mean squared error of the modified RRI decreased to 31% in PVC and 73% in PAC from the original RRI data by DAEM. In addition, the proposed framework was validated through application to a clinical epileptic seizure problem, which showed that it correctly suppressed the false positives caused by PVC. Thus, the proposed framework can contribute to realizing accurate HRV-based health monitoring and medical sensing systems.

Scale-Free Soft Sensor for Monitoring of Water Content in Fluid Bed Granulation Process.

In-line monitoring of granule water content during fluid bed granulation is important to control drug product qualities. In this study, a practical scale-free soft sensor to predict water content was proposed to cope with the manufacturing scale changes in drug product development. The proposed method exploits two key ideas to construct a scale-free soft sensor. First, to accommodate the changes in the manufacturing scale, the process parameters (PPs) that are critical to water content at different manufacturing scales were selected as input variables. Second, to construct an accurate statistical model, locally weighted partial least squares regression (LW-PLSR), which can cope with collinearity and nonlinearity, was utilized. The soft sensor was developed using both laboratory (approx. 4 kg) data and pilot (approx. 25 kg) scale data, and the prediction accuracy in the commercial (approx. 100 kg) scale was evaluated based on the assumption that the process was scaled-up from the pilot scale to the commercial scale. The developed soft sensor exhibited a high prediction accuracy, which was equivalent to the commonly used near-infrared (NIR) spectra-based method. The proposed method requires only standard instruments; therefore, it is expected to be a cost-effective alternative to the NIR spectra-based method.

Wearable Epileptic Seizure Prediction System with Machine-Learning-Based Anomaly Detection of Heart Rate Variability.

A warning prior to seizure onset can help improve the quality of life for epilepsy patients. The feasibility of a wearable system for predicting epileptic seizures using anomaly detection based on machine learning is evaluated. An original telemeter is developed for continuous measurement of R-R intervals derived from an electrocardiogram. A bespoke smartphone app calculates the indices of heart rate variability in real time from the R-R intervals, and the indices are monitored using multivariate statistical process control by the smartphone app. The proposed system was evaluated on seven epilepsy patients. The accuracy and reliability of the R-R interval measurement, which was examined in comparison with the reference electrocardiogram, showed sufficient performance for heart rate variability analysis. The results obtained using the proposed system were compared with those obtained using the existing video and electroencephalogram assessments; it was noted that the proposed method has a sensitivity of 85.7% in detecting heart rate variability change prior to seizures. The false positive rate of 0.62 times/h was not significantly different from the healthy controls. The prediction performance and practical advantages of portability and real-time operation are demonstrated in this study.

Data-Based Prediction and Stochastic Analysis of Entrained Flow Coal Gasification under Uncertainty.

Entrained flow gasification is a commonly used method for conversion of coal into syngas. A stable and efficient operation of entrained flow coal gasification is always desired to reduce consumption of raw materials and utilities, and achieve higher productivity. However, uncertainty in the process hinders the stability and efficiency. In this work, a quantitative analysis of the effect of uncertainty on the conversion efficiency of the entrained flow gasification is performed. A data-driven, i.e., ensemble, model of the process was developed to predict conversion efficiency of the process. Then sensitivity analysis methods, i.e., Sobol and Fourier amplitude sensitivity test, were used to analyze the effect of each individual process variables on conversion efficiency. For analyzing the collective impact of uncertainty in process variables on conversion efficiency, a non-intrusive polynomial chaos expansion (PCE) method was used. The PCE predicts probability distribution of the conversion efficiency. Reliability of the process was determined on the basis of percentage of the probability distribution falling within control limits. Measured data is used to derive the control limits for off-line reliability analysis. For on-line reliability analysis of the process, measured data is not available so a just-in-time method, i.e., k⁻d tree, was used. The k⁻d tree searches the nearest neighbor sample from a database of historical data to determine the control limits.

Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment.

Steam methane reforming (SMR) is a dominant technology for hydrogen production. For the highly energy-efficient operation, robust energy analysis is crucial. In particular, exergy analysis has received the attention of researchers due to its advantage over the conventional energy analysis. In this work, an exergy analysis based on the computational fluid dynamics (CFD)-based method was applied to a monolith microreactor of SMR. Initially, a CFD model of SMR was developed using literature data. Then, the design and operating conditions of the microreactor were optimized based on the developed CFD model to achieve higher conversion efficiency and shorter length. Exergy analysis of the optimized microreactor was performed using the custom field function (CFF) integrated with the CFD environment. The optimized catalytic monolith microreactor of SMR achieved higher conversion efficiency at a smaller consumption of energy, catalyst, and material of construction than the reactor reported in the literature. The exergy analysis algorithm helped in evaluating length-wise profiles of all three types of exergy, namely, physical exergy, chemical exergy, and mixing exergy, in the microreactor.

Missing RRI Interpolation Algorithm based on Locally Weighted Partial Least Squares for Precise Heart Rate Variability Analysis.

The R-R interval (RRI) fluctuation in electrocardiogram (ECG) is called heart rate variability (HRV), which reflects activities of the autonomic nervous system (ANS) and has been used for various health monitoring services. Accurate R wave detection is crucial for success in HRV-based health monitoring services; however, ECG artifacts often cause missing R waves and deteriorate the accuracy of HRV analysis. The present work proposes a new missing RRI interpolation technique based on Just-In-Time (JIT) modeling. In the JIT modeling framework, a local regression model is built by weighing samples stored in the database according to the distance from a query and output is estimated only when an estimate is requested. The proposed method builds a local model and estimates missing RRI only when an RRI detection error is detected. Locally weighted partial least squares (LWPLS) is adopted for local model construction. The proposed method is referred to as LWPLS-based RRI interpolation (LWPLS-RI). The performance of the proposed LWPLS-RI was evaluated through its application to RRI data with artificial missing RRIs. We used the MIT-BIH Normal Sinus Rhythm Database for nominal RRI dataset construction. Missing RRIs were artificially introduced and they were interpolated by the proposed LWPLS-RI. In addition, MEAN that replaces the missing RRI by a mean of the past RRI data was compared as a conventional method. The result showed that the proposed LWPLS-RI improved root mean squared error (RMSE) of RRI by about 70% in comparison with MEAN. In addition, the proposed method realized precise HRV analysis. The proposed method will contribute to the realization of precise HRV-based health monitoring services.

Visualisation of the T cell differentiation programme by Canonical Correspondence Analysis of transcriptomes.

BACKGROUND: Currently, in the era of post-genomics, immunology is facing a challenging problem to translate mutant phenotypes into gene functions based on high-throughput data, while taking into account the classifications and functions of immune cells, which requires new methods. RESULTS: Here we propose a novel application of a multidimensional analysis, Canonical Correspondence Analysis (CCA), to reveal the molecular characteristics of undefined cells in terms of cellular differentiation programmes by analysing two transcriptomic datasets. Using two independent datasets, whether RNA-seq or microarray data, CCA successfully visualised the cross-level relationships between genes, cells, and differentiation programmes, and thereby identified the immunological features of mutant cells (Gata3-KO T cells and Stat3-KO T cells) in a data-oriented manner. With a new concept, differentiation variable, CCA provides an automatic classification of cell samples, which had a high sensitivity and a comparable performance to other classification methods. In addition, we elaborate how CCA results can be interpreted, and reveal the features of CCA in comparison with other visualisation techniques. CONCLUSIONS: CCA is a visualisation tool with a classification ability to reveal the cross-level relationships of genes, cells and differentiation programmes. This can be used for characterising the functional defect of cells of interest (e.g. mutant cells) in the context of cellular differentiation. The proposed approach fits with common hypothesis-oriented studies in immunology, and can be used for a wide range of molecular and genomic studies on cellular differentiation mechanisms.

Feed factor profile prediction model for two-component mixed powder in the twin-screw feeder.

In continuous pharmaceutical manufacturing processes, it is crucial to control the powder flow rate. The feeding process is characterized by the amount of powder delivered per screw rotation, referred to as the feed factor. This study aims to develop models for predicting the feed factor profiles (FFPs) of two-component mixed powders with various formulations, while most previous studies have focused on single-component powders. It further aims to identify the suitable model type and to determine the significance of material properties in enhancing prediction accuracy by using several FFP prediction models with different input variables. Four datasets from the experiment were generated with different ranges of the mass fraction of active pharmaceutical ingredients (API) and the powder weight in the hopper. The candidates for the model inputs are (a) the mass fraction of API, (b) process parameters, and (c) material properties. It is desirable to construct a high-performance prediction model without the material properties because their measurement is laborious. The results show that using (c) as input variables did not improve the prediction accuracy as much, thus there is no need to use them.

Novel predictive approaches for drug-induced convulsions in non-human primates using machine learning and heart rate variability analysis.

Drug-induced convulsions are a major challenge to drug development because of the lack of reliable biomarkers. Using machine learning, our previous research indicated the potential use of an index derived from heart rate variability (HRV) analysis in non-human primates as a biomarker for convulsions induced by GABAA receptor antagonists. The present study aimed to explore the application of this methodology to other convulsants and evaluate its specificity by testing non-convulsants that affect the autonomic nervous system. Telemetry-implanted males were administered various convulsants (4-aminopyridine, bupropion, kainic acid, and ranolazine) at different doses. Electrocardiogram data gathered during the pre-dose period were employed as training data, and the convulsive potential was evaluated using HRV and multivariate statistical process control. Our findings show that the Q-statistic-derived convulsive index for 4-aminopyridine increased at doses lower than that of the convulsive dose. Increases were also observed for kainic acid and ranolazine at convulsive doses, whereas bupropion did not change the index up to the highest dose (1/3 of the convulsive dose). When the same analysis was applied to non-convulsants (atropine, atenolol, and clonidine), an increase in the index was noted. Thus, the index elevation appeared to correlate with or even predict alterations in autonomic nerve activity indices, implying that this method might be regarded as a sensitive index to fluctuations within the autonomic nervous system. Despite potential false positives, this methodology offers valuable insights into predicting drug-induced convulsions when the pharmacological profile is used to carefully choose a compound.

Simple algorithm for judging equivalence of differential-algebraic equation systems.

Mathematical formulas play a prominent role in science, technology, engineering, and mathematics (STEM) documents; understanding STEM documents usually requires knowing the difference between equation groups containing multiple equations. When two equation groups can be transformed into the same form, we call the equation groups equivalent. Existing tools cannot judge the equivalence of two equation groups; thus, we develop an algorithm to judge such an equivalence using a computer algebra system. The proposed algorithm first eliminates variables appearing only in either equation group. It then checks the equivalence of the equations one by one: the equations with identical algebraic solutions for the same variable are judged equivalent. If each equation in one equation group is equivalent to an equation in the other, the equation groups are judged equivalent; otherwise, non-equivalent. We generated 50 pairs of equation groups for evaluation. The proposed method accurately judged the equivalence of all pairs. This method is expected to facilitate comprehension of a large amount of mathematical information in STEM documents. Furthermore, this is a necessary step for machines to understand equations, including process models.

Greedy design space construction based on regression and latent space extraction for pharmaceutical development.

Implementation of the design space (DS) is a scientific concept for ensuring quality to be submitted as a part of the regulatory filing of a drug product for approval to market. An empirical approach is constructing the DS based on the regression model whose inputs are process parameters and material attributes over the different unit operations, i.e., a high-dimensional statistical model. While the high-dimensional model assures quality and process flexibility through a comprehensive process understanding, it has difficulty visualizing the feasible range of input parameters, i.e., DS. Therefore, this study proposes a greedy approach to constructing the extensive and flexible low-dimensional DS based on the high-dimensional statistical model and the observed internal representations that satisfies both comprehensive process understanding and the DS visualization capability. Introducing the observed correlation structure enabled the dimensionality reduction of the DS. The non-critical controllable parameters were fixed to the target values in visualizing the low-dimensional DS as a function of critical parameters. The expected variation of non-critical non-controllable parameters was considered the source of variation in prediction. The case study demonstrated the proposed approach's usefulness for developing the pharmaceutical manufacturing process.

Evaluation of infrared-reflection absorption spectroscopy measurement and locally weighted partial least-squares for rapid analysis of residual drug substances in cleaning processes.

The usefulness of infrared-reflection absorption spectroscopy (IR-RAS) for the rapid measurement of residual drug substances without sampling was evaluated. In order to realize the highly accurate rapid measurement, locally weighted partial least-squares (LW-PLS) with a new weighting technique was developed. LW-PLS is an adaptive method that builds a calibration model on demand by using a database whenever prediction is required. By adding more weight to samples closer to a query, LW-PLS can achieve higher prediction accuracy than PLS. In this study, a new weighting technique is proposed to further improve the prediction accuracy of LW-PLS. The root-mean-square error of prediction (RMSEP) of the IR-RAS spectra analyzed by LW-PLS with the new weighting technique was compared with that analyzed by PLS and locally weighted regression (LWR). The RMSEP of LW-PLS with the proposed weighting technique was about 36% and 14% smaller than that of PLS and LWR, respectively, when ibuprofen was a residual drug substance. Similarly, LW-PLS with the weighting technique was about 39% and 24% better than PLS and LWR in RMSEP, respectively, when magnesium stearate was a residual excipient. The combination of IR-RAS and LW-PLS with the proposed weighting technique is a very useful rapid measurement technique of the residual drug substances.

R-R interval-based sleep apnea screening by a recurrent neural network in a large clinical polysomnography dataset.

OBJECTIVE: Easily detecting patients with undiagnosed sleep apnea syndrome (SAS) requires a home-use SAS screening system. In this study, we validate a previously developed SAS screening methodology using a large clinical polysomnography (PSG) dataset (N = 938). METHODS: We combined R-R interval (RRI) and long short-term memory (LSTM), a type of recurrent neural networks, and created a model to discriminate respiratory conditions using the training dataset (N = 468). Its performance was validated using the validation dataset (N = 470). RESULTS: Our method screened patients with severe SAS (apnea hypopnea index; AHI ≥ 30) with an area under the curve (AUC) of 0.92, a sensitivity of 0.80, and a specificity of 0.84. In addition, the model screened patients with moderate/severe SAS (AHI ≥ 15) with an AUC of 0.89, a sensitivity of 0.75, and a specificity of 0.87. CONCLUSIONS: Our method achieved high screening performance when applied to a large clinical dataset. SIGNIFICANCE: Our method can help realize an easy-to-use SAS screening system because RRI data can be easily measured with a wearable heart rate sensor. It has been validated on a large dataset including subjects with various backgrounds and is expected to perform well in real-world clinical practice.

Medical checkup data analysis method based on LiNGAM and its application to nonalcoholic fatty liver disease.

Although medical checkup data would be useful for identifying unknown factors of disease progression, a causal relationship between checkup items should be taken into account for precise analysis. Missing values in medical checkup data must be appropriately imputed because checkup items vary from person to person, and items that have not been tested include missing values. In addition, the patients with target diseases or disorders are small in comparison with the total number of persons recorded in the data, which means medical checkup data is an imbalanced data analysis. We propose a new method for analyzing the causal relationship in medical checkup data to discover disease progression factors based on a linear non-Gaussian acyclic model (LiNGAM), a machine learning technique for causal inference. In the proposed method, specific regression coefficients calculated through LiNGAM were compared to estimate the causal strength of the checkup items on disease progression, which is referred to as LiNGAM-beta. We also propose an analysis framework consisting of LiNGAM-beta, collaborative filtering (CF), and a sampling approach for causal inference of medical checkup data. CF and the sampling approach are useful for missing value imputation and balancing of the data distribution. We applied the proposed analysis framework to medical checkup data for identifying factors of Nonalcoholic fatty liver disease (NAFLD) development. The checkup items related to metabolic syndrome and age showed high causal effects on NAFLD severity. The level of blood urea nitrogen (BUN) would have a negative effect on NAFLD severity. Snoring frequency, which is associated with obstructive sleep apnea, affected NAFLD severity, particularly in the male group. Sleep duration also affected NAFLD severity in persons over fifty years old. These analysis results are consistent with previous reports about the causes of NAFLD; for example, NAFLD and metabolic syndrome are mutual and bi-directionally related, and BUN has a negative effect on NAFLD progression. Thus, our analysis result is plausible. The proposed analysis framework including LiNGAM-beta can be applied to various medical checkup data and will contribute to discovering unknown disease factors.

Work Habit-Related Sleep Debt; Insights From Factor Identification Analysis of Actigraphy Data.

The present study investigates the factors of "Weekday sleep debt (WSD)" by comparing activity data collected from persons with and without WSD. Since it has been reported that the amount of sleep debt as well the difference between the social clock and the biological clock is associated with WSD, specifying the factors of WSD other than chronotype may contribute to sleep debt prevention. We recruited 324 healthy male employees working at the same company and collected their 1-week wrist actigraphy data and answers to questionnaires. Because 106 participants were excluded due to measurement failure of the actigraphy data, the remaining 218 participants were included in the analysis. All participants were classified into WSD or non-WSD groups, in which persons had WDS if the difference between their weekend sleep duration and the mean weekday sleep duration was more than 120 min. We evaluated multiple measurements derived from the collected actigraphy data and trained a classifier that predicts the presence of WSD using these measurements. A support vector machine (SVM) was adopted as the classifier. In addition, to evaluate the contribution of each indicator to WSD, permutation feature importance was calculated based on the trained classifier. Our analysis results showed significant importance of the following three out of the tested 32 factors: (1) WSD was significantly related to persons with evening tendency. (2) Daily activity rhythms and sleep were less stable in the WSD group than in the non-WSD group. (3) A specific day of the week had the highest importance in our data, suggesting that work habit contributes to WSD. These findings indicate some WSD factors: evening chronotype, instability of the daily activity rhythm, and differences in work habits on the specific day of the week. Thus, it is necessary to evaluate the rhythms of diurnal activities as well as sleep conditions to identify the WSD factors. In particular, the diurnal activity rhythm influences WSD. It is suggested that proper management of activity rhythm may contribute to the prevention of sleep debt.

Model-Based Quality, Exergy, and Economic Analysis of Fluidized Bed Membrane Reactors.

In petroleum refineries, naphtha reforming units produce reformate streams and as a by-product, hydrogen (H2). Naphtha reforming units traditionally deployed are designed as packed bed reactors (PBR). However, they are restrained by a high-pressure drop, diffusion limitations in the catalyst, and radial and axial gradients of temperature and concentration. A new design using the fluidized bed reactor (FBR) surpasses the issues of the PBR, whereby the incorporation of the membrane can improve the yield of products by selectively removing hydrogen from the reaction side. In this work, a sequential modular simulation (SMS) approach is adopted to simulate the hydrodynamics of a fluidized bed membrane reactor (FBMR) for catalytic reforming of naphtha in Aspen Plus. The reformer reactor is divided into five sections of plug flow reactors and a continuous stirrer tank reactor with the membrane module to simulate the overall FBMR. Similarly, a fluidized bed reactor (FBR), without membrane permeation phenomenon, is also modelled in the Aspen Plus environment for a comparative study with FBMR. In FBMR, the continuous elimination of permeated hydrogen enhanced the production of aromatics compound in the reformate stream. Moreover, the exergy and economic analyses were carried out for both FBR and FBMR.