Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization methods”

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

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

At least 1,423 records · Page 79

Uncertainty Assessment of Space-Borne Passive Soil Moisture Retrievals

The uncertainty associated with passive soil moisture retrieval is hard to quantify, and known to be underlain by various, diverse, and complex causes. Factors affecting space-borne retrieved soil moisture estimation include: (i) the optimization or inversion method applied to the radiative transfer model (RTM), such as e.g. the Single Channel Algorithm (SCA), or the Land Parameter Retrieval Model (LPRM), (ii) the selection of the observed brightness temperatures (Tbs), e.g. polarization and incidence angle, (iii) the definition of the cost function and the impact of prior information in it, and (iv) the RTM parameterization (e.g. parameterizations officially used by the SMOS L2 and SMAP L2 retrieval products, ECMWF-based SMOS assimilation product, SMAP L4 assimilation product, and perturbations from those configurations). This study aims at disentangling the relative importance of the above-mentioned sources of uncertainty, by carrying out soil moisture retrieval experiments, using SMOS Tb observations in different settings, of which some are mentioned above. The ensemble uncertainties are evaluated at 11 reference CalVal sites, over a time period of more than 5 years. These experimental retrievals were inter-compared, and further confronted with in situ soil moisture measurements and operational SMOS L2 retrievals, using commonly used skill metrics to quantify the temporal uncertainty in the retrievals.

Quets, Jan↗

Preliminary Design of the Superconducting Rotor for NASA's High-Efficiency Megawatt Motor

NASA Glenn Research Center is developing a 1.4~MW high-efficiency electric machine for future electrified aircraft to reduce energy consumption, emissions, and noise. This wound-field, synchronous machine employs a self-cooled, superconducting rotor to achieve excellent specific power and efficiency. The design of the superconducting rotor and the optimization of its electromagnetic and structural responses are analyzed in this paper. Candidate designs are evaluated in terms of absolute performance, specific performance, and performance per cost. It is found that optimizing the electromagnetic response for absolute or specific performance yields designs that are similar to each other, but opposite of the cost-optimized design. A method to define the thermal requirements of the superconducting coils and integrated cryocooler is also presented.

electric machines↗

Preliminary Design of the Superconducting Rotor for NASA's High-Efficiency Megawatt Motor

NASA Glenn Research Center is developing a 1.4 MW high-efficiency electric machine for future electrified aircraft to reduce energy consumption, emissions, and noise. This wound-field, synchronous machine employs a self-cooled, superconducting rotor to achieve excellent specific power and efficiency. The design of the superconducting rotor and the optimization of its electromagnetic and structural responses are analyzed in this paper. Candidate designs are evaluated in terms of absolute performance, specific performance, and performance per cost. It is found that optimizing the electromagnetic response for absolute or specific performance yields designs that are similar to each other, but opposite of the cost-optimized design. A method to define the thermal requirements of the superconducting coils and integrated cryocooler is also presented.

Scheidler, Justin↗

Dawn Spacecraft Performance at Ceres: Results of Hybrid Control for Ceres Mapping

Dawn is a low-thrust interplanetary spacecraft currently orbiting the dwarf planet Ceres, to better understand the early creation of the solar system. Launched in September 2007, Dawn arrived at Vesta in July 2011. After a 16-month successful science campaign at Vesta, Dawn departed for Ceres, arriving in early 2015. The Dawn spacecraft uses both reaction wheel assemblies (RWA) and a reaction control system (RCS) to provide 3-axis attitude control for the spacecraft. Reaction wheels were designed to be the primary system for attitude control, however two wheels have shown high friction anomalies and have been removed from service. The project has implemented a hybrid control algorithm using two reaction wheels and RCS thrusters. This hybrid control capability enabled Dawn to achieve very high science return, while significantly conserving remaining hydrazine propellant. With only two remaining healthy RWA’s, hybrid control became part of the baseline plan for Ceres science operations. The Dawn team developed specific operational approaches in which sequences were developed with careful consideration of science versus resource trades. Commanding and sequence planning also incorporated contingency planning, in the event that another reaction wheel may fail. Despite the differences in operational approach between Vesta and Ceres, both campaigns achieved very rich scientific data return. This paper highlights Dawn’s recent flight experience with hybrid attitude control during Ceres orbit operations. The discussion includes the approaches utilized by the Dawn team to address unique operational challenges presented by the hybrid approach, and reviews spacecraft performance under hybrid control in low orbit at Ceres. Additionally, methods used to optimize hydrazine use and thereby extend the science campaign will be presented. Finally, a preliminary assessment of an orbit transfer with two reaction wheels, during extended mission operations, is discussed.

Smith, Brett A.↗

Liftoff and Transition Database Generation for Launch Vehicles Using Data-Fusion-Based Modeling

A data fusion technique for merging multiple data sources with differing fidelity and resolution was developed to support the production of aerodynamic line load databases for the Liftoff and Transition (LOT) flight phase of the Space Launch System (SLS). The technique uses a reduced order model based on a high-fidelity line load data set from Computational Fluid Dynamics (CFD) to predict solutions for a much larger solution space. Even higher-fidelity force and moment information (from wind-tunnel tests) is then used to adjust the model. The adjustment uses constrained optimization through the method of Lagrange multipliers in order to minimize the deviation of the line load distribution from the spatially-dense CFD solution, while ensuring that the integrated force and moment values match those observed in physical wind tunnel measurements. Though the wind-tunnel data are operationally-dense (available at many flow conditions), they are spatially coarse (as only the overall forces and moments are available). Conversely, CFD for such complex configurations is expensive, and thus operationally sparse. Data fusion techniques are necessary to make the most efficient use of available information, delivering accurate results within time and resource constraints.

Wignall, T. J.↗

Validation of CrIS Ozone and Carbon Oxide Products from CrIS Single-Field-of-View Retrievals

Single Field-of-view Sounder Atmospheric Products (SiFSAP) from Cross-track Infrared Sounder (CrIS) on S-NPP have been developed at NASA Langley Research center, and recently a significant improvement in its trace gases products has been made. These products are derived based on the Principal Component (PC)-based Radiative Transfer Model (PCRTM) and an optimal estimation retrieval method (PCRTM-RA). Use of the PCRTM in this method allows for fast and accurate calculations for thousands of spectral channels under clear and cloudy sky conditions. The use of principal component to compress information from all spectral channels into the PC domain enables to reduce the observational noise and use the information from all hyperspectral infrared channels. The retrieval is made on PC-domain, and temperature, water vapor, trace gases, cloud and surface parameters are retrieved simultaneously. This presentation will present some validation to the SiFSAP ozone (O3) and carbon monoxide (CO) products. Since the SiFSAP has a spatial resolution of about 14 km at nadir, which is better than most global weather and climate models and enables to make process-oriented analysis of the transport of these gases, this study will focus on validation of CO and O3 products under some extreme weather and/or climate events, such as CO emission and transport from large biomass burning, and O3 variation impacted by fire emission, stratospheric intrusion, and tropical cyclones. Data used in our validation include CO and O3 from NASA aircraft campaigns, plus other satellite products such as CO and O3 from AIRS, IASI and TROPOMI, MOPPIT CO and OMPS O3. Some comparison of SiFSAP O3 with the O3 from the fifth-generation ECMWF reanalysis (ERA5) data the NASA Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) will also present. Such a validation not only helps to quantify the performance of the retrievals in these extreme conditions, but also displays the capability of satellite observation to capture the fine-scale spatial pattern of CO and O3 variation.

Xiaozhen Xiong↗

Global Ionospheric Electron Density from GNSS-POD Limb Measurements

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in studying and monitoring the global ionospheric electron density (Ne). In this study we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles. The hmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006-2019), COSMIC-2 (2019-present), and Spire (2019-present) constellations, showing a consistent ionospheric Ne morphology. A detailed analysis of the frequency-wavenumber spectra is made for the Ne variability at different heights. In the lower ionosphere (~150 km) we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2 and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4 and SW4. The co-exist of eastward and westward propagating WN4 components implies the presence of a stationary planetary wave4 (sPW4), as suggested by other earlier studies.

GNSS-POD↗

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems.

Tank Chilldown↗

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

Abstract. In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems. Keywords: tank chilldown, charge-hold-vent, no-vent fill

Tank Chilldown↗

Modeling of Multicomponent Trace Contaminant Adsorption on Activated Carbon

Modeling of multicomponent trace contaminant adsorption on activated carbon for spacecraft life support applications presents challenges due to the very low contaminant concentrations in cabin air coupled with very strong adsorption of heavier components such as cyclic polydimethylsiloxanes. Both gas-phase and solid-phase micropore mass transfer resistances can be important, and competitive interactions between adsorbed components can result in behavior such as roll-up where lighter components displaced by heavier components have a higher concentration at the outlet than at the inlet. This paper describes the further development of a dynamic trace contaminant adsorption model based on Ideal Adsorbed Solution Theory (IAST). Alternative 1-dimensional and 2-dimensional versions of the model are described with different mass transfer resistance and driving force assumptions. Optimization of solution methods for speed and stability is also described. Mass balance equations are scaled to provide sensitivity over the wide range of component concentrations in cabin air. Model parameters are estimated from isotherm and breakthrough test data on single adsorbate and binary adsorbate systems. Multicomponent results are compared against predictions of a heritage model used extensively by NASA for trace contaminant control system (TCCS) carbon bed sizing. Finally, efforts to implement the model in a second software platform for wider access are described.

Kevin E Lange↗

Singular perturbation analysis of AOTV-related trajectory optimization problems

The problem of real time guidance and optimal control of Aeroassisted Orbit Transfer Vehicles (AOTV's) was addressed using singular perturbation theory as an underlying method of analysis. Trajectories were optimized with the objective of minimum energy expenditure in the atmospheric phase of the maneuver. Two major problem areas were addressed: optimal reentry, and synergetic plane change with aeroglide. For the reentry problem, several reduced order models were analyzed with the objective of optimal changes in heading with minimum energy loss. It was demonstrated that a further model order reduction to a single state model is possible through the application of singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy level of the vehicle. A separate boundary layer analysis is used to account for altitude and flight path angle dynamics, and to obtain lift and bank angle control solutions. By considering alternative approximations to solve the boundary layer problem, three guidance laws were derived, each having an analytic feedback form. The guidance laws were evaluated using a Maneuvering Reentry Research Vehicle model and all three laws were found to be near optimal. For the problem of synergetic plane change with aeroglide, a difficult terminal boundary layer control problem arises which to date is found to be analytically intractable. Thus a predictive/corrective solution was developed to satisfy the terminal constraints on altitude and flight path angle. A composite guidance solution was obtained by combining the optimal reentry solution with the predictive/corrective guidance method. Numerical comparisons with the corresponding optimal trajectory solutions show that the resulting performance is very close to optimal. An attempt was made to obtain numerically optimized trajectories for the case where heating rate is constrained. A first order state variable inequality constraint was imposed on the full order AOTV point mass equations of motion, using a simple aerodynamic heating rate model.

Calise, Anthony J.↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗

On the constrained Chebyshev approximation problem on ellipses

Constrained Chebyshev approximation problems of the type with minimum (p is an element of Pi(sub n):p(c)=1) and maximum (z is an element of E) with /p(z)/ are considered. Here Pi(sub n) denotes the set of all complex polynomials of degree at most n, E is any ellipse in the complex plane, and c is an element of C/E. Such approximation problems arise in the context of optimizing semi-iterative methods for the solution of large, sparse systems of linear equations Ax=b with complex non-Hermitian coefficient matrices A. The problem of obtaining optimal polynomial preconditioners for conjugate gradient type methods for Ax=b also leads to problems of this type. A new family of polynomials -- q(sub n)(z;c), n is an element of N, and c is an element of C/E -- are introduced as the polynomials which are optimal for a modified version of the Chebyshev approximation problem with Pi(sub n) replaced by a certain subfamily. Some simple properties of q(sub n) are also listed. A necessary and sufficient condition for q(sub n) to be the extremal polynomial for the approximation problem is then derived. Finally, it is shown that q(sub n) is indeed optimal for the problem for all fixed n whenever the distance between c and E is sufficiently large. Results of some numerical tests are presented.

Freund, Roland↗

Enhancements on the Convex Programming Based Powered Descent Guidance Algorithm for Mars Landing

In this paper, we present enhancements on the powered descent guidance algorithm developed for Mars pinpoint landing. The guidance algorithm solves the powered descent minimum fuel trajectory optimization problem via a direct numerical method. Our main contribution is to formulate the trajectory optimization problem, which has nonconvex control constraints, as a finite dimensional convex optimization problem, specifically as a finite dimensional second order cone programming (SOCP) problem. SOCP is a subclass of convex programming, and there are efficient SOCP solvers with deterministic convergence properties. Hence, the resulting guidance algorithm can potentially be implemented onboard a spacecraft for real-time applications. Particularly, this paper discusses the algorithmic improvements obtained by: (i) Using an efficient approach to choose the optimal time-of-flight; (ii) Using a computationally inexpensive way to detect the feasibility/ infeasibility of the problem due to the thrust-to-weight constraint; (iii) Incorporating the rotation rate of the planet into the problem formulation; (iv) Developing additional constraints on the position and velocity to guarantee no-subsurface flight between the time samples of the temporal discretization; (v) Developing a fuel-limited targeting algorithm; (vi) Initial result on developing an onboard table lookup method to obtain almost fuel optimal solutions in real-time.

Guidance↗

Structural/control system optimization with variable actuator masses

A method is presented to integrate the design space for structural/control system optimization problems in the case of linear state feedback control. Nonstructural lumped masses and control system design variables as well as structural sizing variables are all treated equally as independent design variables in the optimization process. Structural and control design variable linking schemes are used in order to avoid a prohibitively large increase in the total number of independent design variables. When actuator masses are treated as nonstructural lumped mass design variables, special consideration is given to the relation between the transient peak responses and the required actuator masses which is formulated as a behavior constraint form. The original nonlinear mathematical programming problem based on a finite element formulation and linear state feedback is replaced by a sequence of explicit approximate problems exploiting various approximation concepts such as design variable linkings, temporary constraint deletion and first order Taylor series expansion of nonlinear behavior constraints in terms of intermediate design variables. Examples which involve a variety of dynamic behavior constraints (including constraints on closed-loop eigenvalues, peak transient displacements, peak actuator forces, and relations between the peak responses and the actuator masses) are effectively solved by using the method presented.

Jin, Ik M.↗

Multifidelity Optimization with Transonic Flutter Constraints

This work considers static and dynamic aeroelastic optimization of a cantilevered platewing in transonic flow. Low- and high-fidelity aeroelastic predictions are fed into a standard trust region model management scheme in order to efficiently solve this expensive optimization problem with multifidelity methods. Differences in fidelity are entirely driven by different aerodynamic solvers: linear panel methods for the low-fidelity method, and inviscid CFD-based solvers for the high-fidelity response. The optimization problem utilizes shape, sizing, and trim variables to satisfy stress, trim, and flutter constraints, and is solved across a range of subsonic and transonic Mach numbers. The multifidelity method is able to provide a speed-up for all cases considered here, despite sizable inaccuracies in the low-fidelity response for transonic flows.

Bret K Stanford↗

Peak-Seeking Optimization of Spanwise Lift Distribution for Wings in Formation Flight

A method is presented for the optimization of the lift distribution across the wing of an aircraft in formation flight. The usual elliptical distribution is no longer optimal for the trailing wing in the formation due to the asymmetric nature of the encountered flow field. Control surfaces along the trailing edge of the wing can be configured to obtain a non-elliptical profile that is more optimal in terms of minimum drag. Due to the difficult-to-predict nature of formation flight aerodynamics, a Newton-Raphson peak-seeking controller is used to identify in real time the best aileron and flap deployment scheme for minimum total drag. Simulation results show that the peak-seeking controller correctly identifies an optimal trim configuration that provides additional drag savings above those achieved with conventional anti-symmetric aileron trim.

Hanson, Curtis E.↗

Multilevel control optimization using subsystem relative performance index sensitivity

A method is presented for the design of optimal feedback controllers for large multivariable systems with subsystem sensitivity constraints. The weighted sum of subsystem and/or operational mode relative performance index sensitivities is defined as the overall performance index. The method is developed for linear systems with quadratic performance criteria and either full or partial state feedback. An example concerning the design of a stability augmentation system for a VTOL aircraft in the transition mode demonstrates the effectiveness of the design method.

Leininger, G. G.↗