Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fractional operators”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cycle-to-cycle variability in spark-assisted compression ignition engines near optimal mean combustion phasing

Stoichiometric spark-assisted compression ignition (SACI) combustion with exhaust gas recirculation (EGR) dilution has demonstrated higher part-load thermal efficiencies compared to spark-ignited engines, while maintaining ultra-low tailpipe emissions. However, SACI is often characterized by high cyclic variability in heat release or torque output, which poses a challenge to its implementation in light-duty vehicles. This paper presents an experimental investigation of cyclic variability in stoichiometric SACI combustion under EGR-dilute conditions, while maintaining mean combustion phasing (θ 50 ) near optimal timing for thermal efficiency. The present work focuses on SACI conditions where the flame-based heat release fraction is between approximately 10% and 40% of the overall heat release, and therefore contributes significantly to the combustion process. For the SACI conditions examined, the variability in θ 50 was driven by variability in autoignition timing, which in turn correlated with the start of measurable heat release (θ 02 ). High variability in θ 50 caused unstable work output due to very late combustion with poor or no end-gas autoignition. The use of a high ignition energy dual-coil offset ignition system had negligible impact in reducing θ 50 variability. Analysis of experimental data from close to 1000 operating conditions showed that the magnitude of θ 50 variability correlates with the flame-based heat release fraction ([Formula: see text]), for a large range of intake pressures, spark timings and exhaust gas recirculation levels. For the SACI conditions examined, combustion phasing variability was largely determined by flame-based combustion and particularly the initial flame formation (θ IGN-02 ), and minimally by the end-gas autoignition heat release. The analysis also demonstrated that θ 50 variability is amplified as the contribution of the flame to the overall heat release increases.

Engineering↗
DOE OSTI2024

Direct internal recycling fractions approaching unity

Direct internal recycling (DIR) refers to the process of recovering pure hydrogen isotopes (D/T) from helium and other impurities in the fusion plasma exhaust and directing them back to the fuel injection system. Increasing the exhaust fraction purified through DIR significantly reduces the size and cost of the tritium plant and provides additional benefits including reduced requirements for both the tritium startup inventory and tritium breeding ratio. Metal foil pumps (MFPs) are the dominant technology for this separation, relying on the concept of superpermeation. We recently demonstrated that PdCu foils operated at low temperature provide both exceptional flux and resilience to helium absorption as the DIR fraction is increased. Herein we design and demonstrate continuous and semi-batch DIR processes using PdCu MFPs. Under continuous processing, stable performance was observed for DIR fractions up to 92 %. In addition, we demonstrate a semi-batch process capable of extending the DIR fraction to unity (> 99.8 %). Under the experimental conditions described within a PdCu MFP area of ~22 m 2 would be sufficient to process the fusion exhaust with 92 % DIR fraction at expected flowrates of 100 Pa·m 3 ·s -1 for a future fusion power plant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Operability thresholds for thermally damaged EBW detonators

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (~406 K). Furthermore, the temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

42 ENGINEERING↗
DOE OSTI2022

High Yield, Economical and Environmentally Benign Production of Rare Earth Elements from Coal Ash [Abstract]

Fly ash stored in landfills and ponds across the United States is an attractive, abundant domestic resource for the cost-effective recovery of rare earth elements (REE) and other critical minerals (CM). Physical Sciences Inc. (PSI) and its team members, Winner Water Services (WWS) and University of Kentucky/Center for Applied Energy Research (UK/CAER) successfully executed a multiphase program that developed and demonstrated pilot plant operations for cost-effective and environmentally-friendly production of rare earth element oxide (REO) concentrates, and the critical minerals scandium and aluminum (in the forms of salts or oxide products). We constructed and operated a sub-scale (0.5 kg/day) micropilot facility to demonstrate the key physical and chemical processing operations, predict yields, and troubleshoot process bottlenecks. The project team then designed, constructed and operated two decoupled pilot plant operations: (1) an operational pilot plant for physical separation processes with a capacity of 0.4 metric tons per day (tpd) where we optimized processes to produce selected ash fractions as the feedstock for chemical processing as well as valuable byproducts such as cement substitute, cenospheres, magnetic ash and secondary fuel carbon, and (2) an operational pilot plant for chemical processing with a capacity of 0.5 tpd that developed optimized processes for the production of: (a) REO concentrates, (b) critical minerals recovery (Sc, Al), and (c) beneficiated ash as valuable byproduct suitable for cement applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

An analysis of fluff formation in metallic fuel via data analyzes from EBR-II experiments and BISON fuel code modeling

During the operation of EBR-II, it was found that a highly porous structure (over 40% area fraction) formed at the top of several fuel columns. Previous work has shown that this structure, designated fluff in this paper, contains a significant fraction of fuel elements (e.g., U and Pu) which could potentially impact neutronics. This work aims in analyzing the formation mechanism of this microstructure so its impact can be incorporated into future metallic fuel modeling codes and algorithms. This paper details a preliminary examination into the formation mechanisms of fluff by performing qualitative and statistical analysis of EBR-II experimental data. Additionally, the operating conditions that have the greatest impact on fluff formation were determined based on this data set. Also, BISON fuel code simulations were used to help postulate potential fluff formation mechanisms. From this analysis it was found that the largest contributors to fluff formation were fuel burnup and composition, with fluff formation exhibiting a roughly linear positive correlation with increasing burnup and a negative correlation with increasing Pu content. It was also found that higher pin operating temperatures decreased fluff formation but only for U-10Zr fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Power handling in a highly-radiative negative triangularity pilot plant

Abstract This work explores detailed power handling solutions for a class of high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept (Rutherford et al 2024 Plasma Phys. Control. Fusion ). The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only ∼50 cm (∼40% of the minor radius). The UEDGE code package is used for the boundary plasma solution to track plasma temperatures and fluxes to the divertor targets. It is found that for P SOL = 25 MW and n sep = 0.96 × 10 20 m −3 , conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For a fixed impurity fraction of just 0.13% Ne in the plasma, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW m −2 , while the electron temperature T e remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, P SOL can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1% neon (Ne) impurity content, the divertor still experiences less than 10 MW m −2 at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550 ∘ C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01% Ne content.

Miller, M. A. (ORCID:0000000265406533)↗
DOE OSTI2025

Mesofluidic Inline Separation for Produced Water Treatment

Mesofluidic inline separation developed by PNNL represents an opportunity to remove a key barrier in the treatment of produced water, namely suspended solids that clog downstream operations. The US alone produces over 800 billion gallons of produced water each year, most of which is reinjected underground (but not into the aquifers) as waste. The impact from treating and reusing even a fraction of this wastewater is immense as aquifers in the Midwest and elsewhere dry. The work described herein is essential to leveraging the opportunity to improve produced water quality (to allow for beneficial use) by combining mesofluidic inline separators with reverse osmosis systems. Providing water for agricultural and industrial uses in the American Southwest and Midwest by allowing the reuse of petroleum produced water may be critical to the long-term economics of the region. This is especially true as drought conditions are rapidly lowering aquifer levels in these regions. The incumbent technology for desalination is reverse osmosis (RO) due to its ability to treat a wide range of feedwaters and technological maturity. However, suspended solids cause RO (and other dissolved solids removal technologies) to become clogged and loose performance. A common strategy is to use prefilters upstream of RO systems, but these membrane-based filters easily clog and require regular maintenance. Unlike membrane filters, mesofluidic inline separators provide removal of suspended solids with much lower pressure drops than conventional systems filtration systems, permitting substantially higher flowrates. Although the amount of produced water from petroleum operations is vast, only a small fraction of it is reused or turned into potable water, for example, because of the lack of technologies to remove both dissolved and suspended solids at a significant throughput. This project shows that mesofluidic inline separators coupled with a commercial dissolved solids removal technology (RO) are positioned to do exactly that. Mesofluidic inline separators fit within commercial piping and are tunable for particle sizes of interest as described below. Indeed, mesofluidic inline separators have remarkably smaller footprints than competing filtration technologies and are easily transportable from jobsite to jobsite. This technology has no moving parts so that solids removal can be accomplished at much lower operating cost.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Performance Comparison of Helium and sCO 2 as Coolants for Modular Finger-Type Divertors

Several studies at the Georgia Institute of Technology have evaluated the thermal and fluids performance of the helium-cooled modular divertor with multiple jets (HEMJ) over the past decade. This finger-type divertor was studied both experimentally at nearly prototypical conditions and numerically at fully prototypical operating conditions using experimentally validated simulations. Recently, supercritical carbon dioxide (sCO 2 ) has been studied as the primary coolant in power cycles and other applications in various systems, in part because CO 2 achieves the high densities typical of supercritical fluids at relatively low temperatures and pressures, with a critical point of (7.38 MPa, 31°C). This density makes it possible to realize very compact and efficient sCO 2 power cycles. The feasibility of sCO 2 as a coolant for plasma-facing components, specifically the divertor, was therefore evaluated as part of the Fusion Energy System Studies design study activities. This work compares the thermal-fluid performance of helium and sCO 2 in the HEMJ divertor geometry using numerical simulations at prototypical conditions: inlet temperatures T i = 600°C to 700°C, pressures p ≈ 10 MPa, and steady-state incident heat fluxes on the tile q″ < 17 MW/m 2 . The performance is quantified here as the maximum heat flux that can be accommodated by the plasma-facing tile, the pumping power fraction, defined as the ratio of the coolant pumping power to the incident thermal power, and the operating stress limits based on ASME pressure vessel criteria. As expected, helium requires lower mass flow rates and pumping power fractions within imposed maximum temperature limits for the HEMJ pressure boundary. However, it also appears that neither helium nor sCO 2 can remove 10 MW/m 2 of incident heat flux while meeting ASME pressure vessel criteria. Finally, the numerical modeling reveals that sCO 2 may remove slightly higher incident heat fluxes than helium due to the imposed stress limits due to the sCO 2 coolant resulting in smaller local temperature gradients, albeit at a considerably higher pumping power fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Reliability of double probe measurements in nanodusty plasmas

Abstract Nonthermal plasmas are attractive sources for nanoparticles synthesis, however, their plasma properties are notoriously difficult to assess due to the chemically reactive environment and high nanoparticle concentrations. Here, we are using a floating double probe to measure the plasma properties of a nanoparticle-forming argon:silane plasma. We demonstrate good stability of current–voltage characteristics over several minutes of operation. However, unexpectedly larger electron temperatures are measured with increasing the silane mole fraction. To test the validity of these results, we developed a zero-dimensional global model to investigate the effect of the presence of nanoparticles on the plasma properties. Using this model, we show that increasing particle concentration leads to an increasing electronegativity of the plasma, causing an increase of the reduced electric field. However, this causes only a moderate increase in mean electron energy, in contrast to the much larger increase measured by the double probe. We argue that these large electron temperatures are based on the fact that a double probe measures an ‘apparent’ electron temperature, which is defined by the negative inverse slope of the logarithm of the electron energy probability function (EEPF) at an energy corresponding to the probe’s floating potential. As the silane mole fraction is increased, the plasma becomes more electronegative and the probe’s floating potential moves closer to the plasma potential. Combined with the strong non-Maxwellian EEPF, this leads to the large apparent electron temperatures obtained by the probe. Thus, the apparent electron temperatures measured with the double probe do not follow the trends in mean electron energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Skill of Mesoscale Models in Forecasting Springtime Macrophysical Cloud Properties at the Savannah River Site in the Southeastern US

Predicting boundary layer clouds is important for the accurate modeling of pollutant dispersion. Higher resolution mesoscale models would be expected to produce better forecasts of cloud properties that affect dispersion. Using ceilometer observations, we assess the skill of two operational mesoscale models (RAMS and WRF) to forecast cloud base altitude and cloud fraction at the Savannah River Site in the southeastern US during the springtime. Verifications were performed at small spatial and temporal scales necessary for dispersion modeling. Both models were unreliable with a 50% (RAMS) and a 46% (WRF) rate of predicting clouds observed by the ceilometer which led to low cloud fraction predictions. Results indicated that WRF better predicted daytime cloud bases from convection that occurred frequently later in the period and RAMS better predicted nighttime cloud bases. Using root mean squared error (RMSE) to score the forecast periods also highlighted this diurnal dichotomy, with WRF scores better during the day and RAMS scores better at night. Analysis of forecast errors revealed divergent model cloud base biases—WRF low and RAMS high. A hybrid solution which weighs more heavily the RAMS nighttime forecasts and WRF daytime forecasts will likely provide the best prediction of cloud properties for dispersion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Organizing Large Data Sets for Efficient Analyses on HPC Systems

Upcoming exascale applications could introduce significant data management challenges due to their large sizes, dynamic work distribution, and involvement of accelerators such as graphical processing units, GPUs. In this work, we explore the performance of reading and writing operations involving one such scientific application on two different supercomputers. Our tests showed that the Adaptable Input and Output System, ADIOS, was able to achieve speeds over 1TB/s, a significant fraction of the peak I/O performance on Summit. We also demonstrated the querying functionality in ADIOS could effectively support common selective data analysis operations, such as conditional histograms. In tests, this query mechanism was able to reduce the execution time by a factor of five. More importantly, ADIOS data management framework allows us to achieve these performance improvements with only a minimal amount of coding effort.

Gu, Junmin↗
DOE OSTI2025

A Data-Driven Framework for Predicting the Sorting and Screening Performance of an Integrated Biomass Feedstock Preprocessing System

The characteristics of mechanically sorted and screened lignocellulosic biomass, such as the mass contents of corn stover anatomical fractions (leaves, husks, stalks, cobs, etc.), can be used to calculate the intermediate feedstock quality attributes “yield” and “purity” that indicate the conversion efficiency of biocrude. No prior study has investigated the correlations from the characteristics of raw biomass and preprocessing unit operation parameters to those intermediate feedstock quality attributes. This work presents a data-driven framework for assessing and predicting the intermediate feedstock quality attributes in an integrated biomass feedstock preprocessing system. Our study used corn stover as a typical type of herbaceous biomass because of its abundance in the U.S. It began with data acquisition of moisture content, particle size distribution, and anatomical fractions of the materials after each unit operation in the system. The objective of this preprocessing system is to minimize husks and leaves and maximizing cobs and stalks by mechanically separating the materials into three streams via disc screen and air separator. Prototype neural network models were then developed to evaluate the feasibility of predicting process outcomes based on measurable parameters. It is found that incorporating physical constraints into these prediction models significantly enhances the accuracy of the predicted yield and purity against the ground truth data. The experimental data and model predictions indicate that decreasing throughput increases purity, while higher throughput results in lower purity. Finally, an optimization problem was introduced to search optimal combinations of feed material properties and preprocessing unit operation parameters, as the intermediate feedstock quality attributes – yield and purity, appeared to be competing factors. The study also suggests the continual need to improve the data-driven framework’s predictability by incorporating more accurate physical models to describe the dynamics in the preprocessing units such as the air separator.

09 - BIOMASS FUELS↗
DOE OSTI2022

HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi

Abstract Background Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing—SIP—remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance—the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling. Results Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to 13 C-AMF hyphosphere DNA from a 13 CO 2 plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10–33 atom% 13 C), even though the soils’ overall enrichment was low (1.8 atom% 13 C). We assembled 212 13 C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived 13 C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera . Conclusions Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP—fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat—generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs’ phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Two-phase operation of a Terry steam turbine using air and water mixtures as working fluids

Terry steam turbines are employed in the safety systems of many nuclear Boiling Water Reactors to drive pumps and provide cooling water to the nuclear reactor core. While the turbine efficiency is low, the more important feature is high reliability under off-normal conditions. An important aspect of reliability is the ability to function with two-phase steam-water injection into the turbine, as most likely occurred in the Fukushima Dai-ichi nuclear accidents. This study investigates the characteristics of a Terry turbine during air-water injection with gas mass fractions ranging from 1 (dry gas) to 0.05 (wet gas), to better understand the Terry turbine’s true operational capabilities and provide justification for extended Terry turbine use for reactor safety. Other parameters investigated are the inlet pressure, the exhaust backpressure and the turbine’s rotational speed. The turbine performance is presented in terms of dynamometer loading and pump performance change as functions of the gas mass fraction.

42 ENGINEERING↗
DOE OSTI2022

Heavy-Duty Vehicle Activity Updates for MOVES Using NREL Fleet DNA and CE-CERT Data

The U.S. Environmental Protection Agency's (EPA's) Motor Vehicle Emission Simulator (MOVES) is a publicly available tool used by researchers and policymakers to help understand motor vehicle emission sources at a national, county, and project level. Estimates of heavy-duty activity in the most recent version of the model at the time this work was conducted, MOVES2014, was identified as an area in need of improvement. The start activity in MOVES2014 is based on a limited and dated data set. In addition, MOVES2014 relies on drive cycles that represent on-network activity but do not account for idling activity that occurs on off-network roads, such as at a distribution center, while the truck is queuing or during loading and unloading. As a result, MOVES2014 may currently underestimate the number of starts and idle and soak time for heavy-duty trucks in real-world operation. The National Renewable Energy Laboratory (NREL) has previously leveraged its expansive Fleet DNA database of heavy-duty vehicles to idle and start activity for six of the nine heavy-duty vehicle source types of classes in the MOVES model. The data available in Fleet DNA from 416 conventional, diesel-powered vehicles provided activity estimates from more than 120,000 hours of operation throughout 14,682 vehicle days between October 2006 and January 2016. NREL calculated start fraction, starts per day, soak fraction, and idle fraction by hour of the day for each vehicle type, state, and vocation, and provided results in .CSV files that can be translated to MOVES table inputs. The idle and start activity from this initial analysis of Fleet DNA data was used to develop default idle and start data for heavy-duty vehicles in MOVES3. Satisfied with the results from the Fleet DNA data used for MOVES3, the EPA asked NREL to extend this start/soak/idle analysis using additional data from a larger number of vehicles for a potential future update to the MOVES model. Such a data set was achieved from a project led by the University of California at Riverside, College of Engineering, Center for Environmental Research & Technology (CE-CERT) and funded by California Air Resources Board. Specifically, this data set consists of 90 heavy-duty vehicles operated mainly in California, which can be separated into five of the nine heavy-duty vehicle classes in the MOVES model. In addition, the heavy-duty activity database collected by CE-CERT provided activity estimates from more than 44,000 hours of operation throughout 4,724 vehicle days between November 2014 and September 2016. This report details the analysis of the heavy-duty activity database collected from the University of California at Riverside by providing graphical analysis and context for the start, soak, and idle distributions. The comparison of the related results from both the Fleet DNA and CE-CERT data sets are documented as well.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Lessons from commissioning and operations of the ICARUS cryogenic system at Fermilab

Results and lessons from three years of commissioning and operations of the cryogenic and purificationsystems of the ICARUS neutrino detector installed at the Fermilab are described. The detector uses liquid argon stored in two separate containers of about 275 m 3 each. The liquid argon must be kept in very stable thermal conditions, and the contamination of electronegative impurities must be consistently kept at small fractions of parts per billion. Since its fill in early 2020, the detector has been commissioned and then operated in neutrino beams at Fermilab for the search of sterile neutrinos and measurements of neutrino-argon cross sections. The results of operating the cryogenic and purification systems and lessons learned, including challenges in reaching the design electron lifetime in argon bath, operating and servicing liquid nitrogen and argon circulation pumps, argon condensers and purifiers, are described. Plans for future operations and list of recommendations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗
DOE OSTI2023

A New Proposal Generalized Predictive Control Algorithm With Polynomial Reference Tracking Applied for Sodium Fast Reactors

This paper proposes a generalized predictive control (GPC) with constraints and orthonormal Laguerre functions using the simplified model of the primary system (reactor core and intermediate heat exchanger (IHX)) of a prototypical sodium fast reactor (SFR). This paper develops a multiple-input multiple-output (MIMO) GPC with input constraints able to track polynomial references of any degree applied in coolant temperature difference across the core and fractional power. The manipulated variables of the GPC-SFR are the reactivity and the sodium flow rate of the primary and secondary pipes. Moreover, orthonormal Laguerre functions and step down condition number techniques were also applied to avoid the numerical ill-conditioning issue in quadratic programming of large systems. Thus, a GPC type-2 was designed to control fractional power, coolant temperature difference across the core and sodium tank temperature of the SFR primary system when temperature references change according to a linear ramp after reaching their steady-state operation, sustaining 100% power operation on the reactor. In order to analyze the load tracking capability of the GPC-SFR type-2, the load following from 100% fractional power (FP) to 60% FP at 0.8% FP/min rate is simulated. Constraints on the rate of coolant temperature difference across the core and reactivity were applied for the design safety. For comparison criteria, this paper compares the GPC-SFR type-2 with the GPC-SFR type-1, i.e., standard model predictive control (MPC), to verify the viability and superior performance of the proposal regarding: (a) ramp-tracking capability of temperature and load; (b) the rejections of a reactivity disturbance of -1 cent and a secondary sodium inlet temperature disturbance of +10°F; and (c) a simulation with uncertainty in reactor design. The simulations show that the GPC-SFR type-2 overcome the GPC-SFR type-1 robustness and performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗