Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

The potential impact of scatterometry on oceanography - A wave forecasting case

A series of observing system simulation experiments have been performed in order to assess the potential impact of marine surface wind data on numerical weather prediction. In addition to conventional data, the experiments simulated the time-continuous assimilation of remotely sensed marine surface wind or temperature sounding data. The wind data were fabricated directly for model grid points intercepted by a Seasat-1 scatterometer swath and were assimilated into the lowest active level (945 mb) of the model using a localized successive correction method. It is shown that Seasat wind data can greatly improve numerical weather forecasts due to better definition of specific features. The case of the QE II storm is examined.

Cane, M. A.↗
NASA NTRS1993

Calculation methods for rotor wake formation in forward flight

The primary focus of this research centered on improving the vorticity-capturing properties of Eulerian flow solvers in convection-dominated flows. The work was motivated by the fact that excessive numerical diffusion observed in current Euler solvers on grids of feasible size is a major obstacle to accurately computing the flow about helicopter rotors, especially in the forward flight case. The work was originally supported by NASA Grant NAG2-421 for two years. The research then continued for four more years under an ONR contract. The research goal of reduction of numerical diffusion in Euler solvers remained the same, although the target application was now the prediction of internal vortical secondary flows. This was part of a larger ONR-sponsored research program at MIT aimed at reducing the noise of turbopumps. The research work culminated in the PhD thesis titled 'A Coupled Eulerian/Lagrangian Method for the Solution of Three-Dimensional Vortical Flows'.

Drela, Mark↗
NASA NTRS2001

Nonlinear Thermomechanical Response of Composite Panels with Continuous and Terminated Stiffeners

A two-phase approach and a computational procedure are used for predicting the variability of the response of stiffened composite panels associated with variations in the geometric and material parameters of the structures. In the first phase, hierarchical sensitivity analysis is used to identify the major parameters that have the most effect on the response quantities of interest. In the second phase, the major parameters are taken to be fuzzy parameters, and a fuzzy set analysis is used to determine the range of variation of the response, associated with preselected variations in the major parameters. Numerical results are presented showing the variability of the response of panels with both continuous and terminated stiffeners associated with variations in the micro mechanical and geometric parameters. Both flat and curved panels are considered.

Noor, Ahmed K.↗
NASA NTRS1993

A global analytical representation of the magnetic field produced by the region 2 Birkeland currents and the partial ring current

A quantitative model is developed of the magnetic field produced by the electric current system of region 2 Birkeland currents, closed via the partial ring current. The distribution of j-perpendicular is computed from a given axially asymmetric spatial distribution of hot isotropic magnetospheric plasma over an infinitely thin L shell in an axisymmetric purely dipolar geomagnetic field, while the field-aligned current density is found from the continuity of the net electric current. The magnetic field distribution is derived by a Biot-Savart integral over the electric current system. An assumed cosine dependence of the plasma pressure on local time makes it possible to reduce the problem of analytical representation of the B field to two dimensions. The obtained numerical fits for the partial ring current/region 2 Birkeland current magnetic field are relatively simple, continuous, and valid throughout the whole extraterrestrial space from ionospheric heights up to tens of earth radii.

Tsyganenko, N. A.↗
NASA NTRS1995

An Experimental and Theoretical Study of Radiative Extinction of Diffusion Flames

The objective of this research was to experimentally and theoretically investigate the radiation-induced extinction of gaseous diffusion flames in microgravity. The microgravity conditions were required because radiation-induced extinction is generally not possible in 1-g but is highly likely in microgravity. In 1-g, the flame-generated particulates (e.g. soot) and gaseous combustion products that are responsible for flame radiation, are swept away from the high temperature reaction zone by the buoyancy-induced flow and a steady state is developed. In microgravity, however, the absence of buoyancy-induced flow which transports the fuel and the oxidizer to the combustion zone and removes the hot combustion products from it enhances the flame radiation due to: (1) transient build-up of the combustion products in the flame zone which increases the gas radiation, and (2) longer residence time makes conditions appropriate for substantial amounts of soot to form which is usually responsible for most of the radiative heat loss. Numerical calculations conducted during the course of this work show that even non-radiative flames continue to become "weaker" (diminished burning rate per unit flame area) due to reduced rates of convective and diffusive transport. Thus, it was anticipated that radiative heat loss may eventually extinguish the already "weak" microgravity diffusion flame. While this hypothesis appears convincing and our numerical calculations support it, experiments for a long enough microgravity time could not be conducted during the course of this research to provide an experimental proof. Space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in microgravity will burn indefinitely. It was hoped that radiative extinction can be experimentally shown by the aerodynamically stabilized gaseous diffusion flames where the fuel supply rate was externally controlled. While substantial progress toward this goal was made during this project, identifying the experimental conditions for which radiative extinction occurs for various fuels requires further study. Details concerning this research which are discussed in published articles are included in the appendices.

Atreya, Arvind↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗
NASA NTRS2013

Numerical Simulation of the Anomalous Transport of High-Energy Cosmic Rays in Galactic Superbubble

A continuous-time random-walk (CTRW) model to simulate the transport and acceleration of high-energy cosmic rays in galactic superbubbles has recently been put forward (Barghouty & Schnee 2102). The new model has been developed to simulate and highlight signatures of anomalous transport on particles' evolution and their spectra in a multi-shock, collective acceleration context. The superbubble is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. This work concentrates on the effects of the bubble's assumed astrophysical characteristics (cf. geometry and roughness) on the particles' spectra.

Barghouty, A. F.↗
NASA NTRS1973

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗
NASA NTRS2017

Numerical Experiments in Error Control for Sound Propagation Using a Damping Layer Boundary Treatment

This paper presents results from numerical experiments for controlling the error caused by a damping layer boundary treatment when simulating the propagation of an acoustic signal from a continuous pressure source. The computations are with the 2D Linearized Euler Equations (LEE) for both a uniform mean flow and a steady parallel jet. The numerical experiments are with algorithms that are third, fifth, seventh and ninth order accurate in space and time. The numerical domain is enclosed in a damping layer boundary treatment. The damping is implemented in a time accurate manner, with simple polynomial damping profiles of second, fourth, sixth and eighth power. At the outer boundaries of the damping layer the propagating solution is uniformly set to zero. The complete boundary treatment is remarkably simple and intrinsically independant from the dimension of the spatial domain. The reported results show the relative effect on the error from the boundary treatment by varying the damping layer width, damping profile power, damping amplitude, propagtion time, grid resolution and algorithm order. The issue that is being addressed is not the accuracy of the numerical solution when compared to a mathematical solution, but the effect of the complete boundary treatment on the numerical solution, and to what degree the error in the numerical solution from the complete boundary treatment can be controlled. We report maximum relative absolute errors from just the boundary treatment that range from O[10-2] to O[10-7].

Algorithm↗
NASA NTRS2017

Numerical Experiments in Error Control for Sound Propagation Using a Damping Layer Boundary Treatment

This paper presents results from numerical experiments for controlling the error caused by a damping layer boundary treatment when simulating the propagation of an acoustic signal from a continuous pressure source. The computations are with the 2D Linearized Euler Equations (LEE) for both a uniform mean flow and a steady parallel jet. The numerical experiments are with algorithms that are third, fifth, seventh and ninth order accurate in space and time. The numerical domain is enclosed in a damping layer boundary treatment. The damping is implemented in a time accurate manner, with simple polynomial damping profiles of second, fourth, sixth and eighth power. At the outer boundaries of the damping layer the propagating solution is uniformly set to zero. The complete boundary treatment is remarkably simple and intrinsically independant from the dimension of the spatial domain. The reported results show the relative effect on the error from the boundary treatment by varying the damping layer width, damping profile power, damping amplitude, propagtion time, grid resolution and algorithm order. The issue that is being addressed is not the accuracy of the numerical solution when compared to a mathematical solution, but the effect of the complete boundary treatment on the numerical solution, and to what degree the error in the numerical solution from the complete boundary treatment can be controlled. We report maximum relative absolute errors from just the boundary treatment that range from O[10-2] to O[10-7].

Computational Aeroacoustics↗
NASA NTRS1985

Hydrodynamic stability problem formulated by numerical solutions of the Navier-Stokes equations

The problem of hydrodynamic stability and the transition from laminar to turbulent flows are reformulated by seeking numerical solutions of the full, unapproximated Navier-Stokes equations. This method differs significantly from the well known Orr-Sommerfeld equation approach. The oncoming laminar flow is disturbed by forced, time dependent perturbations. The magnitudes of these perturbations are arbitrary. Then, the ensuing spatial and temporal development of the imposed perturbations on the basic flow is calculated by direct numerical solutions of the time dependent Navier-Stokes equations. One of the main advantages of this method is its ability to simulate nonlinear processes. As a specific application of this technique to SSME (Space Shuttle Main Engine) flow configurations, computer programs have been written for the two dimensional flow over a backward facing step. This numerical code will be tested for operational use as part of a continued research collaboration effort with the NASA/MSFC counterparts.

Hyun, J. M.↗
NASA NTRS1987

Self-consistent inclusion of space-charge in the traveling wave tube

It is shown how the complete field of the electron beam may be incorporated into the transmission line model theory of the traveling wave tube (TWT). The fact that the longitudinal component of the field due to the bunched beam is not used when formulating the beam-to-circuit coupling equation is not well-known. The fundamental partial differential equation for the traveling wave field is developed and compared with the older (now standard) one. The equation can be solved numerically using the same algorithms, but now the coefficients can be updated continuously as the calculation proceeds down the tube. The coefficients in the older equations are primarily derived from preliminary measurements and some trial and error. The newer coefficients can be found by a recursive method, since each has a well defined physical interpretation and can be calculated once a reasonable first trial solution is postulated. The results of the new expression were compared with those of the older forms, as well as to a field theory model to show the ease in which a reasonable fit to the field prediction is obtained. A complete summary of the existing transmission line modeling of the TWT is given to explain the somewhat vague ideas and techniques in the general area of drifting carrier-traveling circuit wave interactions. The basic assumptions and inconsistencies of the existing theory and areas of confusion in the general literature are examined and hopefully cleared up.

Freeman, Jon C.↗
NASA NTRS1989

Analysis and modeling of summertime convective cloud and precipitation structure over the Southeastern United States

Data analysis and numerical modeling efforts that are related to the 1986 Satellite Precipitation and Cloud Experiment (SPACE) are discussed. Progress during this period includes the following: further testing and development of the RAMS numerical modeling system on the Alabama CRAY X-MP/24; a continuation of the observational analysis of the 13 July 1986 mesoscale convective system (MCS); and an initial investigation of a small MCS that formed over the COHMEX region on 15 July 1986. Details for each of these individual tasks are given.

Knupp, Kevin R.↗
NASA NTRS1989

The QBO and weak external forcing by solar activity: A three dimensional model study

A better understanding is attempted of the physical mechanisms leading to significant correlations between oscillations in the lower and middle stratosphere and solar variability associated with the sun's rotation. A global 3-d mechanistic model of the middle atmosphere is employed to investigate the effects of minor artificially induced perturbations. The aim is to explore the physical mechanisms of the dynamical response especially of the stratosphere to weak external forcing as it may result from UV flux changes due to solar rotation. First results of numerical experiments dealing about the external forcing of the middle atmosphere by solar activity were presented elsewhere. Different numerical studies regarding the excitation and propagation of weak perturbations have been continued since then. The model calculations presented are made to investigate the influence of the quasi-biennial oscillation (QBO) on the dynamical response of the middle atmosphere to weak perturbations by employing different initial wind fields which represent the west and east phase of the QBO.

Dameris, M.↗
NASA NTRS1992

Fractional-Step, Finite-Volume Computation Of Flow

Method of solving Navier-Stokes equations of incompressible flow in general nonorthogonal curvilinear coordinates incorporates fractional-step and finite-volume approaches. Developed in continuing effort to achieve accuracy without need for excessive computation time in numerical simulation of time-dependent, three-dimensional flows bounded by surfaces of complicated shape.

Kwak, Dochan↗
NASA NTRS1992

Ceramics and ceramic matrix composites - Aerospace potential and status

Thermostructural ceramics and ceramic-matrix composites are attractive in numerous aerospace applications; the noncatastrophic fracture behavior and flaw-insensitivity of continuous fiber-reinforced CMCs renders them especially desirable. The present development status evaluation notes that, for most highly-loaded high-temperature applications, the requisite fiber-technology base is at present insufficient. In addition to materials processing techniques, the life prediction and NDE methods are immature and require a projection of 15-20 years for the maturity of CMC turbine rotors. More lightly loaded, moderate temperature aircraft engine applications are approaching maturity.

Levine, Stanley R.↗
NASA NTRS1993

Trajectory optimization for the National aerospace plane

While continuing the application of the inverse dynamics approach in obtaining the optimal numerical solutions, the research during the past six months has been focused on the formulation and derivation of closed-form solutions for constrained hypersonic flight trajectories. Since it was found in the research of the first year that a dominant portion of the optimal ascent trajectory of the aerospace plane is constrained by dynamic pressure and heating constraints, the application of the analytical solutions significantly enhances the efficiency in trajectory optimization, provides a better insight to understanding of the trajectory and conceivably has great potential in guidance of the vehicle. Work of this period has been reported in four technical papers. Two of the papers were presented in the AIAA Guidance, Navigation, and Control Conference (Hilton Head, SC, August, 1992) and Fourth International Aerospace Planes Conference (Orlando, FL, December, 1992). The other two papers have been accepted for publication by Journal of Guidance, Control, and Dynamics, and will appear in 1993. This report briefly summarizes the work done in the past six months and work currently underway.

Lu, Ping↗