Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear approximation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A new approach to approximating the linear quadratic optimal control law for hereditary systems with control delays

A factorization approach is presented for deriving approximations to the optimal feedback gain for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the feedback kernels.

Milman, M. H.↗
NASA NTRS1987

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗
NASA NTRS1988

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗
NASA NTRS1992

Efficient yet accurate solution of the linear transport equation in the presence of internal sources - The exponential-linear-in-depth approximation

The present solutions to the linear transport equation pertain to monoenergetic particles' interaction with a multiple scattering/absorbing layered medium with a general anisotropic internal source term. Attention is given to a novel exponential-linear approximation to the internal source, as a function of scattering depth, which furnishes an at-once efficient and accurate solution to the linear transport equation through its reduction of the spatial mesh size. The great superiority of the proposed method is demonstrated by the numerical results obtained in the illustrative cases of (1) an embedded thermal source and (2) a rapidly varying beam pseudosource.

Kylling, Arve↗
NASA NTRS2009

Bounded Linear Stability Analysis - A Time Delay Margin Estimation Approach for Adaptive Control

This paper presents a method for estimating time delay margin for model-reference adaptive control of systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The time delay margin of this system represents a local stability measure and is computed analytically by a matrix measure method, which provides a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation.

Nguyen, Nhan T.↗
NASA NTRS2024

Global Optimization via Quadratic Disjunctive Programming for Water Networks Design with Energy Recovery

Generalized disjunctive programming (GDP) models with bilinear and concave constraints, often seen in water network design, are challenging optimization problems. This work proposes quadratic and piecewise linear approximations for nonlinear terms to reformulate GDP models into quadratic GDP (QGDP) models that suitable solvers may solve more efficiently. We illustrate the benefits of the quadratic reformulation with a water treatment network design problem in which nonconvexities arise from bilinear terms in the mixers’ mass balances and concave investment cost functions of treatment units. Given the similarities with water network design problems, we suggest quadratic approximation for the GDP model for the optimal design of a large-scale reverse electrodialysis (RED) process. This power technology can recover energy from salinity differences between by-product streams of the water sector, such as desalination brine mixed with regenerated wastewater effluents. The solver Gurobi excels in handling QGDP problems, but weighing the problem’s precision and tractability balance is crucial. The piecewise linear approximation yields more accurate, yet larger QGDP models that may require longer optimization times in large-scale process synthesis problems.

Water Networks↗
NASA NTRS1986

Output time response approximation

The approximation of the output response of a nonlinear system by the output response of a linear system to a desired order irrespective of the admissible input applied should prove useful for purposes of control generation and simulation. Given a nonlinear system, an integer k, and an open subset of state space, sufficient conditions are stated that such a linear approximation exists to order k for every point in the set. In addition, a method for finding the approximating linear systems is presented.

Hunt, L. R.↗
NASA NTRS1992

Stacking sequence optimization of simply supported laminates with stability and strain constraints

An integer programming formulation for the design of symmetric and balanced rectangular composite laminates with simply supported boundary conditions subject to buckling and strain constraints is presented. The design variables that define the stacking sequence of the laminate are ply-identity zero-one integers. The buckling constraint is linear in terms of the ply-identity design variables, but strains are nonlinear functions of these variables. A linear approximation is developed for the strain constraints so that the problem can be solved by sequential linearization using the branch and bound algorithm. Examples of graphite-epoxy plates under biaxial compression are presented. Optimum stacking sequences obtained using the linear approximation are compared with global optimum designs obtained using a genetic search procedure.

Nagendra, S.↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2023

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗
NASA NTRS1988

Improving cover type identification in speckled SAR images by prefiltering and sequential classification

Synthetic aperture radar utilizes coherent microwaves to produce images of the earth's surface. Due to the interference of coherent wavelets, the images appear speckled. This reduces the performance of per-pel classifiers. One way to increase the performance is to filter the image first, then classify the filtered image. For this purpose, several novel filters that have been reported in the literature are investigated. These are the geometric filter, adaptive LMMSE filter, and linear approximation filter. For comparison, conventional mean and median filters are also considered. It is found that the mean filter with seven iterations gives the best result. The overall performance increased from 65.2 to 88.9 percent. The capability of these filters to preserve edges in the original image are also assessed. It is seen that the geometric and median filters are the best in preserving edges, and that the linear approximation and adaptive LMMSE filters are the best in discriminating roads. Prefiltering the image effectively provides contextual information to the per-pel classifier. An alternate approach is to directly design a contextual classifier. A new contextual classifier based on sequential decision theory is proposed. With this classifier, it is found that the overall performance increases to 89.5 percent.

Lin, Qian↗
DOE OSTI2026

Two-stage formation-energy correction (NbZr, TaZr, VZr)

This bundle contains the scripts, the raw and corrected per-structure data, and the manuscript plots for the NbZr / TaZr / VZr BCC binary formation energies and the associated RMSDs. Why a two-stage correction is necessary: The "raw" formation energy of every relaxed VASP configuration is computed in the usual way, FE_raw(c) = E_alloy(c) - sum_i x_i * E_pure_i , where E_pure_i are the per-atom total energies of the pure-element reference structures (Nb, Ta, V, Zr in the same BCC supercell, with identical INCAR / KPOINTS / PAW choices). With perfectly consistent reference runs the raw FE should vanish at the two pure-element endpoints (x = 0 and x = 1) by construction. In practice this does not hold for two reasons that are present in our dataset: 1. Reference-energy inconsistency (composition-dependent bias). Even with identical input parameters, the pure-element runs (stored in `corrected_DFT_pure_element_runs/`) differ slightly from the values that would be implied by the alloy runs at near-pure compositions (a few meV/atom). This bias is approximately linear in concentration, because the residual error in E_pure_Nb (or E_pure_Ta / E_pure_V) propagates into FE_raw(c) as (1 - x) * dE_pure_1, and the corresponding error in E_pure_Zr propagates as x * dE_pure_2. Left uncorrected, this produces a non-physical "tilt" of FE_raw(x) and shifts the entire FE-vs-x cloud away from zero at the endpoints. 2. Endpoint anchoring against the audited true endpoints. The strict endpoint values (FE_x0_meVatom, FE_x1_meVatom in `corrected_fe_strict_endpoints_20260518/strict_endpoint_check_20260518.csv`) were re-derived from an independent cross-check of the pure-element runs. After stage 1 removes the linear bias, the near-pure compositions in the alloy dataset still extrapolate to values that differ slightly from these audited endpoints — because stage 1 is fit from a few near-end alloy bins, not from the audited pure-element references themselves. The README.txt file discusses how these issues are addressed by the two-stage correction, and describes folder layout, pipeline summary, and how to re-run.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Combining global and local approximations

A method based on a linear approximation to a scaling factor, designated the 'global-local approximation' (GLA) method, is presented and shown capable of extending the range of usefulness of derivative-based approximations to a more refined model. The GLA approach refines the conventional scaling factor by means of a linearly varying, rather than constant, scaling factor. The capabilities of the method are demonstrated for a simple beam example with a crude and more refined FEM model.

Haftka, Raphael T.↗
NASA NTRS1972

Weighted least squares stationary approximations to linear systems.

Investigation of the problem of replacing a certain time-varying linear system by a stationary one. Several quadratic criteria are proposed to aid in determining suitable candidate systems. One criterion for choosing the matrix B (in the stationary system) is initial-condition dependent, and another bounds the 'worst case' homogeneous system performance. Both of these criteria produce weighted least square fits.

Bierman, G. J.↗
NASA NTRS1993

Time accuracy and the use of implicit methods

Some of the approximations used to make implicit methods more efficient and practical for the solution of the Euler and Navier-Stokes equations are addressed. In particular, approximate factorizations, diagonalizations, and linearization approximations are reviewed and categorized. A subiteration correction scheme commonly in use at present is introduced, improved, demonstrated, and analyzed. This scheme is used to produce a second-order accurate, more robust implicit method for unsteady flow computations. The subiteration approach can be employed to recover time accuracy without increasing computational time (in most cases producing substantial savings).

Pulliam, Thomas H.↗
NASA NTRS1969

Numerical results in the statistical theory of thermal turbulence

Several statistical theories of turbulence are applied to thermal convection problem, and are compared in a series of numerical calculations. Results are presented for both the case of very large and very small Prandtl numbers, and for moderate Rayleigh numbers. The theories compared are the direct-interaction approximation, the self-consistent field approximation, the quasi-linear approximation, the quasi-normal approximation, and a 'Markovianized' quasi-normal approximation. Results for these theories are then compared with some simple exact solutions to the statistical problem. The results indicate that the direct-interaction and the self-consistent field method give satisfactory results at all Rayleigh numbers investigated, whereas the procedures based on the quasi-normal approximation give sensible results only at small Rayleigh numbers.

Herring, J. R.↗
NASA NTRS1996

Droplet Deformation in an Extensional Flow: The Role of Surfactant Physical Chemistry

Surfactant-induced Marangoni effects strongly alter the stresses exerted along fluid particle interfaces. In low gravity processes, these stresses can dictate the system behavior. The dependence of Marangoni effects on surfactant physical chemistry is not understood, severely impacting our ability to predict and control fluid particle flows. A droplet in an extensional flow allows the controlled study of stretching and deforming interfaces. The deformations of the drop allow both Marangoni stresses, which resist tangential shear, and Marangoni elasticities, which resist surface dilatation, to develop. This flow presents an ideal model system for studying these effects. Prior surfactant-related work in this flow considered a linear dependence of the surface tension on the surface concentration, valid only at dilute surface concentrations, or a non-linear framework at concentrations sufficiently dilute that the linear approximation was valid. The linear framework becomes inadequate for several reasons. The finite dimensions of surfactant molecules must be taken into account with a model that includes surfaces saturation. Nonideal interactions between adsorbed surfactant molecules alter the partitioning of surfactant between the bulk and the interface, the dynamics of surfactant adsorptive/desorptive exchange, and the sensitivity of the surface tension to adsorbed surfactant. For example, cohesion between hydrocarbon chains favors strong adsorption. Cohesion also slows the rate of desorption from interfaces, and decreases the sensitivity of the surface tension to adsorbed surfactant. Strong cohesive interactions result in first order surface phase changes with a plateau in the surface tension vs surface concentration. Within this surface concentration range, the surface tension is decoupled from surface concentration gradients. We are engaged in the study of the role of surfactant physical chemistry in determining the Marangoni stresses on a drop in an extensional flow in a numerical and experimental program. Using surfactants whose dynamics and equilibrium behavior have been characterized in our laboratory, drop deformation will be studied in ground-based experiment. In an accompanying numerical study, predictive drop deformations will be determined based on the isotherm and equation of state determined in our laboratory. This work will improve our abilities to predict and control all fluid particle flows.

Stebe, Kathleen J.↗