Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 649 records · Page 36

A bilevel multistage stochastic self-scheduling model with indivisibilities for trading in the continuous intraday electricity market

In this paper, we study the profit maximization problem of a virtual power plant trading in the continuous intraday electricity market. Our virtual power plant model is compatible with renewable, and thermal assets, covering a range of virtual power plants currently participating in energy markets. We model the trading problem as a bilevel multistage stochastic program. The upper level of the problem accounts for the profit maximization of the virtual power plant with explicit modeling of the technical constraints of the operational status of the thermal power plant including minimum start-up and shut-down times, ramp-up and ramp-down rates, and minimum generation level. The upper level also decides which continuous and indivisible (fill-or-kill) orders are submitted to the market. The lower-level problem accounts for the clearing of the continuous intraday market, i.e., matching of buy and sell orders. Because of the presence of fill-or-kill orders, the lower-level problem is mixed-integer, which prevents its direct conversion to a single-level problem using duality. In order to solve this challenging problem, we develop a convex-hull extended formulation for the lower-level problem, apply duality theory to obtain a single-level stochastic equivalent formulation, and employ McCormick envelopes to turn the problem into a multistage stochastic mixed-integer linear problem, which we solve using the stochastic dual dynamic integer programming algorithm. We conduct numerical experiments and analyze the optimal trading behavior of a virtual power plant trading in an ideal continuous market without arbitrage.

Bilevel multistage stochastic programming problem↗

Zigzag flow reactor for weekly thermochemical energy storage

This paper describes theoretical models and experimental performance of a novel Zigzag Flow Reactor (ZFR) for weekly thermochemical energy storage. The ZFR reduces redox-active metal oxide (MO x ) particles at high temperature (up to ~1100 °C) under inert gas sweep. A physical model demonstrates the approach to process equilibrium by minimizing the associated exergy destruction in a finite number of reaction steps, establishing the thermodynamic requirements for a practical reactor. The model results show several cost-relevant parameter tradeoffs, and the tradeoff analysis implies a cost-optimized set of boundary conditions. Numerical models and prototypes show that the ZFR enables significant gas phase homogenization while simultaneously enabling a customizable MO x residence time in the reactor, both key requirements for approaching an equilibrium process. A scaling model demonstrates the simplicity and affordability of sizing the ZFR to grid-scale levels, with fabrication costs at least five times lower than previously proposed scalable reactor concepts. As a result, a laboratory ZFR prototype achieved an energy storage density of ~90 Wh/kg with CaAl 0.2 Mn 0.8 O 3-δ as the MO x , at temperatures of ~850 °C in >10 h of total runtime.

Thermochemical energy storage↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New technique for optimal smoothing of data

Recursive method for the optimal smoothing of data has numerical superiority and is more easily understood in terms of physical reasoning than earlier methods. Using a Kalman filter, the smoothing technique, applied to a nonlinear parameter identification problem, is useful in those situations where linearization about a reference solution is valid.

Fraser, D. C.↗

A linear stochastic model of the human operator

A linear stochastic model of the human operator is developed and applied to the problem of piloted control of an aircraft. The pilot and aircraft are modeled as linear time-invariant systems containing both process and measurement noise. The loop closure by the pilot is determined by formulating the problem as an optimal stochastic control problem. The solution to the optimal control problem yields not only the pilot's optimal control output which he uses to control the vehicle, but also the optimal combination of his observations of the vehicle states upon which the pilot bases his control. A method is presented so that, using experimental pilot vehicle data, the cost functional which is minimized in the optimal control problem will be numerically equal to the pilot rating that the pilot would associate with the given vehicle and task.

Durrett, J. C.↗

Optimum quantisers for a Gaussian input probability density and for the magnitude-error distortion measure

The parameters of non-uniform and uniform quantizers up to ten bits of quantization, optimum for a Gaussian input probability and for the magnitude-error distortion criterion are computed. Optimum quantizers must be understood as quantizers with minimum distortion. The numerical method used for the optimization converges relatively rapidly. The comparison between optimum non-uniform quantizers and optimum uniform quantizers is made.

Demaret, J. C.↗

Propeller Study. Part 2: the Design of Propellers for Minimum Noise

The design of propellers which are efficient and yet produce minimum noise requires accurate determinations of both the flow over the propeller. Topics discussed in relating aerodynamic propeller design and propeller acoustics include the necessary approximations and assumptions involved, the coordinate systems and their transformations, the geometry of the propeller blade, and the problem formulations including the induced velocity, required in the determination of mean lines of blade sections, and the optimization of propeller noise. The numerical formulation for the lifting-line model are given. Some applications and numerical results are included.

Ormsbee, A. I.↗

Mathematical model for preventive maintenance scheduling

A model is formulated to describe the effect of the time interval chosen for preventive maintenance upon the frequency of failure and frequency of total maintenance. Trade-offs between these two frequencies are determined by computation of an optimal interval in the case where the failure distribution is known. For unknown distributions, an adaptive statistical technique is developed that converges to an optimal preventive maintenance interval. A numerical illustration is given.

Lorden, G.↗

Hybrid thermoelectric solar collector design and analysis

A flat-plate solar collector is conceived where energy cascades through thermoelectric power modules generating direct-current electricity. The intent of this work was to choose a collector configuration and to perform a steady-state thermal performance assessment. A set of energy balance equations were written and solved numerically for the purpose of optimizing collector thermal and electrical performance. The collector design involves finned columns of thermoelectric modules imbedded in the absorber plate (hot junction) over a parallel array of vertical tubes. The thermoelectric power output is limited by the small hot-junction/cold-junction temperature difference which can be maintained under steady-state conditions. The electric power per unit tube pass area is found to have a maximum as a function of a geometric parameter, while electric power is maximized with respect to an electric resistance ratio. Although the electric power efficiency is small, results indicate that there is sufficient electric power production to drive a coolant circulator, suggesting the potential for a stand-alone system.

Roberts, A. S., Jr.↗

Rapid torque-limited line-of-sight pointing of SCOLE (Spacecraft Control Laboratory Experiment) configuration

The design concept of a control for rapid torque-limited slewing of a rigid-mast version of the NASA SCOLE configuration is presented and demonstrated by means of numerical simulation. The time-optimal control problem for the system is decomposed into separate single-axis problems, expanding analytically the implicit nonlinear transcendental expression for the SCOLE line-of-sight error, and the final Euler attitude angles and slew angles are determined. The simulation results are presented in tables and graphs, and it is found that bang-bang or bang-pause-bang slew maneuvers with control moment applied to the Shuttle and control force applied to the reflector, and with a 5-deg/s slew-rate limit, produce the best pointing accuracy and the shortest slew times, although the specified line-of-sight error of 0.02 deg cannot be achieved using such open-loop single-axis maneuvers.

Lin, J. G.↗

Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation

Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.

Goorjian, Peter M.↗

A General Approach to the Geostationary Transfer Orbit Mission Recovery

This paper discusses recovery scenarios for geosynchronous satellites injected in a non-nominal orbit due to a launcher underperformance. The theory on minimum-fuel orbital transfers is applied to develop an operational tool capable to design a recovery mission. To obtain promising initial guesses for the recovery three complementary techniques are used: p-optimized impulse function contouring, a numerical impulse function minimization and the solutions to the switching equations. The tool evaluates the feasibility of a recovery with the on-board propellant of the spacecraft and performs the complete mission design. This design takes into account for various mission operational constraints such as e.g., the requirement of multiple finite-duration burns, third-body orbital perturbations, spacecraft attitude constraints and ground station visibility. In a final case study, we analyze the consequences of a premature breakdown of an upper rocket stage engine during injection on a geostationary transfer orbit, as well as the possible recovery solution with the satellite on-board propellant.

Faber, Nicolas↗

Two-Phase Flow in Filters and Orifices: In preparation for Measurements on ISS Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS)

Understanding the hydrodynamics of adiabatic two-phase flows in packed bed reactors offers numerous benefits. It enables the optimization of chemical reactions rates and products, crucial for pharmaceuticals and energy production industries. Additionally, this understanding aids in designing more efficient and compact reactors, reducing space and resource requirements. Insights gained from studying such flows in microgravity contribute to advancements of space technologies and the enhancement of our capabilities for undertaking future long duration safe, and sustainable space exploration missions. The Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS) is a flight experiment planned for operation on the Microgravity Science Glove box (MSG) facility of the International Space Station (ISS). The objective of the experiment is to validate hydrodynamic models at a range of gas and liquid flow rates. The microgravity environment allows for measurement and observation of aspects of fluid dynamics that are unique when compared to observations made in full or partial gravity environments. The experiment consists of testing two-phase flow hydrodynamics in three different filters, four (4) orifices and one check valve test articles. The experiment utilizes the hardware of the previously flown PBRE and PBRE-2 with packed beds of different packings materials and sizes. The fluid system of the PBRE-WRS consists of a nitrogen gas delivery and water delivery subsystems. The gas delivery subsystem can accommodate two ranges of gas flow rates, just as for the water delivery system. The nitrogen gas loop is open whereas the water loop is closed. Gas is separated from the water using a gas-liquid phase separator which is located in the water delivery module. Preliminary results from ground testing show the pressure drop increasing at different rates with the liquid and gas flow rates in the Brine filter. In this work, a detailed system description is presented along with a summary of results from ground performance testing of different test sections in an end-to-end preflight testing campaign.

Packed Bed Reactor Experiment↗

Two-Phase Flow in Filters and Orifices: ISS Packed Bed Reactor Experiment Water Recovery Series (PBRE-WRS)

Understanding the hydrodynamics of adiabatic two-phase flows in packed bed reactors offers numerous benefits. It enables the optimization of chemical reactions rates and products which results in more efficient and compact reactors, thus reducing space and resource requirements, crucial for long duration space missions. The Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS) is a flight experiment planned for operation on the Microgravity Science Glove box (MSG) facility of the International Space Station (ISS). The objective of the experiment is to validate hydrodynamic models at a range of gas and liquid flow rates, where these models are used in the design of processes for life support systems in microgravity. The microgravity environment allows for measurement and observation of aspects of fluid dynamics that are unique when compared to observations made in full or partial gravity environments. The experiment consists of testing two-phase flow hydrodynamics in three different filters, four (4) orifices and one check valve test articles. The experiment utilizes the hardware of the previously flown PBRE and PBRE-2 with packed beds of different packings materials and sizes. The fluid system of the PBRE-WRS consists of a nitrogen gas delivery and water delivery subsystems. The gas delivery subsystem can accommodate two ranges of gas flow rates, just as for the water delivery system. The nitrogen gas loop is open whereas the water loop is closed. Gas is separated from the water using a gas-liquid phase separator which is located in the water delivery module. Preliminary results from ground testing show the pressure drop across the filters increasing at different rates with the liquid and gas flow rates. In this work, a detailed system description is presented along with a summary of results from ground performance testing of different test sections in an end-to-end preflight testing campaign.

ISS↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

A forward method for optimal stochastic nonlinear and adaptive control

A computational approach is taken to solve the optimal nonlinear stochastic control problem. The approach is to systematically solve the stochastic dynamic programming equations forward in time, using a nested stochastic approximation technique. Although computationally intensive, this provides a straightforward numerical solution for this class of problems and provides an alternative to the usual dimensionality problem associated with solving the dynamic programming equations backward in time. It is shown that the cost degrades monotonically as the complexity of the algorithm is reduced. This provides a strategy for suboptimal control with clear performance/computation tradeoffs. A numerical study focusing on a generic optimal stochastic adaptive control example is included to demonstrate the feasibility of the method.

Bayard, David S.↗