Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional convergence”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

An input adaptive, pursuit tracking model of the human opertor

Developed and evaluated is a simple model of the input adaptive behavior of the human operator (HO) in a pursuit tracking task in which the plant controlled consists of a pure gain. If it is assumed that the HO is approximately an optimal predictor using only position and velocity information, then there is a simple method of computing the values of the model parameters in terms of the autocorrelation function of the input signal. Experimental evidence indicates that the ability of the HO to use velocity information decreases with increasing signal velocity indicating that a biased estimator of the velocity weighting should be used. A suitable approximation is derived which has rapid convergence and low variance. The model thus derived is compared to actual subject transfer functions and is found to be in close agreement. In addition to tracking random processes the model can adapt to and track deterministic signals, such as sine waves, up to approximately the frequency at which human operators begin to track precognitively.

Ware, J. R.↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
NASA NTRS2012

P-Wave Electron-Hydrogen Scattering

A variational wave function incorporating short range correlations via Hylleraas type functions plus long-range polarization terms of the polarized orbital type but with smooth cut-off factors has been used to calculate P-wave phase shifts for electron-hydrogen scattering. This approach gives the direct r(exp -4) potential and a non-local optical potential which is definite. The resulting phase shifts have rigorous lower bounds and the convergence is much faster than those obtained without the modification of the target function. Final results will be presented at the conference.

Bhtia, Anand↗
NASA NTRS1991

Optimal design of solidification processes

An optimal design algorithm is presented for the analysis of general solidification processes, and is demonstrated for the growth of GaAs crystals in a Bridgman furnace. The system is optimal in the sense that the prespecified temperature distribution in the solidifying materials is obtained to maximize product quality. The optimization uses traditional numerical programming techniques which require the evaluation of cost and constraint functions and their sensitivities. The finite element method is incorporated to analyze the crystal solidification problem, evaluate the cost and constraint functions, and compute the sensitivities. These techniques are demonstrated in the crystal growth application by determining an optimal furnace wall temperature distribution to obtain the desired temperature profile in the crystal, and hence to maximize the crystal's quality. Several numerical optimization algorithms are studied to determine the proper convergence criteria, effective 1-D search strategies, appropriate forms of the cost and constraint functions, etc. In particular, we incorporate the conjugate gradient and quasi-Newton methods for unconstrained problems. The efficiency and effectiveness of each algorithm is presented in the example problem.

Dantzig, Jonathan A.↗
DOE OSTI2026

ZEUS: An Efficient GPU Optimization Method Integrating PSO, BFGS, and Automatic Differentiation

We introduce a novel, efficient computational method, ZEUS, for numerical optimization, and provide an open-source implementation. It has four key ingredients: (1) particle swarm optimization (PSO), (2) the use of the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method, (3) automatic differentiation (AD), and (4) GPUs. Our approach addresses the computational challenges inherent in high-dimensional, non-convex optimization problems. In the first phase of the algorithm, we get a potentially good set of starting points using PSO. Thereafter, we run BFGS independently in parallel from these starting points. BFGS is one of the best-performing algorithms for numerical optimization. However, it requires the gradient of the function being optimized. ZEUS integrates automatic differentiation into BFGS thus avoiding the need for the user to calculate derivatives explicitly. The use of GPUs allows ZEUS to speed up the calculations substantially. We carry out systematic studies to explore the trade-offs between the number of PSO iterations taken, starting points, and BFGS iteration depth. We show that a handful of iterations of PSO can improve global convergence when combined with BFGS. We also present performance studies using common test functions. The source code can be found at https://github.com/fnal-numerics/global-optimizer-gpu.

Soos, Dominik [Old Dominion U.]↗
NASA NTRS1975

A parallel Jacobson-Oksman optimization algorithm

A gradient-dependent optimization technique which exploits the vector-streaming or parallel-computing capabilities of some modern computers is presented. The algorithm, derived by assuming that the function to be minimized is homogeneous, is a modification of the Jacobson-Oksman serial minimization method. In addition to describing the algorithm, conditions insuring the convergence of the iterates of the algorithm and the results of numerical experiments on a group of sample test functions are presented. The results of these experiments indicate that this algorithm will solve optimization problems in less computing time than conventional serial methods on machines having vector-streaming or parallel-computing capabilities.

Straeter, T. A.↗
NASA NTRS1974

Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations

This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.

Gupta, N. K.↗
NASA NTRS1978

Global analysis of a buck regulator

Sufficient conditions for global stability of a buck regulator using a discrete control law are found. The method of paired systems and Liapunov functions are used to establish global stability and to study the convergence of the regulator. A heuristic argument is given that the optimal switching curves associated with the paired continuous systems approximate the optimal switching curves of the discrete systems.

Edwards, D. B.↗
NASA NTRS1979

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗
NASA NTRS1982

A penalty finite element algorithm for parabolic flow problems

The thin-layer simplification of the two-dimensional Navier-Stokes equations for steady viscous flow are developed using an order of magnitude analysis. A space marching finite element solution algorithm is developed, wherein the first order continuity effects are enforced as a penalty function differential constraint. Numerical results are presented to document accuracy and convergence features of the numerical solution algorithm.

Baker, A. J.↗
NASA NTRS1984

Estimation of discontinuous coefficients in parabolic systems: Applications to reservoir simulation

Spline based techniques for estimating spatially varying parameters that appear in parabolic distributed systems (typical of those found in reservoir simulation problems) are presented. The problem of determining discontinuous coefficients, estimating both the functional shape and points of discontinuity for such parameters is discussed. Convergence results and a summary of numerical performance of the resulting algorithms are given.

Lamm, P. D.↗
NASA NTRS1990

Asymptotic boundary conditions for dissipative waves: General theory

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗
NASA NTRS1990

Measured forces and moments on a delta wing during pitch-up

A series of low-speed wind tunnel tests on a 70-deg, sharp, leading-edge delta wing undergoing ramp pitching motion of high amplitude were performed to investigate the aerodynamic forces and moments. Forces and moments were obtained from a six-component interanl balance. Large amplitude oscillatory motion was produced by sinusoidally oscillating the model over a range of reduced frequencies. Ramp motion was produced by pitching the model through a half cycle of sinusoidal motion at a root chord Reynolds number of 1.54 million. The effect of ramp and oscillatory motions on the forces and moments are almost identical at matched pitch rates. Pitch rate had strong effect on the magnitude of the aerodynamic forces and moments. Upon completion of the model motion, some time is required for the forces and moments to decay to their static values. This convergence of the dynamic values to the static ones was a function of the pitch rate.

Bragg, M. B.↗
NASA NTRS1991

Asymptotic boundary conditions for dissipative waves - General theory

An outstanding issue in computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗
NASA NTRS1992

Flap-lag damping in hover and forward flight with a three-dimensional wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in-plane direction. Moreover, induced drag is sensitive to dynamic wake, both shed and trailing components, and thus its prediction requires adequate unsteady-wake representation. Accordingly, rigid-blade flap-lag equations are coupled with a three-dimensional finite-state wake model; three isolated rotor configurations with three, four and five blades are treated over a range of thrust levels, Lock numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of shape functions and harmonics of the wake model for multiblade modes of low frequency (less than 1/rev.) to high frequency (greater than l/rev.). Predicted flap and lag damping levels are then compared with similar predictions with (1) rigid wake (no unsteady induced flow), (2) Loewy lift deficiency, and (3) dynamic inflow.

Manjunath, A. R.↗
NASA NTRS1996

Study of Analytic Statistical Model for Decay of Light and Medium Mass Nuclei in Nuclear Fragmentation

The angular momentum independent statistical decay model is often applied using a Monte-Carlo simulation to describe the decay of prefragment nuclei in heavy ion reactions. This paper presents an analytical approach to the decay problem of nuclei with mass number less than 60, which is important for galactic cosmic ray (GCR) studies. This decay problem of nuclei with mass number less than 60 incorporates well-known levels of the lightest nuclei (A less than 11) to improve convergence and accuracy. A sensitivity study of the model level density function is used to determine the impact on mass and charge distributions in nuclear fragmentation. This angular momentum independent statistical decay model also describes the momentum and energy distribution of emitted particles (n, p, d, t, h, and a) from a prefragment nucleus.

Cucinotta, Francis A.↗