SEARCH · Search NASA
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.
Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.
Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.
The NASA Langley Research Center in 1984 initiated a rotorcraft structural dynamics program, designated DAMVIBS (Design Analysis Methods for VIBrationS), with the objective of establishing the technology base needed by the rotorcraft industry for developing an advanced finite-element-based dynamics design analysis capability for vibrations. An assessment of the program showed that the DAMVIBS Program has resulted in notable technical achievements and major changes in industrial design practice, all of which have significantly advanced the industry's capability to use and rely on finite-element-based dynamics analyses during the design process.
The motions of throw away detectors (TAD's) are analyzed using the linearized relative motion equations. The TAD's are to be used in the amps program as diagnostic instruments for making various measurements near the shuttle. The TAD's are ejected from the shuttle in arbitrary directions with small relative velocities (0.1 to 1.0 m/s) their subsequent trajectories relative to the shuttle are analyzed. Initial conditions that are likely to result in recontact between the TAD and the shuttle are identified. The sensitivity of the motion to variations in the initial conditions, possibly resulting from inaccuracy in the ejection mechanism, are analyzed as are effects of atmospheric drag. A targeting method, a method of giving the TAD correct initial conditions such that it will pass through a given point relative to the shuttle at a given time, is developed. The results of many specific cases are presented in graphical form.
An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.
As computer power continues to increase, mesoscale models are initialized at all hours of the day and continue to be run at higher and higher spatial resolutions. As a result, initializing land surface temperature can be problematic. The majority of research-based models are initialized at 00 and 12 UTC when upper air observations and reanalysis fields are available. The landsea surface temperatures are then set equal to the two-meter air temperature produced by the preprocessor analysis system. This particular procedure might be valid in the early morning hours just prior to sunrise, but it becomes less valid during the remainder of the diurnal cycle. Operational models, such as the Rapid Update Cycle run at the National Centers for Environmental Prediction (NCEP), are initialized every hour on a daily basis. This presents a unique challenge to the initial specification of the land temperature, especially during the first several hours of the solar heating cycle when land and overlying air temperatures are far from being the same. Another issue that needs to be addressed is the spatial variability of land surface temperature. By early next year, the NCEP operational Eta model will be employed at 8 km resolution. Methods to accurately specify the initial land surface temperature at such high resolution need to be explored. This paper presents the results of using data from the NASA Moderate Imaging Sensor aboard the TERRA Satellite to initialize land and sea surface temperatures within the Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) 5'th generation Mesoscale Model (MM5). We have simulated a northern Gulf Coast sea breeze case to demonstrate the utility of using the MODIS data to initialize both the land and sea surface temperature fields. Model grid configurations of 12-, 4-, and l-km are employed.
Characterization of the three-dimensional structure of solar transients using incomplete plane of sky data is a difficult problem whose solutions have potential for societal benefit in terms of space weather applications. In this paper transients are characterized in three dimensions by means of conic coronal mass ejection (CME) approximation. A novel method for the automatic determination of cone model parameters from observed halo CMEs is introduced. The method uses both standard image processing techniques to extract the CME mass from white-light coronagraph images and a novel inversion routine providing the final cone parameters. A bootstrap technique is used to provide model parameter distributions. When combined with heliospheric modeling, the cone model parameter distributions will provide direct means for ensemble predictions of transient propagation in the heliosphere. An initial validation of the automatic method is carried by comparison to manually determined cone model parameters. It is shown using 14 halo CME events that there is reasonable agreement, especially between the heliocentric locations of the cones derived with the two methods. It is argued that both the heliocentric locations and the opening half-angles of the automatically determined cones may be more realistic than those obtained from the manual analysis
In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.
Anion exchange membranes (AEMs) offer a cost-effective alternative to proton exchange membranes as alkaline fuel cells and electrolyzers permit the use of non-platinum group electrodes and components. Despite continued progress, the operational lifetime and stability of these membranes limit the widespread adoption of AEM-based electrochemical technologies. This study presents a flexible and easily implemented ultraviolet (UV)-initiated nitrene-based cross-linking method which uses a small, facile organic azide precursor. As a proof of concept, we demonstrate this approach on the well-studied poly(2,6-dimethyl-1,4-phenylene oxide) quaternary ammonium AEM (QPPO) polyelectrolyte. Further, a survey of cross-linker density (2.5–10 mol %) found that the addition of 10 mol % results in a 59% reduction in water uptake, a 58.8% decrease in the swell ratio, and a 31% increase in tensile strength vs the un-cross-linked material. Nitrene cross-linking also enhanced the membrane’s durability, enabling 1000 h of stable performance under electrochemical load. This UV-initiated cross-linking method may be easily integrated into production processes, allowing chemical cross-linking at any stage, including posthydration.
Acoustic tomography (AT) is explored as a remote sensing technique to obtain instantaneous snapshots of temperature and velocity fluctuations for wind energy applications. This study integrates Large Eddy Simulation (LES) with the Stochastic Inversion (SI) method to validate the algorithm’s capacity for accurate reconstruction of atmospheric fluctuations. The initial findings demonstrate the efficacy of the method in accurately capturing the predominant flow structures. Normalized L2 error evaluations further inform the algorithm’s precision, with errors accentuated in less sampled peripheral regions. The results underscore the method’s promise as a non-intrusive observational tool, with ongoing development poised to improve its precision and reliability.
The problem of calculating the stability of steady state solutions of differential equations is addressed. Leading eigenvalues of large matrices that arise from discretization are calculated, and an efficient multigrid method for solving these problems is presented. The resulting grid functions are used as initial approximations for appropriate eigenvalue problems. The method employs local relaxation on all levels together with a global change on the coarsest level only, which is designed to separate the different eigenfunctions as well as to update their corresponding eigenvalues. Coarsening is done using the FAS formulation in a nonstandard way in which the right-hand side of the coarse grid equations involves unknown parameters to be solved on the coarse grid. This leads to a new multigrid method for calculating the eigenvalues of symmetric problems. Numerical experiments with a model problem are presented which demonstrate the effectiveness of the method.
Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.
Fault detection, diagnosis, and prognosis are essential tasks in the operation of autonomous spacecraft, instruments, and in situ platforms. One of NASA s key mission requirements is robust state estimation. Sensing, using a wide range of sensors and sensor fusion approaches, plays a central role in robust state estimation, and there is a need to diagnose sensor failure as well as component failure. Sensor validation can be considered to be part of the larger effort of improving reliability and safety. The standard methods for solving the sensor validation problem are based on probabilistic analysis of the system, from which the method based on Bayesian networks is most popular. Therefore, these methods can only predict the most probable faulty sensors, which are subject to the initial probabilities defined for the failures. The method developed in this work is based on a model-based approach and provides the faulty sensors (if any), which can be logically inferred from the model of the system and the sensor readings (observations). The method is also more suitable for the systems when it is hard, or even impossible, to find the probability functions of the system. The method starts by a new mathematical description of the problem and develops a very efficient and systematic algorithm for its solution. The method builds on the concepts of analytical redundant relations (ARRs).
OBJECTIVE: A new computerized semi-automatic method for left ventricular (LV) chamber segmentation is presented. METHODS: The LV is imaged by real-time three-dimensional echocardiography (RT3DE). The surface detection model, based on level set techniques, is applied to RT3DE data for image analysis. The modified level set partial differential equation we use is solved by applying numerical methods for conservation laws. The initial conditions are manually established on some slices of the entire volume. The solution obtained for each slice is a contour line corresponding with the boundary between LV cavity and LV endocardium. RESULTS: The mathematical model has been applied to sequences of frames of human hearts (volume range: 34-109 ml) imaged by 2D and reconstructed off-line and RT3DE data. Volume estimation obtained by this new semi-automatic method shows an excellent correlation with those obtained by manual tracing (r = 0.992). Dynamic change of LV volume during the cardiac cycle is also obtained. CONCLUSION: The volume estimation method is accurate; edge based segmentation, image completion and volume reconstruction can be accomplished. The visualization technique also allows to navigate into the reconstructed volume and to display any section of the volume.
Intake flow for a dual intake valved diesel engine is modeled using moving valves and realistic geometries. The objectives are to obtain accurate initial conditions for combustion calculations and to provide a tool for studying intake processes. Global simulation parameters are compared with experimental results and show good agreement. The intake process shows a 30 percent difference in mass flows and average swirl in opposite directions across the two intake valves. The effect of the intake process on the flow field at the end of compression is examined. Modeling the intake flow results in swirl and turbulence characteristics that are quite different from those obtained by conventional methods in which compression stroke initial conditions are assumed.
An approach was developed for incorporating flow and geometric constraints into the Direct Iterative Surface Curvature (DISC) design method. In this approach, an initial target pressure distribution is developed using a set of control points. The chordwise locations and pressure levels of these points are initially estimated either from empirical relationships and observed characteristics of pressure distributions for a given class of airfoils or by fitting the points to an existing pressure distribution. These values are then automatically adjusted during the design process to satisfy the flow and geometric constraints. The flow constraints currently available are lift, wave drag, pitching moment, pressure gradient, and local pressure levels. The geometric constraint options include maximum thickness, local thickness, leading-edge radius, and a 'glove' constraint involving inner and outer bounding surfaces. This design method was also extended to include the successive constraint release (SCR) approach to constrained minimization.
Visual inspection (VI) of composite structures provides an effective, wide field survey to ensure design and material compliance is maintained for the entire service life of the manufactured component. Visual inspection is one of the most commonly used nondestructive inspection methods to assess surface defects on composites. By applying visual inspection elements, mechanical damage that could affect component residual strength can be identified and a disposition reached before a potentially catastrophic event occurs. This technique is non-contact and applied at all stages of the composite structure s processing and use. Additionally, VI is required to be performed in all service environments until decommission. By following Level I and Level II damage accept/reject criteria set forth in this document, the material review (MR) process can be initiated. This nondestructive evaluation (NDE) method should be complemented with additional NDE methods to best understand the nature of the observed indications. When VI identification is followed by interrogation using effective complimentary NDE methods, the final Material Review Board (MRB) disposition of the component can be effectively achieved. Although elements of the visual inspection method are discussed and expressly required by the various range standards (KNPR 8715.3 and AFSPCMAN 91-710), emphasis should be placed on supporting the interpretation of recorded visual observations. Sound engineering practices and procedures should be applied with the interpretation of nondestructive testing (NDT) results when the residual strength data specific to a vessel design is incomplete or absent.
The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.