On a class of numerical methods with an adaptive integration subprocedure for optimal control problems
Explore the source record for details and available documents.
Polak, E.↗
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.
Comparative study of numerical schemes of TVD3, UNO3-ACM and optimized compact scheme
Three different schemes are employed to solve the benchmark problem. The first one is a conventional TVD-MUSCL (Monotone Upwind Schemes for Conservation Laws) scheme. The second scheme is a UNO3-ACM (Uniformly Non-Oscillatory Artificial Compression Method) scheme. The third scheme is an optimized compact finite difference scheme modified by us: the 4th order Runge Kutta time stepping, the 4th order pentadiagonal compact spatial discretization with the maximum resolution characteristics. The problems of category 1 are solved by using the second (UNO3-ACM) and third (Optimized Compact) schemes. The problems of category 2 are solved by using the first (TVD3) and second (UNO3-ACM) schemes. The problem of category 5 is solved by using the first (TVD3) scheme. It can be concluded from the present calculations that the Optimized Compact scheme and the UN03-ACM show good resolutions for category 1 and category 2 respectively.
Lee, Duck-Joo↗
A wavelet-optimized, very high order adaptive grid and order numerical method
Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.
Jameson, Leland↗
Annual Research Briefs - 2000: Center for Turbulence Research
This report contains the 2000 annual progress reports of the postdoctoral Fellows and visiting scholars of the Center for Turbulence Research (CTR). It summarizes the research efforts undertaken under the core CTR program. Last year, CTR sponsored sixteen resident Postdoctoral Fellows, nine Research Associates, and two Senior Research Fellows, hosted seven short term visitors, and supported four doctoral students. The Research Associates are supported by the Departments of Defense and Energy. The reports in this volume are divided into five groups. The first group largely consists of the new areas of interest at CTR. It includes efficient algorithms for molecular dynamics, stability in protoplanetary disks, and experimental and numerical applications of evolutionary optimization algorithms for jet flow control. The next group of reports is in experimental, theoretical, and numerical modeling efforts in turbulent combustion. As more challenging computations are attempted, the need for additional theoretical and experimental studies in combustion has emerged. A pacing item for computation of nonpremixed combustion is the prediction of extinction and re-ignition phenomena, which is currently being addressed at CTR. The third group of reports is in the development of accurate and efficient numerical methods, which has always been an important part of CTR's work. This is the tool development part of the program which supports our high fidelity numerical simulations in such areas as turbulence in complex geometries, hypersonics, and acoustics. The final two groups of reports are concerned with LES and RANS prediction methods. There has been significant progress in wall modeling for LES of high Reynolds number turbulence and in validation of the v(exp 2) - f model for industrial applications.
Source record↗
Robust eigenstructure assignment by a projection method - Applications using multiple optimization criteria
A methodology for robust eigenstructure assignment for multivariable feedback systems is presented. The algorithm is based upon a pole placement technique using projections onto subspaces of admissible eigenvectors. New ideas are introduced to generate target (desired) sets of unitary eigenvectors and determine optimal feasible eigenvectors in a least-square sense. Useful connections are established between the pole-placement by independent modal space control and the method introduced in this paper. A multicriterion optimization algorithm is also presented, which takes efficient advantage of the present eigenstructure assignment method. These developments show significant improvement over an earlier version of this algorithm in both computational cost and accuracy. This optimization process appears to be numerically robust and suitable for high-dimensional multicriterion optimizations; it is especially attractive for computer-aided design of control systems.
Rew, D. W.↗
A least-squares finite element method for incompressible Navier-Stokes problem
Most finite element schemes for solving the Navier-Stokes equations can be categorized into the Galerkin mixed method and the penalty method. The mixed method leads to a saddle-point problem. In order to guarantee the existence of a solution, the combination of velocity and pressure interpolations requires satisfaction of the Ladyzhenskaya Babuska Brezzi (LBB) consistency condition which precludes the use of equal order interpolations and many seemingly natural pairs of velocity and pressure elements. In a previous paper a least-squares finite element method based on the first order velocity-pressure-vorticity formulation for the Stokes problem was proposed. This method leads to a minimization problem. The choice of combination of elements is thus not subject to the LBB condition. The numerical experiments exhibit the optimal rate of convergence for all variables with equal order interpolations. A theoretical error analysis supports the numerical results. In this paper the least-squares finite element method is extended to solving the incompressible Navier-Stokes problem.
Jiang, Bonan↗
Enhancement of disposal efficiency for deep geological repositories based on three design factors - Decay heat optimization, increased thermal limit of the buffer and double-layer concept
This study investigates the enhancement of disposal efficiency for deep geological repositories (DGRs) based on three design factors: decay heat optimization, increased thermal limit of the buffer, and double-layer concept using coupled thermo-hydro-mechanical (THM) numerical simulations. Decay heat optimization is achieved by iteratively emplacing spent nuclear fuels having the maximum and minimum decay heat in a canister. Disposal areas can be reduced by 20 % to 40 % compared to the current reference disposal system in Korea (KRS+) in accordance with the combinations of the three design factors, alleviating challenges in site selection for the DGR. This study additionally identifies an optimal layer spacing of 500 m for the double-layer concept in the viewpoint of the buffer temperature, where thermal interaction between the upper and lower layers nearly disappears. However, determining the ultimate disposal and layer spacing requires engineering judgement, considering not only the thermal performance of the DGR but also various factors such as cost and difficulties of the construction and rock mass stability. DGRs designed with an increased thermal limit of the buffer poses a greater probability of rock mass failure around disposal tunnels and deposition holes due to elevated thermal stresses. Densely arranged heat sources for the DGRs with enhanced disposal efficiency lead to larger temperature increase even at the far-field scale, raising a possibility of thermally driven fracture shear activation with associated hydraulic, mechanical, and seismic changes.
58 GEOSCIENCES↗
Deterministic optimal maneuver strategy for multi-target missions
This paper presents an optimal strategy for making impulsive correction to a multi-target trajectory by a single maneuver. The concept of an optimal maneuver time is introduced. The choice of suitable weighting functions is explored to enable one to properly translate the subjective desire of mission success into an objective cost function whose minimization yields the optimal strategy. It is shown that a number of strategies previously formulated are derivable from one general expression. A number of other interesting properties of the optimal strategy are described. Numerical results are presented for a typical two-target mission. It is shown that the strategy formulated is optimal. For some perturbations, there exists an optimal maneuver time different from the time of initiation of the perturbation. That is, the physical properties of the trajectory can be exploited to select the optimal time of making a corrective maneuver.
Dwivedi, N. P.↗
Applications of generalized inverses to estimation in dynamic systems.
Generalized matrix inverses are used to obtain an estimation procedure for estimation of the state vector of a dynamic system. This sequential procedure is studied analytically with respect to the choice of an arbitrary vector. The covariance matrix of the estimator is determined and compared to the optimal Kalman type procedure. A numerical example illustrates the procedure and compares it to the optimal one.
Prochaska, B. J.↗
Robust Neighboring Optimal Guidance for the Advanced Launch System
In recent years, optimization has become an engineering tool through the availability of numerous successful nonlinear programming codes. Optimal control problems are converted into parameter optimization (nonlinear programming) problems by assuming the control to be piecewise linear, making the unknowns the nodes or junction points of the linear control segments. Once the optimal piecewise linear control (suboptimal) control is known, a guidance law for operating near the suboptimal path is the neighboring optimal piecewise linear control (neighboring suboptimal control). Research conducted under this grant has been directed toward the investigation of neighboring suboptimal control as a guidance scheme for an advanced launch system.
Hull, David G.↗
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↗
An Initial Assessment of Variable Depth Liner Optimization for Ducted Proprotor Applications
The rise of the Urban Air Mobility market has spurred the design of a new generation of novel aircraft. To aid industry and researchers interested in these types of aircraft, the Revolutionary Vertical Lift Technology project at NASA has developed a fleet of reference vehicles for system studies. A new six-passenger reference vehicle has recently been added to the research portfolio that has ducted proprotors for propulsors. The ducts present an opportunity to apply acoustic treatment to the interior of the duct that could target both tonal and broadband noise, representative of the sound produced by this type of propulsor. In this paper, design methodologies are presented to design a variable depth liner for this application. An optimizer is used to design multiple liners with variable chamber depths. Experimental results from normal impedance testing are compared to numerical predictions using the optimizer model and a finite element model. Results show that an optimizer can be used to design an acoustic liner with favorable performance for a broad range of frequencies, which could be appropriate for a ducted proprotor application.
Matthew B Galles↗
Locally purified maximally mixed states at scale: Entanglement pruning and symmetries
Locally Purified Density Operators (LPDOs) are state-of-the-art tensor network ansatze candidates that efficiently represent mixed quantum states at scale. However, given their non-uniqueness, their representational complexity is generally sub-optimal in practical computations. Here, in this work we perform a comprehensive numerical and analytical analysis and resolve this issue in the experimentally relevant limit where noise depolarizes the density operator into a maximally mixed state. To resolve the sub-optimality issue, we analyze two numerical tools, one analytic method, and detail the relations between them. The numerical tools used are fidelity-preserving truncations and isometric gauge transformations leveraging Riemannian optimizations over entropic objective functions. In addition, by invoking the injectivity and symmetry constraints of the maximally mixed LPDO, we also present analytical closed-form expressions for the disentangler and discuss their relation to numerical optimizers. Further, away from the maximally mixed state, our simulations highlight how the truncation threshold smoothly interpolate, as a function of depolarization, between established matrix product results and our new results. Our work shows how, by minimizing the resources required to represent key states of practical interest in experiment, the efficiency of tensor network algorithms can be substantially increased. This paves the path for uncovering tensor network’s fundamental scalability limits and latent potential in representing the wide locus of mixed quantum states that are accessible on near-term quantum devices.
Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
Numerical Modeling of a Vortex Stabilized Arcjet
Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.
Pawlas, Gary E.↗
Analysis of Photothermal Characterization of Layered Materials: Design of Optimal Experiments
In this paper numerical calculations are presented for the steady-periodic temperature in layered materials and functionally-graded materials to simulate photothermal methods for the measurement of thermal properties. No laboratory experiments were performed. The temperature is found from a new Green s function formulation which is particularly well-suited to machine calculation. The simulation method is verified by comparison with literature data for a layered material. The method is applied to a class of two-component functionally-graded materials and results for temperature and sensitivity coefficients are presented. An optimality criterion, based on the sensitivity coefficients, is used for choosing what experimental conditions will be needed for photothermal measurements to determine the spatial distribution of thermal properties. This method for optimal experiment design is completely general and may be applied to any photothermal technique and to any functionally-graded material.
Cole, Kevin D.↗
Shape optimization governed by the Euler equations using an adjoint method
A numerical approach for the treatment of optimal shape problems governed by the Euler equations is discussed. Focus is on flows with embedded shocks. A very simple problem is considered: the design of a quasi-one-dimensional Laval nozzle. A cost function and a set of Lagrange multipliers are introduced to achieve the minimum. The nature of the resulting costate equations is discussed. A theoretical difficulty that arises for cases with embedded shocks is pointed out and solved. Finally, some results are given to illustrate the effectiveness of the method.
Iollo, Angelo↗
Optimal selection of passes
Preliminary numerical results obtained from the application of a linear feature selection technique to the determination of combinations of passes which best discriminate between a given set of crops in a given area of interest, are reported. The results obtained are not purported to hold in a general situation, but only for the given set of crops and the given, but unknown, levels of several factors-such as soil type, and fertilizer practice, holding in the area of interest. However, by identifying the various factors affecting the spectral signatures, and by formulating a regression model one could use the feature selection technique to determine the regression coefficients for predicting optimal passes for a given set of crops. Another use of the feature selection technique as applied to multiple pass registered data is the generation of enhanced grey scale displays by using a single linear combination of all channels of all designated passes as opposed to a single channel within a single pass.
Guseman, L. F., Jr.↗
Heat Transfer Analysis of a Closed Brayton Cycle Space Radiator
This paper presents a mathematical analysis of the heat transfer processes taking place in a radiator for a closed cycle gas turbine (CCGT), also referred to as a Closed Brayton Cycle (CBC) space power system. The resulting equations and relationships have been incorporated into a radiator sub-routine of a numerical triple objective CCGT optimization program to determine operating conditions yielding maximum cycle efficiency, minimum radiator area and minimum overall systems mass. Study results should be of interest to numerical modeling of closed cycle Brayton space power systems and to the design of fluid cooled radiators in general.
Juhasz, Albert J.↗