Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Predictive Control”

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 1,081 records · Page 60

A smart alarm for particle accelerator beamline operations

We present the initial results of a proof-of-concept 'smart alarm' for the Continuous Electron Beam Accelerator Facility injector beamline at Jefferson Lab. To minimize machine downtime and improve operational efficiency, an autonomous alarm system able to identify and diagnose unusual machine states is needed. Our approach leverages a trained neural network capable of alerting operators (a) when an anomalous condition exists in the beamline and (b) identifying the element setting that is the root cause. The tool is based on an inverse model that maps beamline readings (diagnostic readbacks) to settings (beamline attributes operators can modify). The model takes as input readings from the machine and computes machine settings which are compared to control setpoints. Instances where predictions differ from setpoints by a user-defined threshold are flagged as anomalous. Given data corresponding to 354 anomalous injector configurations, the model can narrow the root cause of an anomalous condition to three potential candidates with 94.6% accuracy. Furthermore, compared to the current method of identifying anomalous conditions which raises an alarm when machine parameters drift outside their normal tolerances, the data-driven model can identify 83% more anomalous conditions.

43 PARTICLE ACCELERATORS↗

Self-diffusion of Ti interstitial based point defects and complexes in TiC

To date, the mechanism of Ti atom self-diffusion is unproven. Prior theoretical work mostly focused on Ti vacancy based mediators, but these do not reproduce the experimental activation energy or entropy. In this work, in density functional theory calculations, Ti interstitials and related defect complexes are systematically considered as possible mediators of Ti self-diffusion. Among these defects, the defect complex of two C vacancies tightly bound to a Ti dumbbell is found to have the lowest formation energy. A sustainable migration of the complex, in a translational or rotational fashion, is enabled in the presence of another (free) carbon vacancy nearby the complex, and thus the rate of Ti self-diffusion by this mechanism is dependent on the concentration of carbon vacancies. The calculated activation energy of the complex agrees well with the experimental value in TiC 0.97 . Similar analyses of the Ti self-diffusion mechanisms mediated by Ti interstitials or dumbbells yield much higher activation energies, but the corresponding migration energies are evaluated to be less than 1 eV, which suggests they can be possible mediators of the radiation-enhanced Ti self-diffusion in TiC. To fully enable the comparison with experiments that are typically conducted at temperatures as high as 2500 K, we also consider the temperature dependent vibrational contribution to the activation energy of the defect complex. The vibrational contribution imposes an additive effect on the defect formation energy, while the migration energies are lowered due to the thermal expansion of the lattice. When combined, these factors give an excellent agreement with the experiments. Furthermore, this work gives strong support to the concept that Ti interstitial based defect complexes are likely diffusion mediators for Ti atom self-diffusion in TiC, further establishes a solid basis for large-scale modeling, and may eventually pave the way to accurately predicting defect-controlled diffusional processes.

36 MATERIALS SCIENCE↗

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗

Airborne measurements of isoprene and monoterpene emissions from southeastern U.S. forests

Isoprene and monoterpene emission rates are essential inputs for atmospheric chemistry models that simulate atmospheric oxidant and particle distributions. Process studies of the biochemical and physiological mechanisms controlling these emissions are advancing our understanding and the accuracy of model predictions but efforts to quantify regional emissions have been limited by a lack of constraints on regional distributions of ecosystem emission capacities. We used an airborne wavelet-based eddy covariance measurement technique to characterize isoprene and monoterpene fluxes with high spatial resolution during the 2013 SAS (Southeast Atmosphere Study) in the southeastern United States. The fluxes measured by direct eddy covariance were comparable to emissions independently estimated using an indirect inverse modeling approach. Isoprene emission factors based on the aircraft wavelet flux estimates for high isoprene chemotypes (e.g., oaks) were similar to the MEGAN2.1 biogenic emission model estimates for landscapes dominated by oaks. Aircraft flux measurement estimates for landscapes with fewer isoprene emitting trees (e.g., pine plantations), were about a factor of two lower than MEGAN2.1 model estimates. The tendency for high isoprene emitters in these landscapes to occur in the shaded understory, where light dependent isoprene emissions are diminished, may explain the lower than expected emissions. This result demonstrates the importance of accurately representing the vertical profile of isoprene emitting biomass in biogenic emission models. Airborne measurement-based emission factors for high monoterpene chemotypes agreed with MEGAN2.1 in landscapes dominated by pine (high monoterpene chemotype) trees but were more than a factor of three higher than model estimates for landscapes dominated by oak (relatively low monoterpene emitting) trees. This results suggests that unaccounted processes, such as floral emissions or light dependent monoterpene emissions, or vegetation other than high monoterpene emitting trees may be an important source of monoterpene emissions in those landscapes and should be identified and included in biogenic emission models.

54 ENVIRONMENTAL SCIENCES↗

Aerosol-Cloud Interactions From Aviation Soot Emissions

Current models estimate global aviation contributes approximately 5% to the total anthropogenic climate forcing, with aerosol-cloud interactions having the greatest effect. However, radiative forcing estimates from aviation aerosol-cloud interactions remain undetermined. There is an expected significant increase in aircraft emissions with aviation demand expected to rise by over 4% per year. Soot may play an important role in the ice nucleation of aircraft-induced cirrus formation due to a high emission rate, but the ice nucleating properties are poorly constrained. Understanding the microphysical processes leading to atmospheric ice crystal formation is crucial for the reliable parameterization of aerosol-cloud interactions in climate models due to their impact on precipitation and cloud radiative properties. Ice nucleation of aircraft-emitted soot is potentially affected by particle morphology with condensation of supercooled water occurring in pores followed by ice nucleation. However, soot has heterogeneous properties and undergoes atmospheric aging and oxidation that could change surface properties and contribute to complex ice nucleation processes. Further, this review synthesizes current knowledge of ice nucleation catalyzed by aviation in the cirrus regime and its effects on global radiative forcing. Further research is required to determine the ice nucleation and microphysical processes of cirrus cloud formation from aviation emissions in both controlled laboratory and field investigations to inform models for more accurate climate predictions and to provide efficient mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

Theory of highly efficient multiexciton generation in type-II nanorods

Multiexciton generation, by which more than a single electron–hole pair is generated on optical excitation, is a promising paradigm for pushing the efficiency of solar cells beyond the Shockley–Queisser limit of 31%. Utilizing this paradigm, however, requires the onset energy of multiexciton generation to be close to twice the band gap energy and the efficiency to increase rapidly above this onset. This challenge remains unattainable even using confined nanocrystals, nanorods or nanowires. Here, we show how both goals can be achieved in a nanorod heterostructure with type-II band offsets. Using pseudopotential atomistic calculation on a model type-II semiconductor heterostructure we predict the optimal conditions for controlling multiexciton generation efficiencies at twice the band gap energy. For a finite band offset, this requires a sharp interface along with a reduction of the exciton cooling and may enable a route for breaking the Shockley–Queisser limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of resonant magnetic perturbations on secondary structures in drift-wave turbulence

Recent experiments showed a decrease of long range correlations during the application of resonant magnetic perturbations (RMPs) [Y. Xu et al., Nucl. Fusion 51, 063020 (2011)]. This finding suggests that RMPs damp zonal flows. To elucidate the effect of the RMPs on zonal structures in drift wave turbulence, we construct a generalized Hasegawa-Wakatani model including RMP fields. The effect of the RMPs is to induce a linear coupling between the zonal electric field and the zonal density gradient, which drives the system to a state of electron radial force balance for large RMP amplitude. A predator-prey model coupling the primary drift wave dynamics to the zonal modes evolution is derived. This model has both turbulence drive and RMP amplitude as control parameters and predicts a novel type of transport bifurcation in the presence of RMPs. The novel regime has a power threshold which increases with RMP amplitude as {gamma}{sub c}{approx}[({delta}B{sub r}/B)]{sup 2}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

De novo Nanopore read quality improvement using deep learning

BACKGROUND: Long read sequencing technologies such as Oxford Nanopore can greatly decrease the complexity of de novo genome assembly and large structural variation identification. Currently Nanopore reads have high error rates, and the errors often cluster into low-quality segments within the reads. The limited sensitivity of existing read-based error correction methods can cause large-scale mis-assemblies in the assembled genomes, motivating further innovation in this area. RESULTS: Here we developed a Convolutional Neural Network (CNN) based method, called MiniScrub, for identification and subsequent "scrubbing" (removal) of low-quality Nanopore read segments to minimize their interference in downstream assembly process. MiniScrub first generates read-to-read overlaps via MiniMap2, then encodes the overlaps into images, and finally builds CNN models to predict low-quality segments. Applying MiniScrub to real world control datasets under several different parameters, we show that it robustly improves read quality, and improves read error correction in the metagenome setting. Compared to raw reads, de novo genome assembly with scrubbed reads produces many fewer mis-assemblies and large indel errors. CONCLUSIONS: MiniScrub is able to robustly improve read quality of Oxford Nanopore reads, especially in the metagenome setting, making it useful for downstream applications such as de novo assembly. We propose MiniScrub as a tool for preprocessing Nanopore reads for downstream analyses. MiniScrub is open-source software and is available at https://bitbucket.org/berkeleylab/jgi-miniscrub .

59 BASIC BIOLOGICAL SCIENCES↗

Influence of optics vibration on synchrotron beam coherence

Vibration of optics is one of the major limiting factors of the performance of state-of-the-art beamlines at low-emittance synchrotron facilities. In this work, we present a theoretical model with experimental data describing vibration-induced effects on X-ray beam coherence. Owing to the Van Cittert-Zernike theorem, the coherence function of an X-ray beam perturbed by optics vibrations can be adequately characterized by modeling the spatial profile of the beam effective source. The measurements were carried out by using grating interferometry and a refractive lens with controlled vibration as a test optics. The experimental results confirm the model prediction and reveal the dependency of the measured beam spatial coherence on the acquisition time. The proposed method can be used to identify the eigenfrequency of the optical system as well as to optimize beamline operation and experimental conditions for coherence-related techniques.

36 MATERIALS SCIENCE↗

Technical Evaluation of Accelerated Basin De-Inventory Material Addition to Sludge Batch 11 (Rev. 1)

The Accelerated Basin De-inventory (ABD) program involves discarding spent nuclear fuel that is currently stored in L-Basin to the Defense Waste Processing Facility (DWPF) for vitrification. The first ABD discards will occur during the preparation of Sludge Batch (SB) 11. Savannah River Mission Completion has requested that the Savannah River National Laboratory assess the technical gaps related to the increased gadolinium poisoning requirement and the impacts of performing the Low Temperature Aluminum Dissolution (LTAD) process in Tank 51 with H-Canyon discards present. The following summarizes the evaluation of the impacts of increasing the quantity of gadolinium (and related topics) from what was previously evaluated in the SRNL studies of gadolinium-poisoned ABD material solubility, the overall ABD flowsheet review, and increasing the fissile mass loading in glass: 1) Based on literature surveys, there is no indication that organic interactions with gadolinium will be significant at the high pH (typically >13) conditions of the Concentration, Storage, and Transfer Facilities. Any interactions of gadolinium with organics in DWPF are not expected to adversely impact DWPF or downstream facilities. Thus, there is little-to-no residual risk from organic interactions with gadolinium [Gap closed]; 2) Adding depleted uranium to ABD material, targeting 235 U enrichment of 4.90% within each transfer window, will mitigate potential impacts from an increase in soluble 235 U enrichment during sludge washing and LTAD. The plan to take advantage of previous transfers and allow 235 U enrichment of >5% during the final transfer window carries a risk that Tank 51 supernate will have a 235 U enrichment of >5%, which should be evaluated for acceptance; 3) Increasing the gadolinium mass ratio to 3.0:1 Gd: 235 U(eq SLU ) should lead to the same or higher partitioning of gadolinium into the solid phase within the DWPF Chemical Process Cell, resulting in both liquid and solid phases with expected partitioning of Gd consistent with the prior solubility study [Gap closed for SB11]; 4) There are no expected impacts on DWPF melt temperature and melter operations due to the minimal ~0.2 weight percent (wt%) increase in Gd concentration relative to previous sludge batches [Gap closed for SB11]; 5) As observed previously, Gd is expected to enter the off-gas system via physical entrainment, but at a slightly higher concentration than what was observed for SB9 melter off-gas pluggage deposits (0.07 wt%) [Gap closed for SB11] ; 6) There are no expected impacts on DWPF recycle or the Recycle Collection Tank glycolate destruction process. [Gap closed for SB11]; 7) Gd is projected to be a trace component in the SB11 glass (<0.5 wt%) and can be ignored for process control. Trace components do not significantly impact glass durability, thus the conclusions of the previous Product Consistency Test evaluation at a fissile mass loading of 2,500 g fissile/m3 glass still applies to SB11. The ~0.1 wt% increase in Gd2O3 concentration relative to the previous study will not impact the predictability of SB11 glass with the DWPF Product Composition Control System (PCCS) models for durability or the acceptability of glass according to the Waste Acceptance Product Specifications (WAPS) criterion for product consistency [Gap closed for SB11]; 8) No additional Toxicity Characteristic Leaching Procedure testing is necessary for SB11 and the hazardous waste specification of the SB11 DWPF waste form is unchanged after the addition of the ABD stream [Gap closed for SB11]. The following summarizes the evaluation of the impacts of adding two-thirds of the ABD material to Tank 51 prior to LTAD: 1) The addition of two-thirds of the ABD increases overall aluminum mass from 1.39×10 4 kg to 1.64×10 4 kg (15.5% ABD Al). The form of the insoluble portion of the Al resulting from ABD addition should be the more readily dissolved Al(OH) 3 and amorphous forms. The portion of the ABD aluminum that is processed by LTAD is expected to be completely soluble, thus requiring that less of the boehmite in the sludge be dissolved to reach the same Al target in the SB. [Gap closed for SB11]; The expected LTAD impact on other components, as related primarily to the components in ABD, are discussed. Gd is expected to remain insoluble during LTAD and not impact the solubility of other components. [Gap closed for SB11]; The addition of two-thirds of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 13,100 kg (63% ABD U) and the projected plutonium mass from 86.0 kg to 89.5 kg (3.9% ABD Pu). The addition of all of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 16,100 kg (70% ABD U) and the projected plutonium mass from 86.0 kg to 90.4 kg (5.3% ABD Pu). The 235 U enrichment will be ≤5%. The fissile uranium will be adequately poisoned by Gd and the fissile Pu will be adequately poisoned by Fe from the sludge. [Gap closed for SB11]; There is a low risk that ABD addition will impact the rheology or pumpability of the slurry. There is a low but higher risk of ABD addition prior to LTAD impacting the settling rate; Based on the evaluation of adding two-thirds of the ABD material and all of the ABD material prior to the LTAD process, there is no volume or mass limit that would need to be imposed on ABD additions prior to LTAD. [Gap closed for SB11]. Revision 1 of this report addresses a variation on the ABD additions and LTAD strategy where sodium hydroxide additions for LTAD may be performed intermittently or concurrently with an ABD addition window. The proposed change does not alter the conclusions of this evaluation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Feasibility of Adding Twitter Data to Aid Drought Depiction: Case Study in Colorado

The use of social media, such as Twitter, has changed the information landscape for citizens’ participation in crisis response and recovery activities. Given that drought progression is slow and also spatially extensive, an interesting set of questions arise, such as how the usage of Twitter by a large population may change during the development of a major drought alongside how the changing usage facilitates drought detection. For this reason, contemporary analysis of how social media data, in conjunction with meteorological records, was conducted towards improvement in the detection of drought and its progression. The research utilized machine learning techniques applied over satellite-derived drought conditions in Colorado. Three different machine learning techniques were examined: the generalized linear model, support vector machines and deep learning, each applied to test the integration of Twitter data with meteorological records as a predictor of drought development. It is found that the integration of data resources is viable given that the Twitter-based model outperformed the control run which did not include social media input. Eight of the ten models tested showed quantifiable improvements in the performance over the control run model, suggesting that the Twitter-based model was superior in predicting drought severity. Future work lies in expanding this method to depict drought in the western U.S.

54 ENVIRONMENTAL SCIENCES↗

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Benchmark models and experimental data for a U(20) polyethylene-moderated critical system

This work involves the analysis of recent experiments performed on the Aerojet General Nucleonics (AGN)-201M (AGN) polyethylene-moderated research reactor at the University of New Mexico (UNM). The experiments include 36 delayed critical (DC) configurations and 11 positive-period and rod-drop measurements (transient sequences). The Even Parity Neutron Transport (EVENT) radiation transport code was chosen to analyze these steady state and time-dependent experimental configurations. The UNM AGN specifications provided in a benchmark calculation report (2007) were used to initiate AGN EVENT model development and to test the EVENT AGN calculation methodology. The results of the EVENT DC experimental analyses compared well with the experimental data; the average AGN EVENT calculation bias in the k eff is –0.0048% for the Legrendre Flux Expansion Order of 11 (P 11 ) cases and +0.0119% for the P 13 cases. The EVENT transient analysis also compared well with the AGN experimental data with respect to predicting the reactor period and control rod worth values. This paper discusses the benchmark models used, the recent experimental configurations, and the EVENT experimental analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗

Benchmark Models and experimental data for a U(20)-Polyethylene-moderated critical system - 14650

This work involves the analysis of recent experiments performed on the Aerojet General Nucleonics (AGN)-201M (AGN) polyethylene-moderated research reactor at the University of New Mexico (UNM). The experiments include 36 delayed critical (DC) configurations and 11 positive-period and rod-drop measurements (transient sequences). The Even Parity Neutron Transport (EVENT) radiation transport code was chosen to analyze these steady state and time-dependent experimental configurations. The UNM AGN specifications provided in a benchmark calculation report (2007) were used to initiate AGN EVENT model development and to test the EVENT AGN calculation methodology. The results of the EVENT DC experimental analyses compared well with the experimental data; the average AGN EVENT calculation bias in the k{sub eff} is -0.0048% for the Legendre Flux Expansion Order of 11 (P{sub 11}) cases and +0.0119% for the P{sub 13} cases. The EVENT transient analysis also compared well with the AGN experimental data with respect to predicting the reactor period and control rod worth values. This paper discusses the benchmark models used, the recent experimental configurations, and the EVENT experimental analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Prediction of Combustion Stability and Flashback in Turbines with High-Hydrogen Fuel

During the duration of this sponsorship, we broadened our understanding of combustion instabilities through both analytical and experimental work. Predictive models were developed for flame response to transverse acoustic instabilities and for quantifying how a turbulent flame responds to velocity and fuel/air ratio forcing. Analysis was performed on the key instability mechanisms controlling heat release response for flames over a wide range of instability frequencies. Importantly, work was done closely with industrial partners to transition existing models into internal instability prediction codes. Experimentally, the forced response of hydrogen-enriched natural gas/air premixed and partially premixed flames were measured. The response of a lean premixed flame was investigated, subjected to velocity, equivalence ratio, and both forcing mechanisms simultaneously. In addition, important physical mechanisms controlling the response of partially premixed flames to inlet velocity and equivalence ratio oscillations were analyzed. This final technical report summarizes our findings and major publications stemming from this program.

03 NATURAL GAS↗

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗