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 145 records · Page 8

Long-time asymptotics of a system for plasma diffusion

A system of parabolic nonlinear equations that describe the diffusion of a fully collisional plasma across a strong magnetic field is discussed. It is demonstrated that the solution to this system tends to a time asymptotic state which is of space-time separable form, theta(t)f(x). Furthermore, f(x) is independent of the initial conditions and theta(t) depends slightly on the initial conditions. The rate of decay of the temporal part is governed by a nonlinear eigenvalue problem. Since the equations are considered in a bounded domain we are able to analyze the effect of boundary conditions on the evolution of the system. Additional effects on radiation, heating, and particle injection can also be accounted for. Essential differences between the behavior of a fully-coupled system and a scalar equation are observed.

Rosenau, P.↗

Evolution of the Orszag-Tang vortex system in a compressible medium. II - Supersonic flow

A study is presented on the effect of embedded supersonic flows and the resulting emerging shock waves on phenomena associated with MHD turbulence, including reconnection, the formation of current sheets and vortex structures, and the evolution of spatial and temporal correlations among physical variables. A two-dimensional model problem, the Orszag-Tang (1979) vortex system, is chosen, which involves decay from nonrandom initial conditions. The system is doubly periodic, and the initial conditions consist of single-mode solenoidal velocity and magnetic fields, each containing X points and O points. The initial mass density is flat, and the initial pressure fluctuations are incompressible, balancing the local forces for a magnetofluid of unit mass density. Results on the evolution of the local structure of the flow field, the global properties of the system, and spectral correlations are presented. The important dynamical properties and observational consequences of embedded supersonic regions and emerging shocks in the Orszag-Tang model of an MHD system undergoing reconnection are discussed. Conclusions are drawn regarding the effects of local supersonic regions on MHD turbulence.

Picone, J. Michael↗

Efficient parametric analysis of performance measures for communication networks

Efficient techniques for estimating performance measures in communication networks in a steady-state or transient setting are developed. These techniques may be used in a simulation environment or in connection with real-time observations. For Markov chain models, the recently proposed standard clock approach is extended, and a class of real-time algorithms for simultaneously generating multiple sample paths under different parameter sets is presented. Attention is focused on the link crash time estimation problem, where a 'crash' is defined as the first time a buffer overflows, given some initial conditions. An algorithm for estimating crash times under various traffic shocks is derived, where all estimates are obtained in parallel to an actual network's normal operation. An algorithm for crash time estimation of a G/D/1 link model is also derived using a different (perturbation-analysis-based) approach. Finally, extensive simulation results are provided to validate the proposed algorithms and compare them to brute-force simulation.

Cassandras, Christos G.↗

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

Backscattering and Nonparaxiality Arrest Collapse of Damped Nonlinear Waves

The critical nonlinear Schrodinger equation (NLS) models the propagation of intense laser light in Kerr media. This equation is derived from the more comprehensive nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. It is known that if the input power of the laser beam (i.e., L(sub 2) norm of the initial solution) is sufficiently high, then the NLS model predicts that the beam will self-focus to a point (i.e.. collapse) at a finite propagation distance. Mathematically, this behavior corresponds to the formation of a singularity in the solution of the NLS. A key question which has been open for many years is whether the solution to the NLH, i.e., the 'parent' equation, may nonetheless exist and remain regular everywhere, in particular for those initial conditions (input powers) that lead to blowup in the NLS. In the current study, we address this question by introducing linear damping into both models and subsequently comparing the numerical solutions of the damped NLH (boundary-value problem) with the corresponding solutions of the damped NLS (initial-value problem). Linear damping is introduced in much the same way as done when analyzing the classical constant-coefficient Helmholtz equation using the limiting absorption principle. Numerically, we have found that it provides a very efficient tool for controlling the solutions of both the NLH and NHS. In particular, we have been able to identify initial conditions for which the NLS solution does become singular. whereas the NLH solution still remains regular everywhere. We believe that our finding of a larger domain of existence for the NLH than that for the NLS is accounted for by precisely those mechanisms, that have been neglected when deriving the NLS from the NLH, i.e., nonparaxiality and backscattering.

Fibich, G.↗

A Computational Intelligence (CI) Approach to the Precision Mars Lander Problem

A Mars precision landing requires a landed footprint of no more than 100 meters. Obstacles to reducing the landed footprint include trajectory dispersions due to initial atmospheric entry conditions such as entry angle, parachute deployment height, environment parameters such as wind, atmospheric density, parachute deployment dynamics, unavoidable injection error or propagated error from launch, etc. Computational Intelligence (CI) techniques such as Artificial Neural Nets and Particle Swarm Optimization have been shown to have great success with other control problems. The research period extended previous work on investigating applicability of the computational intelligent approaches. The focus of this investigation was on Particle Swarm Optimization and basic Neural Net architectures. The research investigating these issues was performed for the grant cycle from 5/15/01 to 5/15/02. Matlab 5.1 and 6.0 along with NASA's POST were the primary computational tools.

Birge, Brian↗

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗

Techniques for trajectory optimization using a hybrid computer

The use of a hybrid computer in the solution of trajectory optimization problems is described. The solution technique utilizes the indirect method and requires iterative computation of the initial condition vector of the co-state variables. Convergence of the iteration is assisted by feedback switching and contour modification. A simulation of the method in an on-line updating scheme is presented.

Neely, P. L.↗

Three dimensional thermal pollution models. Volume 1: Review of mathematical formulations

A mathematical model package for thermal pollution analyses and prediction is presented. These models, intended as user's manuals, are three dimensional and time dependent using the primitive equation approach. Although they have sufficient generality for application at sites with diverse topographical features; they also present specific instructions regarding data preparation for program execution and sample problems. The mathematical formulation of these models is presented including assumptions, approximations, governing equations, boundary and initial conditions, numerical method of solution, and same results.

Lee, S. S.↗

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.

Gottlieb, D.↗

The dynamics and optimal control of spinning spacecraft and movable telescoping appendages, part A

The problem of optimal control with a minimum time criterion as applied to a single boom system for achieving two axis control is discussed. The special case where the initial conditions are such that the system can be driven to the equilibrium state with only a single switching maneuver in the bang-bang optimal sequence is analyzed. The system responses are presented. Application of the linear regulator problem for the optimal control of the telescoping system is extended to consider the effects of measurement and plant noises. The noise uncertainties are included with an application of the estimator - Kalman filter problem. Different schemes for measuring the components of the angular velocity are considered. Analytical results are obtained for special cases, and numerical results are presented for the general case.

Bainum, P. M.↗

The Effect of Lateral- and Longitudinal- Range Control on Allowable Entry Conditions for a Point Return from Space

The problem of return to a specified landing point on the earth from flight in space is considered by studying the interaction between an assumed control over the lateral and longitudinal range and the initial conditions of approach to the earth, given by orbital-plane inclination, vacuum perigee location, and time of arrival. The maneuvering capability in the atmosphere permits a point return for a range of entry conditions. A lateral-range capability of +/- 500 miles from the center line of an entry trajectory can allow a variation in the time of arrival of over 3.5 hours. Variation in the orbital-plane inclination angle can be as much as +/- 13 deg.

Boissevain, Alfred G.↗

Decision-problem state analysis methodology

A methodology for analyzing a decision-problem state is presented. The methodology is based on the analysis of an incident in terms of the set of decision-problem conditions encountered. By decomposing the events that preceded an unwanted outcome, such as an accident, into the set of decision-problem conditions that were resolved, a more comprehensive understanding is possible. All human-error accidents are not caused by faulty decision-problem resolutions, but it appears to be one of the major areas of accidents cited in the literature. A three-phase methodology is presented which accommodates a wide spectrum of events. It allows for a systems content analysis of the available data to establish: (1) the resolutions made, (2) alternatives not considered, (3) resolutions missed, and (4) possible conditions not considered. The product is a map of the decision-problem conditions that were encountered as well as a projected, assumed set of conditions that should have been considered. The application of this methodology introduces a systematic approach to decomposing the events that transpired prior to the accident. The initial emphasis is on decision and problem resolution. The technique allows for a standardized method of accident into a scenario which may used for review or the development of a training simulation.

Dieterly, D. L.↗

Benchmark problems in computational aeroacoustics

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Porter-Locklear, Freda↗

A study and analysis of the MSFC lunar roving vehicle dust profile test program

The dust problem and fender design for the LRV were studied under reduced gravity with a lunar soil simulant. The test equipment, soil characteristics of the lunar soil simulant, and the test procedures are described. It is concluded: (1) The fender plus flap design is adequate. (2) Vacuum conditions tend to eliminate or reduce suspended dust clouds. (3) Reduced gravity conditions tend to increase the dust problems. (4) Slow starting speeds are necessary to minimize slip and reduce initial dust generation.

Mullis, C. H.↗

Transmitted sound field due to an impulsive line acoustic source bounded by a plate followed by a vortex sheet

The propagation of sound due to a line acoustic source in the moving stream across a semiinfinite vortex sheet which trails from a rigid plate is examined in a linear theory for the subsonic case. A solution for the transmitted sound field is obtained with the aid of multiple integral transforms and the Wiener-Hopf technique for both the steady state (time harmonic) and initial value (impulsive source) situations. The contour of inverse transform and hence the decomposition of the functions are determined through causality and radiation conditions. The solution obtained satisfies causality and the full Kutta conditions. The transmitted sound field is composed of two waves in both the stady state and initial value problems. One is the wave scattered from the edge of the plate which is associated with the bow wave and the instability wave. These waves exist in the downstream sectors. The other is the wave transmitted through the vortex sheet which is also associated with the instability wave. Regional divisions of the transmitted sound field are identified.

Miura, T.↗

Calculation of two-dimensional inlet flow fields by an implicit method including viscous effects: User's manual

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time marching finite difference procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. The mathematical formulation of the problem and the solution algorithm are given. Numerical stability and accuracy as well as the initial and boundary conditions used are discussed. Instructions for program use and operation along with the overall program logic are also given.

Biringen, S.↗