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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Integrated Demand Management: Minimizing Unanticipated Excessive Departure Delay while Ensuring Fairness from a Traffic Management Initiative

This paper introduces NASA's Integrated Demand Management (IDM) concept and presents the results from an early proof-of-concept evaluation and an exploratory experiment. An initial development of the concept was focused on integrating two systems - i.e. the FAA's newly deployed Traffic Flow Management System (TFMS) tool called the Collaborative Trajectory Options Program (CTOP) and the Time-Based Flow Management (TBFM) system with Extended Metering (XM) capabilities to manage projected heavy traffic demand into a capacity-constrained airport. A human-in-the-loop (HITL) simulation experiment was conducted to demonstrate the feasibility of the initial development of the concept by adapting it to an arrival traffic problem at Newark Liberty International Airport (EWR) during clear weather conditions. In this study, the CTOP was utilized to strategically plan the arrival traffic demand by controlling take-off times of both short- and long-haul flights (long-hauls specify aircraft outside TBFM regions and short-hauls specify aircraft within TBFM regions) in a way that results in equitable delays among the groups. Such strategic planning allows less airborne delay to occur within TBFM by feeding manageable long-haul traffic demand while reserving sufficient slots in the overhead streams for the short-haul departures. The manageable traffic demand indicates the TBFM scheduler assigns no more airborne delay than its assigned airspace is capable of absorbing. TBFM then uses its time-based metering capabilities to deliver the desirable throughput by tactically rescheduling the TBFM entered long-haul flights and short-haul departures. Additional research was also performed to explore use of Required Time of Arrival (RTA) capabilities as a potential control mechanism for the airborne flights to improve arrival traffic delivery accuracy of scheduled long-haul traffic demand. The study results show that both short- and long-haul flights received similar ground delays. In addition, there was a noticeable reduction in the total amount of excessive unanticipated last-minute ground delays, i.e. delays that are frequently imposed on the short-haul flight in current day operations due to saturation in the overhead stream, commonly referred to as 'double penalty'. Furthermore, the concept achieved the target throughput while minimizing the expected cost associated with overall delays in arrival traffic. Assessment of the RTA capabilities showed that there was indeed improvement of the scheduled entry times into TBFM regions by using RTA capabilities. However, with respect to reduction in delays incurred within TBFM, there was no observable benefit of improving the precision of long-haul flights entry times.

Collaborative Trajectory Options Program (CTOP)↗

NASA NTRS1977

Calculations of the effects of angular momentum on the early evolution of Jupiter

Early phases in the evolutionary development of a rotating cloud of gas and dust are investigated to determine whether physically reasonable initial conditions are likely to produce a planet-satellite system with the values of specific orbital angular momentum currently observed for Jupiter and its inner satellites. Spherically symmetric computations of Jupiter's evolutionary history are reviewed, and transport of angular momentum is considered along with evolutionary sequences of individual differentially rotating models with decreasing values of equatorial radius. The problem of accounting for the present specific orbital angular momenta of Jupiter and its regular satellites is examined qualitatively, and a range of possible angular-momentum distributions in an initially tidally unstable cloud is assessed. The results show that particular physically plausible angular-momentum distributions, which could either be established as initial conditions or be caused by angular-momentum transport, will lead to favorable situations for the formation of a central planet and a surrounding rotating cloud with specific angular momentum comparable to that of the regular satellites.

Bodenheimer, P.↗

NASA NTRS1986

A practical assessment of spectral accuracy for hyperbolic problems with discontinuities

Physical and transform space filtering has been applied to the Fourier spectral collocation solution of the constant coefficient scalar wave equation with a discontinuous initial condition. High order accuracy can be extracted from the unfiltered solution. Smooth, high order Fourier space filtering gives expected polynomial order solutions away from the discontinuity. Spectral accuracy is observed with the physical space filter of Gottlieb and Tadmor.

Kopriva, David A.↗

NASA NTRS1987

A practical assessment of spectral accuracy for hyperbolic problems with discontinuities

Kopriva, David A.↗

NASA NTRS1984

STARS: A general-purpose finite element computer program for analysis of engineering structures

STARS (Structural Analysis Routines) is primarily an interactive, graphics-oriented, finite-element computer program for analyzing the static, stability, free vibration, and dynamic responses of damped and undamped structures, including rotating systems. The element library consists of one-dimensional (1-D) line elements, two-dimensional (2-D) triangular and quadrilateral shell elements, and three-dimensional (3-D) tetrahedral and hexahedral solid elements. These elements enable the solution of structural problems that include truss, beam, space frame, plane, plate, shell, and solid structures, or any combination thereof. Zero, finite, and interdependent deflection boundary conditions can be implemented by the program. The associated dynamic response analysis capability provides for initial deformation and velocity inputs, whereas the transient excitation may be either forces or accelerations. An effective in-core or out-of-core solution strategy is automatically employed by the program, depending on the size of the problem. Data input may be at random within a data set, and the program offers certain automatic data-generation features. Input data are formatted as an optimal combination of free and fixed formats. Interactive graphics capabilities enable convenient display of nodal deformations, mode shapes, and element stresses.

Gupta, K. K.↗

DOE OSTI2023

Fast Bayesian optimization of Needle-in-a-Haystack problems using zooming memory-based initialization (ZoMBI)

Abstract Needle-in-a-Haystack problems exist across a wide range of applications including rare disease prediction, ecological resource management, fraud detection, and material property optimization. A Needle-in-a-Haystack problem arises when there is an extreme imbalance of optimum conditions relative to the size of the dataset. However, current state-of-the-art optimization algorithms are not designed with the capabilities to find solutions to these challenging multidimensional Needle-in-a-Haystack problems, resulting in slow convergence or pigeonholing into a local minimum. In this paper, we present a Zooming Memory-Based Initialization algorithm, entitled ZoMBI, that builds on conventional Bayesian optimization principles to quickly and efficiently optimize Needle-in-a-Haystack problems in both less time and fewer experiments. The ZoMBI algorithm demonstrates compute time speed-ups of 400× compared to traditional Bayesian optimization as well as efficiently discovering optima in under 100 experiments that are up to 3× more highly optimized than those discovered by similar methods.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Computer analysis of multicircuit shells of revolution by the field method

The field method, presented previously for the solution of even-order linear boundary value problems defined on one-dimensional open branch domains, is extended to boundary value problems defined on one-dimensional domains containing circuits. This method converts the boundary value problem into two successive numerically stable initial value problems, which may be solved by standard forward integration techniques. In addition, a new method for the treatment of singular boundary conditions is presented. This method, which amounts to a partial interchange of the roles of force and displacement variables, is problem independent with respect to both accuracy and speed of execution. This method was implemented in a computer program to calculate the static response of ring stiffened orthotropic multicircuit shells of revolution to asymmetric loads. Solutions are presented for sample problems which illustrate the accuracy and efficiency of the method.

Cohen, G. A.↗

NASA NTRS1974

COMOC: Three dimensional boundary region variant, programmer's manual

The three-dimensional boundary region variant of the COMOC computer program system solves the partial differential equation system governing certain three-dimensional flows of a viscous, heat conducting, multiple-species, compressible fluid including combustion. The solution is established in physical variables, using a finite element algorithm for the boundary value portion of the problem description in combination with an explicit marching technique for the initial value character. The computational lattice may be arbitrarily nonregular, and boundary condition constraints are readily applied. The theoretical foundation of the algorithm, a detailed description on the construction and operation of the program, and instructions on utilization of the many features of the code are presented.

Orzechowski, J. A.↗

NASA NTRS1987

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗

NASA NTRS1986

Turbulent relaxation processes in magnetohydrodynamics

Competing processes previously called 'selective decay' and 'dynamic alignment' are studied numerically for two-dimensional magnetohydrodynamic turbulence. In selective decay, the energy decays relatively to mean-square vector potential, and in dynamic alignment, the energy decays relatively to cross helicity. In the former case, the kinetic (fluid) energy decays to zero and the magnetic energy occupies the largest scales allowed by the boundary conditions. In the latter case, the velocity field and magnetic field become aligned and energetically equipartitioned. An extensive study of the initial value problem, with viscous and resistive dissipation, indicates that four distinct regimes of behavior are possible: a magnetically dominated regime, a velocity dominated regime, a dynamic alignment dominated regime, and a transition regime separating the others in parameter space. An analytical variational problem predicts several features seen in the computations, including geometrical alignment of velocity and magnetic fields, constraints on global quantities in the long-time limit, and in some cases, condensation to long wavelength modes.

Ting, Antonio C.↗

NASA NTRS1991

In-phase and out-of-phase axial-torsional fatigue behavior of Haynes 188 at 760 C

Isothermal, in-phase and out-of-phase axial-torsional fatigue experiments have been conducted at 760 C on uniform gage section, thin-walled tubular specimens of a wrought cobalt-base superalloy, Haynes 188. Test-control and data acquisition were accomplished with a minicomputer. Fatigue lives of the in- and out-of-phase axial-torsional fatigue tests have been estimated with four different multiaxial fatigue life prediction models that were developed primarly for predicting axial-torsional fatigue lives at room temperature. The models investigated were: (1) the von Mises equivalent strain range; (2) the Modified Multiaxiality Factor Approach; (3) the Modified Smith-Watson-Topper Parameter; and (4) the critical shear plane method of Fatemi, Socie, and Kurath. In general, life predictions by the von Mises equivalent strain range model were within a factor of 2 for a majority of the tests and the predictions by the Modified Multiaxiality Factor Approach were within a factor of 2, while predictions of the Modified Smith-Watson-Topper Parameter and of the critical shear plane method of Fatemi, Socie, and Kurath were unconservative and conservative, respectively, by up to factors of 4. In some of the specimens tested under combined axial-torsional loading conditions, fatigue cracks initiated near extensometer indentations. Two design modifications have been proposed to the thin-walled tubular specimen to overcome this problem.

Kalluri, Sreeramesh↗

NASA NTRS1994

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗

NASA NTRSDate not supplied

Comprehensive Analysis Aeroelastic Stability Predictions for the Mars Sample Recovery Helicopters and Ingenuity

A rotorcraft comprehensive analysis was used to generate aeroelastic stability predictions for both Ingenuity and the Mars Sample Recovery Helicopters. The analysis was performed multiple times for each rotor system with sequential increases to the level of modeling complexity. Initial cases used an independent blade analysis, trim inflow distribution, no airframe degrees of freedom, and a hover condition (a time-invariant problem). The final model used all blades, dynamic inflow, airframe degrees of freedom, and a Floquet theory approach to the flutter analysis. The approach of sequentially increasing modeling complexity was employed to ensure the causes of any observed instabilities could be identified. For each case, stability results were determined across a range of air densities, tip Mach numbers, and rotor collective pitch settings. For the Sample Recovery Helicopter, a forward flight condition analysis with the highest level of modeling complexity was also performed. The aeroelastic stability predictions showed both the Ingenuity and Sample Recovery Helicopter rotor systems to be stable for all analyzed conditions. Parameter excursions of structural properties were performed to explore conditions resulting in aeroelastic instabilities.

Aeroelastic↗

NASA NTRS1997

Three-Dimensional Viscous Alternating Direction Implicit Algorithm and Strategies for Shape Optimization

A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.

Pandya, Mohagna J.↗

NASA NTRS1980

Space Shuttle stability and control flight test techniques

A unique approach for obtaining vehicle aerodynamic characteristics during entry has been developed for the Space Shuttle. This is due to the high cost of Shuttle testing, the need to open constraints for operational flights, and the fact that all flight regimes are flown starting with the first flight. Because of uncertainties associated with predicted aerodynamic coefficients, nine flight conditions have been identified at which control problems could occur. A detailed test plan has been developed for testing at these conditions and is presented. Due to limited testing, precise computer initiated maneuvers are implemented. These maneuvers are designed to optimize the vehicle motion for determining aerodynamic coefficients. Special sensors and atmospheric measurements are required to provide stability and control flight data during an entire entry. The techniques employed in data reduction are proven programs developed and used at NASA/DFRC.

Cooke, D. R.↗

NASA NTRS1980

Computer code simulations of the formation of Meteor Crater, Arizona - Calculations MC-1 and MC-2

It has been widely accepted that hypervelocity impact processes play a major role in the evolution of the terrestrial planets and satellites. In connection with the development of quantitative methods for the description of impact cratering, it was found that the results provided by two-dimensional finite difference, computer codes is greatly improved when initial impact conditions can be defined and when the numerical results can be tested against field and laboratory data. In order to address this problem, a numerical code study of the formation of Meteor (Barringer) Crater, Arizona, has been undertaken. A description is presented of the major results from the first two code calculations, MC-1 and MC-2, that have been completed for Meteor Crater. Both calculations used an iron meteorite with a kinetic energy of 3.8 Megatons. Calculation MC-1 had an impact velocity of 25 km/sec and MC-2 had an impact velocity of 15 km/sec.

Roddy, D. J.↗

NASA NTRS1982

Utilization of normal mode initial conditions for detecting errors in the dynamics part of primitive equation global models

When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Chao, W. C.↗

NASA NTRS1991

A multi-point inverse airfoil design method based on conformal mapping

An exact method of multipoint inverse airfoil design for incompressible flow is presented. Multipoint design is handled by dividing the airfoil into a number of desired segments. For each segment the velocity distribution is prescribed together with an angle of attack at which the prescribed velocity distribution is to be achieved. In this manner, multipoint design objectives can be taken into account in the initial specification of the velocity distribution. In order for the multipoint inverse airfoil design problem to be well posed, three integral constraints and several conditions arise which must be satisfied. Further restrictions are imposed if the airfoil is to have a specified pitching moment, thickness ratio, or other constraints. The system of equations is solved partly as a linear system and partly through multidimensional Newton iteration. Since the velocity distribution is prescribed about the circle in the angular coordinate, specification of the velocity in terms of arc length is handled through the multidimensional Newton iteration as well. The current formulation sets the stage for a more general multipoint inverse airfoil design method in which it will be possible to specify the velocity distribution, some boundary-layer development, or the surface geometry along a segment.

Selig, Michael S.↗