Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained 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 685 records · Page 38

Optimizing infrasound observations for sample return capsule re-entry: Insights from OSIRIS-REx and Hayabusa2

The atmospheric entry of meteoroids presents a rare and unpredictable phenomenon, posing challenges for systematic observation and detailed characterization. Such events are nonetheless critical for advancing understanding of acoustic wave propagation, atmospheric structure, and entry dynamics. In contrast, sample return capsules (SRCs) from space missions follow well constrained re-entry trajectories, enabling planned observations of shock wave generation and propagation under controlled conditions. This study compares two SRC atmospheric entries, Hayabusa2 in 2020 and OSIRIS-REx in 2023, to assess how different infrasound array configurations influence shock wave detection and trajectory validation. Hayabusa2’s re-entry was monitored using a distributed network of 28 portable infrasound sensors across seven arrays in Woomera, Australia, permitting three-dimensional reconstruction of the trajectory and analysis of wave characteristics. For OSIRIS-REx, a compact four-sensor array deployed near Eureka Airport provided trajectory confirmation through arrival-time differences and back-azimuth estimates. Spectral and waveform analyses revealed differences in signal properties associated with variations in entry angle and velocity. The results illustrate both the strengths and the limitations of the deployed array configurations. The dense and distributed arrays during Hayabusa2’s re-entry enabled detailed trajectory reconstruction, whereas the compact array at Eureka primarily provided confirmation of signal coherence and back-azimuth consistency. These case studies highlight that even relatively small arrays, if located close to the predicted ground track, can still capture useful information on arrival direction and timing. Such insights provide practical guidance for planning future observational campaigns of SRC returns and other controlled atmospheric entries, and may also inform approaches to opportunistic observations of meteoroid events and other atmospheric acoustic phenomena.

celestial mechanics: orbit determination↗

The future of the NPT: REVCON and beyond

The Treaty on the Nonproliferation of Nuclear Weapons (NPT), with its quinquennial review conferences (RevCons), has become the foundation of the international nonproliferation regime over the course of its 50-year history. With more states parties than any other arms control agreement (currently 191), the Treaty provides a legally binding framework for all states parties, constraining proliferation of nuclear weapons and committing the five recognized states who possess nuclear weapons to work toward disarmament. The NPT has enjoyed wide success in constraining the spread of nuclear weapons. Only one state party has left the agreement and gone on to acquire these weapons, and the overall track record of compliance with IAEA safeguards has been robust. However, for all its effectiveness in the area of nonproliferation, its achievements on the disarmament front have been relatively unimpressive. While the NWS – most notably the United States and the Soviet Union/Russia – have decreased the overall size of their nuclear arsenals significantly since 1970, these arsenals have remained an important feature of NWS security policies. The limited progress towards disarmament has created friction between the NNWS and NWS and has been the source of mounting tension within the NPT regime. The 2020 RevCon, the 50th anniversary of the NPT and 25th anniversary of its indefinite extension, faces a range of contentious challenges, old and new, which call into question the Treaty’s future viability.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An advanced generation land mobile satellite system and its critical technologies

A conceptual design for a Land Mobile Satellite System (LMSS) for the 1990s is presented. LMSS involves small tranceivers accessing satellites directly, with ground reception through small car-top antennas. The satellite would have a large antenna and blanket coverage areas in the UHF. The call may originate from a home, be carried by wire to a gateway, transmitted to satellite on the S-band, converted to UHF on the satellite, and transmitted to the vehicle. The system design is constrained by the number of users in an area during the busiest hours, Shuttle storage, controllability factors, and the total area served. A 55-m antenna has been selected, with 87 spot beams and two 10 MHz UHF bands in the 806-890 MHz band. A 17 dB interbeam isolation level is required, implying that sufficient sub-bands can be generated to assure 8265 total channels. The mobile satellite (MSAT) would have an 83 m mast lower segment, a 34 m upper segment, and a second, 10 m antenna made of a deployable mesh. Various antenna function modes are considered.

Naderi, F.↗

Frequency domain identification experiment on a large flexible structure

Recent experiences in the field of flexible structure control in space have indicated a need for on-orbit system identification to support robust control redesign to avoid in-flight instabilities and maintain high spacecraft performance. The authors highlight an automated frequency domain system identification methodology recently developed to fill this need. The methodology supports (1) the estimation of system quantities useful for robust control analysis and design, (2) experiment design tailored to performing system identification in a typically constrained on-orbit environment, and (3) the automation of operations to reduce human-in-the-loop requirements. A basic overview of the methodology is presented first, followed by an experimental verification of the approach performed on the JPL/AFAL testbed facility.

Bayard, D. S.↗

Automated on-orbit frequency domain identification for large space structures

Recent experiences in the field of flexible structure control in space have indicated a need for on-orbit system identification to support robust control redesign to avoid in-flight instabilities and maintain high spacecraft performance. This paper highlights an automated frequency domain system identification methodology recently developed to fulfill this need. The methodology is focused to support (1) the estimation of system quantities useful for robust control analysis and design; (2) experiment design tailored to performing system identification in a typically constrained on-orbit environment; and (3) the automation of operations to reduce 'human in the loop' requirements.

Bayard, D. S.↗

Development of homotopy algorithms for fixed-order mixed H2/H(infinity) controller synthesis

A major difficulty associated with H-infinity and mu-synthesis methods is the order of the resulting compensator. Whereas model and/or controller reduction techniques are sometimes applied, performance and robustness properties are not preserved. By directly constraining compensator order during the optimization process, these properties are better preserved, albeit at the expense of computational complexity. This paper presents a novel homotopy algorithm to synthesize fixed-order mixed H2/H-infinity compensators. Numerical results are presented for a four-disk flexible structure to evaluate the efficiency of the algorithm.

Whorton, M.↗

Europa Clipper Mission: System Integration Review Report

Jupiter’s icy moon Europa is a prime target in our exploration of potentially habitable worlds beyond Earth. The combination of a subsurface liquid water layer in contact with a rocky seafloor may yield an ocean rich in the elements and energy needed for the emergence of life, and for potentially sustaining life through time. Europa may hold the clues to one of NASA’s long-standing quests – to determine whether or not we are alone in the universe. The Europa Clipper mission will characterize Europa’s habitability as the first step in the search for potential life at Europa by conducting approximately four dozen flybys. The project is preparing for a System Integration Review in November of 2021. This paper will summarize changes from the mission’s Preliminary Design Review (PDR) baseline, including science objectives, launch vehicle, interplanetary trajectory, science tour, flight system, payload and mission phases. Since PDR several of the Level 1 science requirements were descoped and the science requirement flow was realigned, in order to control cost growth. The realignment deemphasized science redundancy, which had resulted in an over-constrained science tour, in favor of a more targeted linkage to the primary responding instrument. Additionally, changes to the PDR flight system design baseline include changes to the sounding radar antenna design, the magnetometer architecture, and the heat redistribution system pumps. All changes have been incorporated into a reference science tour used to demonstrate that the requirements and design of the flight system and payload are mature enough to support flight builds.

Cook, Kendra↗

Switched low-dropout voltage regulator

High-resolution switched digital regulators are disclosed having fast cross corner and variable temperature response, with constrained ripple. The strength of the power transistors utilized by the regulator are adjusted to control the current delivered to the load. The regulators utilize a slow control loop in parallel with a primary fast switching loop. The slow loop uses the switching signal of the primary loop to estimate the load current and set the power transistor size accordingly.

Kudva, Sudhir Shrikantha↗

Organic sulfur fluxes and geomorphic control of sulfur isotope ratios in rivers

Pyrite oxidation plays a critical role in the relationship between weathering and climate, and its impact on the global carbon cycle has previously been constrained through inversion models utilizing observations of river sulfate ($SO^{2–}_{4}$) and its 34 S/ 32 S isotope ratio (δ 34 S SO4 ). However, measurements from some rivers have suggested that SSO4 can be substantially impacted by processes such as microbial sulfate reduction and/or sulfur assimilation and cycling, rather than simply reflecting a weighted mixture of lithologic sulfur sources. To study the prevalence and controls on $SO^{2–}_{4}$ transformations, in this study we measured dissolved major element concentrations and δ 34 S SO4 in river water samples from throughout western Iceland. Our analyses focused on samples from a small catchment hosting the Efri Haukadalsá river, a system with relatively uniform and isotopically constrained basaltic bedrock. We also measured sediment δ 34 S and sulfur speciation using sulfur K-edge X-ray absorption spectroscopy on sediment and vegetation samples from this catchment. Values of dissolved δ 34 S SO4 in the Efri Haukadalsá ranged from 2.5‰ to 23.7‰ and had a linear relationship with Cl – /$SO^{2–}_{4}$ ratios, indicating that $SO^{2–}_{4}$ predominantly derived from basalt weathering and meteoric precipitation. The lower δ 34 S SO4 values were found in fluvial valleys with V-shaped cross sections, while higher values of δ 34 S SO4 occurred in U-shaped, glacially eroded valleys with thick alluvial fills blanketing the valley floor. Spectroscopic observations identified organic sulfur phases in suspended river sediment, floodplain deposits, and vegetation. Mass balance calculations quantified the organic sulfur flux as less than 10% of $SO^{2–}_{4}$ export, and sediment δ 34 S values were comparable to river δ 34 S SO4 . We interpreted these isotopic and chemical patterns as reflecting differences in the availability of unweathered bedrock across the Efri Haukadalsá catchment, with V-shaped valleys having greater access to fresh sulfide-bearing minerals than alluviated U-shaped valleys; this interpretation is in contrast to one in which the elevated δ 34 S SO4 values reflect fractionation during sulfur transformations along alluvial reaches. These results validated the application of river inversion models for constraining weathering fluxes and affirmed that pyrite oxidation globally, even in the presence of river sulfur cycling, modulates the abundance of atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗

Fast model-based scenario optimization in NSTX-U enabled by analytic gradient computation

Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

Interconnection-wide hour-ahead scheduling in the presence of intermittent renewables and demand response: A surplus maximizing approach

This study describes a new approach for solving the multi-area electricity resource allocation problem when considering both intermittent renewables and demand response. The method determines the hourly inter-area export/import set that maximizes the interconnection (global) surplus satisfying transmission, generation and load constraints. The optimal inter-area transfer set effectively makes the electricity price uniform over the interconnection apart from constrained areas, which overall increases the consumer surplus more than it decreases the producer surplus. The method is computationally efficient and suitable for use in simulations that depend on optimal scheduling models. The method is demonstrated on a system that represents North America Western Interconnection for the planning year of 2024. Simulation results indicate that effective use of interties reduces the system operation cost substantially. Excluding demand response, both the unconstrained and the constrained scheduling solutions decrease the global production cost (and equivalently increase the global economic surplus) by 12.30B and 10.67B per year, respectively, when compared to the standalone case in which each control area relies only on its local supply resources. This cost saving is equal to 25% and 22% of the annual production cost. Including 5% demand response, the constrained solution decreases the annual production cost by 10.70B, while increases the annual surplus by 9.32B in comparison to the standalone case.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Stochastic Model Predictive Control With Gaussian Wind Direction Preview for Wake Steering

This article addresses the problem of wake steering control for wind farms that explicitly consider the tradeoff between farm-level power generation and yaw duty cycle under variable and uncertain wind conditions. A novel stochastic model predictive control (MPC) algorithm is presented, which utilizes a stochastic model of the freestream wind field components in a receding horizon framework to compute optimal yaw set points that maximize the expected value of the farm power while constraining the yaw actuation. Different configurations of the algorithm are evaluated using a steady-state wind farm simulator. The proposed stochastic MPC algorithm can plan control actions over a future prediction horizon based on probabilistic estimates of the incoming wind magnitude and direction.

17 WIND ENERGY↗

SMA Hybrid Composites for Dynamic Response Abatement Applications

A recently developed constitutive model and a finite element formulation for predicting the thermomechanical response of Shape Memory Alloy (SMA) hybrid composite (SMAHC) structures is briefly described. Attention is focused on constrained recovery behavior in this study, but the constitutive formulation is also capable of modeling restrained or free recovery. Numerical results are shown for glass/epoxy panel specimens with embedded Nitinol actuators subjected to thermal and acoustic loads. Control of thermal buckling, random response, sonic fatigue, and transmission loss are demonstrated and compared to conventional approaches including addition of conventional composite layers and a constrained layer damping treatment. Embedded SMA actuators are shown to be significantly more effective in dynamic response abatement applications than the conventional approaches and are attractive for combination with other passive and/or active approaches.

Turner, Travis L.↗

Advancements in 3D Structural Analysis of Geothermal Systems

Robust geothermal activity in the Great Basin, USA is a product of both anomalously high regional heat flow and active fault-controlled extension. Elevated permeability associated with some fault systems provides pathways for circulation of geothermal fluids. Constraining the local-scale 3D geometry of these structures and their roles as fluid flow conduits is crucial in order to mitigate both the costs and risks of geothermal exploration and to identify blind (no surface expression) geothermal resources. Ongoing studies have indicated that much of the robust geothermal activity in the Great Basin is associated with high density faulting at structurally complex fault intersection/interaction areas, such as accommodation/transfer zones between discrete fault systems, step-overs or relay ramps in fault systems, intersection zones between faults with different strikes or different senses of slip, and horse-tailing fault terminations. These conceptualized models are crucial for locating and characterizing geothermal systems in a regional context. At the local scale, however, pinpointing drilling targets and characterizing resource potential within known or probable geothermal areas requires precise 3D characterization of the system. Employing a variety of surface and subsurface data sets, we have conducted detailed 3D geologic analyses of two Great Basin geothermal systems. Using EarthVision (Dynamic Graphics Inc., Alameda, CA) we constructed 3D geologic models of both the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These 3D models allow spatial comparison of disparate data sets in 3D and are the basis for quantitative structural analyses that can aid geothermal resource assessment and be used to pinpoint discrete drilling targets. The relatively abundant data set at Brady’s, ~80 km NE of Reno, NV, includes 24 wells with lithologies interpreted from careful analysis of cuttings and core, a 1:24,000 scale detailed geologic map and cross-sections, 2D seismic reflection profiles and other geophysical data, and downhole temperature data. The 3D geologic model based on these data consists of 61 fault planes, 25 distinct stratigraphic units, and 2 intrusive bodies. Geothermal fluids are produced from a left step-over/relay ramp within the Brady’s Fault Zone (BFZ). Under local stress conditions, fault segments that strike NNE-to-NE are most likely to slip and/or dilate, and therefore transmit geothermal fluids. The 3D model defines the locations of discrete fault intersections within the BFZ and indicates that the densest zones of structurally controlled fracture permeability are ~10-to-10s of meters in diameter and plunge ~55° NW-NNW beneath the heart of the BFZ step over. The locations of high intersection density, high fault slip and dilation tendency, high subsurface temperature, and lithologies known to support high fracture permeability are combined to produce 3D ‘fairway’ maps useful in both assessments of geothermal resource potential and for defining drilling targets. Astor Pass is located on the Pyramid Lake Paiute Reservation, ~80 km north of Reno, NV. It is a prospective ‘greenfield’ geothermal area, and thus subsurface data are relatively sparse. Available data include: two relatively deep wells (~1400 m) and one shallower well (~500 m) with lithologies interpreted from drill cuttings, several 2D seismic reflection profiles, a 1:24,000 scale geologic map and cross-section, a shallow temperature survey, and downhole temperature data. 3D modeling based on these data has defined 19 distinct fault planes and 16 stratigraphic units. Based on the stress field calculated from borehole breakouts, drilling induced tensile cracks and petal-centerline cracks in the two relatively deep wells, 3D slip and dilation tendency analysis indicates that northerly striking fault segments are most likely to slip and/or dilate, and therefore transmit geothermal fluids. Analysis of fault intersection density indicates that the highest density of structurally controlled permeability within the field lies in a narrow (10-to-10s of m) zone plunging moderately (~35°) to the NNW beneath Pleistocene tufa deposits. This zone of increased fracture density, which we interpret as the primary upflow zone, is controlled by the intersection of N-to-NNW striking normal faults and a WNW striking dextral fault zone and represents the most promising target for future drilling. Construction of a 3D geologic model involves integration of a variety of data into an internally consistent framework. A robust model allows for spatial comparison between the various types of data (structural, stratigraphic, geophysical, temperature, etc.) that are commonly used independently to site geothermal wells. Furthermore, highly detailed 3D geologic models provide the basis for additional quantitative analysis, including 3D fault slip and dilation tendency analysis and the precise location of structurally controlled permeability pathways. These analyses provide detailed information relating to the internal dynamics of geothermal systems and can mitigate the costs and risks of geothermal exploration and development by contributing to better well targeting and more accurate evaluations of resource potential.

15 GEOTHERMAL ENERGY↗

Constraining Remote Oxidation Capacity with ATom Observations

The global oxidation capacity, defined as the tropospheric mean concentration of the hydroxyl radical (OH), controls the lifetime of reactive trace gases in the atmosphere such as methane and carbon monoxide (CO). Models tend to underestimate the methane lifetime and CO concentrations throughout the troposphere, which is consistent with excessive OH. Approximately half the oxidation of methane and non-methane volatile organic compounds (VOCs) is thought to occur over the oceans where oxidant chemistry has received little validation due to a lack of observational constraints. We use observations from the first two deployments of the NASA ATom aircraft campaign during July–August 2016 and January–February 2017 to evaluate the oxidation capacity over the remote oceans and its representation in the GEOS-Chem chemical transport model. The model successfully simulates the magnitude and vertical profile of remote OH within the measurement uncertainties. Comparisons against the drivers of OH production (water vapor, ozone, and NOy concentrations, ozone photolysis frequencies) also show minimal bias with the exception of wintertime NOy, for which a model overestimate may indicate insufficient wet scavenging and/or missing loss on seasalt aerosol but large uncertainties remain that require further studies of NOy partitioning and removal in the troposphere. During the ATom-1 deployment, OH reactivity (OHR) below 3 km is significantly enhanced, and this is not captured by the sum of its measured components (cOHRobs) or by the model (cOHRmod). This enhancement could suggest missing reactive VOCs but cannot be explained by new estimates of ocean VOC sources and additional modeled reactivity in this region would be difficult to reconcile with the full suite of ATom measurement constraints. The model generally reproduces the magnitude and seasonality of cOHRobs but underestimates the contribution of oxygenated VOC, mainly acetaldehyde, which is severely underestimated throughout the troposphere despite its calculated lifetime of less than a day. Missing model acetaldehyde in previous studies was attributed to measurement uncertainties that have been largely resolved. Observations of peroxyacetic acid (PAA) provide new support for remote levels of acetaldehyde. The underestimate in modeled acetaldehyde and PAA is present throughout the year in both hemispheres and peaks during Northern Hemisphere summer. The addition of ocean VOC sources in the model increases annual surface cOHRmod by 10 % and improves model-measurement agreement for acetaldehyde particularly in winter but cannot resolve the model summertime bias. Doing so would require a 100 Tg yr−1 source of a long-lived unknown precursor throughout the year with significant additional emissions in the Northern Hemisphere summer. Improving the model bias for remote acetaldehyde and PAA is unlikely to fully resolve previously reported model global biases in OH and methane lifetime, suggesting that future work should examine the sources and sinks of OH over land.

Travis, Katherine R.↗

Characteristics of constrained optimum trajectories with specified range

Necessary conditions of optimality are derived for trajectories whose structure is limited to climb, steady cruise, and descent segments. The performance function consists of the sum of fuel and time costs, referred to as direct operating cost (DOC). The state variable is range to go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (500 n. mi. and less) a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur. Computer calculations of typical short-range trajectories obtained about a 1% cost penalty for constraining the thrust.

Erzberger, H.↗