SEARCH · Search NASA

Results for “local solution bounds”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Self adaptive solution strategies: Locally bound constrained Newton Raphson solution algorithms

A summary is given of strategies which enable the automatic adjustment of the constraint surfaces recently used to extend the range and numerical stability/efficiency of nonlinear finite element equation solvers. In addition to handling kinematic and material induced nonlinearity, both pre-and postbuckling behavior can be treated. The scheme employs localized bounds on various hierarchical partitions of the field variables. These are used to resize, shape, and orient the global constraint surface, thereby enabling essentially automatic load/deflection incrementation. Due to the generality of the approach taken, it can be implemented in conjunction with the constraints of an arbitrary functional type. To benchmark the method, several numerical experiments are presented. These include problems involving kinematic and material nonlinearity, as well as pre- and postbuckling characteristics. Also included is a list of papers published in the course of the work.

Padovan, Joe↗

NASA NTRS1986

Locally bound constrained Newton-Raphson solution algorithms

This paper develops strategies which enable the automatic adjustment of the constraint surfaces recently used to extend the range and numerical stability/efficiency of nonlinear finite-element equation solvers. In addition to handling kinematic and material induced nonlinearity, both pre- and postbuckling behavior can be treated. The scheme developed employs localized bounds on various hierarchical partitions of the field variables. These are used to resize, shape, and orient the global constraint surface, thereby enabling essentially automatic load/deflection incrementation. Due to the generality of the approach taken, it can be implemented in conjunction with constraints of arbitrary functional type. To benchmark the method, several numerical experiments are presented. These include problems involving kinematic and material nonlinearity, as well as, pre- and postbuckling characteristics.

Padovan, J.↗

NASA NTRS1998

A Posteriori Bounds for Linear-Functional Outputs of Crouzeix-Raviart Finite Element Discretizations of the Incompressible Stokes Problem

We present a finite element technique for the efficient generation of lower and upper bounds to outputs which are linear functionals of the solutions to the incompressible Stokes equations in two space dimensions; the finite element discretization is effected by Crouzeix-Raviart elements, the discontinuous pressure approximation of which is central to our approach. The bounds are based upon the construction of an augmented Lagrangian: the objective is a quadratic "energy" reformulation of the desired output; the constraints are the finite element equilibrium equations (including the incompressibility constraint), and the intersubdomain continuity conditions on velocity. Appeal to the dual max-min problem for appropriately chosen candidate Lagrange multipliers then yields inexpensive bounds for the output associated with a fine-mesh discretization; the Lagrange multipliers are generated by exploiting an associated coarse-mesh approximation. In addition to the requisite coarse-mesh calculations, the bound technique requires solution only of local subdomain Stokes problems on the fine-mesh. The method is illustrated for the Stokes equations, in which the outputs of interest are the flowrate past, and the lift force on, a body immersed in a channel.

Patera, Anthony T.↗

NASA NTRS1995

An hp-adaptivity and error estimation for hyperbolic conservation laws

This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Bey, Kim S.↗

NASA NTRS1989

A survey on the structured singular value

The structured singular value, U, is an important linear algebra tool to study a class of matrix perturbation problems. It is useful for analyzing the robustness of stability and performance of uncertain, (nominally) linear systems. Computation of (M) is difficult, and usually, upper and lower bounds are all that can be reliably computed. Upper bounds give conservative estimates of the sizes of allowable perturbations. The maximum singular value of a matrix M is an upper bound for (M). As an upper bound, it can be improved by finding a transformations to the data (i.e. M) which do not change the structured singular value, but do reduce the maximum singular value. Typically, upper bound algorithms involve searches over sets of transformations to yield the tightest bound. Lower bound algorithms are intelligent searches for minimum-norm solutions to multivariable polynomial equations, and are based on various optimality conditions that hold at the global (and, unfortunately, some local) minima. The current methods to compute both of these types of bounds are reviewed. Theoretical justification and extensive numerical experience with the various algorithms are covered.

Packard, Andy↗

NASA NTRS1993

An integrated collision prediction and avoidance scheme for mobile robots in non-stationary environments

A formulation that makes possible the integration of collision prediction and avoidance stages for mobile robots moving in general terrains containing moving obstacles is presented. A dynamic model of the mobile robot and the dynamic constraints are derived. Collision avoidance is guaranteed if the distance between the robot and a moving obstacle is nonzero. A nominal trajectory is assumed to be known from off-line planning. The main idea is to change the velocity along the nominal trajectory so that collisions are avoided. A feedback control is developed and local asymptotic stability is proved if the velocity of the moving obstacle is bounded. Furthermore, a solution to the problem of inverse dynamics for the mobile robot is given. Simulation results verify the value of the proposed strategy.

Kyriakopoulos, K. J.↗

NASA NTRS1973

Relaxation factors for supercritical flows.

Relaxation procedures for solution of steady supercritical transonic flows are investigated. Von Neumann (Fourier-mode) stability analysis is used to find bounds of relaxation factors. The bounds depend on local Mach number and local mesh aspect ratio. Long wave instability of Murman-Cole implicit method is indicated. Two new relaxation procedures are introduced. Both employ central differencing exclusively. Group velocities of Fourier modes are used to study signal propagation. It was found necessary to avoid or to damp out signals propagating upstream in supersonic zones in order to obtain physically meaningful transonic solutions. It appears that requirements of high rate of convergence, of stability and accuracy, are in conflict, and that a combination of relaxation methods must be used in order not to compromise the requirements.

Kentzer, C. P.↗

NASA NTRS1987

A note on the generation of Tollmien-Schlichting waves by sudden surface-curvature change

This note is primarily concerned with the generation of spatially growing Tollmien-Schlichting waves by the interaction of very long-wavelength free-stream disturbances with a discontinuity in the curvature of a bounding surface (whose slope may or may not be continuous). The theory is combined with a numerical solution of the local Orr-Sommerfeld equation, and the result is used to predict the Tollmien-Schlichting amplitude in a relevant experiment carried out by Leehey and Shapiro (1980). The calculated results are in satisfactory agreement with their observations.

Goldstein, M. E.↗

NASA NTRS1977

Test particle propagation in magnetostatic turbulence. 3: The approach to equilibrium

The asymptotic behavior, for large time, of the quasi-linear diabatic solutions and their local approximations is considered. A time averaging procedure is introduced which yields the averages of these solutions over time intervals which contain only large time values. A discussion of the quasi-linear diabatic solutions which is limited to those solutions that are bounded from below as functions of time is given. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity the time averaged quasi-linear diabatic solutions must approach isotropy (mu-independence). The first derivative with respect to mu of these solutions is also considered. This discussion is limited to first derivatives which are bounded functions of time. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity, the time averaged first derivative must approach zero everywhere in mu except at mu = 0 where it must approach a large value which is calculated. The impact of this large derivative on the quasi-linear expansion scheme is discussed. An H-theorem for the first local approximation to the quasi-linear diabatic solutions is constructed. Without time averaging, the H-theorem is used to determine sufficient conditions for the first local approximate solutions to asymptote, with increasing time, to exactly the same final state which the time averaged quasi-linear diabatic solutions must approach as discussed above.

Klimas, A. J.↗

NASA NTRS1990

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

NASA NTRS1991

On computing the global time-optimal motions of robotic manipulators in the presence of obstacles

A method for computing the time-optimal motions of robotic manipulators is presented that considers the nonlinear manipulator dynamics, actuator constraints, joint limits, and obstacles. The optimization problem is reduced to a search for the time-optimal path in the n-dimensional position space. A small set of near-optimal paths is first efficiently selected from a grid, using a branch and bound search and a series of lower bound estimates on the traveling time along a given path. These paths are further optimized with a local path optimization to yield the global optimal solution. Obstacles are considered by eliminating the collision points from the tessellated space and by adding a penalty function to the motion time in the local optimization. The computational efficiency of the method stems from the reduced dimensionality of the searched spaced and from combining the grid search with a local optimization. The method is demonstrated in several examples for two- and six-degree-of-freedom manipulators with obstacles.

Shiller, Zvi↗

NASA NTRS1998

Fluorescence Studies of Lysozyme Nucleation

Fluorescence is one of the most powerful tools available for the study of macromolecules. For example, fluorescence can be used to study self association through methods such as anisotropy (the rotational rate of the molecule in solution), quenching (the accessibility of a bound probe to the bulk solution), and resonance energy transfer (measurement of the distance between two species). Fluorescence can also be used to study the local environment of the probe molecules, and the changes in that environment which accompany crystal nucleation and growth. However fluorescent techniques have been very much underutilized in macromolecular growth studies. One major advantage is that the fluorescent species generally must be at low concentration, typically ca 10-5 to 10-6 M. Thus one can study a very wide range of solution conditions, ranging from very high to very low protein concentration, he latter of which are not readily accessible to scattering techniques. We have prepared a number of fluorescent derivatives of chicken egg white lysozyme (CEWL). Fluorescent probes have been attached to two different sites, ASP 101 and the N-terrninal amine, with a sought for use in different lines of study. Preliminary resonance energy transfer studies have been -carried out using pyrene acetic acid (Ex 340 mn, Em 376 nm) lysozyme as a donor and cascade blue (Ex 377 run, Em 423 nm) labeled lysozyme as an acceptor. The emission of both the pyrene and cascade blue probes was followed as a function of the salt protein concentrations. The data show an increase in cascade blue and a concomitant decrease in the pyrene fluorescence as either the salt or protein concentrations are increased, suggesting that the two species are approaching each other close enough for resonance energy transfer to occur. This data can be analyzed to measure the distance between the probe molecules and, knowing their locations on the protein molecule their distances from and orientations with respect to each other. The results of these and other studies will be discussed.

Pusey, Marc L.↗

NASA NTRS2019

Using Gaia Data 2 to Constrain Local Dark Matter Density and Thin Dark Disk

We use stellar kinematics from the latest Gaia data release (DR2) to measure the local dark matter (DM) density ρDM in a heliocentric cylinder of radius R = 150 pc and half-height z = 200 pc. We also explore the prospect of using our analysis to estimate the DM density in local substructure by setting constraints on the surface density and scale height of a thin dark disk aligned with the baryonic disk and formed due to dissipative dark matter self-interactions. Performing the statistical analysis within a Bayesian framework for three types of tracers, we obtain ρDM = 0.016 ± 0.010 Mꙩ/p c 3 for A stars; early G stars give a similar result, while F stars yield a significantly higher value. For a thin dark disk, A stars set the strongest constraint: excluding surface densities (5–12) Mꙩ/pc 2 for scale heights below 100 pc with 95% confidence. The upper bound of this constraint implies . 1% of the Milky Way DM mass is present in a dissipative dark sector. Comparing our results with those derived using Tycho-Gaia Astrometric Solution (TGAS) data, we find that the uncertainty in our measurements of the local DM content is dominated by systematic errors that arise from assumptions of our dynamical analysis in the low z region. Furthermore, there will only be a marginal reduction in these uncertainties with more data in the Gaia era. We comment on the robustness of our method and discuss potential improvements for future work.

dark matter theory↗

NASA NTRS2004

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗

NASA NTRS1991

An approximation based global optimization strategy for structural synthesis

A global optimization strategy for structural synthesis based on approximation concepts is presented. The methodology involves the solution of a sequence of highly accurate approximate problems using a global optimization algorithm. The global optimization algorithm implemented consists of a branch and bound strategy based on the interval evaluation of the objective function and constraint functions, combined with a local feasible directions algorithm. The approximate design optimization problems are constructed using first order approximations of selected intermediate response quantities in terms of intermediate design variables. Some numerical results for example problems are presented to illustrate the efficacy of the design procedure setforth.

Sepulveda, A. E.↗

NASA NTRS1994

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The steady state solution of the system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been obtained using a multigrid strategy of unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time-stepping scheme with a stability-bound local time step, while turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positivity. Low-Reynolds-number modifications to the original two-equation model are incorporated in a manner which results in well-behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved, initializing all quantities with uniform freestream values. Rapid and uniform convergence rates for the flow and turbulence equations are observed.

MULTIGRID SOLUTION↗

NASA NTRS1993

The development of a mixing layer under the action of weak streamwise vortices

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Goldstein, M. E.↗

NASA NTRS1993

The development of a mixing layer under the action of weak streamwise vortices

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds-number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon(sup -1)), where epsilon is much less than 1 is the cross flow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by the appearance of a singularity. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the singular location.

Goldstein, Marvin E.↗