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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Numerical method and FORTRAN program for the solution of an axisymmetric electrostatic collector design problem

The numerical calculation is described of the steady-state flow of electrons in an axisymmetric, spherical, electrostatic collector for a range of boundary conditions. The trajectory equations of motion are solved alternately with Poisson's equation for the potential field until convergence is achieved. A direct (noniterative) numerical technique is used to obtain the solution to Poisson's equation. Space charge effects are included for initial current densities as large as 100 A/sq cm. Ways of dealing successfully with the difficulties associated with these high densities are discussed. A description of the mathematical model, a discussion of numerical techniques, results from two typical runs, and the FORTRAN computer program are included.

Reese, O. W.↗
NASA NTRS1980

A new algorithm for horizontal capture trajectories

An algorithm which transfers an aircraft from an initial position and heading to a final position and heading was developed for onboard synthesis of horizontal flight paths. The algorithm finds all solutions possible, and selects the one with minimum path length. Degenerate conditions in which one or more of the basic segments is missing are handled without difficulty. The solution to this problem is derived, and a FORTRAN listing of the algorithm is provided.

Mclean, J. D.↗
NASA NTRS1981

Minimum fuel horizontal flightpaths in the terminal area

The problem of minimum fuel airplane trajectories from arbitrary initial states to be fixed final state is considered. There are four state variables (two position coordinates, heading, and constrained velocity) and two constrained controls (thrust and bank angle). The fuel optimality of circular and straight line flightpaths is examined. Representative extremals (trajectories satisfying the necessary conditions of the minimum principle) of various types are computed and used to evaluate trajectories generated by an on line algorithm. Attention is paid to the existence of Darboux points (beyond which an extremal ceases to be globally optimal). One fuel flow rate model includes a term quadratic in thrust; hence, the optimal thrust is continuous and nonsingular. The other fuel flow rate model is linear in thrust, and consequently the optimal thrust is discontinuous and singular.

Kreindler, E.↗
NASA NTRS1985

A sequential linear optimization approach for controller design

A linear optimization approach with a simple real arithmetic algorithm is presented for reliable controller design and vibration suppression of flexible structures. Using first order sensitivity of the system eigenvalues with respect to the design parameters in conjunction with a continuation procedure, the method converts a nonlinear optimization problem into a maximization problem with linear inequality constraints. The method of linear programming is then applied to solve the converted linear optimization problem. The general efficiency of the linear programming approach allows the method to handle structural optimization problems with a large number of inequality constraints on the design vector. The method is demonstrated using a truss beam finite element model for the optimal sizing and placement of active/passive-structural members for damping augmentation. Results using both the sequential linear optimization approach and nonlinear optimization are presented and compared. The insensitivity to initial conditions of the linear optimization approach is also demonstrated.

Horta, L. G.↗
NASA NTRS1986

Fast algorithm for calculating chemical kinetics in turbulent reacting flow

This paper addresses the need for a fast batch chemistry solver to perform the kinetics part of a split operator formulation of turbulent reacting flows, with special attention focused on the solution of the ordinary differential equations governing a homogeneous gas-phase chemical reaction. For this purpose, a two-part predictor-corrector algorithm which incorporates an exponentially fitted trapezoidal method was developed. The algorithm performs filtering of ill-posed initial conditions, automatic step-size selection, and automatic selection of Jacobi-Newton or Newton-Raphson iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm, termed CREK1D (combustion reaction kinetics, one-dimensional), compared favorably with the code LSODE when tested on two representative problems drawn from combustion kinetics, and is faster than LSODE.

Radhakrishnan, K.↗
NASA NTRS1992

Thermally-induced, geometrically nonlinear response of symmetrically laminated composite plates

This paper discusses the thermally-induced geometrically nonlinear response of symmetrically laminated composite plates. The plate response is due to a temperature increase that is uniform in the plane of the plate but has a slight gradient through the thickness. The case of a completely uniform temperature increase but with an initial out-of-plane imperfection in the plate is also considered. Because they are closely allied problems, thermal buckling and postbuckling are discussed. Using variational methods in conjunction with a Rayleigh-Ritz formulation, these responses are investigated for two laminates, a (+/- 45/0/90)s and a (+/- 45/02)s, under two different simple support conditions.

Meyers, C. A.↗
NASA NTRS1993

Computation of optimal low- and medium-thrust orbit transfers

This paper presents the formulation of the optimal low- and medium-thrust orbit transfer control problem, numerical methods for solution, and numerical solutions of the problem. The problem formulation is for final mass maximization and allows for second-harmonic oblateness, atmospheric drag, and 3D noncoplanar nonaligned elliptic terminal orbits. We set up examples to demonstrate the ability of two indirect methods to solve the resulting two point boundary value problems (TPBVP). The methods demonstrated are the multiple point shooting method as formulated in Oberle's (1987) subroutine BOUNDSCO, and the minimizing boundary-condition method (MBCM). We find that although both methods can converge solutions, there are tradeoffs to using either method. We present numerical solutions of planar transfers in which both the initial orbit exit and final orbit entry points have been optimized. These solutions include two- and three-burn transfers. The methods used show an ability to handle thrust down to at least T/W(o) = O(10 exp -3). They also show similar convergence abilities with or without the oblateness and drag terms. We discuss the issue of maximizing with respect to the final time and provide evidence that implies a local optimum at a maximum final time for a given number of burns.

Chuang, C.-H.↗
NASA NTRS1998

The Origin and Time Dependence of the Amount and Composition of Non-Constituent Gases Present in Crystal Growth Systems

Presence of different, non-constituent gases may be a critical factor in crystal growth systems. In Physical Vapor Transport processes the cras(es) can be used intentionally (to prevent excessively high, unstable growth conditions), or can evolve unintentionally during the course of the process (which may lead to undesired reduction in the -rowth rate). In melt growth, particularly under low gravity conditions (reduced hydrostatic pressure) the gas present in the system may contribute to formation of voids in the growing crystals and even to a separation of the crystal and the liquid phase [1]. On the other hand, some amount of gas may facilitate 'contactless' crystal growth particularly under reduced gravity conditions [2 - 6]. Different non-constituent gases may be present in growth ampoules, and their amount and composition may change during the crystallization process. Some gases can appear even in empty ampoules sealed originally under high vacuum: they may diffuse in from the outside, and/or desorb from the ampoule walls. Residual gases can also be generated by the source materials: even very high purity commercial elements and compounds may contain trace amounts of impurities, particularly oxides. The oxides may have low volatilities themselves but their reaction with other species, particularly carbon and hydrogen, may produce volatile compounds like water or carbon oxides. The non-constituent gases, either added initially to the system or evolved during the material processing, may diffuse out of the ampoule during the course of the experiment. Gases present outside (e.g. as a protective atmosphere or thermal conductor) may diffuse into the ampoule. In either case the growth conditions and the quality of the crystals may be affected. The problem is of a particular importance in sealed systems where the amount of the gases cannot be directly controlled. Therefore a reasonable knowledge and understanding of the origin, composition, magnitude, and change with time of gases present in sealed ampoules may be important for a meaningful control and interpretation of crystal growth processes. This problem is of a particular importance for processing of electronic materials in space because (i) safety considerations require using sealed systems only, and (ii) high cost of crystal growth experiments in microgravity calls for a throughout, accurate description of the processing conditions necessary for a meaningful, efficient, and conclusive interpretation of the space results. In this paper we present the results of our extensive studies on gases in closed crystal growth systems which include: (a) Degassing properties of fused silica; (b) Generation of inert gases by source materials (CdTe, ZnTe, CdZnTe, ZnSe, PbTe, PbSe, PbSeTe); (c) Diffusive cas losses from silica glass ampoules.

Palosz, Witold↗
DOE OSTI2020

FY2018 Fusion Theory and Simulation Milestone on Plasma Materials Interaction (PMI) (Final Report)

The interaction of the boundary plasma with the material surfaces in magnetically confined plasmas is among the most critical problems in fusion energy science. In FY2018, perform high performance computing simulations with coupled boundary plasma physics and materials surface models to predict the fuel recycling and tritium retention of the ITER divertor for D-T burning plasma conditions, accounting for erosion, re-deposition and impurity transport in the plasma boundary, and an initial evaluation of the influence of material deposition on the recycling and retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Orbit transfer rocket engine technology program: Advanced engine study

In Task D.6 of the Advanced Engine Study, three primary subtasks were accomplished: (1) design of parametric data; (2) engine requirement variation studies; and (3) vehicle study/engine study coordination. Parametric data were generated for vacuum thrusts ranging from 7500 lbf to 50,000 lbf, nozzle expansion ratios from 600 to 1200, and engine mixture ratios from 5:1 to 7:1. Failure Modes and Effects Analysis (FMEA) was used as a departure point for these parametric analyses. These data are intended to assist in definition and trade studies. In the Engine Requirements Variation Studies, the individual effects of increasing the throttling ratio from 10:1 to 20:1 and requiring the engine to operate at a maximum mixture ratio of 12:1 were determined. Off design engine balances were generated at these extreme conditions and individual component operating requirements analyzed in detail. Potential problems were identified and possible solutions generated. In the Vehicle Study/Engine Study coordination subtask, vehicle contractor support was provided as needed, addressing a variety of issues uncovered during vehicle trade studies. This support was primarily provided during Technical Interchange Meetings (TIM) in which Space Exploration Initiative (SEI) studies were addressed.

Erickson, C. M.↗
NASA NTRS1996

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

Boiling is an effective mode of heat transfer since high heat flux levels are possible driven by relatively small temperature differences. The high heat transfer coefficients associated with boiling have made the use of these processes increasingly attractive to aerospace engineering. Applications of this type include compact evaporators in the thermal control of aircraft avionics and spacecraft environments, heat pipes, and use of boiling to cool electronic equipment. In spite of its efficiency, cooling based on liquid-vapor phase change processes has not yet found wide application in aerospace engineering due to specific problems associated with the low gravity environment. After a heated surface has reached the superheat required for the initiation of nucleate boiling, bubbles will start forming at nucleation sites along the solid interface by evaporation of the liquid. Bubbles in contact with the wall will continue growing by this mechanism until they detach. In terrestrial conditions, bubble detachment is determined by the competition between body forces (e.g. buoyancy) and surface tension forces that act to anchor the bubble along the three phase contact line. For a given body force potential and a balance of tensions along the three phase contact line, bubbles must reach a critical size before the body force can cause them to detach from the wall. In a low gravity environment the critical bubble size for detachment is much larger than under terrestrial conditions, since buoyancy is a less effective means of bubble removal. Active techniques of heat transfer enhancement in single phase and phase change processes by utilizing electric fields have been the subject of intensive research during recent years. The field of electrohydrodynamics (EHD) deals with the interactions between electric fields, flow fields and temperature fields. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50 as compared to values obtained for the same system without electric fields. Imposing an external electric field holds the promise to improve pool boiling heat transfer in low gravity, since a phase separation force other than gravity is introduced. The goal of our research is to experimentally investigate the potential of EHD and the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions.

Herman, Cila↗
NASA NTRS2003

Scheduling Earth Observing Satellites with Evolutionary Algorithms

We hypothesize that evolutionary algorithms can effectively schedule coordinated fleets of Earth observing satellites. The constraints are complex and the bottlenecks are not well understood, a condition where evolutionary algorithms are often effective. This is, in part, because evolutionary algorithms require only that one can represent solutions, modify solutions, and evaluate solution fitness. To test the hypothesis we have developed a representative set of problems, produced optimization software (in Java) to solve them, and run experiments comparing techniques. This paper presents initial results of a comparison of several evolutionary and other optimization techniques; namely the genetic algorithm, simulated annealing, squeaky wheel optimization, and stochastic hill climbing. We also compare separate satellite vs. integrated scheduling of a two satellite constellation. While the results are not definitive, tests to date suggest that simulated annealing is the best search technique and integrated scheduling is superior.

Globus, Al↗
NASA NTRS1993

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗
DOE OSTI2026

Instant folded strings, dark energy and a cyclic bouncing universe

We present a wholly self-consistent, complete cyclic bouncing cosmology based on components drawn from string theory and constructed in a way that is under perturbative control throughout (e.g., with temperature much less than the string scale and string coupling g s ≪ 1 at all times). The cyclic evolution is governed by standard dilaton-gravity in (3+1)-dimensions with a perturbatively generated potential and a coupling between the dilaton and a second field that becomes massless at ϕ = ϕ ESP , resulting in an enhanced symmetry point (ESP) that prevents the dilaton from running all the way to zero coupling. A central role is played by instant folded strings (IFSs) — fundamental strings with the unusual property of being much lighter than the string mass while extending far beyond the string length, and violating the Null Energy Condition (NEC). IFSs are produced classically when the string coupling grows with time, which occurs at two critical points in each cycle. In turn, they fulfill a dual function: Enabling cosmological bounces and initiating transient epochs of dark-energy domination that naturally transition into slow contraction. The resulting cosmology eliminates the cosmic singularity and multiverse problems of big bang inflationary models and robustly predicts time-varying IFS-induced dark energy and the absence of primordial B-mode polarization in the cosmic microwave background.

alternatives to inflation↗
DOE OSTI2025

Instant Folded Strings, Dark Energy and a Cyclic Bouncing Universe

We present a wholly self-consistent, complete cyclic bouncing cosmology based on components drawn from string theory and constructed in a way that is under perturbative control throughout (e.g., with temperature much less than the string scale and string coupling $g_s \ll 1$ at all times). The cyclic evolution is governed by standard dilaton-gravity in $(3+1)$-dimensions with a perturbatively generated potential and a coupling between the dilaton and a second field that becomes massless at $ϕ= ϕ_{ESP}$, resulting in an enhanced symmetry point (ESP) that prevents the dilaton from running all the way to zero coupling. A central role is played by instant folded strings (IFSs) - fundamental strings with the unusual property of being much lighter than the string mass while extending far beyond the string length, and violating the Null Energy Condition (NEC). IFSs are produced classically when the string coupling grows with time, which occurs at two critical points in each cycle. In turn, they fulfill a dual function: enabling cosmological bounces and initiating transient epochs of dark-energy domination that naturally transition into slow contraction. The resulting cosmology eliminates the cosmic singularity and multiverse problems of big bang inflationary models and robustly predicts time-varying IFS-induced dark energy and the absence of primordial B-mode polarization in the cosmic microwave background.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
NASA NTRS1972

Calculation of free fall trajectories based on numerical optimization techniques

The development of a means of computing free-fall (nonthrusting) trajectories from one specified point in the solar system to another specified point in the solar system in a given amount of time was studied. The problem is that of solving a two-point boundary value problem for which the initial slope is unknown. Two standard methods of attack exist for solving two-point boundary value problems. The first method is known as the initial value or shooting method. The second method of attack for two-point boundary value problems is to approximate the nonlinear differential equations by an appropriate linearized set. Parts of both boundary value problem solution techniques described above are used. A complete velocity history is guessed such that the corresponding position history satisfies the given boundary conditions at the appropriate times. An iterative procedure is then followed until the last guessed velocity history and the velocity history obtained from integrating the acceleration history agree to some specified tolerance everywhere along the trajectory.

Source record↗
NASA NTRS1984

Hovering model helicopter rotor testing

A relatively large test chamber, chosen to minimize room interference effects during helicopter hovering rotor tests, was found in initial operations to exhibit large, random variations of thrust, torque, and blade flapping. The source of these problems, unsteady recirculation of test chamber flow, was addressed by the addition of a 'flow straightener' that reduced these parameter variations to 1.5 percent of mean values. Rotor performance was then measured outside the chamber in calm, free air conditions and compared with chamber results obtained before and after modifications. Attention is given to wake flow visualization conducted in the test chamber, as well as outside.

Piziali, R. A.↗
NASA NTRS1988

A comparison of experimental and calculated thin-shell leading-edge buckling due to thermal stresses

High-temperature thin-shell leading-edge buckling test data are analyzed using NASA structural analysis (NASTRAN) as a finite element tool for predicting thermal buckling characteristics. Buckling points are predicted for several combinations of edge boundary conditions. The problem of relating the appropriate plate area to the edge stress distribution and the stress gradient is addressed in terms of analysis assumptions. Local plasticity was found to occur on the specimen analyzed, and this tended to simplify the basic problem since it effectively equalized the stress gradient from loaded edge to loaded edge. The initial loading was found to be difficult to select for the buckling analysis because of the transient nature of thermal stress. Multiple initial model loadings are likely required for complicated thermal stress time histories before a pertinent finite element buckling analysis can be achieved. The basic mode shapes determined from experimentation were correctly identified from computation.

Jenkins, Jerald M.↗