Deep Reinforcement Learning for Engineering Design through Topology Optimization of Elementally Discretized Design Domains.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
A design optimization campaign was conducted to search for improved combustion profiles that enhance gasoline compression ignition in a heavy-duty diesel engine with a geometric compression ratio of 17.3. A large-scale design of experiments approach was used for the optimization, employing three-dimensional computational fluid dynamics simulations. The main parameters explored include geometric features, injector specifications, and swirl motion. Both stepped-lip and re-entrant bowls were included in order to assess their respective performance implications. A total of 256 design candidates were prepared using the software package CAESES for automated and simultaneous geometry generation and combustion recipe perturbation. The design optimization was conducted for three engine loads representing light to medium load conditions. The design candidates were evaluated for fuel efficiency, emissions, fuel-air mixing, and global combustion behavior. Simulation results showed that the optimum designs were all stepped-lip bowls, due to improvements in fuel-air mixing, as well as reduced heat loss and emissions formation. Improvements in indicated specific fuel consumption of up to 3.2% were achieved while meeting engine-out NOx emission targets of 1-1.5 g/kW center dot h. Re-entrant bowls performed worse compared to the baseline design, and significant performance variations occurred across the load points. Specifically, the re-entrant bowls were on par with the stepped-lip bowls under light load conditions, but significant deteriorations occurred under higher load conditions. As a final task, selected optimized designs were then evaluated under full-load conditions.
VTMOP is a Fortran 2008 software package containing two Fortran modules for solving computationally expensive bound-constrained blackbox multiobjective optimization problems. VTMOP implements the algorithm of [32], which handles two or more objectives, does not require any derivatives, and produces well-distributed points over the Pareto front. The first module contains a general framework for solving multiobjective optimization problems by combining response surface methodology, trust region methodology, and an adaptive weighting scheme. The second module features a driver subroutine that implements this framework when the objective functions can be wrapped as a Fortran subroutine. Lastly, support is provided for both serial and parallel execution paradigms, and VTMOP is demonstrated on several test problems as well as one real-world problem in the area of particle accelerator optimization.
Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.
In engineering system design, minimizing the variations of the quality measurements while guaranteeing their overall quality up to certain levels, namely the robust parameter design (RPD), is crucial. Recent works have dealt with the design of a system whose response-control variables relationship is a deterministic function with a complex shape and function evaluation is expensive. In this work, we propose a Bayesian optimization method for the RPD of stochastic functions. Dual stochastic response models are carefully designed for stochastic functions. The heterogeneous variance of the sample mean is addressed by the predictive mean of the log variance surrogate model in a two-step approach. We establish an acquisition function that favors exploration across the feasible and optimality-improvable regions to effectively and efficiently solve the stochastic constrained optimization problem. Further, the performance of our proposed method is demonstrated by the extensive numerical and case studies. Note to Practitioners-Many manufacturing processes involve undesirable variations, which create variations in the final products. For example, many emerging manufacturing processes, such as nanomanufacturing, involve complex physical and chemical dynamics and transformation, creating variations in the manufacturing output. In such processes, it is crucial to design the manufacturing processes or products so that they have minimum variations in their quality. Meanwhile, it is also important to maintain the overall quality of the designed processes or products. Furthermore, acquiring data from many advanced manufacturing processes is often very costly, especially in the designing stage. In this work, we propose a data-driven method that automatically finds the best setting of manufacturing processes or products with the minimum variations of quality and a given constraint on the average quality satisfied. Our proposed method is used before conducting every experiment; It analyzes the historical data from previous experiments and provides a setting to be used in the next experiment. Our proposed method efficiently utilizes the historical data, and thus finds the best robust setting by conducting only a small number of experiments.
Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.
A method for optimization of a high-stroke, high-frequency flexural bearing was completed. The main goal was to design a flexural bearing system as the energy restoration component of a one kW free piston linear engine alternator (LEA) operating at 90 Hz and maximum targeted displacement of 22 mm. Defined electricity output of 1 kW, with more than 30% efficiency, required engine operation under high speed and high stroke conditions. Under such challenging stroke and frequency requirements, a method of design was required to select between a wide range of parameters in flexural bearing design. A parametric CAD model was prepared and used in the optimization package of a finite element analysis (FEA) software. Four main parameters of flexural bearing geometry were identified and included the outer diameter, thickness, number of arms, and sweep angle of spiral cuts. Other design parameters including spiral gap width, shape factor and spiral start points angle were shown to be less important compared to these four main parameters in primary design and were examined separately. A neural network algorithm was used to investigate the interactions and effects of design parameters. The sensitivity of each parameter on output criteria such as maximum stress and natural frequency was examined and design charts developed. Methods for further modification of the spiral curve at end points were introduced to reduce the maximum stress level by up to 20%. Several steel alloys and a titanium-based alloy were identified as proper material candidates for high cycle applications and the final designed spring was manufactured and tested to validate the FEA results of axial stiffness and strain level throughout the flexure arms. The percent of the flexure’s moving mass was estimated to be 29% with FEA analysis and compared well with the experimental value of 31%. The stress/strain results of the FEA analysis matched with the experimental results – within ±5%.
This report was prepared as an account of work sponsored by the Office of Fossil Energy and Carbon Management of U.S. Department of Energy under Funding Opportunity Announcement Number DE-FOA-0002332 “Energy Storage for Fossil Power Generation”. The work aimed to explore and advance an innovative hydrogen energy storage system – the synergistically integrated hydrogen energy storage system (SIHES) – that has the following characteristics: • Compatible with existing or new coal and gas fuel electricity generation units, • best suited for intermediate to long duration energy storage, from 12 hours to weeks even months, and • capable of storing energy at the utility scale – hundreds of MWh to GWh energy storage with power output level in tens to hundreds of MW. Preliminary front-end engineering design (Pre-FEED) studies was carried out to develop and refine a site-specific SIHES as peaking power generation units (so named as HyPeaker) as the first market entry point, to demonstrate both the technical feasibility and the economic viability to integrate the HyPeaker “within the fence” of a fossil power plant. This specific site was TVA’s Johnsonville Combustion Turbine Plant. The HyPeaker was designed and engineered to integrate with a 60MW aeroderivative gas turbine unit already available at TVA’s Johnsonville site. This site-specific HyPeaker consists of an alkaline electrolyzer to produce hydrogen from CO 2 free electricity sources, an innovative low-cost high-pressure hydrogen storage system (Big-Ton) and the aero gas turbine to generate electricity using blend of hydrogen and natural gas. A holistic system level technoeconomic analysis tool specific to HyPeaker was developed to optimize the engineering design of the Johnsonville site-specific HyPeaker for cost and performance. The optimal design and specification of the Johnsonville site-specific HyPeaker are the following: • Alkaline electrolyzer: 3MW • Big-Ton storage vessel: 11,000kg H 2 at 3000psi. • 4-stage diaphragm hydrogen compressor: 55kg-H 2 /hr from 150psi to 3000psi. The HyPeaker is designed to provide sufficient hydrogen for 90% continuous operation of the HyPeaker. All major components have design life of 30 years. The capex of HyPeaker is estimated at $\$$7.1M. This included $\$$1.5M for the electrolyzer, $\$$5.6M for the storage vessel and compressor. The cost of aero gas turbine was included as it is already available at the site. Key findings are: • HyPeaker can be designed, manufactured, installed and integrated with the fossil power plants, with sub-systems and components commercially available on the market today, even when it is scaled up to an order of magnitude larger than the one at the Johnsonville site. HyPeaker is a technologically viable solution to cover a wide range of energy storage duration needs, from daily peaking operation to seasonal shifting for fossil fueled assets. • The cost advantage of SCCV based Big-Ton H 2 storage vessel made it possible to “oversize” the H 2 storage subsystem to achieve overall system level cost optimization. The benefits are two-fold. First, it allows to significantly reduce the capacity and cost of electrolyzer by spreading H 2 production over a much longer period of time when the fuel cost for electricity production is low. Second, it allows to balance the hydrogen production and usage shift over weeks to months to meet the peak demands. As such, the capital cost of HyPeaker system using the Big-Ton was less than half of the cost of a system with today’s steel tube based H 2 storage system. The HyPeaker has even better cost advantage Li-ion battery based energy storage system. The estimated capital cost of Li-Ion battery system would be at $\$$38M, under the same projected 20-year electricity generation profile of the Johnsonville site. This is over 5 times more expensive than the HyPeaker system. • Since industry scale energy storage systems do not have 100% energy conversion and storage efficiency, energy storage systems using fossil fuel generated electricity would increase the CO 2 emission. This is particularly the case for HyPeaker due to its low round trip efficiency. Therefore, a more sensible solution would be to the excessive or curtailed electricity from CO 2 emission free sources such as solar farms, wind farms or nuclear power plants, to produce hydrogen, and integrate them with the HyPeaker. Electricity from TVA’s nuclear power plants was used for the Johnsonville HyPeaker. • The economic viability of HyPeaker is expected to be further improved when global supply chains are taken into consideration. For the same Johnsonville site specific HyPeaker, the capex would be reduced to ~$\$$3.6M from ~$\$$7.1M, and the added LCOE is reduced to ~$\$$85/MWh. With the bipartisan Infrastructure Investment and Jobs Act, the cost of domestically produced HyPeaker sub-systems would be at the level of today’s global suppliers. Since the peaking units generally operate at peak usage period, thereby demanding higher price, the projected $\$$85/MWh LCOE would be within the realm of financial viability for utility operators.
In this work, a novel design optimization technique based on active learning, which involves dynamic exploration and exploitation of the design space of interest using an ensemble of machine learning algorithms, is presented. In this approach, a hybrid methodology incorporating an explorative weak learner (regularized basis function model) that fits high-level information about the response surface and an exploitative strong learner (based on committee machine) that fits finer details around promising regions identified by the weak learner is employed. For each design iteration, an aristocratic approach is used to select a set of nominees, where points that meet a threshold merit value as predicted by the weak learner are selected for evaluation. In addition to these points, the global optimum as predicted by the strong learner is also evaluated to enable rapid convergence to the actual global optimum once the most promising region has been identified by the optimizer. Additionally, this methodology is first tested by applying it to the optimization of a two-dimensional multi-modal surface and, subsequently, to a complex internal combustion (IC) engine combustion optimization case with nine control parameters related to fuel injection, initial thermodynamic conditions, and in-cylinder flow. It is found that the new approach significantly lowers the number of function evaluations that are needed to reach the optimum design configuration (by up to 80%) when compared to conventional optimization techniques, such as particle swarm and genetic algorithm-based optimization techniques.
Engine knock remains one of the major barriers to further improve thermal efficiency of Spark Ignition (SI) engines. Knock can be suppressed by lowering the compression ratio, or retarding the spark ignition timing, however, at an expense of efficiency penalty. SI engine is usually operated at knock-limited spark advance (KLSA) to achieve possibly maximum efficiency with given engine hardware and fuel properties, such as Research Octane Number (RON), Motor Octane Number (MON), and heat of vaporization, etc. Co-optimization of engine design and fuel properties is promising to improve the engine efficiency and predictive CFD models can be used to facilitate this optimization process. However, difficulties exist in predicting KLSA in CFD simulations. First, cyclic variability of SI engine demands that multi-cycle results are required to capture the extreme conditions. Secondly, Mach Courant-Friedrichs-Lewy (CFL) number of 1 is desired to accurately predict the knock intensity (KI), resulting in unaffordable computational cost, especially for multi-cycle simulations. In this study, a new approach to numerically predict KLSA using large Mach CFL number of 50 is proposed. This approach is validated against experimental data for a boosted Direct Injection Spark Ignition (DISI) engine at multiple loads and spark timings. G-equation combustion model coupled with well-mixed chemical kinetic model are used to predict the turbulent flame propagation and end-gas auto-ignition, respectively. Simulations run for 10 consecutive engine cycles at each condition. The results show good agreement between model predictions and experiments in terms of cylinder pressure, combustion phasing and cyclic variation. Engine knock is predicted with early spark ignition timing, indicated by significant pressure wave oscillation and end-gas heat release. Maximum Amplitude of Pressure Oscillation (MAPO) analysis is performed to quantify the KI, and the slope change point in KI extrema is used to indicate the KLSA accurately. Using a smaller Mach CFL number of 5 also results in the same conclusions thus demonstrating that this approach is insensitive to the Mach CFL number. The use of large Mach CFL number allows us to achieve fast turn-around time for multi-cycle engine CFD simulations.
Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).
Advanced reactor concepts currently under development in the United States [e.g., the High-Temperature Gas Reactor (HTGR) of X-Energy and the Molten Chloride Fast Reactor (MCFR) of TerraPower] are striving to improve the economics of nuclear power primarily by using passively and ‘walk-away’ safe technologies in the reactor. The recently published report, ‘Future of Nuclear Energy in a Carbon-Constrained World’ by the Massachusetts Institute of Technology Energy Institute (MIT EI) has shown that a significant contributor to the capital cost of NPPs is the ‘civil works’ i.e., the construction of the balance of plant (BoP: all SSCs except those involved in power generation, e.g., reactor vessel) such as the buildings, containment dome, foundation, etc., which amount to almost half of the overnight capital cost. The report suggests that using advanced construction technologies such as seismic isolation and high-performance concrete to significantly reduce the capital costs of advanced reactors. This report demonstrates the application of two cost-cutting solutions on advanced reactor designs that are currently under development. These solutions include, (1) seismic base isolation, and (2) risk- and cost-based seismic design optimization including seismic isolation of individual components.
The work in this paper describes neural network surrogate models for calculating the effective mechanical properties of a periodic composites. The models achieve good accuracy even when only provided with training data sampling a small portion of the design space. As an example, the surrogate models are applied to solving the inverse design problem of finding structures with optimal mechanical properties. The surrogate models are sufficiently accurate to recover optimal solutions in general agreement with established topology optimization methods. However, improvements will be required to develop robust, efficient neural network-based surrogate models and several directions for future research are highlighted here.
An analytical and experimental investigation of air-cooled condensers for small-scale ammonia-water absorption heat pumps is presented. Conventional round-tube corrugated fin design and a novel multi-pass tube array design are considered for a condenser operating at severe ambient conditions within a compact footprint. Segmental heat and mass transfer models are developed to gain a deeper understanding of the characteristics of each design. Predicted heat duties are within 5% of the measured values. The limitations of correlations to predict air-side pressure drop in corrugated fins are identified. The performance of the novel multi-pass tube array condenser was comparable to that of conventional condensers, but is determined to be air-side limited throughout length of the heat exchanger. An improved multi-pass condenser design with fins for increased air-side area is proposed using the insights gained from the modeling and experiments. In closing this design meets the heat duty requirements at severe ambient conditions, while decreasing the size of the heat exchanger by 33%.
Multiobjective simulation optimization (MOSO) problems are optimization problems with multiple conflicting objectives, where evaluation of at least one of the objectives depends on a black-box numerical code or real-world experiment, which we refer to as a simulation. Whereas an extensive body of research is dedicated to developing new algorithms and methods for solving these and related problems, it is challenging and time-consuming to integrate these techniques into real-world production-ready solvers. This is partly because of the diversity and complexity of modern state-of-the-art MOSO algorithms and methods and partly because of the complexity and specificity of many real-world problems and their corresponding computing environments. The complexity of this problem is only compounded when introducing potentially complex and/or domain-specific surrogate-modeling techniques, problem formulations, design spaces, and data acquisition functions. Here, this paper carefully surveys the current state of the art in MOSO algorithms, techniques, and solvers, as well as problem types and computational environments where MOSO is commonly applied. We then present several key challenges in the design of a parallel multiobjective simulation optimization framework (ParMOO) and how they have been addressed. Finally, we provide two case studies demonstrating how customized ParMOO solvers can be quickly built and deployed to solve real-world MOSO problems.
This paper formulates laws for scaling wind turbine rotors. Although the analysis is general, the article primarily focuses on the subscaling problem, i.e., on the design of a smaller-sized model that mimics a full-scale machine. The present study considers both the steady-state and transient response cases, including the effects of aerodynamic, elastic, inertial, and gravitational forces. The analysis reveals the changes to physical characteristics induced by a generic change of scale, indicates which characteristics can be matched faithfully by a subscaled model, and states the conditions that must be fulfilled for desired matchings to hold. Based on the scaling laws formulated here, the article continues by considering the problem of designing scaled rotors that match desired indicators of a full-scale reference. To better illustrate the challenges implicit in scaling and the necessary tradeoffs and approximations, two different approaches are contrasted. The first consists in a straightforward geometric zooming. An analysis of the consequences of zooming reveals that, although apparently simple, this method is often not applicable in practice, because of physical and manufacturing limitations. This motivates the formulation of scaling as a constrained optimal aerodynamic and structural matching problem of wide applicability. Practical illustrations are given considering the scaling of a large reference 10 MW wind turbine of about 180 m in diameter down to three different sizes of 54, 27, and 2.8 m. Results indicate that, with the proper choices, even models characterized by very significant scaling factors can accurately match several key performance indicators. Additionally, when an exact match is not possible, relevant trends can at least be captured.
Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.
This study delineates a meticulous exploration of technologies to enhance the energy efficiency of rooftop air conditioning units, employing the DOE/ORNL heat pump design model for comprehensive engineering design and optimization. A baseline rooftop air conditioning unit, featuring a 13 ton (45.7 kW) cooling capacity and a 17.9 integrated energy efficiency ratio, served as the point of departure for substantive efficiency enhancements. Key modifications included the consolidation of two refrigerant circuits into one, integrating three parallel 2-stage (dual-speed) compressors, fan replacements with high-efficiency substitutes. Notably, a lower global warming potential refrigerant, R452B, was evaluated as a substitute for R-410A, demonstrating better performance in the lab prototype. Further, the achieved measured integrated energy efficiency ratio of 21.4 in the lab prototype surpassed the baseline integrated energy efficiency ratio. Comparative evaluations between R410A and R452B indicated heightened efficiency with the latter, showcasing a lab-demonstrated integrated energy efficiency ratio of 22.4 at the rated capacity of 13.8 ton (48.5 kW) and 23.9 integrated energy efficiency ratio at the rated capacity of 10 ton (35.2 kW). This research underscores the successful development of a rigorous, energy efficient rooftop air conditioning unit prototype with noteworthy environmental and economic implications.