Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complexity bounding”

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 325 records · Page 18

Martian Polar Impact Craters: A Preliminary Assessment Using Mars Orbiter Laser Altimeter (MOLA)

Our knowledge of the age of the layered polar deposits and their activity in the volatile cycling and climate history of Mars is based to a large extent on their apparent ages as determined from crater counts. Interpretation of the polar stratigraphy (in terms of climate change) is complicated by reported differences in the ages of the northern and southern layered deposits. The north polar residual ice deposits are thought to be relatively young, based on the reported lack of any fresh impact craters in Viking Orbiter Images. Herkenhoff et al., report no craters at all on the North polar layered deposits or ice cap, and placed an upper bound on the surface age (or, alternatively, the vertical resurfacing rate) of 100 thousand years to 10 million years, suggesting that the north polar region is an active resurfacing site. In contrast, the southern polar region was found to have at least 15 impact craters in the layered deposits and cap. Plaut et al, concluded that the surface was less than or = 120 million years old. This reported age difference factor of 100 to 1000 increases complexity in climate and volatile modeling. Recent MOLA results for the topography of the northern polar cap document a handful or more of possible craters, which could result in revised age or resurfacing estimates for the northern cap. This study is a preliminary look at putative craters in both polar caps. Additional information is contained in the original extended abstract.

Sakimoto, S. E. H.↗

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗

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↗

Formal Safety Certification of Aerospace Software

In principle, formal methods offer many advantages for aerospace software development: they can help to achieve ultra-high reliability, and they can be used to provide evidence of the reliability claims which can then be subjected to external scrutiny. However, despite years of research and many advances in the underlying formalisms of specification, semantics, and logic, formal methods are not much used in practice. In our opinion this is related to three major shortcomings. First, the application of formal methods is still expensive because they are labor- and knowledge-intensive. Second, they are difficult to scale up to complex systems because they are based on deep mathematical insights about the behavior of the systems (t.e., they rely on the "heroic proof"). Third, the proofs can be difficult to interpret, and typically stand in isolation from the original code. In this paper, we describe a tool for formally demonstrating safety-relevant aspects of aerospace software, which largely circumvents these problems. We focus on safely properties because it has been observed that safety violations such as out-of-bounds memory accesses or use of uninitialized variables constitute the majority of the errors found in the aerospace domain. In our approach, safety means that the program will not violate a set of rules that can range for the simple memory access rules to high-level flight rules. These different safety properties are formalized as different safety policies in Hoare logic, which are then used by a verification condition generator along with the code and logical annotations in order to derive formal safety conditions; these are then proven using an automated theorem prover. Our certification system is currently integrated into a model-based code generation toolset that generates the annotations together with the code. However, this automated formal certification technology is not exclusively constrained to our code generator and could, in principle, also be integrated with other code generators such as RealTime Workshop or even applied to legacy code. Our approach circumvents the historical problems with formal methods by increasing the degree of automation on all levels. The restriction to safety policies (as opposed to arbitrary functional behavior) results in simpler proof problems that can generally be solved by fully automatic theorem proves. An automated linking mechanism between the safety conditions and the code provides some of the traceability mandated by process standards such as DO-178B. An automated explanation mechanism uses semantic markup added by the verification condition generator to produce natural-language explanations of the safety conditions and thus supports their interpretation in relation to the code. It shows an automatically generated certification browser that lets users inspect the (generated) code along with the safety conditions (including textual explanations), and uses hyperlinks to automate tracing between the two levels. Here, the explanations reflect the logical structure of the safety obligation but the mechanism can in principle be customized using different sets of domain concepts. The interface also provides some limited control over the certification process itself. Our long-term goal is a seamless integration of certification, code generation, and manual coding that results in a "certified pipeline" in which specifications are automatically transformed into executable code, together with the supporting artifacts necessary for achieving and demonstrating the high level of assurance needed in the aerospace domain.

Denney, Ewen↗

Standards for Radiation Effects Testing: Ensuring Scientific Rigor in the Face of Budget Realities and Modern Device Challenges

An overview is presented of the space radiation environment and its effects on electrical, electronic, and electromechanical parts. Relevant test standards and guidelines are listed. Test standards and guidelines are necessary to ensure best practices, minimize and bound systematic and random errors, and to ensure comparable results from different testers and vendors. Test standards are by their nature static but exist in a dynamic environment of advancing technology and radiation effects research. New technologies, failure mechanisms, and advancement in our understanding of known failure mechanisms drive the revision or development of test standards. Changes to standards must be weighed against their impact on cost and existing part qualifications. There must be consensus on new best practices. The complexity of some new technologies exceeds the scope of existing test standards and may require development of a guideline specific to the technology. Examples are given to illuminate the value and limitations of key radiation test standards as well as the challenges in keeping these standards up to date.

Test Standard↗

Spectral Derivatives of Optical Depth for Partitioning Aerosol Type and Loading

Quantifying aerosol compositions (e.g., type, loading) from remotely sensed measurements by spaceborne, suborbital and ground-based platforms is a challenging task. In this study, the first and second-order spectral derivatives of aerosol optical depth (AOD) with respect to wavelength are explored to determine the partitions of the major components of aerosols based on the spectral dependence of their particle optical size and complex refractive index. With theoretical simulations from the Second Simulation of a Satellite Signal in the Solar Spectrum (6S) model, AOD spectral derivatives are characterized for collective models of aerosol types, such as mineral dust (DS) particles, biomass-burning (BB) aerosols and anthropogenic pollutants (AP), as well as stretching out to the mixtures among them. Based on the intrinsic values from normalized spectral derivatives, referenced as the Normalized Derivative Aerosol Index (NDAI), a unique pattern is clearly exhibited for bounding the major aerosol components; in turn, fractions of the total AOD (fAOD) for major aerosol components can be extracted. The subtlety of this NDAI method is examined by using measurements of typical aerosol cases identified carefully by the ground-based Aerosol Robotic Network (AERONET) sun–sky spectroradiometer. The results may be highly practicable for quantifying fAOD among mixed-type aerosols by means of the normalized AOD spectral derivatives.

aerosol partition↗

The Sensitivity of the Equatorial Pacific ODZ to Particulate Organic Matter Remineralization in A Climate Model Under Pre-Industrial Conditions

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ O2 ol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Oxygen-deficient zone↗

Evaluation Applied to Reliability Analysis of Reconfigurable, Highly Reliable, Fault-Tolerant, Computing Systems for Avionics

Emulation techniques are proposed as a solution to a difficulty arising in the analysis of the reliability of highly reliable computer systems for future commercial aircraft. The difficulty, viz., the lack of credible precision in reliability estimates obtained by analytical modeling techniques are established. The difficulty is shown to be an unavoidable consequence of: (1) a high reliability requirement so demanding as to make system evaluation by use testing infeasible, (2) a complex system design technique, fault tolerance, (3) system reliability dominated by errors due to flaws in the system definition, and (4) elaborate analytical modeling techniques whose precision outputs are quite sensitive to errors of approximation in their input data. The technique of emulation is described, indicating how its input is a simple description of the logical structure of a system and its output is the consequent behavior. The use of emulation techniques is discussed for pseudo-testing systems to evaluate bounds on the parameter values needed for the analytical techniques.

Migneault, G. E.↗

Emulation applied to reliability analysis of reconfigurable, highly reliable, fault-tolerant computing systems

Emulation techniques applied to the analysis of the reliability of highly reliable computer systems for future commercial aircraft are described. The lack of credible precision in reliability estimates obtained by analytical modeling techniques is first established. The difficulty is shown to be an unavoidable consequence of: (1) a high reliability requirement so demanding as to make system evaluation by use testing infeasible; (2) a complex system design technique, fault tolerance; (3) system reliability dominated by errors due to flaws in the system definition; and (4) elaborate analytical modeling techniques whose precision outputs are quite sensitive to errors of approximation in their input data. Next, the technique of emulation is described, indicating how its input is a simple description of the logical structure of a system and its output is the consequent behavior. Use of emulation techniques is discussed for pseudo-testing systems to evaluate bounds on the parameter values needed for the analytical techniques. Finally an illustrative example is presented to demonstrate from actual use the promise of the proposed application of emulation.

Migneault, G. E.↗

Generalized Orbital Projections of a Sublimating Ice Particle

The issue of orbital debris resulting from human activities in space is a growing concern for the space users' community. Waste generated in Low Earth Orbit (LEO) can stay in orbit for a long time, creating significant hazards for other spacecraft flying at lower intercepting orbits. Many spacecraft, especially crewed vehicles, are required to vent fluids into space. These fluids include propellant, wastewater, excess condensate, and others. It is important to analyze the behavior of particles that result from these activities, since each individual particle is capable of damaging or destroying a spacecraft in a lower, crossing orbit, and such particles are invisible to tracking radar systems on the ground. The deorbit trajectory of an ice particle is complex. It depends on factors including attitude of the vehicle during vent, initial velocities of particles, altitude at which the vent occurred, and numerous evaporation and sublimation factors. These include contamination within the vented water, evolution of bubbles within the clear water, and sun flux factors such as time of the year and current beta angle. The purpose of this study is to examine the influences of these factors on the trajectories of ice particles resulting from condensate water dumps, and to bound the safe trajectories of spacecraft that lie below the venting spacecraft.

Menkin, Evgeny↗

On the masses of giant molecular cloud complexes

A method of mass estimation for molecular clouds is presented which is based on approximate balance in the outer cloud layers between the cloud's gravitation, the galactic tide, and internal pressure. The largest observed clouds, which have greatest linear extents of 100 pc, are found to have masses of at least 200,000 solar masses. The cloud masses cannot exceed this lower limit by more than a factor of 3, or the velocity distributions of disk stars would be more relaxed than is actually observed. This implied upper limit to cloud masses combined with the galactic tide may be related to the absence of clouds at galactocentric radii less than 4 kpc. If Sagittarius B2 is bound, its mass must be more than 50 million solar masses.

Stark, A. A.↗

Transmission-zero bounds for large space structures, with applications

Many large space structure control problems lead quite naturally to the application of an optimal regulator, so the transmission zeros of the open-loop system give fundamental information about the speed of response achievable by the closed-loop system. Despite the importance of this and other well-known zeros properties, little attention has been given to the transmission zeros of large space structures, except for the special case of a rigid spacecraft with flexible appendages. The object of this paper is to remedy this deficiency. In particular, it is proved that the zeros of a structure with colocated sensors and actuators must lie in a region of the complex plane that is defined by its natural frequencies and damping ratios. This generic result, a consequence of the special form of the equations of motion of structural dynamics, admits a very simple graphical interpretation: it is the generalization of the classical pole-zero interlacing property of undamped single-input/single-output structures. The number of sensor/actuator pairs, and their locations, specify where in the permissible region transmission zeros actually lie, thus quantifying the effect of sensor/actuator placement on closed-loop system performance. These points are illustrated by simple examples.

Williams, Trevor↗

Experimental study of adaptive pointing and tracking for large flexible space structures

This paper describes an experimental study of adaptive pointing and tracking control for flexible spacecraft conducted on a complex ground experiment facility. The algorithm used in this study is based on a multivariable direct model reference adaptive control law. Several experimental validation studies were performed earlier using this algorithm for vibration damping and robust regulation, with excellent results. The current work extends previous studies by addressing the pointing and tracking problem. As is consistent with an adaptive control framework, the plant is assumed to be poorly known to the extent that only system level knowledge of its dynamics is available. Explicit bounds on the steady-state pointing error are derived as functions of the adaptive controller design parameters. It is shown that good tracking performance can be achieved in an experimental setting by adjusting adaptive controller design weightings according to the guidelines indicated by the analytical expressions for the error.

Boussalis, D.↗

Standard test case for a low-speed, turbulent junction vortex flow

The mean flow structure upstream, around within, and in the near wake of a turbulent junction or horseshoe vortex is reported for an incompressible subsonic flow. Measurements of the primitive variables of velocity and pressure are reported on all surfaces bounding the region of the vortex flow and on three transverse and one streamwise plane within the flowfield itself for comparisons between the measured and any calculated flow variables. Detailed surface flow visualizations and some direct force measurements of surface shear stress are also available. The data is a highly detailed, coherent, self-consistent set offered to the computational fluid mechanics community as a standard test case for the evaluation of the capability of numerical solvers intended for predicting the flowfield in such a complex, separated, three-dimensional turbulent flow. This data base is available for copying to user supplied tapes or for transmission via BITNET, as well as in two National Technical Information Service (NTIS) reports.

Pierce, F. J.↗

Extraction-Separation Performance and Dynamic Modeling of Orion Test Vehicles with Adams Simulation: 2nd Edition

NASA's Orion Capsule Parachute Assembly System (CPAS) project has advanced into the third generation of its parachute test campaign and requires technically comprehensive modeling capabilities to simulate multi-body dynamics (MBD) of test articles released from a C-17. Safely extracting a 30,000 lbm mated test article from a C-17 and performing stable mid-air separation maneuvers requires an understanding of the interaction between elements in the test configuration and how they are influenced by extraction parachute performance, aircraft dynamics, aerodynamics, separation dynamics, and kinetic energy experienced by the system. During the real-time extraction and deployment sequences, these influences can be highly unsteady and difficult to bound. An avionics logic window based on time, pitch, and pitch rate is used to account for these effects and target a favorable separation state in real time. The Adams simulation has been employed to fine-tune this window, as well as predict and reconstruct the coupled dynamics of the Parachute Test Vehicle (PTV) and Cradle Platform Separation System (CPSS) from aircraft extraction through the mid-air separation event. The test-technique for the extraction of CPAS test articles has evolved with increased complexity and requires new modeling concepts to ensure the test article is delivered to a stable test condition for the programmer phase. Prompted by unexpected dynamics and hardware malfunctions in drop tests, these modeling improvements provide a more accurate loads prediction by incorporating a spring-damper line-model derived from the material properties. The qualification phase of CPAS testing is on the horizon and modeling increasingly complex test-techniques with Adams is vital to successfully qualify the Orion parachute system for human spaceflight.

Fraire, Usbaldo, Jr.↗

A Mixed Integer Efficient Global Optimization Algorithm with Multiple Infill Strategy - Applied to a Wing Topology Optimization Problem

With the advancement in high performance computing and numerical optimization techniques,engineering design optimization problems are becoming more complex, larger scale,higher fidelity, and computationally more demanding, requiring longer run times than ever before. There exists methodologies and techniques that can address some of these challenges but very few can address all, and most are limited in the extent that these concerns can be addressed. With the goal of addressing such challenging engineering problems, we developed anew optimization framework, named AMIEGO, that combines concepts from surrogate-based optimization approaches, gradient-based numerical methods, Partial Least Squares, evolutionary algorithms, and Branch-and-Bound, providing newer capabilities that were not previouslyperceived. However, the original version of this framework, in the process of adaptive samplingto explore and exploit the design space, finds only a single sample point per iteration. The efforthere builds upon this previously developed optimization framework to include multiple infillsampling capability that combines the concept of generalized expected improvement function,unsupervised learning, and multi-objective evolutionary technique. To demonstrate, AMIEGOwith the multiple infill capability (called AMIEGO-MIMOS) solves a series of increasingly difficultengineering design optimization problems. The results reveal the performance of the newapproach is problem dependent. When applied to a ten-bar truss problem, the newly proposedmultiple infill strategy consistently leads to a better design solutions when compared to theexisting CPTV method (implemented with the context of the AMIEGO framework). On theother hand, when applied to a mixed-integer high fidelity wing topology optimization problem- MIMOS, despite showing a steeper convergence at the start, eventually leads to an inferiorsolution as compared to CPTV approach. These results also reveal that a small number ofstarting points, in general, are sufficient to lead to a good overall solution.

Mixed-integer optimization↗

Search for resonances in positron-atom systems

No one likes to see a scattering cross section curve that is too smooth; it is much more interesting to find bumps and wiggles and most interesting if it is possible to understand their cause. Several types of resonances were clearly established in positron-containing systems: those lying just below a degenerate threshold (like 2s-2p in hydrogenic atoms or ions) and those representing Coulomb bound states in a rearranged channel (like Ps + H reversing positron +H(-)). Recently, two new sorts of resonances have been reported for which the resonant mechanism is not clear. The first is a very low lying resonance in the e-Ps system (obtained by an adiabatic expansion method), and the second is a similarly low lying two channel resonance in the positron -H system (obtained by a close-coupling technique.) These developments encouraged the examination of such systems using the standard methods of stabilization and complex rotation. Most of the results are negative; the low lying resonances in either system are verified. Some indication of new resonances in the e(+)-He(+) system is found; this may be caused by the attraction between Ps in the n=2 state and the He(++) nucleus.

Bhatia, A. K.↗

Lunar Pole Illumination and Communications Maps Computed from GSSR Elevation Data

A Digital Elevation Model of the lunar south pole was produced using Goldstone Solar System RADAR (GSSR) data obtained in 2006.12 This model has 40-meter horizontal resolution and about 5-meter relative vertical accuracy. This Digital Elevation Model was used to compute average solar illumination and Earth visibility with 100 kilometers of the lunar south pole. The elevation data were converted into local terrain horizon masks, then converted into lunar-centric latitude and longitude coordinates. The horizon masks were compared to latitude, longitude regions bounding the maximum Sun and Earth motions relative to the moon. Estimates of Earth visibility were computed by integrating the area of the region bounding the Earth's motion that was below the horizon mask. Solar illumination and other metrics were computed similarly. Proposed lunar south pole base sites were examined in detail, with the best site showing yearly solar power availability of 92 percent and Direct-To-Earth (DTE) communication availability of about 50 percent. Similar analysis of the lunar south pole used an older GSSR Digital Elevation Model with 600-meter horizontal resolution. The paper also explores using a heliostat to reduce the photovoltaic power system mass and complexity.

lunar base↗