Search NASA⌕ Search

SEARCH · Search NASA

Results for “riser flows”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗
DOE OSTI2025

Design Approach for Circulating Fluidized Bed Riser Reactors such as NETL ASURE: A Review for Clustering Flow

Circulating fluidized bed riser reactors have been a commercial reality for some 80 to 90 years. Through all this time and hundreds if not thousands of research papers on the topic, no clear understanding exists on what is required to design and build a process reactor. There are numerous reasons for this, but the most significant reason is likely that riser reactors are chaotic in nature and unless this chaotic behavior is accounted for, there will be wide discrepancies in the predicted performance and in attempts to scale these reactors without significant trial and error. This work puts forth a new methodology to design riser reactors based upon the desired operating range using the chaotic microstructure behavior. This is accomplished through the use of a nondimensional regime map that was developed using chaotic parameters of memory and order/disorder along with higher moments of skewness and kurtosis. Validation of the approach is provided with circulating fluidized bed combustors and applied to the design of the NETL ASURE facility.

42 ENGINEERING↗
NASA NTRS2013

Electrical Capacitance Volume Tomography for the Packed Bed Reactor ISS Flight Experiment

Fixed packed bed reactors are compact, require minimum power and maintenance to operate, and are highly reliable. These features make this technology a highly desirable unit operation for long duration life support systems in space. NASA is developing an ISS experiment to address this technology with particular focus on water reclamation and air revitalization. Earlier research and development efforts funded by NASA have resulted in two hydrodynamic models which require validation with appropriate instrumentation in an extended microgravity environment. To validate these models, the instantaneous distribution of the gas and liquid phases must be measured.Electrical Capacitance Volume Tomography (ECVT) is a non-invasive imaging technology recently developed for multi-phase flow applications. It is based on distributing flexible capacitance plates on the peripheral of a flow column and collecting real-time measurements of inter-electrode capacitances. Capacitance measurements here are directly related to dielectric constant distribution, a physical property that is also related to material distribution in the imaging domain. Reconstruction algorithms are employed to map volume images of dielectric distribution in the imaging domain, which is in turn related to phase distribution. ECVT is suitable for imaging interacting materials of different dielectric constants, typical in multi-phase flow systems. ECVT is being used extensively for measuring flow variables in various gas-liquid and gas-solid flow systems. Recent application of ECVT include flows in risers and exit regions of circulating fluidized beds, gas-liquid and gas-solid bubble columns, trickle beds, and slurry bubble columns. ECVT is also used to validate flow models and CFD simulations. The technology is uniquely qualified for imaging phase concentrations in packed bed reactors for the ISS flight experiments as it exhibits favorable features of compact size, low profile sensors, high imaging speed, and flexibility to fit around columns of various shapes and sizes. ECVT is also safer than other commonly used imaging modalities as it operates in the range of low frequencies (1 MHz) and does not radiate radioactive energy. In this effort, ECVT is being used to image flow parameters in a packed bed reactor for an ISS flight experiment.

Tomography↗
DOE OSTI2021

Saturation Carrying Capacity for Group A Particles in a Circulating Fluidized Bed

Empirical models continue to play a significant role in the design process of multiphase chemical reactors, particularly riser reactors in circulating fluidized bed (CFB) processes. It is imperative that accurate, industrial relevant correlations are developed to aid these design efforts. Using poor correlations could result in startup issues and significant redesign work. In this work, a new correlation is proposed to predict the saturation carrying capacity of Geldart Group A particles. This new correlation improves upon the currently available correlations for these materials and covers a broad range of Geldart Group A particles (particle diameters from 52 to 70 µm, and Archimedes numbers ranging from 5 to 20), superficial gas velocities (1 to 4 m/s), and riser diameters (0.066 to 0.3048 m). The new correlation has an Absolute Average Percent Deviation of only 17.6%, making it the most accurate correlation for Geldart Group A particles in the current literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

FY23 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY23. As a continuation of the modeling effort from FY22, this year’s work focuses on improving the RELAP5-3D model developed previously for two-phase flow simulations. The RELAP5-3D model is updated to more accurately capture the heat loss experienced by the facility. The updated model is compared against experimental data for benchmarking purposes of the RELAP5-3D input model. By correctly accounting for heat loss, the updated RELAP5-3D model can now predict the two-phase baseline case more accurately. The onset and the duration of instability are captured well by the model. Furthermore, analyses are performed to better understand the instability mechanism experienced by the flow where the expansion of the boiling boundary in the chimney is studied in details and the fundamental frequencies of the oscillations are obtained. The updated RELAP5-3D model is further compared against four fault conditions, namely the reduction of riser header inlet flow area, depletion of system inventory, blocked riser channels, and static boiling scenario. For each fault condition, minor modifications and tuning are performed to improve the predictions of the model. The purpose of the analyses is to investigate the capability of RELAP5-3D in predicting complex two-phase flows in possible accident scenarios in actual Reactor Cavity Cooling System (RCCS). Overall, the model is able to capture the behaviors and trends of these fault conditions relatively well. Some discrepancies remain between the experimental data and the predictions, many of which are likely due to the differences in the predicted and experimental vapor generation rate. Future work will focus on the continued development of the current RELAP5-3D input model of the NSTF to both improve the accuracy of the model’s predictive capability and continue supporting the experimental program needs. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1989

The role of particle collisions in pneumatic transport

A model of dilute gas-solid flow in vertical risers is developed in which the particle phase is treated as a granular material, the balance equations for rapid granular flow are modified to incorporate the drag force from the gas, and boundary conditions, based on collisional exchanges of momentum and energy at the wall, are employed. In this model, it is assumed that the particle fluctuations are determined by inter-particle collisions only and that the turbulence of the gas is unaffected by the presence of the particles. The model is developed in the context of, but not limited to, steady, fully developed flow. A numerical solution of the resulting governing equations provides concentration profiles generally observed in dilute pneumatic flow, velocity profiles in good agreement with the measurements of Tsuji, et al. (1984), and an explanation for the enhancement of turbulence that they observed.

Mastorakos, E.↗
DOE OSTI2020

Effects of Riser Diameter on Solids Holdup and Particle Velocity Profiles in Circulating Fluidized Bed Riser Systems

Tests were performed in a 0.1-m diameter small circulating fluidized bed (SCFB) and 0.3 m diameter cold flow circulating fluidized bed (CFCFB) riser systems located at the National Energy Technology Laboratory (NETL) to study the effects of riser diameter on the riser hydrodynamics. These tests were performed at solids circulation rates of G s = 20 and 75 kg/m 2 s and superficial gas velocities of U g = 5.8 and 6.5 m/s using high-density polyethylene (HDPE) pellets with a density of 0.863 g/cm 3 , particle size range of 600–1400 µm (with a Sauter mean diameter of 871 µm, placing them in the Geldart B classification). Comparisons of riser axial pressure and solids fraction profiles, radial particle velocity profiles, and radial profiles of higher statistical moments and select chaos analysis parameters were considered. The results showed that for a given U g and G s , the smaller diameter riser exhibited characteristics associated with more dilute solids flow than that observed in the larger diameter riser. Additionally, the larger diameter riser exhibited a downward flow of solids near the wall under all test conditions, whereas the smaller diameter riser data exhibited little or no indications of solids downflow near the wall. Overall, these findings suggest that, from an industrial standpoint, a direct scaleup of small-scale tests cannot readily be accomplished as the solids holdup and the solids velocity profiles in small units (those normally tested in the laboratory) are not similar to those of large units and the performance of large units can therefore not be predicted from small-scale tests.

42 ENGINEERING↗
DOE OSTI2021

Comparative study with two reduced test facilities for air-cooled RCCS scaling law

Korea Atomic Energy Research Institute (KAERI) and Argonne National Laboratory (ANL) conducted steady-state tests at constant heat flux conditions to investigate the scale effect on heat removal behavior of the air-cooled Reactor Cavity Cooling System (RCCS). Two differently scaled-down test facilities were used in this study, KAERI’s ¼ scale with 4.0-m cavity height, and ANL’s ½ scale with 6.8-m cavity height. Two scaling laws were proposed to simulate the radiation across the cavity and the buoyancy-driven duct flow in the riser, respectively. The test matrix for KAERI and ANL test facilities focused on the scaling effect of the heated riser length with constant heat flux at the PMR200 RCCS design. The test results showed that mixed convection in the riser duct is an important factor for accurately extrapolating the thermo-fluid behavior in the prototype from the test results in the scale-down facilities. The system analysis code, GAMMA+, with improved heat transfer models predicted fairly well the air-cooled RCCS test data from two facilities. GAMMA+ analysis results showed that the predicted radiation fractions on the heated plate in the scale-down test conditions were larger than those on the reactor vessel in the prototype. The scaling law for air-cooled RCCS was improved by considering the mixed convection in the riser duct and the same radiation fraction on the heated plate. The mixed convection effect was calculated by the height ratio. The same radiation fraction provided a conservative extrapolation from the scale-down test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Passive techniques for long-term reactor cooling

In a pressurized water reactor (PWR), emergency core cooling (ECC) responds to depressurization due to a vessel penetration break at the top of the pressure vessel by draining water from a body of water through an injection line into the pressure vessel. A barrier operates concurrently with the ECC to suppress flow of liquid water from the pressure vessel out the vessel penetration break. The barrier may comprise one or more of: (1) an injection line extension passing through the central riser to drain water into the central riser; (2) openings in a lower portion of a central riser to shunt some upward flow from the central riser into a lower portion of the downcomer annulus; and (3) a surge line providing fluid communication between a pressurizer volume at the top of the pressure vessel and the remainder of the pressure vessel which directs water outboard toward the downcomer annulus.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Passive techniques for long-term reactor cooling

In a pressurized water reactor (PWR), emergency core cooling (ECC) responds to depressurization due to a vessel penetration break at the top of the pressure vessel by draining water from a body of water through an injection line into the pressure vessel. A barrier operates concurrently with the ECC to suppress flow of liquid water from the pressure vessel out the vessel penetration break. The barrier may comprise one or more of: (1) an injection line extension passing through the central riser to drain water into the central riser; (2) openings in a lower portion of a central riser to shunt some upward flow from the central riser into a lower portion of the downcomer annulus; and (3) a surge line providing fluid communication between a pressurizer volume at the top of the pressure vessel and the remainder of the pressure vessel which directs water outboard toward the downcomer annulus.

Malloy, III, John D.↗
DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗
DOE OSTI2022

EXPERIMENTAL AND MODELLING ANALYSIS OF A LARGE-SCALE TWO-PHASE LOOP THERMOSYPHON

Liquid pumping requires large quantities of electrical energy, including about 7% of the energy of building heating, ventilation, and air conditioning (HVAC) systems. To replace pumped condenser-cooling tower water loops with a passive alternative system, we implemented a commercial-scale two-phase loop thermosyphon (TPLT). The unit consists of a 13 m riser integrated with a commercially available cooling tower and circulation water heater that simulates heat loads up to 25 kW. In addition to providing passive cooling capabilities, the cooling tower unit is also maintenance free, reliable, and can operate in both dry and wet modes. This study characterizes the performance (temperature difference between the evaporator and condenser) and the two-phase flow behavior of the loop under various refrigerant charges. Sight glasses installed throughout the loop are used to identify the operating flow regimes in the riser and downcomer. Over the range of operating conditions, we identified that there is an optimal refrigerant charge range for a specific heat load at which near-isothermal operation can be achieved. We further developed a model to predict the operating flow rate and gravitational height in the downcomer and compared it with the experimental data. The results show that the model agrees with the experimental data, in particular the threshold gravitational head height at which there will be subcooled liquid at the condenser exit, indicating that it can be used as a predictive tool for identifying the optimal loop charge for a given heat load.

Loop thermosyphon, passive two phase↗
DOE OSTI2024

Fluidization Regime Mapping Analysis for 200 μm Glass Beads in a Circulating Fluidized Bed Riser

Powerpoint slides for an oral presentation at the 2024 Clearwater Clean Energy Conference. It includes a brief discussion of a recently developed non-dimensional map for scaling circulating fluidized bed risers using velocity and solids circulation rate ratios. Data is presented for a set of experiments conducted in NETL’s 12-inch cold flow circulating fluidized bed (CFCFB) riser with 200 micron diameter glass beads, showing a strong dependance between solids volume fraction profiles and the solids circulation rate ratio.

Rowan, Steven↗
DOE OSTI2024

Fluidization Regime Mapping Analysis for 200 μm Glass Beads in a Circulating Fluidized Bed Riser

This paper was an executive summary for an oral presentation to be included in the conference proceedings of the 2024 Clearwater Clean Energy Conference. It includes a brief discussion of a recently developed non-dimensional map for scaling circulating fluidized bed risers using velocity and solids circulation rate ratios. Data is presented for a set of experiments conducted in NETL’s 12-inch cold flow circulating fluidized bed (CFCFB) riser with 200 micron diameter glass beads, showing a strong dependance between solids volume fraction profiles and the solids circulation rate ratio.

Rowan, Steven↗
DOE OSTI2020

Two-phase flow measurements and observations in a cooling panel of the reactor cavity cooling system

The analysis of the flow behavior within a cooling panel of a scaled water-cooled reactor cavity cooling system has been performed supported by high-resolution measurements of the void fraction at the risers’ outlet. The results have shown a stable, symmetric distribution of the flow through the risers during subcooled boiling conditions, characterized by small, slow bubbles with average void fraction lower than 0.3 for all risers. As the temperature increased to saturation, flow distribution appeared strongly asymmetrical and unstable, with peaks of void fraction as high as 0.9. Large, faster bubbles were observed during these peaks, immediately followed by flow of subcooled liquid water. Flow stagnation and inversion has also been observed. The high resolution void fraction results, complemented by the measurements of the coolant flow rate and temperatures, not only have further described the behavior of the system under hypothetical accident conditions, but represent a unique set of data to support the validation of system-level codes, and more advanced computational tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗