Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 415 records · Page 23

Laser-heated rotating specimen autoignition test

Specimens of 440 C steel were rotated in a chamber pressurized with oxygen gas and heated with a 5-kW CO2 laser to determine the temperature required for autoignition to occur. Tests included exposures of static and rotating (25,000 rpm) specimens in oxygen pressurized to 5.51 MPa, and with focused laser fluences of more than 3.5 billion W/sq m. Specimen surface temperatures were monitored with a scanning infrared camera. Temperature measurement difficulties were experienced due to a problem with internal reflection inside the test chamber; however, posttest specimen examinations confirmed that surface melt (1371 C) was achieved in several tests. No sustained combustion was initiated in any rotating specimen. One static specimen was ignited. Results indicated that conditions necessary for autoignition of 440 C steel are more dependent on specimen geometry and available heat removal mechanisms. Sustained combustion occurred in the ignited static specimen with an estimated 130 C/sec cooling rate due to conduction. The rotating specimens could not sustain combustion due to a greater conductive/convective cooling rate of about 4000 C/sec and ejection of molten material. These results were applied to the Space Shuttle Main Engine (SSME) oxygen turbopump bearings to conclude that the LOX-cooled 440 C steel bearings cannot sustain combustion initiated by skidding friction.

Au, A. C.↗

Initial Results of an MDO Method Evaluation Study

The NASA Langley MDO method evaluation study seeks to arrive at a set of guidelines for using promising MDO methods by accumulating and analyzing computational data for such methods. The data are collected by conducting a series of re- producible experiments. In the first phase of the study, three MDO methods were implemented in the SIGHT: framework and used to solve a set of ten relatively simple problems. In this paper, we comment on the general considerations for conducting method evaluation studies and report some initial results obtained to date. In particular, although the results are not conclusive because of the small initial test set, other formulations, optimality conditions, and sensitivity of solutions to various perturbations. Optimization algorithms are used to solve a particular MDO formulation. It is then appropriate to speak of local convergence rates and of global convergence properties of an optimization algorithm applied to a specific formulation. An analogous distinction exists in the field of partial differential equations. On the one hand, equations are analyzed in terms of regularity, well-posedness, and the existence and unique- ness of solutions. On the other, one considers numerous algorithms for solving differential equations. The area of MDO methods studies MDO formulations combined with optimization algorithms, although at times the distinction is blurred. It is important to

Alexandrov, Natalia M.↗

A Module for Assimilating Hyperspectral Infrared Retrieved Profiles into the Gridpoint Statistical Interpolation System for Unique Forecasting Applications

Hyperspectral infrared sounder radiance data are assimilated into operational modeling systems however the process is computationally expensive and only approximately 1% of available data are assimilated due to data thinning as well as the fact that radiances are restricted to cloud-free fields of view. In contrast, the number of hyperspectral infrared profiles assimilated is much higher since the retrieved profiles can be assimilated in some partly cloudy scenes due to profile coupling other data, such as microwave or neural networks, as first guesses to the retrieval process. As the operational data assimilation community attempts to assimilate cloud-affected radiances, it is possible that the use of retrieved profiles might offer an alternative methodology that is less complex and more computationally efficient to solve this problem. The NASA Short-term Prediction Research and Transition (SPoRT) Center has assimilated hyperspectral infrared retrieved profiles into Weather Research and Forecasting Model (WRF) simulations using the Gridpoint Statistical Interpolation (GSI) System. Early research at SPoRT demonstrated improved initial conditions when assimilating Atmospheric Infrared Sounder (AIRS) thermodynamic profiles into WRF (using WRF-Var and assigning more appropriate error weighting to the profiles) to improve regional analysis and heavy precipitation forecasts. Successful early work has led to more recent research utilizing WRF and GSI for applications including the assimilation of AIRS profiles to improve WRF forecasts of atmospheric rivers and assimilation of AIRS, Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI) profiles to improve model representation of tropopause folds and associated non-convective wind events. Although more hyperspectral infrared retrieved profiles can be assimilated into model forecasts, one disadvantage is the retrieved profiles have traditionally been assigned the same error values as the rawinsonde observations when assimilated with GSI. Typically, satellitederived profile errors are larger and more difficult to quantify than traditional rawinsonde observations (especially in the boundary layer), so it is important to appropriately assign observation errors within GSI to eliminate potential spurious innovations and analysis increments that can sometimes arise when using retrieved profiles. The goal of this study is to describe modifications to the GSI source code to more appropriately assimilate hyperspectral infrared retrieved profiles and outline preliminary results that show the differences between a model simulation that assimilated the profiles as rawinsonde observations and one that assimilated the profiles in a module with the appropriate error values.

Berndt, Emily↗

The Nexus for Exoplanet System Science

NExSS is a research coordination network dedicated to the study of planetary habitability. A NASA cross-division initiative bringing astrophysicists, planetary scientists, Earth scientists, and heliophysicists together to bring a systems science approach to this problem. NExSS's goals being to investigate the diversity of planets, understanding how planet history, geology, and climate interact to create the conditions for life. Also, to put planets into an architectural context as stellar systems built over time by dynamical processes and sculpted by stars. Use experience from solar system (including Earth) history to identify where habitable niches are most likely to occur and which planets are most likely to be habitable. Leverage NASA investments in research and missions to accelerate discovery and characterization of potential life-bearing worlds.

coordination↗

Trajectory Design Employing Convex Optimization for Landing on Irregularly Shaped Asteroids

Mission proposals that land on asteroids are becoming popular. However, in order to have a successful mission the spacecraft must reliably and softly land at the intended landing site. The problem under investigation is how to design a fuel-optimal powered descent trajectory that can be quickly computed on- board the spacecraft, without interaction from ground control. An optimal trajectory designed immediately prior to the descent burn has many advantages. These advantages include the ability to use the actual vehicle starting state as the initial condition in the trajectory design and the ease of updating the landing target site if the original landing site is no longer viable. For long trajectories, the trajectory can be updated periodically by a redesign of the optimal trajectory based on current vehicle conditions to improve the guidance performance. One of the key drivers for being completely autonomous is the infrequent and delayed communication between ground control and the vehicle. Challenges that arise from designing an asteroid powered descent trajectory include complicated nonlinear gravity fields, small rotating bodies and low thrust vehicles. There are two previous studies that form the background to the current investigation. The first set looked in-depth at applying convex optimization to a powered descent trajectory on Mars with promising results.1, 2 This showed that the powered descent equations of motion can be relaxed and formed into a convex optimization problem and that the optimal solution of the relaxed problem is indeed a feasible solution to the original problem. This analysis used a constant gravity field. The second area applied a successive solution process to formulate a second order cone program that designs rendezvous and proximity operations trajectories.3, 4 These trajectories included a Newtonian gravity model. The equivalence of the solutions between the relaxed and the original problem is theoretically established. The proposed solution for designing the asteroid powered descent trajectory is to use convex optimization, a gravity model with higher fidelity than Newtonian, and an iterative solution process to design the fuel optimal trajectory. The solution to the convex optimization problem is the thrust profile, magnitude and direction, that will yield the minimum fuel trajectory for a soft landing at the target site, subject to various mission and operational constraints. The equations of motion are formulated in a rotating coordinate system and includes a high fidelity gravity model. The vehicle's thrust magnitude can vary between maximum and minimum bounds during the burn. Also, constraints are included to ensure that the vehicle does not run out of propellant, or go below the asteroid's surface, and any vehicle pointing requirements. The equations of motion are discretized and propagated with the trapezoidal rule in order to produce equality constraints for the optimization problem. These equality constraints allow the optimization algorithm to solve the entire problem, without including a propagator inside the optimization algorithm.

Pinson, Robin M.↗

Active Control Analysis for Aeroelastic Instabilities in Turbomachines

Turbomachines onboard aircraft operate in a highly complex and harsh environment. The unsteady flowfield inherent to turbomachines leads to several problems associated with safety, stability, performance and noise. In-flight surge or flutter incidents could be catastrophic and impact the safety and reliability of the aircraft. High-Cycle-Fatigue (HCF), on the other hand, can significantly impact safety, readiness and maintenance costs. To avoid or minimize these problems generally a more conservative design method must be initiated which results in thicker blades and a loss of performance. Actively controlled turbomachines have the potential to reduce or even eliminate the instabilities by impacting the unsteady aerodynamic characteristics. By modifying the unsteady aerodynamics, active control may significantly improve the safety and performance especially at off-design conditions, reduce noise, and increase the range of operation of the turbomachine. Active control can also help improve reliability for mission critical applications such as the Mars Flyer. In recent years, HCF has become one of the major issues concerning the cost of operation for current turbomachines. HCF alone accounts for roughly 30% of maintenance cost for the United States Air-Force. Other instabilities (flutter, surge, rotating-stall, etc.) are generally identified during the design and testing phase. Usually a redesign overcomes these problems, often reducing performance and range of operation, and resulting in an increase in the development cost and time. Despite a redesign, the engines do not have the capabilities or means to cope with in-flight unforeseen vibration, stall, flutter or surge related instabilities. This could require the entire fleet worldwide to be stood down for expensive modifications. These problems can be largely overcome by incorporating active control within the turbomachine and its design. Active control can help in maintaining the integrity of the system in unforeseen events and provide for more aggressive designs to reduce the weight and improve efficiency of the turbomachine. Another area where active control can be useful is in controlling and suppressing rotating stall and surge in compressors, thereby increasing its operating range. Although some of these benefits will be offset by the added cost and weight penalty of the control system, the potential benefits in safety, reliability, performance, and noise characteristics are significant enough to warrant research in the area of active control of turbomachines. There is renewed interest within industry to understand unsteady aerodynamic behavior. This improved understanding not only leads to better design of turbomachines, which avoid instabilities but also which helps in understanding the controllability of the instabilities. The proliferation of micro-electro-mechanical-system (MEMS) devices has made available new tools to designers for employing feedback controls at reasonable costs. MEMS have also made the control devices small and unobtrusive enough to be implemented within the turbomachine without significant obstruction to the flow path. This has made active-control very attractive especially for systems requiring extreme confidence.

Keith, Theo G., Jr.↗

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Digital system identification and its application to digital flight control

On-line system identification of linear discrete systems for implementation in a digital adaptive flight controller is considered by the conventional extended Kalman filter and a decoupling process in which the linear state estimation problem and the linear parameter identification problem are each treated separately and alternately. Input requirements for parameter identifiability are established using the standard conditions of observability for a time variant system. Experimental results for simulated linearized lateral aircraft motion are included along with the effect of different initialization and updating procedures for the priming trajectory used by the filter.

Kotob, S.↗

Determination of instrumentation errors from measured data using maximum likelihood method

The maximum likelihood method is used for estimation of unknown initial conditions, constant bias and scale factor errors in measured flight data. The model for the system to be identified consists of the airplane six-degree-of-freedom kinematic equations, and the output equations specifying the measured variables. The estimation problem is formulated in a general way and then, for practical use, simplified by ignoring the effect of process noise. The algorithm developed is first applied to computer generated data having different levels of process noise for the demonstration of the robustness of the method. Then the real flight data are analyzed and the results compared with those obtained by the extended Kalman filter algorithm.

Keskar, D. A.↗

Stability and bifurcation of quasiparallel Alfven solitons

The inverse scattering transformation (IST) is used to study the one-parameter and two-parameter soliton families of the derivative nonlinear Schroedinger (DNLS) equation. The two-parameter soliton family is determined by the discrete complex eigenvalue spectrum of the Kaup-Newell scattering problem and the one-parameter soliton family corresponds to the discrete real eigenvalue spectrum. The structure of the IST is exploited to discuss the existence of discrete real eigenvalues and to prove their structural stability to perturbations of the initial conditions. Also, though the two-parameter soliton is structurally stable in general, it is shown that a perturbation of the initial conditions may change the two-parameter soliton into a degenerate soliton which, in turn, is structurally unstable. This degenerate, or double pole, soliton may bifurcate due to a perturbation of the initial conditions into a pair of one-parameter solitons. If the initial profile is on compact support, then this pair of one-parameter solitons must be compressive and rarefactive respectively. Finally, the Gelfand-Levitan equations appropriate for the double pole soliton are solved.

Hamilton, R. L.↗

Residual Stress Reversal in Highly Strained Shot Peened Structural Elements

The purpose of this research was to further the understanding of a crack initiation problem in a highly strained pressure containment housing. Finite Element Analysis methods were used to model the behavior of shot peened materials undergoing plastic deformation. Analytical results are in agreement with laboratory tensile tests that simulated the actual housing load conditions. These results further validate the original investigation finding that the shot peened residual stress had reversed, changing from compressive to tensile, and demonstrate that analytical finite element methods can be used to predict this behavior.

Mitchell, William S.↗

Understanding metal propagation in solid electrolytes due to mixed ionic-electronic conduction

Metal penetration into a solid electrolyte (SE) is one of the critical problems impeding the practical application of solid-state batteries. In this study, we investigate the conditions under which electronic conductivity of the SE can lead to metal deposition and fracture within the SE. Three different stages for void filling (metal plating initiation, metal growth, and metal compression) in the SE are identified. We show that a micron-size isolated void in the SE near the anode can be quickly filled in by metal and fractured when the developed pressure in the void grows larger than the maximum pressure the SE material can sustain. We find that the anode voltage and applied current density play a significant role in determining the vulnerability to metal deposition. We discuss several strategies to prevent electronic conductivity-driven metal propagation in electrolytes that are not fully dense, including the densified layers between the anode and SE.

25 ENERGY STORAGE↗

Optimal horizontal guidance law for aircraft in the terminal area.

Study of the horizontal guidance of aircraft in and near the terminal area. The problem of guiding an aircraft in minimum time from an arbitrary point to the outer marker is formulated as a nonlinear optimal control problem, and the control law solution is obtained by the application of the maximum principle. It is found that for some initial states the problem is singular. Furthermore, the extremal controls for this problem are not unique. Consequently, the optimal controls must be obtained on the basis of the value of the performance index. The control law is implemented in the form of a digital computer program which computes the optimal trajectory for arbitrary initial conditions.

Pecsvaradi, T.↗

Turbulence processes and simple closure schemes

The closure problem in turbulence is reviewed, and some simple closure schemes are introduced. Processes occurring in turbulent flow are discussed on the basis of solutions for some elementary flows such as homogeneous turbulence with and without uniform shear. Closure by specification of initial conditions, and finally practical closure schemes for more complicated flows are considered. The latter include Reynolds stress, eddy viscosity, and mixing-length closures.

Deissler, R. G.↗

Data presentation techniques for rotating machinery malfunction diagnosis

Baseline steady state data is excellent for documentation of vibration signals at normal operating conditions. Assuming that a set of initial data was acquired with the machinery in a good state of repair, any future changes or deterioration in mechanical condition can be easily compared to the baseline information. Often this type of comparison will yield sufficient information for evaluation of the problem. However, many malfunctions require the analysis of transient data in order to identify the malfunction. Steady-state data formats consist of: Time Base Waveform, Orbit, Spectrum. Transient data formats consist of: Polar, Bode, Cascade. Our objective is to demonstrate the use of the above formats to diagnose a machine malfunction. A turbine-driven compressor train is chosen as an example. The machine train outline drawing is shown.

Spettel, T.↗

An approximate loop transfer recovery method for designing fixed-order compensators

A method for designing fixed-order dynamic compensators for multivariable time invariant linear systems is presented which is based on the minimization of a linear quadratic performance index. The present formulation is performed in an output feedback setting which uses an observer cononical form to represent the compensator dynamics. Techniques for penalizing the plant and compensator states and for selecting the distribution on initial conditions such that the loop transfer matrix approximates that of a full-state feedback design have been developed. The effectiveness of the method is demonstrated using the examples of the pointing of a flexible structure and a helicopter flight control problem.

Calise, Anthony J.↗

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity were included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary (CMB) and the sub-adjacent fluid motion; it also features Gauss' method for solving the non-linear inverse problem associated with steady motional induction. The tradeoff between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models (DGRF's) is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition, a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Coerte V.↗