Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak estimation”

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 469 records · Page 26

Cool Star Binaries with ALEXIS

We proposed to search for high-temperature, flare-produced Fe XXIII line emission from active cool star binary systems using the ALEXIS all-sky survey. Previous X-ray transient searches with ARIEL V and HEAO-1, and subsequent shorter duration monitoring with the GINGA and EXOSAT satellites demonstrated that active binaries can produce large (EM approximately equals 10(exp 55-56/cu cm) X-ray flares lasting several hours or longer. Hot plasma from these flares at temperatures of 10(exp 7)K or more should produce Fe XXIII line emission at lambda = 132.8 A, very near the peak response of ALEXIS telescopes 1A and 2A. Our primary goals were to estimate flare frequency for the largest flares in the active binary systems, and, if the data permitted, to derive a distribution of flare energy vs. frequency for the sample as a whole. After a long delay due to the initial problems with the ALEXIS attitude control, the heroic efforts on the part of the ALEXIS satellite team enabled us to carry out this survey. However, the combination of the higher than expected and variable background in the ALEXIS detectors, and the lower throughput of the ALEXIS telescopes resulted in no convincing detections of large flares from the active binary systems. In addition, vignetting-corrected effective exposure times from the ALEXIS aspect solution were not available prior to the end of this contract; therefore, we were unable to convert upper limits measured in ALEXIS counts to the equivalent L(sub EUV).

Stern, Robert A.↗

Photospheric Radius Expansion and a Double-peaked Type-I X-Ray Burst from GRS 1741.9-2853

We present an analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star (NS) low-mass X-ray binary GRS 1741.9–2853 during a period of outburst in 2020 May. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe Kα emission, which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of L(bol) = 7.03(+0.04, -0.05) x 10^(36) erg/s, assuming a distance of 7 kpc, corresponding to an Eddington ratio of 4.5%. This persistent luminosity combined with light-curve analysis leads us to infer that the bursts were the result of pure He burning rather than mixed H/He burning. Time-resolved spectroscopy reveals that the bolometric flux of the first burst exhibits a double-peaked structure, placing the source within a small population of accreting NSs that exhibit multiple-peaked type-I X-ray bursts. We find that the second, brighter burst shows evidence for photospheric radius expansion (PRE) and that at its peak, this PRE event had an unabsorbed bolometric flux of F(peak)= 2.94 (+0.28, -0.26) 10^(-8) erg/sq.cm s. This yields a new distance estimate of d = 9.0 ± 0.5 kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical NS. Additionally, we performed a detailed timing analysis that failed to find evidence for quasi-periodic oscillations or burst oscillations, and we place an upper limit of 16% on the rms variability around 589 Hz, the frequency at which oscillations have previously been reported.

Sean N. Pike↗

The Diurnal Dynamics of Gross Primary Productivity Using Observations From the Advanced Baseline Imager on the Geostationary Operational Environmental Satellite‐R Series at an Oak Savanna Ecosystem

Gross primary productivity (GPP) is the largest flux in the global carbon cycle and satellite-based GPP estimates have long been used to study the trends and interannual variability of GPP. With recent updates to geostationary satellites, we can now explore the diurnal variability of GPP at a comparable spatial resolution to polar-orbiting satellites and at temporal frequencies comparable to eddy covariance (EC) tower sites. We used observations from the Advanced Baseline Imager on the Geostationary Operational Environmental Satellite-R series (GOES-R) to test the ability of subdaily satellite data to capture the shifts in the diurnal course of GPP at an oak savanna EC site in California, USA that is subject to seasonal soil moisture declines. We compared three methods to estimate GPP: (a) a light-use efficiency model, (b) a linear relationship between the product of near-infrared reflectance of vegetation and photosynthetically active radiation (LIN-NIR v P) and EC tower GPP, and (c) a light response curve (LRC-NIR v P) between NIR v P and EC GPP. The LRC-NIR v P achieved the lowest mean absolute error for winter (2 µmol CO 2 m –2 s –1 ), spring (2.51 µmol CO 2 m –2 s –1 ), summer (1.43 µmol CO 2 m –2 s –1 ), and fall (1.35 µmol CO 2 m –2 s –1 ). The ecosystem experienced the largest shift in daily peak GPP in relation to the peak of incoming solar radiation toward the morning hours during the dry summers. The LRC-NIRvP and the light-use efficiency model were in agreement with these patterns of a shift in peak daily GPP toward the morning hours during summer. Our results can help develop diurnal estimates of GPP from geostationary satellites that are sensitive to fluctuating environmental conditions during the day.

54 ENVIRONMENTAL SCIENCES↗

Histogrammatic Method for Determining Relative Abundance of Input Gas Pulse

To satisfy the Major Constituents Analysis (MCA) requirements for the Vehicle Cabin Atmosphere Monitor (VCAM), this software analyzes the relative abundance ratios for N2, O2, Ar, and CO2 as a function of time and constructs their best-estimate mean. A histogram is first built of all abundance ratios for each of the species vs time. The abundance peaks corresponding to the intended measurement and any obfuscating background are then separated via standard peak-finding techniques in histogram space. A voting scheme is then used to include/exclude this particular time sample in the final average based on its membership to the intended measurement or the background population. This results in a robust and reasonable estimate of the abundance of trace components such as CO2 and Ar even in the presence of obfuscating backgrounds internal to the VCAM device. VCAM can provide a means for monitoring the air within the enclosed environments, such as the ISS (International Space Station), Crew Exploration Vehicle (CEV), a Lunar Habitat, or another vehicle traveling to Mars. Its miniature pre-concentrator, gas chromatograph (GC), and mass spectrometer can provide unbiased detection of a large number of organic species as well as MCA analysis. VCAM s software can identify the concentration of trace chemicals and whether the chemicals are on a targeted list of hazardous compounds. This innovation s performance and reliability on orbit, along with the ground team s assessment of its raw data and analysis results, will validate its technology for future use and development.

Mandrake, Lukas↗

Analytical Predictions of Thermal Stress in the Stardust PICA Heatshield Under Reentry Flight Conditions

We performed finite element analyses on a model of the Phenolic Impregnated Carbon Ablator (PICA) heatshield from the Stardust sample return capsule (SRC) to predict the thermal stresses in the PICA material during reentry. The heatshield on the Stardust SRC was a 0.83 m sphere cone, fabricated from a single piece of 5.82 cm-thick PICA. The heatshield performed successfully during Earth reentry of the SRC in January 2006. Material response analyses of the full, axisymmetric PICA heatshield were run using the Two-Dimensional Implicit Ablation, Pyrolysis, and Thermal Response Program (TITAN). Peak surface temperatures were predicted to be 3385K, while the temperature at the PICA backface remained at the estimated initial cold-soak temperature of 278K. Surface recession and temperature distribution results from TITAN, at several points in the reentry trajectory, were mapped onto an axisymmetric finite element model of the heatshield. We used the finite element model to predict the thermal stresses in the PICA from differential thermal expansion. The predicted peak compressive stress in the PICA heatshield was 1.38 MPa. Although this level of stress exceeded the chosen design limit for compressive stresses in PICA tiles for the design of the Orion crew exploration vehicle heatshield, the Stardust heatshield exhibited no obvious mechanical failures from thermal stress. The analyses of the Stardust heatshield were used to assess and adjust the level of conservatism in the finite element analyses in support of the Orion heatshield design.

Squire, Thomas↗

Arctic sea ice signatures for passive microwave algorithms

Two new sets of global pure type signatures are defined, and the precision of the concentration estimates made using the Team algorithm and the new signatures is estimated. One set varies temporally but not spatially and is defined by the modes of major peaks in histograms of brightness temperature. The other set is defined as the overall means of the temporally varying signatures; these are called 'new constant signatures'. Both sets of new signatures have essentially the same precision: 0.05 (fractional area) for winter ice concentration, 0.10 for summer ice concentration, and 0.15 to 0.3 for winter multilayer ice concentrations. Both sets result in about 0.6 to 0.7 winter multiyear ice for the Arctic Ocean, substantially more than the 0.3 to 0.4 computed using the original team algorithm signatures, and in the same range as several independent estimates of the multiyear ice, 0.54 to 0.75.

Thomas, D. R.↗

Cloud top height variability of strong convective cells

Cloud top height variability of cloud domes above the cirrostratus anvil was investigated for two severe local storm occurrences over Texas on May 12 and 13, 1972. Twenty-three domes were monitored and the average maximum height was 1.6 km above the anvil. It was estimated that a dome would take a maximum of about 6 min to grow from the anvil level to peak altitude based on measurements over the last 2 min of the growth period when the upward vertical velocity averaged 4 m/sec. The extreme vertical growth rate was 27 m/sec for a 30-sec period. This vertical growth rate compared with a range of 11.4-19.4 m sec determined from a method of estimating vertical motions based on laboratory measurements. Average dome collapse was at the same rate as the average ascent for the first 90 sec after the maximum altitude was reached.

Shenk, W. E.↗

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗

Electronic structures and pseudogap nature of the strongly correlated Kondo semimetals and insulators CeNiSn and CeRhX(X=As, Sb)

Employing temperature (T)-dependent angle-resolved photoemission spectroscopy (ARPES) near the Ce 4d absorption edge, we have investigated the electronic structures of the strongly correlated Kondo insulators (KIs) and/or Kondo semimetals of CeNiSn and CeRhX (X=As, Sb), which also belong to potential topological KIs (TKIs). Good agreement is found between the Fermi surfaces (FSs) measured at the Ce 4d→4f on resonance and those calculated and unfolded into the Ce-only Brillouin zone, supporting the crucial Ce 4f contribution to the Fermi-edge states in them. While the FSs of CeRhSb and CeNiSn are similar to each other, they are quite different from those of CeRhAs, which is understood to originate from the semimetallic Kondo ground states of CeNiSn and CeRhSb in contrast to the insulating Kondo ground state of CeRhAs. T-dependent ARPES demonstrates the existence of the Kondo resonance above the Fermi level (EF) and its T-driven decoherence in CeNiSn and CeRhSb, and the existence of the Ce 4f Kondo-like peak well below EF with a pseudogap in CeRhAs. With increasing T, the tail of the Kondo resonance in CeNiSn and CeRhSb loses its coherence, whereas the Ce 4f Kondo-like peak in CeRhAs persists up to T≥200 K. The Kondo temperatures (TKs) are estimated via the analysis of T-dependent ARPES, yielding TK≈80 K (CeNiSn), TK≈180 K (CeRhSb), and TK≫200 K (CeRhAs). High-resolution, low-T ARPES confirms the Kondo resonance just above EF in CeNiSn and CeRhSb Kondo semimetals and the KI ground state of CeRhAs with a pseudogap of Δ≈30 meV, demonstrating the importance of the coherent Kondo states in determining their potential topological properties. The differences in the T-dependent ARPES of CeNiSn and CeRhSb suggest the weaker f-c hybridization in CeNiSn than in CeRhSb (c: conduction electron), which is likely due to the more localized Ni 3d orbitals than the delocalized Rh 4d orbitals. The differences between the T-dependent ARPES of CeNiSn/CeRhSb Kondo semimetals and those of CeRhAs Kondo insulator are likely due to the smaller volume of CeRhAs than those of CeNiSn/CeRhSb.

Kang, J-S↗

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) Version 4.1 Technical Reference Manual

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 4.1. The two most significant changes for HyRAM+ version 4.1 from HyRAM+ version 4.0 are direct incorporation of unconfined overpressure into the QRA calculations and modification of the models for cryogenic liquid flow through an orifice. In QRA mode, the user no longer needs to input peak overpressure and impulse values that were calculated separately; rather, the unconfined overpressure is estimated for the given system inputs, leak size, and occupant location. The orifice flow model now solves for the maximum mass flux through the orifice at constant entropy while conserving energy, which does not require a direct speed of sound calculation. This does not affect the mass flow for all-gaseous releases; the method results in the same speed of sound for choked flow. However, this method does result in a higher (and more realistic) mass flow rate for a given leak size for liquid releases than was previously calculated.

08 HYDROGEN↗

Estimation of the vortex length scale and intensity from two-dimensional samples

A method is proposed for estimating flow features that influence flame wrinkling in reciprocating internal combustion engines, where traditional statistical measures of turbulence are suspect. Candidate methods were tested in a computed channel flow where traditional turbulence measures are valid and performance can be rationally evaluated. Two concepts are tested. First, spatial filtering is applied to the two-dimensional velocity distribution and found to reveal structures corresponding to the vorticity field. Decreasing the spatial-frequency cutoff of the filter locally changes the character and size of the flow structures that are revealed by the filter. Second, vortex length scale and intensity is estimated by computing the ensemble-average velocity distribution conditionally sampled on the vorticity peaks. The resulting conditionally sampled 'average vortex' has a peak velocity less than half the rms velocity and a size approximately equal to the two-point-correlation integral-length scale.

Reuss, D. L.↗

Luminous radio-quiet sources in the W3(MAIN) cloud core

We have resolved 450 micrometer and 800 micrometer emission from the W3(Main) star forming region into three major peaks, using 8 inch - 14 inch beams with the James Clerk Maxwell Telescope on the summit of Mauna Kea. One of the submillimeter sources is identified with W3 - IRS5, a well-known candidate protostar. However, to our surprise, we find that none of the submillimeter peaks coincides with any of the prominent compact HII regions in the area. We estimate that the three submillimeter sources together contribute 30-50 percent of the total bolometric luminosity of the region, and speculate that the contribution of luminous radio-quiet sources to the total luminosity of HII region/molecular cloud complexes may be larger than is often assumed.

Wynnwilliams, C. Gareth↗

Venus O2 visible and IR nightglow: Implications for lower thermosphere dynamics and chemistry

The National Center for Atmospheric Research thermospheric general circulation model for the Venus thermosphere is modified to examine two observed night airglow features, both of which serve as sensitive tracers of the thermospheric circulation. New O2 nightglow data from the Pioneer Venus Orbiter (PVO) star tracker (O2 Herzberg II at 400-800 nm) and ground-based telescopes (O2 IR at 1.27 microns) yield additional model constraints for estimating Venus winds over 100-130 km. Atomic oxygen, produced by dayside CO2 photolysis peaking near 110 km, and transported to the nightside by the global wind system, is partially destroyed through three-body recombination, yielding the O2 Herzberg II visible nightglow. This emission is very sensitive to horizontal winds at altitudes between 100 and 130 km. Other trace species catalytic reactions also contribute to the production of the very strong nightside infrared (1.27 microns) emission. This paper examines the dynamical and chemical implications of these new data using the Venus thermospheric general circulation model (VTGCM) as an analysis tool. Three-dimensional calculations are presented for both solar maximum and solar medium conditions, corresponding to early PVO (1979-1981) and PVO entry (mid-1992) time periods. Very distinct periods are identified in which zonal winds are alternately weak and strong in the Venus lower thermosphere. VTGCM sensitivity studies are conducted to assess the impacts of potential changes in thermospheric zonal and day-to-night winds, and eddy diffusion on the corresponding nightglow intensities. It appears that cyclostrophic balance extends above 80 km periodically, owing to a reversal of the upper mesosphere latitudinal temperature gradient, and thereby producing strong zonal winds and correspondingly modified O2 nightglow distributions that are observed.

Bougher, S. W.↗

Comparison of Postural Recovery Following Short and Long Duration Spaceflights

INTRODUCTION: Post-flight postural ataxia reflects adaptive changes to vestibulo-spinal reflexes and control strategies adopted for movement in weightlessness. Quantitative measures obtained during computerized dynamic posturography (CDP) from US and Russian programs provide insight into the effect of spaceflight duration in terms of both the initial decrements and recovery of postural stability. METHODS: CDP was obtained on 117 crewmembers following Shuttle flights lasting 4-17 days, and on 64 crewmembers following long-duration missions lasting 48-380 days. Although the number and timing of sessions varied, the goal was to characterize postural recovery pooling similar measures from different research and flight medicine programs. This report focuses on eyes closed, head erect conditions with either a fixed or sway-referenced base of support. A smaller subset of subjects repeated the sway-referenced condition while making pitch head movements (+/- 20deg at 0.33Hz). Equilibrium scores were derived from peak-to-peak anterior-posterior sway. Fall probability was modeled using Bayesian statistical methods to estimate parameters of a logit function. RESULTS: The standard Romberg condition was the least sensitive. Longer duration flights led to larger decrements in stability with sway-reference support during the first 1-2 days, although the timecourse of recovery was similar across flight duration with head erect. Head movements led to increased incidence of falls during the first week, with a significantly longer recovery following long duration flights. CONCLUSIONS: The diagnostic assessment of postural instability, and differences in the timecourse of postural recovery between short and long flight durations, are more pronounced during unstable support conditions requiring active head movements.

Wood, S. J.↗

CGRO/BATSE Data Support the New Paradigm For GRB Prompt Emission and the New L-i(nTh)-E-peak,i(nTh,rest) Relation

The paradigm for gamma-ray burst (GRB) prompt emission is changing. Since early in the Compton Gamma RayObservatory (CGRO) era, the empirical Band function has been considered a good description of the keV-MeV-gamma-ray prompt emission spectra despite the fact that its shape was very often inconsistent with the theoretical predictions, especially those expected in pure synchrotron emission scenarios. We have recently established a new observational model analyzing data of the NASA Fermi Gamma-ray Space Telescope. In this model, GRB prompt emission would be a combination of three main emission components: (i) a thermal-like component that we have interpreted so far as emission from the jet photosphere, (ii) a non-thermal component that we have interpreted so far as either synchrotron radiation from the propagating and accelerated charged particles within the jet or reprocessed jet photospheric emission, and (iii) an additional non-thermal (cutoff) power law (PL) extending from low to high energies in gamma-rays and most likely of inverse Compton origin. In this article we reanalyze some of the bright GRBs, namely GRBs 941017, 970111, and 990123, observed with the Burst And Transient Source Experiment (BATSE) on board CGRO with the new model. We conclude that BATSE data for these three GRBs are fully consistent with the recent results obtained with Fermi: some bright BATSE GRBs exhibit three separate components during the prompt phase with similar spectral parameters as those reported from Fermi data. In addition, the analysis of the BATSE GRBs with the new prompt emission model results in a relation between the time-resolved energy flux of the non-thermal component, F(in)(Th), and its corresponding nuFnu spectral peak energy,Epeak,inTh (i.e., FinThEpeak,inTh ), which has a similar index when fitted to a PL as the one initially derived from Fermi data. For GRBs with known redshifts (z) this results in a possible universal relation between the luminosity of the non-thermal component, LinTh, and its corresponding nuFnu spectral peak energy in the rest frame, Epeak,iNT,rest (i.e.,LinThEpeak,iNT,rest). We estimated the redshifts of GRBs 941017 and 970111 using GRB 990123 with z = 1.61as a reference. The estimated redshift for GRB 941017 is typical for long GRBs and the estimated redshift for GRB970111 is right in the range of the expected values for this burst.

acceleration of particles – black hole physics â↗

A Spaceborne Demonstration of P-Band Signals-of-Opportunity (SoOp) Reflectometry

Land-reflected signals from a geosynchronous communication satellite broadcasting in P-band (367.5 MHz) were captured in low Earth orbit (LEO) using a simple dipole antenna. A delay-Doppler map (DDM) was generated through autocorrelation. Estimates of the specular point delay were obtained from the lag of the second peak in the DDM with a bias of 239.4 m and a standard deviation of 44 m (12 m over a frozen lake) with respect to a predicted orbit model. Relative magnitudes of the first and second DDM peaks fell within the range of values predicted using dielectric models for the frozen ground and lake. Finally, retrievals of surface reflection coefficient were generated using a range of realistic values for the transmitter link budget G/T , and these also fell within the range of possible values for the antenna gain pattern. Given the lack of calibration and the large uncertainties in the receiver orbit and attitude, this agreement is sufficient to conclude a successful demonstration of the fundamental principle of single-antenna reflectometry in P-band. P-band reflectometry may offer a new approach to remote sensing of sub-canopy and root-zone soil moisture.

Delays↗

A verification study of fatigue-based methods in API RP 1183 for estimating fatigue life of pipeline dents

Pipeline dents are mechanical damage caused by a third party in buried pipelines during construction or maintenance repair. Dents are recognized as a leading cause of pipeline failures for both liquid and gas transmission lines, and their assessment and management are critical to pipeline integrity. To meet the industry need, American Institute of Petroleum (API) in 2020 published a recommended practice (RP) 1183 – Assessment and Management of Pipeline Dents. This API code provides the pipeline industry “standard methods” for evaluating severity and fatigue life of pipeline dents with a single peak. This paper performed an initial verification study of fatigue-based methods prescribed in API RP 1183 for estimating fatigue life of pipeline dents. A brief review was first given to the dent fatigue-based methods, including three dent screening methods for estimating dent fatigue life in Sections 7.4.1, 7.4.2, and 7.4.3 of API 1183 and one dent assessment method for predicting fatigue life in Section 8.3.4. In order to verify these fatigue life prediction methods, a set of “standard examples” of dents provided in Annexes A.1 to A.5 of API RP 1183 were utilized in this study to calculate fatigue life, to determine restraint condition, and to evaluate both screening and assessment methods for single-peak dents. Through detailed calculations and comparisons, this study confirmed some specifications of API RP 1183, but also found some self-inconsistences and errors existing in this new API code. Furthermore, this includes the self-inconsistences between three screening methods and the self-inconsistence between the screening and assessment methods. On this basis, recommendations were made on how the operators should use API RP 1183, and what actions API should take to improve this code.

42 ENGINEERING↗

Model predictive combustion control of a Gasoline Compression Ignition engine

Gasoline Compression Ignition is a novel combustion concept that derives its superiority from the high compression ratio of a compression ignition engine as well as the properties of gasoline fuel, such as longer ignition delay and higher volatility compared to diesel. Here, this combustion concept was experimentally tested on a 12.4L Class 8 truck engine. Based on these experimental data, prior efforts by the authors focused on the development of an engine model for a heavy-duty engine operating on a low-reactivity fuel. This engine model was leveraged within this study to investigate a combustion control strategy at different engine conditions and injection methods and was augmented to incorporate cycle-to-cycle combustion variations. State estimation is performed by means of a Kalman filter which feeds into a model predictive controller. The model predictive controller chooses control actions based on a predefined cost function under consideration of bounds reflecting physical constraints. The engine model was utilized to establish a state-space model that serves the Kalman filter and model predictive controller for estimation and prediction. A comparative study investigating control actions and engine behavior was performed with and without limiting in-cylinder peak pressure as well as combustion noise, which is of particular interest for early pilot injection strategies. In addition, the proposed control architecture was investigated at two different levels of cycle-to-cycle variations and compared to the performance of a control structure with input disturbance rejection. For increased cycle-to-cycle variations, disturbance estimation reduces state fluctuations and control effort. In general, this investigation highlights control aspects specific to a compression-ignited combustion regime with low-reactivity fuel. The control algorithm is able to maintain the desired references for brake mean effective pressure and combustion phasing while controlling peak in-cylinder pressure and combustion noise.

42 ENGINEERING↗