Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

An Evaluation of Architectural Platforms for Parallel Navier-Stokes Computations

We study the computational, communication, and scalability characteristics of a computational fluid dynamics application, which solves the time accurate flow field of a jet using the compressible Navier-Stokes equations, on a variety of parallel architecture platforms. The platforms chosen for this study are a cluster of workstations (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), and distributed memory multiprocessors with different topologies - the IBM SP and the Cray T3D. We investigate the impact of various networks connecting the cluster of workstations on the performance of the application and the overheads induced by popular message passing libraries used for parallelization. The work also highlights the importance of matching the memory bandwidth to the processor speed for good single processor performance. By studying the performance of an application on a variety of architectures, we are able to point out the strengths and weaknesses of each of the example computing platforms.

Jayasimha, D. N.↗
NASA NTRS1997

Parallelizing Navier-Stokes Computations on a Variety of Architectural Platforms

We study the computational, communication, and scalability characteristics of a Computational Fluid Dynamics application, which solves the time accurate flow field of a jet using the compressible Navier-Stokes equations, on a variety of parallel architectural platforms. The platforms chosen for this study are a cluster of workstations (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), distributed memory multiprocessors with different topologies-the IBM SP and the Cray T3D. We investigate the impact of various networks, connecting the cluster of workstations, on the performance of the application and the overheads induced by popular message passing libraries used for parallelization. The work also highlights the importance of matching the memory bandwidth to the processor speed for good single processor performance. By studying the performance of an application on a variety of architectures, we are able to point out the strengths and weaknesses of each of the example computing platforms.

Jayasimha, D. N.↗
NASA NTRS1999

Development of Formation Deployment and Intialization Concepts

NASA's Cross-Cutting Technology Development Program identified formation flying as a key enabler for the next generation Earth and Sciences campaign. It is hoped that this technology will allow a distributed network of autonomous satellites to act collaboratively as a single collective unit paving the way for extensive co-observing campaigns, coordinated multi-point observing programs, improved space-based interferometry, and entirely new approaches to conducting science. APL as a team member with GSFC, funded by the Earth Sciences and Technology Organization (ESTO), investigated formation deployment and initialization concepts which is central to the formation flying concept. This paper presents the analytical approach and preliminary results of the study. The study investigated a simple mission involving the deployment of six micro-satellites, one at a time, from a bus. At the initialization state, the satellites fly in an along-track trajectory separated by nominal spacing. The study entailed the development of a two-body (bus and satellite) relative motion propagator based on Clohessy-Wiltshire (C-W) equations with drag from which the relative motion of the micro-satellites is deduced. This code was used to investigate cluster development characteristics subject to "tip-off' (ejection) conditions. Results indicate that cluster development is very sensitive to the ballistic coefficients of the bus and satellites, and to relative ejection velocity. This information can be used to identify optimum deployment parameters, along with accuracy bounds for a particular mission, and to develop a cluster control strategy minimizing global fuel and cost. A suitable control strategy concept has been identified, however, it needs to be developed further.

Badesha, Surjit S.↗
NASA NTRS1999

Premixed-Gas Flame Propagation in Hele-Shaw Cells

It is well known that buoyancy and thermal expansion affect the propagation ra and shapes of premixed gas flames. The understanding of such effects is complicated by the large density ratio between the reactants and products, which induces a baroclinic production of vorticity due to misalignment of density and pressure gradients at the front, which in turn leads to a complicated multi-dimensional flame/flow interaction. The Hele-Shaw cell, i.e., the region between closely-spaced flat parallel plates, is probably the simplest system in which multi-dimensional convection is presents consequently, the behavior of fluids in this system has been studied extensively (Homsy, 1987). Probably the most important characteristic of Hele-Shaw flows is that when the Reynolds number based on gap width is sufficiently small, the Navier-Stokes equations averaged over the gap reduce to a linear relation, namely a Laplace equation for pressure (Darcy's law). In this work, flame propagation in Hele-Shaw cells is studied to obtain a better understanding of buoyancy and thermal expansion effects on premixed flames. This work is also relevant to the study of unburned hydrocarbon emissions produced by internal combustion engines since these emissions are largely a result of the partial burning or complete flame quenching in the narrow, annular gap called the "crevice volume" between the piston and cylinder walls (Heywood, 1988). A better understanding of how flames propagate in these volumes through experiments using Hele-Shaw cells could lead to identification of means to reduce these emissions.

Sharif, J.↗
NASA NTRS2003

Universal Responses of Cyclic-Oxidation Models Studied

Oxidation is an important degradation process for materials operating in the high-temperature air or oxygen environments typical of jet turbine or rocket engines. Reaction of the combustion gases with the component material forms surface layer scales during these oxidative exposures. Typically, the instantaneous rate of reaction is inversely proportional to the existing scale thickness, giving rise to parabolic kinetics. However, more realistic applications entail periodic startup and shutdown. Some scale spallation may occur upon cooling, resulting in loss of the protective diffusion barrier provided by a fully intact scale. Upon reheating, the component will experience accelerated oxidation due to this spallation. Cyclic-oxidation testing has, therefore, been a mainstay of characterization and performance ranking for high-temperature materials. Models simulate this process by calculating how a scale spalls upon cooling and regrows upon heating (refs. 1 to 3). Recently released NASA software (COSP for Windows) allows researchers to specify a uniform layer or discrete segments of spallation (ref. 4). Families of model curves exhibit consistent regularity and trends with input parameters, and characteristic features have been empirically described in terms of these parameters. Although much insight has been gained from experimental and model curves, no equation has been derived that can describe this behavior explicitly as functions of the key oxidation parameters.

Smialek, James L.↗
NASA NTRS1997

SYNTHESIS OF NOVEL ALL-DIELECTRIC GRATING FILTERS USING GENETIC ALGORITHMS

We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.

Zuffada, Cinzia↗
NASA NTRS2012

Eclipse-Free-Time Assessment Tool for IRIS

IRIS_EFT is a scientific simulation that can be used to perform an Eclipse-Free- Time (EFT) assessment of IRIS (Infrared Imaging Surveyor) mission orbits. EFT is defined to be those time intervals longer than one day during which the IRIS spacecraft is not in the Earth s shadow. Program IRIS_EFT implements a special perturbation of orbital motion to numerically integrate Cowell's form of the system of differential equations. Shadow conditions are predicted by embedding this integrator within Brent s method for finding the root of a nonlinear equation. The IRIS_EFT software models the effects of the following types of orbit perturbations on the long-term evolution and shadow characteristics of IRIS mission orbits. (1) Non-spherical Earth gravity, (2) Atmospheric drag, (3) Point-mass gravity of the Sun, and (4) Point-mass gravity of the Moon. The objective of this effort was to create an in-house computer program that would perform eclipse-free-time analysis. of candidate IRIS spacecraft mission orbits in an accurate and timely fashion. The software is a suite of Fortran subroutines and data files organized as a "computational" engine that is used to accurately predict the long-term orbit evolution of IRIS mission orbits while searching for Earth shadow conditions.

Eagle, David↗
NASA NTRS2016

Determination of the Earth's Plasmapause Location from the CE-3 EUVC Images

The Moon-based Extreme Ultraviolet Camera (EUVC) aboard China's Chang'e-3 (CE-3) mission has successfully imaged the entire Earth's plasmasphere for the first time from the side views on lunar surface. An EUVC image on 21 April 2014 is used in this study to demonstrate the characteristics and configurations of the Moon-based EUV imaging and to illustrate the determination algorithm of the plasmapause locations on the magnetic equator. The plasmapause locations determined from all the available EUVC images with the Minimum L Algorithm are quantitatively compared with those extracted from insitu observations (Defense Meteorological Satellite Program, Time History of Events and Macroscale Interactions during Substorms, and Radiation Belt Storm Probes). Excellent agreement between the determined plasmapauses seen by EUVC and the extracted ones from other satellites indicates the reliability of the Moon-based EUVC images as well as the determination algorithm. This preliminary study provides an important basis for future investigation of the dynamics of the plasmasphere with the Moon-based EUVC imaging.

He, Fei↗
NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. III: The Vector Radiative Transfer Equation

A vector radiative transfer equation with an additional source term typical of dense media is obtained. The analysis includes (i) the derivation of an integral equation for the correlation matrix of the exciting field coefficients accounting for the correlation between the particles, (ii) the derivation of an integral representation for the specific coherency dyadicin terms of this matrix, and (iii) the simplification of the integral equation for the correlation matrix and of the integral representation for the specific coherency dyadic by employing a series of approximations which are characteristic of sparse media.

Adrian Doicu↗
NASA NTRS1977

Dynamic characteristics of rotor blades with pendulum absorbers

The point transmission matrix for a vertical plane pendulum on a rotating blade undergoing combined flapwise bending, and chordwise bending and torsion is derived. The equilibrium equation of the pendulum is linearized for small oscillations about the steady state. A FORTRAN program was written for the case of a vertical plane pendulum attached to a uniform blade with flapwise bending degree of freedom for cantilever boundary conditions. The frequency has a singular value right at the uncoupled pendulum natural frequency and thus introduces two frequencies corresponding to the nearest natural frequency of the blade without pendulum. In both of these modes it was observed that the pendulum deflection is large. One frequency can be thought of as a coupled pendulum frequency and the other as a coupled bending and pendulum frequency.

Murthy, V. R.↗
DOE OSTI2022

Size Effect of Local Current-Voltage Characteristics of MX 2 Nanoflakes: Local Density of States Reconstruction from Scanning Tunneling Microscopy Experiments

Local current-voltage characteristics for low-dimensional transition-metal dichalcogenides (LDTMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments, are of fundamental interest and can be useful for advanced applications. Most of the existing models either have limited applicability for complex-shaped LDTMDs (e.g., those based on the Simmons approach) or require solving of an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on the Tersoff approach). Using a serial expansion of the Tersoff formulas, we propose a flexible method to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We establish a set of key physical parameters, which characterize the tunneling current of a STM-probe–sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different numbers of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around the Fermi level, which allows us to gain insight into the local band structure of LDTMDs. The reconstructed LDOS reveal pronounced size effects for the single-layer, two-layer, and three-layer MoS 2 nanoflakes, which we relate to the low dimensionality and strong bending or corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach, allowing LDOS reconstruction, will be of critical interest for the quantitative description of STM experiments and also of use to better understand the microscopic physical aspects of the surface, strain, and bending contributions to the LDTMDs’ electronic properties.

42 ENGINEERING↗
NASA NTRS2006

Dynamic Response Testing in an Electrically Heated Reactor Test Facility

Non-nuclear testing can be a valuable tool in development of a space nuclear power or propulsion system. In a non-nuclear test bed, electric heaters are used to simulate the heat from nuclear fuel. Standard testing allows one to fully assess thermal, heat transfer, and stress related attributes of a given system, but fails to demonstrate the dynamic response that would be present in an integrated, fueled reactor system. The integration of thermal hydraulic hardware tests with simulated neutronic response provides a bridge between electrically heated testing and full nuclear testing. By implementing a neutronic response model to simulate the dynamic response that would be expected in a fueled reactor system, one can better understand system integration issues, characterize integrated system response times and response characteristics, and assess potential design improvements at a relatively small fiscal investment. Initial system dynamic response testing was demonstrated on the integrated SAFE-100a heat pipe cooled, electrically heated reactor and heat exchanger hardware, utilizing a one-group solution to the point kinetics equations to simulate the expected neutronic response of the system (Bragg-Sitton, 2005). The current paper applies the same testing methodology to a direct drive gas cooled reactor system, demonstrating the applicability of the testing methodology to any reactor type and demonstrating the variation in system response characteristics in different reactor concepts. In each testing application, core power transients were controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. Although both system designs utilize a fast spectrum reactor, the method of cooling the reactor differs significantly, leading to a variable system response that can be demonstrated and assessed in a non-nuclear test facility.

Bragg-Sitton, Shannon M.↗
NASA NTRS1987

The effect of along-front temperature variation in a two-dimensional primitive equation model of a surface frontogenesis

The effects of along-front temperature variations on the structure and evolution of surface fronts are investigated using a two-dimensional frontogenesis model. The formulation and characteristics of the confluence-horizontal shear model are described. The structures of three frontal zones forced by confluence and characterized by zero, negative, and positive along-front potential temperature variations are analyzed; particular attention is given to the ageostrophic vertical circulation and their forcings. The temporal evolution of the model frontal zones is examined in terms of the low-level maxima of vorticity, cross-front potential temperature gradient, and convergence. The kinematic and dynamic processes influencing the frontogenesis are discussed. It is noted that the two-dimensional model is applicable to the study of the dynamical properties of frontogenesis.

Keyser, Daniel↗
NASA NTRS1995

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗
NASA NTRS1953

On the application of transonic similarity rules to wings of finite span

The transonic aerodynamic characteristics of wings of finite span are discussed from the point of view of a unified small perturbation theory for subsonic, transonic, and supersonic flows about thin wings. This approach avoids certain ambiguities which appear if one studies transonic flows by means of equations derived under the more restrictive assumption that the local velocities are everywhere close to sonic velocity. The relation between the two methods of analysis of transonic flow is examined, the similarity rules and known solutions of transonic flow theory are reviewed, and the asymptotic behavior of the lift, drag, and pitching-moment characteristics of wings of large and small aspect ratio is discussed. It is shown that certain methods of data presentation are advantageous for the effective display of these characteristics.

Spreiter, John R↗
DOE OSTI2020

Axionic dark matter halos in the gravitational field of baryonic matter

We consider a dark matter halo (DMH) of a spherical galaxy as a Bose–Einstein condensate (BEC) of the ultra-light axions (ULA) interacting with the baryonic matter. In the mean-field (MF) limit, we have derived the integro-differential equation of the Hartree–Fock type for the spherically symmetrical wave function of the DMH component. This equation includes two independent dimensionless parameters: (i) β is the ratio of baryon and axion total mases and (ii) ξ is the ratio of characteristic baryon and axion spatial parameters. We extended our “dissipation algorithm” for studying numerically the ground state of the axion halo in the gravitational field produced by the baryonic component. We calculated the characteristic size, x c of DMH as a function of β and ξ and obtained an analytical approximation for x c .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Combustion of a premixed system in stagnation flow. I - Theoretical.

The theoretical steady-state and transient characteristics of a premixed combustible system in laminar axisymmetric stagnation flow are studied by way of numerical solution of the governing boundary-layer equations. Computations, which are based on the oxidation reaction of carbon monoxide in humid air, invoke a minimum of simplifying assumptions. Steady-state results include prediction of (1) profiles of the state variables, (2) extinction of blow-off conditions, (3) ignition conditions, and (4) the heat flux to the stagnation point. Computer simulations of the unsteady state provide an indication of the stability of a given steady state to perturbations. Comparisons between some of the predictions of the mathematical model and the results of laboratory experiments will be reported in Part II.

Smith, H. W.↗
NASA NTRS1973

Development and application of a local linearization algorithm for the integration of quaternion rate equations in real-time flight simulation problems

High angular rates encountered in real-time flight simulation problems may require a more stable and accurate integration method than the classical methods normally used. A study was made to develop a general local linearization procedure of integrating dynamic system equations when using a digital computer in real-time. The procedure is specifically applied to the integration of the quaternion rate equations. For this application, results are compared to a classical second-order method. The local linearization approach is shown to have desirable stability characteristics and gives significant improvement in accuracy over the classical second-order integration methods.

Barker, L. E., Jr.↗