Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

Single-layered organic photovoltaics with double cascading charge transport pathways: 18% efficiencies

The chemical structure of donors and acceptors limit the power conversion efficiencies achievable with active layers of binary donor-acceptor mixtures. Here, using quaternary blends, double cascading energy level alignment in bulk heterojunction organic photovoltaic active layers are realized, enabling efficient carrier splitting and transport. Numerous avenues to optimize light absorption, carrier transport, and charge-transfer state energy levels are opened by the chemical constitution of the components. Record-breaking PCEs of 18.07% are achieved where, by electronic structure and morphology optimization, simultaneous improvements of the open-circuit voltage, short-circuit current and fill factor occur. The donor and acceptor chemical structures afford control over electronic structure and charge-transfer state energy levels, enabling manipulation of hole-transfer rates, carrier transport, and non-radiative recombination losses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

String instability mitigation of adaptive cruise control without modifying control laws: trajectory shaper and parameter estimation

Vehicle automation technologies equip vehicles with adaptive cruise control (ACC) systems, which relieve driving fatigue. However, recent studies have shown that the current ACC systems are string-unstable (i.e., exacerbate traffic congestion). To achieve string stability, most existing studies directly modify the control algorithms of ACC systems. Alternatively, this study proposes a trajectory shaper (TS)-based method, which only modifies the trajectory information of the predecessor vehicle, so that the ego vehicle driven by a string-unstable ACC system leverages the modified trajectory information to achieve string stability. To devise the TS-based method, an offline-online parameter estimation method integrating batch optimization and an extended Kalman filter is applied to estimate the parameters of an ACC system. The proposed TS-based method is cost-effective during implementation, as it avoids modifying existing ACC control algorithms (which entails a complex analysis of control systems and parameter tuning). In conclusion, the effectiveness of the proposed TS-based method is validated through extensive numerical experiments.

33 ADVANCED PROPULSION SYSTEMS↗

Distributed Aerodynamic Control using Active Trailing-Edge Flaps for Large Wind Turbines

This work presents a numerical framework to investigate distributed aerodynamic control devices for application in large wind turbines. Tool capabilities were extended to facilitate multiple aerodynamic polar tables. The airfoil aerodynamics characteristics were automatically determined, and blade-pitch, generator-torque, and trailing-edge-flap controllers were tuned in-the-loop according to a specific blade design. This automated workflow allows analysis and optimization of trailing-edge flaps, enabling codesign studies. Results targeted reductions of root-flap-bending moment derivatives. The applied trailing-edge-flap control reduced the standard deviation of root-flap-bending moments by more than 6% and benefit related parameters, e.g., reduce blade-tip deflections, by up to 8%. Because of varying thrust distributions along the blade span, different flap designs have nonlinear characteristics in terms of the control objective and show best performance when located at the radial position with maximum thrust. In general, larger flaps provide a greater influence to reduce the target control objective.

17 WIND ENERGY↗

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

Effect of epitope variant co-delivery on the depth of CD8 T cell responses induced by HIV-1 conserved mosaic vaccines

To stop the HIV-1 pandemic, vaccines must induce responses capable of controlling vast HIV-1 variants circulating in the population as well as those evolved in each individual following transmission. Numerous strategies have been proposed, of which the most promising include focusing responses on the vulnerable sites of HIV-1 displaying the least entropy among global isolates and using algorithms that maximize vaccine match to circulating HIV-1 variants by vaccine cocktails of optimized complementing sequences. In this study, we investigated CD8 T cell responses induced by a bi-valent mosaic of highly conserved HIVconsvX regions delivered by a combination of simian adenovirus ChAdOx1 and poxvirus MVA. We compared partially and fully mono- and bi-valent prime-boost regimens and their ability to elicit T cells recognizing natural epitope variants using an interferon-g enzyme-linked immunospot (ELISPOT) assay. We used 11 well-defined CD8 T cell epitopes in two mouse haplotypes and, for each epitope, assessed recognition of the two vaccine forms together with the other most frequent epitope variants in the HIV-1 database. We conclude that for the magnitude and depth of epitope recognition, CD8 T cell responses benefitted in most comparisons from the combined bi-valent mosaic and envisage the main advantage of the bi-valent vaccine during its deployment to diverse populations.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-Agent Safe Policy Learning for Power Management of Networked Microgrids

This paper presents a supervised multi-agent safe policy learning (SMAS-PL) method for optimal power management of networked microgrids (MGs) in distribution systems. While unconstrained reinforcement learning (RL) algorithms are black-box decision models that could fail to satisfy grid operational constraints, our proposed method considers AC power flow equations and other operational limits. Accordingly, the training process employs the gradient information of operational constraints to ensure that the optimal control policy functions generate safe and feasible decisions. Furthermore, we have developed a distributed consensus-based optimization approach to train the agents’ policy functions while maintaining MGs’ privacy and data ownership boundaries. After training, the learned optimal policy functions can be safely used by the MGs to dispatch their local resources, without the need to solve a complex optimization problem from scratch. Lastly, numerical experiments have been devised to verify the performance of the proposed method.

, policy gradient↗

Modular Hydronic Room Conditioning System

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

CUDO: closed-form universal dwell-time optimization for computer-controlled optical surfacing

Precision optical figuring demands fast and accurate dwell time optimization to reach nanometer- and sub-nanometer-level accuracy in next-generation optical systems. We introduce CUDO (closed-form universal dwell-time optimization), the first, to the best of our knowledge, unified closed-form analytical framework that supports both function-form and matrix-form dwell time models in computer-controlled optical surfacing (CCOS). In contrast to traditional methods, which rely on iterative optimization and hyperparameter tuning, our framework derives direct analytical solutions with no adjustable parameters. This approach unifies the solution principles of existing methods within a single mathematical model, delivering three key advantages: (1) accuracy on par with, or superior to, iterative solvers, (2) substantial reduction in computation time, and (3) numerical robustness. Comparative studies with prior art confirm that closed-form solutions achieve equivalent residual error while removing runtime bottlenecks. By simplifying the implementation and enabling real-time, scalable deployment, CUDO establishes a practical foundation for future deterministic fabrication of large-aperture and high-performance optics.

36 MATERIALS SCIENCE↗

Travels, optimizations, and ordering

With all the necessary build pathing in place from the preceding steps, the final task is to organize the build pathing via directionality, ordering, and connectivity. Directionality controls closed-loop pathing by deciding whether construction occurs clockwise or counterclockwise. Ordering determines whether construction occurs outside-in or inside-out, by manipulating what path types will be printed first. Connectivity links all the pathing together, which is achieved through the insertion of travels. Connectivity has numerous approaches, but only the most common will be addressed in this chapter. In addition to these manipulations, other considerations, such as spiral pathing and unique cases for infill, will be discussed.

Borish, Michael↗

Parametric and sensitivity analysis of a PCM-integrated wall for optimal thermal load modulation in lightweight buildings

Load modulation in buildings is becoming increasingly important due to growing disparity in energy demand during peak and off-peak hours. Integrating phase change material (PCM) in building envelopes and using a controlled precooling strategy can provide substantial thermal load modulation; however, it may greatly increase the total energy use. Previous studies have employed PCM in building envelopes primarily for energy savings and, to some extent, peak load shedding and shifting. However, the load modulation capacity of a PCM-integrated envelope has not been well explored in the literature. In this study, we perform an extensive parametric and sensitivity analysis on PCM-integrated lightweight building walls and examine the combinatory effects of various PCM parameters on thermal load modulation and wall-related heat gains in buildings. Using numerical simulations, we investigate eight PCM parameters: PCM location in the wall, transition temperature, thickness, latent heat, transition range, density, specific heat, and thermal conductivity. Here, we evaluate their impact and relative importance to achieve maximum load modulation in buildings without compromising occupants’ thermal comfort or total energy use. The results show that the optimized PCM proposed in this study can completely invert the transient heat gain profile of the wall, providing up to 70% reduction of wall-related heat gain during peak hours without a major increase in the cumulative heat gain.

42 ENGINEERING↗

Investigations into Balancing Peak-to-Average Power Ratio and Mean Power Extraction for a Two-Body Point-Absorber Wave Energy Converter

The power harnessed by wave energy converters (WECs) in oceans is highly variable and, therefore, has a high peak-to-average power (PTAP) ratio. To minimize the cost of a WEC power take off (PTO) system, it is desirable to reduce the PTAP ratio while maximizing the mean power extracted by WECs. The important issue of how PTAP ratio reduction measures (such as adding an inertia element) can affect the mean power extracted in a reference model has not been thoroughly addressed in the literature. To investigate this correlation, this study focuses on the integration of the U.S. Department of Energy’s Reference Model 3, a two-body point absorber, with a slider-crank WEC for linear-to-rotational conversion. In the first phase of this study, a full-scale numerical model was developed that predicts how PTO system parameters, along with an advanced control algorithm, can potentially affect the proposed WEC’s PTAP ratio as well as the mean power extracted. In the second phase, an appropriate scaled-down model was developed, and extracted power results were successfully validated against the full-scale model. Finally, numerical and hardware-in-the-loop (HIL) simulations based on the scaled-down model were designed and conducted to optimize or make trade-offs between the operational performance and PTAP ratio. The initial results with numerical and HIL simulations reveal that gear ratio, crank radius, and generator parameters substantially impact the PTAP ratio and mean power extracted.

16 TIDAL AND WAVE POWER↗

Efficient, direct compilation of SU(N) operations into SNAP & Displacement gates

We present a function which connects the parameter of a previously published short sequence of selective number-dependent arbitrary phase (SNAP) and displacement gates acting on a qudit encoded into the Fock states of a superconducting cavity, $V_k(\alpha)=D(\alpha)R_\pi(k)D(-2\alpha)R_\pi(k)D(\alpha)$ to the angle of the Givens rotation $G(\theta)$ on levels $|k\rangle,|k+1\rangle$ that sequence approximates, namely $\alpha=\Phi(\theta) = \frac{\theta}{4\sqrt{k+1}}$. Previous publications left the determination of an appropriate $\alpha$ to numerical optimization at compile time. The map $\Phi$ gives us the ability to compile directly any $d$-dimensional unitary into a sequence of SNAP and displacement gates in $O(d^3)$ complex floating point operations with low constant prefactor, avoiding the need for numerical optimization. Numerical studies demonstrate that the infidelity of the generated gate sequence $V_k$ per Givens rotation $G$ scales as approximately $O(\theta^6)$. We find numerically that the error on compiled circuits can be made arbitrarily small by breaking each rotation into $m$$\theta/m$ rotations, with the full $d\times d$ unitary infidelity scaling as approximately $O(m^{-4})$. This represents a significant reduction in the computational effort to compile qudit unitaries either to SNAP and displacement gates or to generate them via direct low-level pulse optimization via optimal control.

Job, Joshua↗

Evaluation of the Energy, Hygrothermal, and Thermal Capacity Performance of Cross-Laminated Timber

Cross-laminated timber (CLT) construction is gaining momentum in the US because it offers multiple advantages over traditional construction methods. Benefits that have received the most attention focus on constructability, the environment, and protection (e.g., blast resistance), although CLT construction is likely to offer other benefits, as well. Still, these have not been studied at length because such evaluations are costly, requiring long-term assessments in an actual building and specialized technical knowledge. Among the possible benefits, CLT construction likely provides a higher-performing building envelope. Using CLT panels to enclose a building means fewer joints in the opaque envelope than what is required in traditional stick-framed construction. Fewer joints mean fewer locations where the air- and water-resistive barrier (WRB) could be compromised; thus, a CLT building enclosure may require less maintenance and have a longer lifespan than a traditionally built structure because of fewer air and water leaks. In addition, CLT’s thermal mass moderates indoor temperatures, allowing the heating, ventilation, and air conditioning (HVAC) system to operate more efficiently during peak hours, reducing operational energy consumption throughout the lifetime of the CLT building (Salonvaara et al., 2022). Furthermore, more stable indoor temperatures can increase occupant comfort. The CLT’s thermal mass can also reduce energy costs by adjusting to utility time-of-use pricing without affecting occupant comfort. The ability of CLT buildings to bridge periods without HVAC operation prepares them for future grid interaction and provides a certain level of resilience against power outages. Researchers have attempted to quantify these benefits; however, their work is based on simplified simulations with numerous assumptions. To correctly understand the benefits, an actual building must be monitored. Therefore, information needs to be gathered on indoor and outdoor temperatures, HVAC energy consumption, thermostat setpoints, temperatures, and thermal transport in CLT components to comprehend how these parameters are affected by the CLT’s thermal mass. These data are needed to reduce the number of assumptions and calibrate simulation models to optimize HVAC controls to minimize overall energy consumption, reduce energy use and higher fees during peak demand, and maintain occupant comfort. Additionally, the calibrated simulation model allows the optimization exercise to be repeated in various US climates. Potential benefits can be tailored to buildings in various locations, and decisions can be made on where CLT construction could be most advantageous. Furthermore, monitoring and simulation results are needed to evaluate the durability of the CLT structures in different climates. This project’s researchers gathered information to help understand and quantify the benefits of CLT buildings concerning operational energy, moderated indoor temperatures, and comfort; the dynamic operation to provide grid services; and resilience in times of power outage. Through the corroboration of simulation models with real-world measurements, this study paves the way for extrapolating findings to other climatic zones and building typologies, thereby broadening the understanding of CLT’s multifaceted benefits and reinforcing its position as a material of choice in sustainable construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Examination of a Methane/Diesel RCCI Engine Using Pele: Preprint

Multi-fuel, advanced injection strategies have become increasingly promising as a strategy to mitigate the emissions generated from internal combustion engines. By carefully controlling the combustion phasing in-cylinder, these new multi-pulse, multi-fuel injection strategies are able to burn in the low-temperature combustion regime where both NOx and soot are not readily produced, reducing the need for extensive exhaust gas recirculation systems. In this study, we examine a reactivity-controlled compression ignition (RCCI) strategy that uses an early pre-filled methane-air mixture with low turbulence background as the low-reactivity fuel and a direct injection of four discrete dodecane jets as a surrogate for the high-reactivity diesel fuel. We use the Pele software suite, a highly optimized, exascale-ready, adaptive mesh refinement codebase to perform high-resolution numerical simulations of a scaled down, single cylinder from the RCCI engine. Here, we resolve the ignition kernels down to micrometer scales and present several statistical quantities evaluating the development of the flow and detailing the onset of ignition and subsequent flame development. Particular attention is paid to the conditions surrounding the onset of the first ignition kernels and discussing what led to the development of those conditions.

CFD↗

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Role of the Digital Thread for Security, Resilience, and Adaptability in Manufacturing

The manufacturing sector is dramatically evolving with recent technical and digital advances. Among these, the digital thread is revolutionizing manufacturing operations well beyond the historic downloading of programs to computer numerically controlled (CNC) machine tools and the uploading of edited programs from CNC controllers. The digital thread is “the communication framework that enables a connected data flow and integrated view of the asset’s data throughout its lifecycle across traditionally siloed functional perspectives” (Leiva 2016). It is pervasive around the world and has changed the way society operates; for example, map apps are used to provide time-optimal directions to destinations using real-time traffic feedback. The digital thread makes it possible to digitally verify products, ensure that the latest technologies are deployed across the entire manufacturing ecosystem, and strengthen the workforce by making each individual more efficient and effective. These three foundational advanced manufacturing concepts will ensure a next-generation secure, resilient, and adaptable manufacturing ecosystem and sustainably support current and future needs of society (Lynn et al. 2020).

99 GENERAL AND MISCELLANEOUS↗