Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation optimization”

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,585 records · Page 88

Manufacturability-based optical design optimization for advanced Kirkpatrick–Baez X-ray focusing mirrors

The advanced Kirkpatrick–Baez (AKB) mirror setup is an effective and compelling solution to provide stable X-ray nano-focusing for synchrotron radiation or free-electron laser beamlines. We propose an AKB mirror design optimization approach to mitigate the difficulties associated with mirror fabrication by minimizing the total slope ranges of the four curved mirrors while achieving the expected focusing performance. In the optimization, we have considered geometry constraints to ensure the beam acceptance with the required clear aperture, the diffraction-limited focal size with the adequate numerical aperture, and the desired mirror gaps for adjustment and the necessary working distance for the sample stage. Additionally, practical constraints linked to mirror metrology and fabrication, such as mirror length limits and curvature uncertainty in measurement, are taken into account. Furthermore, progressive objective optimization eliminates the need for any initial guess, fully automating the AKB optimization process. This approach facilitates the development of an elegant Wolter-I or Wolter-III type AKB design solution that satisfies these multiple constraints. In cases where constraints cannot be simultaneously satisfied, the optimization results provide valuable insights into areas where trade-offs need to be considered. Simulations with ray tracing and wavefront propagation validate the optimized AKB design showing high tolerance to the beam incident angle.

36 MATERIALS SCIENCE↗

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗

Computational Engineering of Human Apolipoprotein for Enhanced Interaction with Bacterial Toxins

Sepsis is caused by an over-activation of the immune system in response to the presence of bacterial endotoxins in the bloodstream. LPS exemplifies a potent immune stimulator which can lead to an excessive activation leading to an overly pro-inflammatory reaction and if left untreated results in multi-organ failure and death. Human lipoproteins (HLP) have been shown to neutralize LPS and facilitate its clearance from circulation via the liver. However, host-derived lipoproteins are inhomogeneous in content and size, presenting a challenge for therapeutic applications. In this work, we have used IC facilities in order to computationally engineer HLP for an enhanced interaction with LPS. This work is a continuation of our previous efforts to derive an optimal Nanodisc (ND) lipid mixture for enhanced association with the endotoxin. We have performed fully atomistic Molecular dynamic (MD) simulations in order to carefully identify the most optimal residues in HLP, which can potentially lead to a stronger association to LPS.

42 ENGINEERING↗

Personalized learning via task load optimization

A method for providing task load-optimized computer-generated training experiences to a user of a training system that includes: a display, a training simulator, a prediction program (ML1), and a training optimization program (ML2). In response to receiving a predicted optimal task load, ML2 provides a first training experience recommendation related to the training content and/or training conditions that, if utilized in providing a training experience to the user, is predicted to result in the predicted actual task load of the user equaling the predicted optimal task load. In response to receiving biometric information or performance metric information, ML1 determines the predicted actual task load. If the predicted actual task load does not match the predicted optimal task load, ML2 provides a second training experience recommendation and a second training experience is provided where at least one of the training content or the training conditions is changed.

Bertolli, Michael G.↗

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation and Optimization of Snow Albedo Scheme in Noah-MP Land Surface Model Using In Situ Spectral Observations in the Colorado Rockies

The Biosphere-Atmosphere Transfer Scheme (BATS) ground snow albedo algorithm is commonly used in land-surface models (LSM), weather forecasting and research applications. This study addresses key uncertainties in BATS simulated ground snow albedo within the Noah-MP LSM framework through evaluation and optimization of the Noah-MP BATS ground snow albedo formulation using 2-band (visible and near-infrared (NIR)) in situ albedo observations at Rocky Mountain field stations. The Noah-MP BATS ground snow albedo scheme is extremely sensitive to its input parameters. Namely, an ensemble generated by varying BATS input parameters within potentially plausible ranges provides an average daily range (maximum ensemble member minus minimum ensemble member) of ground snow albedo exceeding 0.45 in visible and NIR bands. Parameter optimization improves agreement between simulated and in situ observed ground snow albedo in visible, NIR and broadband spectrums. Importantly, optimized parameters result in reduced biases relative to observed fresh-snow albedo and better agreement with observed albedo decay. Our analysis across different sites supports that the optimized BATS ground snow albedo parameters are appropriate to transfer in space and time, at least within the region studied (the central-southern Rocky Mountains). The primary error source remaining after parameter optimization is that observed fresh-snow albedo is highly variable, particularly in the NIR spectrum, whereas BATS fresh-snow albedo is constant, an issue which requires further investigation. This study shows significant correlations between observed fresh-snow albedo and surface meteorological conditions (e.g., downward shortwave radiation and temperature) which can support future model development that attempts to include a time-varying formulation for fresh-snow albedo.

54 ENVIRONMENTAL SCIENCES↗

Full-wave simulations on helicon and parasitic excitation of slow waves near the edge plasma

Helicon waves are thought to be promising in various tokamaks, such as DIII-D, because they can penetrate reactor-grade high-density cores and drive the off-axis current with higher efficiency. In the frequency regime ~476 MHz, both slow electrostatic and fast electromagnetic helicon waves can coexist in DIII-D. If the antenna parasitically excites the slow mode, these waves can propagate along the magnetic field line into the scrape-off layer (SOL). Although the importance of the misalignment of the Faraday screen and the electron density in the SOL on the excitation and propagation of slow modes is well known, the conditions for minimizing slow mode excitation have yet to be optimized. Using the Petra-M simulation code in the 2D domain, we analyze the effects of the misalignment of the antenna in the poloidal direction, the misalignment of the Faraday screen in the toroidal direction, and the density in front of the antenna on slow mode generation. Our results suggest that the misalignment of the Faraday screen is a critical factor in reducing the slow mode and that the misalignment angle should be below ~5° to minimize the slow wave excitation. When the electron density is higher than 3.5 x 10 18 m –3 in the SOL, the generation of the slow mode from the antenna is minimized and unaffected by the misalignment of the Faraday screen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗

Designing a boron nitride polyethylene composite for shielding neutrons

Neutrons are encountered in many different fields, including condensed matter physics, space exploration, nuclear power, and healthcare. Neutrons interacting with a biological target produce secondary charged particles that are damaging to human health. The most effective way to shield neutrons is to slow them to thermal energies and then capture the thermalized neutrons. These factors lead us to consider potential materials solutions for neutron shields that maximize the protection of humans while minimizing the shield mass and adapt well to modern additive manufacturing techniques. Using hexagonal boron nitride (hBN) as a capture medium and high-density polyethylene (HDPE) as a thermalization medium, we aim to design the optimal internal structure of h 10 BN/HDPE composites by minimizing the effective dose, which is a measure of the estimated radiation damage exposure for a human. Through Monte Carlo simulations in Geant4, we find that the optimal structure reduces the effective dose up to a factor of 72 over aluminum (Al) and up to a factor of 4 over HDPE; this is a significant improvement in shielding effectiveness that could dramatically reduce the radiation exposure of occupational workers.

36 MATERIALS SCIENCE↗

Baseline Response of Silicon Strip Detectors and Comparison of Cost-Effective Amplifiers

Silicon strip detectors are a proven and more affordable technology. However, they are susceptible to high radiation damage. As part of the characterization of the effects of high radiation on silicon detectors, a baseline performance of these detectors under normal, low radiation conditions needs to be established. We did this by using an optimized setup and an infrared laser to simulate particle incidents. The detectors responded to applied bias voltage as expected. We then compared the performance of both economy amplifiers and an expensive amplifier as part of the optimized setup. For all intents and purposes, the economy amplifier worked just as well as the expensive amplifier.

43 PARTICLE ACCELERATORS↗

Monte-Carlo simulation of noise in hard X-ray Transmission Crystal Spectrometers: Identification of contributors to the background noise and shielding optimization

Transmission crystal spectrometers (TCS) are used on many laser facilities to record hard X-ray spectra. During experiments, signal recorded on imaging plates is often degraded by a background noise. Monte-Carlo simulations made with the code GEANT4 show that this background noise is mainly generated by diffusion of MeV electrons and very hard X-rays. An experiment, carried out at LULI2000, confirmed that the use of magnets in front of the diagnostic, that bent the electron trajectories, reduces significantly this background. The new spectrometer SPECTIX (Spectromètre PETAL à Cristal en TransmIssion X), built for the LMJ/PETAL facility, will include this optimized shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improved quantitative circuit model of realistic patch-based nanoantenna-enabled detectors

Improving the sensitivity of infrared detectors is an essential step for future applications, including satellite- and terrestrial-based systems. We investigate nanoantenna-enabled detectors (NEDs) in the infrared, where the nanoantenna arrays play a fundamental role in enhancing the level of absorption within the active material of a photodetector. The design and optimization of nanoantenna-enabled detectors via full-wave simulations is a challenging task given the large parameter space to be explored. Here, we present a fast and accurate fully analytic circuit model of patch-based NEDs. This model allows for the inclusion of real metals, realistic patch thicknesses, non-absorbing spacer layers, the active detector layer, and absorption due to higher-order evanescent modes of the metallic array. We apply the circuit model to the design of NED devices based on Type II superlattice absorbers, and show that we can achieve absorption of ~70% of the incoming energy in subwavelength (~λ/5) absorber layers. In conclusion, the accuracy of the circuit model is verified against full-wave simulations, establishing this model as an efficient design tool to quickly and accurately optimize NED structures.

47 OTHER INSTRUMENTATION↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Visualizing an Exascale Data Center Digital Twin: Considerations, Challenges and Opportunities

Digital twins are an excellent tool to model, visualize, and simulate complex systems, to understand and optimize their operation. In this work, we present the technical challenges of real-time visualization of a digital twin of the Frontier supercomputer.We show the initial prototype and current state of the twin and highlight technical design challenges of visualizing such a large High Performance Computing (HPC) system. The goal is to understand the use of augmented reality as a primary way to extract information and collaborate on digital twins of complex systems. This leverages the spatio-temporal aspect of a 3D representation of a digital twin, with the ability to view historical and real-time telemetry, triggering simulations of a system state and viewing the results, which can be augmented via dashboards for details. Finally, we discuss considerations and opportunities for augmented reality of digital twins of large-scale, parallel computers.

Maiterth, Matthias↗

Flex Fuel Optimized SI and HCCI Engine

The central objective of the proposed work is to demonstrate an HCCI (homogeneous charge compression ignition) capable SI (spark ignited) engine that is capable of fast and smooth mode transition between SI and HCCI combustion modes. The model-based control technique was used to develop and validate the proposed control strategy for the fast and smooth combustion mode transition based upon the developed control-oriented engine; and an HCCI capable SI engine was designed and constructed using production ready two-step valve-train with electrical variable valve timing actuating system. Finally, smooth combustion mode transition was demonstrated on a metal engine within eight engine cycles. The Chrysler turbocharged 2.0L I4 direct injection engine was selected as the base engine for the project and the engine was modified to fit the two-step valve with electrical variable valve timing actuating system. To develop the model-based control strategy for stable HCCI combustion and smooth combustion mode transition between SI and HCCI combustion, a control-oriented real-time engine model was developed and implemented into the MSU HIL (hardware-in-the-loop) simulation environment. The developed model was used to study the engine actuating system requirement for the smooth and fast combustion mode transition and to develop the proposed mode transition control strategy. Finally, a single cylinder optical engine was designed and fabricated for studying the HCCI combustion characteristics. Optical engine combustion tests were conducted in both SI and HCCI combustion modes and the test results were used to calibrate the developed control-oriented engine model. Intensive GT-Power simulations were conducted to determine the optimal valve lift (high and low) and the cam phasing range. Delphi was selected to be the supplier for the two-step valve-train and Denso to be the electrical variable valve timing system supplier. A test bench was constructed to develop control strategies for the electrical variable valve timing (VVT) actuating system and satisfactory electrical VVT responses were obtained. Target engine control system was designed and fabricated at MSU for both single-cylinder optical and multi-cylinder metal engines. Finally, the developed control-oriented engine model was successfully implemented into the HIL simulation environment. The Chrysler 2.0L I4 DI engine was modified to fit the two-step vale with electrical variable valve timing actuating system. A used prototype engine was used as the base engine and the cylinder head was modified for the two-step valve with electrical VVT actuating system. Engine validation tests indicated that cylinder #3 has very high blow-by and it cannot be reduced with new pistons and rings. Due to the time constraint, it was decided to convert the four-cylinder engine into a single cylinder engine by blocking both intake and exhaust ports of the unused cylinders. The model-based combustion mode transition control algorithm was developed in the MSU HIL simulation environment and the Simulink based control strategy was implemented into the target engine controller. With both single-cylinder metal engine and control strategy ready, stable HCCI combustion was achived with COV of 2.1% Motoring tests were conducted to validate the actuator transient operations including valve lift, electrical variable valve timing, electronic throttle, multiple spark and injection controls. After the actuator operations were confirmed, 15-cycle smooth combustion mode transition from SI to HCCI combustion was achieved; and fast 8-cycle smooth combustion mode transition followed. With a fast electrical variable valve timing actuator, the number of engine cycles required for mode transition can be reduced down to five. It was also found that the combustion mode transition is sensitive to the charge air and engine coolant temperatures and regulating the corresponding temperatures to the target levels during the combustion mode transition is the key for a smooth combustion mode transition. As a summary, the proposed combustion mode transition strategy using the hybrid combustion mode that starts with the SI combustion and ends with the HCCI combustion was experimentally validated on a metal engine. The proposed model-based control approach made it possible to complete the SI-HCCI combustion mode transition within eight engine cycles utilizing the well controlled hybrid combustion mode. Without intensive control-oriented engine modeling and HIL simulation study of using the hybrid combustion mode during the mode transition, it would be impossible to validate the proposed combustion mode transition strategy in a very short period.

33 ADVANCED PROPULSION SYSTEMS↗