Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Quasi-axisymmetric circulation and superrotation in planetary atmospheres

Atmospheric superrotation is ubiquitous in the solar system, and numerous theoretical models have been proposed to describe the effect. However, no model is fully accepted. Problems are primarily related to an understanding of the angular momentum and the energy budgets. In connection with the present investigation, attention is given to hybrid models, where three-dimensional models may describe eddy processes and two-dimensional symmetric models may describe the global scale circulation in the framework of a viscous medium. It is believed that the concept of an axisymmetric circulation is still useful. It is the main purpose of the current investigation to explore some of the properties of the circulation in planetary atmospheres. It is assumed that the eddy transports are downgradient or diffusive in nature. Energy and momentum are perceived to cascade from the global scale circulation driven by the sun to the smaller scale eddies which in turn determine the diffusive properties. A numerical model of the Venusian circulation is presented, and results are discussed for different conditions of static stability.

Mayr, H. G.↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

A defining feature of the discontinuous Galerkin (DG) method for ODE is that the piecewise polynomial solution can have a jump discontinuity at the beginning of each step. Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠 − 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK-DG, collocation, and CG schemes.

numerical methods↗
NASA NTRS2013

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings↗
NASA NTRS1988

Numerical investigations in three-dimensional internal flows

Renewed interest in simulation of high enthalpy flows was brought about by interest in the design and testing of hypersonic propulsion components. This renewed interest has led to refurbishment of high speed ground based wind tunnel facilities. Simulation of the high enthalpy flows occurring in the NASP propulsion system is of particular interest. Studying the combustor presentes some of the more challenging problems in the ground based simulation of the fluid flow. A difficulty in simulating the required conditions at the entrance to the combustor is related to achieveing a sufficiently high static pressure such that combustion would be possible at a high Mach number. An analysis of the flow within the original nozzle was not performed previously using numerical simulations, although it was operated experimentally. A Navier-Stokes code was chosen in the present study to examine analytically the flow within representative selected nozzles for use in combustion experiments. The numerical simulation of these fluid flows is popular and may be able to shed light on the direction that redesign of the shock tube nozzles should take. The numerical codes were modified to account for the axially symmetric flow field, for the inclusion of equilibrium real gas effects and for ideal gases whose thermodynamic properties are frozen at the stagnation conditions. Numerical simulation of flows within the original and two proposed nozzles, including the energy lost to the relatively cold walls, are investigated.

Rose, William C.↗
NASA NTRS1972

The interstellar line spectra of zeta Ophiuchi and zeta Persei and their relation to the short wavelength microwave background radiation

Thirty-one high dispersion Coude spectrograms of zeta Ophiuchi and seven of zeta Persei were numerically synthesized to produce high resolution, low noise spectra in the interval 3650 A to 4350 that yield data on atomic and molecular absorption in well-defined regions of the interstellar medium. The detection threshold is improved by as much as a factor 5 over single plates. Several interstellar lines were discovered in the zeta Oph - 15km/sec cloud and the zeta Per + 13 km/sec cloud.

Bortolot, V. J., Jr.↗
NASA NTRS1972

Investigation of Some Characteristics Related to PCM Thermal Capacitors

In order to provide greater flexibility, shorter computer time, and a parametric study, attention was devoted to an implicit numerical scheme. The experimental work involved outlining and monitoring approximately fourteen tests of an invertable capacitor. Surface and PCM temperatures were measured at several locations and photographs were made of the PCM during melting and solidification. Some typical experimental results are presented and discussed.

Griggs, E. I.↗
NASA NTRS1974

Seasonal and magnetic storm related changes in the thermosphere induced by eddy mixing

The possibility of explaining the seasonal and magnetic storm related changes in the thermosphere and the ionosphere through variation in eddy diffusion coefficient in the lower thermosphere is investigated theoretically by obtaining simultaneous numerical solutions of the relevant continuity, momentum, and energy balance equations in a self-consistent manner. It is shown that various important features of thermospheric seasonal behavior, including the winter helium bulge phenomenon, can be explained by assuming effective reduction in the winter-time eddy diffusion coefficient by about one order of magnitude. Storm related changes in the thermal structure and composition are seen to arise by assuming a downward shift of turbulence. Physical implications and relative significance of variation in thermospheric turbulence and large scale inter-hemispheric circulation are briefly discussed.

Sinha, A. K.↗
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

Stability, accuracy, and efficiency of some underintegrated methods in finite element computations

In an attempt to increase computational efficiency in the numerical solution of highly nonlinear problems in solid and fluid mechanics, underintegrated finite element methods have been employed by many analysts. Underintegration refers to the use of a rule of an order lower than that required to integrate polynomial integrands exactly. The main drawback of this technique is related to the production of rank-deficient stiffness matrices, or equivalently an expanded kernel of the governing linear momentum operators. Such a development can introduce numerical instabilities. In order to overcome this difficulty, artificial stiffness or viscosity methods, or other stabilization methods have been proposed. One approach involves the elimination of spurious modes in a postprocessing operation. The present study is concerned with this a posteriori elimination method, taking into account the results which can be expected from it, and some of its possible extensions.

Jacquotte, O.-P.↗
NASA NTRS2016

An Introduction to Transient Engine Applications Using the Numerical Propulsion System Simulation (NPSS) and MATLAB

This document outlines methodologies designed to improve the interface between the Numerical Propulsion System Simulation framework and various control and dynamic analyses developed in the Matlab and Simulink environment. Although NPSS is most commonly used for steady-state modeling, this paper is intended to supplement the relatively sparse documentation on it's transient analysis functionality. Matlab has become an extremely popular engineering environment, and better methodologies are necessary to develop tools that leverage the benefits of these disparate frameworks. Transient analysis is not a new feature of the Numerical Propulsion System Simulation (NPSS), but transient considerations are becoming more pertinent as multidisciplinary trade-offs begin to play a larger role in advanced engine designs. This paper serves to supplement the relatively sparse documentation on transient modeling and cover the budding convergence between NPSS and Matlab based modeling toolsets. The following sections explore various design patterns to rapidly develop transient models. Each approach starts with a base model built with NPSS, and assumes the reader already has a basic understanding of how to construct a steady-state model. The second half of the paper focuses on further enhancements required to subsequently interface NPSS with Matlab codes. The first method being the simplest and most straightforward but performance constrained, and the last being the most abstract. These methods aren't mutually exclusive and the specific implementation details could vary greatly based on the designer's discretion. Basic recommendations are provided to organize model logic in a format most easily amenable to integration with existing Matlab control toolsets.

Matlab↗
NASA NTRS1990

Observation simulation experiments with regional prediction models

Research efforts in FY 1990 included studies employing regional scale numerical models as aids in evaluating potential contributions of specific satellite observing systems (current and future) to numerical prediction. One study involves Observing System Simulation Experiments (OSSEs) which mimic operational initialization/forecast cycles but incorporate simulated Advanced Microwave Sounding Unit (AMSU) radiances as input data. The objective of this and related studies is to anticipate the potential value of data from these satellite systems, and develop applications of remotely sensed data for the benefit of short range forecasts. Techniques are also being used that rely on numerical model-based synthetic satellite radiances to interpret the information content of various types of remotely sensed image and sounding products. With this approach, evolution of simulated channel radiance image features can be directly interpreted in terms of the atmospheric dynamical processes depicted by a model. Progress is being made in a study using the internal consistency of a regional prediction model to simplify the assessment of forced diabatic heating and moisture initialization in reducing model spinup times. Techniques for model initialization are being examined, with focus on implications for potential applications of remote microwave observations, including AMSU and Special Sensor Microwave Imager (SSM/I), in shortening model spinup time for regional prediction.

Diak, George↗
NASA NTRS1999

Lymphocyte Functions in Microgravity

To understand the mechanism of immunity impairment in space it is important to analyze the direct effects of space-related conditions on different lymphocytes functions. Since 1992, we are investigating the effect of modeled and true microgravity (MG) on numerous lymphocyte functions. We had shown that modeled (MMG) and true microgravity inhibit lymphocyte locomotion through type I collagen. Modeled microgravity also suppresses polyclonal and antigen-specific lymphocyte activation. Polyclonal activation of lymphocytes prior to exposure to MMG abrogates the MG-induced inhibition of lymphocyte locomotion. The relationship between activation deficits and the loss of locomotion in MG was investigated using PKC activation by phorbol ester (PMA) and calcium ionophore (ionomycin). Direct activation of PKC by PMA substantially restored the MMG-inhibited lymphocyte locomotion and PHA-induced lymphocyte activation lonomycin by itself did not restore either locomotion or activation of the lymphocytes, indicating that these changes are not related to the impairment in the calcium flux in MMG. Treatment of lymphocytes with PMA before exposure to MMG prevented the loss of locomotion. It was observed that DNA synthesis is not necessary for restoration of locomotion since mitomicin C treated and untreated cells recovered their locomotion to the same level after PKC activation. Our recent data indicate that microgravity may selectively effect the expression of novel Ca2+ independent isoforms of PKC, in particularly PKC sigma and delta. This provides a new insight in understanding of the mechanisms of MG-sensitive cellular functions.

Pellis, Neal R.↗
NASA NTRS1977

Post-Newtonian gravitational bremsstrahlung

Formulae and numerical results are presented for the gravitational radiation emitted during a low-deflection encounter between two massive bodies. Results are valid through post-Newtonian order within general relativity. The gravitational waveform, the total luminosity and total emitted energy, the angular distribution of emitted energy, and the frequency spectrum are discussed in detail. A method boosting the accuracy of these quantities to post Newtonian order is also presented. A numerical comparison of results with those of Peters, and of Kovacs and Thorne shows that the post Newtonian method is reliable to better than 0.1 percent at v = 0.1 c, to a few percent at v = 0.35 c, and to 10 to 20 percent at v = 0.5 c.

Turner, M.↗
NASA NTRS1987

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗
NASA NTRS1972

Ground effect for V/STOL aircraft configurations and its simulation in the wind tunnel. Part 2: Experimental studies

Wind tunnel tests on a finite, knee-blowing-flapped wing and on a direct jet lift configuration were performed in a 30 x 42 inch wind tunnel. The objectives of these tests were, in experimental order:(1) to obtain definitive fixed-versus-moving ground comparisons on a powered, high-lift finite wing and on a lifting jet configuration as a function of model height and lift, (2) to understand the flow mechanism which lead to observed differences between moving and fixed-ground results, (3) to discover whether the above effects of a moving ground plane can be simulated by tangential-blowing boundary layer control at the wind tunnel floor, and (4) to check the form and the numerical constants for the wall jet blowing equation derived theoretically in Part I of this report. This relates wall jet blowing quantities to model lift coefficient.

Hackett, J. E.↗
NASA NTRS1974

A study of the temporal changes recorded by ERTS and their geological significance

The temporal changes that are recorded by ERTS were evaluated for an area around Bathurst Inlet in the North West Territories. The seasons represented by the images included: early winter, spring, early summer, summer, and fall. Numerous surface characteristics (vegetation, drainage patterns, surface texture, lineament systems and topographic relief, etc.) were used to relate the change in observable features with the different seasons. It was found that the time of year when an observation is made has a strong control over the amount and type of information that can be derived by an experienced interpreter. It was concluded that a detailed study of temporal changes is an important part of any ERTS interpretation for geology.

Moore, H. D.↗
NASA NTRS1975

Interpolation in numerical optimization

The present work discusses the generation of the cubic-spline interpolator in numerical optimization methods which use a variable-step integrator with step size control based on local relative truncation error. An algorithm for generating the cubic spline with successive over-relaxation is presented which represents an improvement over that given by Ralston and Wilf (1967). Rewriting the code reduces the number of N-vectors from eight to one. The algorithm is formulated in such a way that the solution of the linear system set up yields the first derivatives at the nodal points. This method is as accurate as other schemes but requires the minimum amount of storage.

Hall, K. R.↗