Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential time integration”

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 73 records · Page 4

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗

Computation of aerodynamic interference effects on oscillating airfoils with controls in ventilated subsonic wind tunnels

Lift interference effects are discussed based on Bland's (1968) integral equation. A mathematical existence theory is utilized for which convergence of the numerical method has been proved for general (square-integrable) downwashes. Airloads are computed using orthogonal airfoil polynomial pairs in conjunction with a collocation method which is numerically equivalent to Galerkin's method and complex least squares. Convergence exhibits exponentially decreasing error with the number n of collocation points for smooth downwashes, whereas errors are proportional to 1/n for discontinuous downwashes. The latter can be reduced to 1/n to the m+1 power with mth-order Richardson extrapolation (by using m = 2, hundredfold error reductions were obtained with only a 13% increase of computer time). Numerical results are presented showing acoustic resonance, as well as the effect of Mach number, ventilation, height-to-chord ratio, and mode shape on wind-tunnel interference. Excellent agreement with experiment is obtained in steady flow, and good agreement is obtained for unsteady flow.

Fromme, J. A.↗

Probabilistic Processing of Asynchronous Geiger-Mode Avalanche Photodiode Arrays for Background Rate Measurements

A system-level performance evaluation of Geiger-mode avalanche photodiode (GmAPD) arrays requires accurate measurement and prediction of the background rate of the device due to dark counts and other spurious detection events. Since a GmAPD detector reports only a binary value and timestamp associated with an avalanche event, dark count rates are typically measured by averaging thousands of frames to support a statistically significant measurement. For both synchronous and asynchronous detector, the Poisson distributed background rates are referenced to the time each pixel is armed. Unlike for synchronous GmAPD imagers where all the pixels are armed to an array-wide arm signal, an asynchronous pixel operates independently from its neighboring pixels; requiring the background rates to be calculated using an interarrival histogram. For both types of imagers, the background rate is typically evaluated by fitting an exponential distribution to a fixed window within a measured histogram of time intervals between detection events However, if the statistics of the background rate are insufficient – whether that is due to low population sizes, saturation, or a large dynamic range of population size across the array, the pixel, or array-wide, performance metrics may report results with varying accuracy. This paper reports on an implementation of an algorithm that evaluates GmAPD background rates based on statistical metrics rather than fixed windows. The algorithm functions by determining the appropriate integration window within the interarrival time histogram based on a per-pixel count rate set by a predetermined tolerable measurement error. The implementation of the algorithm allows us to characterize GmAPD arrays with orders of magnitude spread in background rates across the detector using common statistical parameters.

Mumolo, J.↗

Some simple properties of stellar pulsation modes

Except for the lowest harmonics, small-amplitude stellar pulsation modes possess many simple properties whose evaluation does not require numerical integration of the fourth-order equations of motion. All antinodes tend to have the same total kinetic energy except for those lying near physical or geometric boundaries. However, when kinetic energy per unit volume is considered, order-of-magnitude enhancements are seen in antinodes lying near the center of the star, and factor-of-2 enhancements occur near the polar axis. The nodes are distributed very regularly along the radius. They follow an exponential law in g-regions, and their separation is proportional to the sound travel time in p-regions. A simple graphical procedure is described for surveying the oscillation frequencies of a new stellar model. A precise condition is derived giving the division of energy between radial and angular motion. Another condition gives the fractional contribution to the velocity field of its two sources, the divergence and the curl. Certain simplifying results of weak coupling among the linear modes are briefly described.

Wolff, C. L.↗

Simulation Propulsion System and Trajectory Optimization

A number of new aircraft concepts have recently been proposed which tightly couple the propulsion system design and operation with the overall vehicle design and performance characteristics. These concepts include propulsion technology such as boundary layer ingestion, hybrid electric propulsion systems, distributed propulsion systems and variable cycle engines. Initial studies examining these concepts have typically used a traditional decoupled approach to aircraft design where the aerodynamics and propulsion designs are done a-priori and tabular data is used to provide inexpensive look ups to the trajectory ana-ysis. However the cost of generating the tabular data begins to grow exponentially when newer aircraft concepts require consideration of additional operational parameters such as multiple throttle settings, angle-of-attack effects on the propulsion system, or propulsion throttle setting effects on aerodynamics. This paper proposes a new modeling approach that eliminated the need to generate tabular data, instead allowing an expensive propulsion or aerodynamic analysis to be directly integrated into the trajectory analysis model and the entire design problem optimized in a fully coupled manner. The new method is demonstrated by implementing a canonical optimal control problem, the F-4 minimum time-to-climb trajectory optimization using three relatively new analysis tools: Open M-DAO, PyCycle and Pointer. Pycycle and Pointer both provide analytic derivatives and Open MDAO enables the two tools to be combined into a coupled model that can be run in an efficient parallel manner that helps to cost the increased cost of the more expensive propulsion analysis. Results generated with this model serve as a validation of the tightly coupled design method and guide future studies to examine aircraft concepts with more complex operational dependencies for the aerodynamic and propulsion models.

Optimization↗

Autonomous Task Management and Decision Support Tools

For some time aircraft manufacturers and researchers have been pursuing mechanisms for reducing crew workload and providing better decision support to the pilots, especially during non-normal situations. Some previous attempts to develop task managers or pilot decision support tools have not resulted in robust and fully functional systems. However, the increasing sophistication of sensors and automated reasoners, and the exponential surge in the amount of digital data that is now available create a ripe environment for the development of a robust, dynamic, task manager and decision support tool that is context sensitive and integrates information from a wide array of on-board and off aircraft sourcesa tool that monitors systems and the overall flight situation, anticipates information needs, prioritizes tasks appropriately, keeps pilots well informed, and is nimble and able to adapt to changing circumstances. This presentation will discuss the many significant challenges and issues associated with the development and functionality of such a system for use on the aircraft flight deck.

decision support↗

Measurements of Cyclotron Features and Pulse Periods in the High-Mass X-Ray Binaries 4U 1538-522 and 4U 1907+09 with the International Gamma-Ray Astrophysics Laboratory

We present a spectral and timing analysis of International Gamma-Ray Astrophysics Laboratory (INTEGRAL) observations of two high-mass X-ray binaries, 4U 1538−522 and 4U 1907+09. Our timing measurements for 4U 1538−522 find the pulse period to have exhibited a spin-up trend until approximately 2009, after which there is evidence for a torque reversal, with the source beginning to spin down to the most recently measured period of 525.407 plus or minus 0.001 seconds. The most recent INTEGRAL observations of 4U 1907+09 are not found to yield statistically significant pulse periods due to the significantly lower flux from the source compared with 4U 1538−522. A spectral model consisting of a power-law continuum with an exponential cutoff and modified by two cyclotron resonance scattering features is found to fit both sources well, with the cyclotron scattering features detected at approximately 22 and approximately 49 kiloelectronvolts for 4U 1538−522 and at approximately 18 and approximately 36 kiloelectronvolts for 4U 1907+09. The spectral parameters of 4U 1538−522 are generally not found to vary significantly with flux and there is little to no variation across the torque reversal. Examining our results in conjunction with previous work, we find no evidence for a correlation between cyclotron line energy and luminosity for 4U 1538−522. 4U 1907+09 shows evidence for a positive correlation between cyclotron line energy and luminosity, which would make it the fourth, and lowest luminosity, cyclotron line source to exhibit this relationship.

stars: magnetic field↗

The Cooling Loop A Anomaly of 2013: A Case Study in Human-Systems Resilience

Throughout the history of human spaceflight, NASA has employed an operational paradigm of 24/7 dependence on experts in Mission Control Center (MCC). In addition to nominal flight control and mission operations, these 85+ experts per shift manage anomaly detection, diagnosis, and response, and support the crew in real-time in performing maintenance and repair, procedure execution, and other complex mission operations. Future long-duration exploration missions (LDEMs) beyond low-Earth orbit (LEO) will not operate successfully using this same Human-Systems Integration Architecture (HSIA) where crew rely on ground controllers, have ready access to resupply, and have a fallback plan of evacuation. As distance from Earth increases and the communication delay grows, crews will need to respond independently and adequately to time-critical vehicle malfunctions. It will not always be sufficient or even possible to ‘safe the system’ and then wait upon ground intervention. A new and radically different HSIA is needed to accommodate the paradigm shift of deep-space travel. Historical International Space Station (ISS) data show that for a 30-day mission, the likelihood of a high-consequence vehicle anomaly of uncertain origin that requires rapid response is greater than 10%. The likelihood of such an event is 50% by the fourth month of the mission, and it grows exponentially with time. Our team has conducted in-depth investigations into these events and their corresponding anomaly resolution activities. Using MCC and Mission Evaluation Room (MER) anomaly resolution artifacts (including meeting summaries, caution and warning data, and ISS daily summaries), we created timelines detailing ground actions and in-orbit events for two significant anomalies. We then mapped these timelines onto Mars transit conditions, introducing a ground-crew communications time delay and shifting immediate response, time-critical task execution, and vehicle commanding to the crew. In detailing successful anomaly resolution in transit to Mars, the timelines highlight where effective resolution requires drastically evolved onboard capabilities. Though this research has yielded a rich data set based on ground response in past missions, there is still insufficient knowledge to assess the potential impact of inflight anomalies on a small autonomous crew on future LDEMs beyond LEO. To begin building an evidence base that will inform future HSIA standards and requirements, we are developing an approach to systematically capture crew anomaly response and procedure execution during early Artemis missions. Being the first human spaceflight beyond LEO since Apollo, early Artemis missions provide a rare and unique opportunity to serve as a testbed for Mars missions. Our work aims to capitalize on planned data collection to derive crew operational responses to anomalous events in real-time. Our team is also researching the level of simulation fidelity required for empirically validating proposed HSIA standards and evaluating HSIA implementations for LDEMs beyond LEO. This work will produce a trade space study of HSIA simulation objectives and fidelity requirements. Ultimately, these research efforts will assist in developing the standards and technologies needed to build a next-generation HSIA for LDEMs beyond LEO.

human-systems integration architecture↗

Variation in light intensity with height and time from subsequent lightning return strokes

Photographic measurements of relative light intensity as a function of height and time have been conducted for seven return strokes in two lightning flashes at 7.8 and 8.7 km ranges, using film which possesses an approximately constant spectral response in the 300-670 nm range. The amplitude of the initial light peak is noted to decrease exponentially with height, with a decay constant of 0.6-0.8 km. The logarithm of the peak light intensity near the ground is found to be approximately proportional to the initial peak electric field intensity, implying that the current decrease with height may be much slower than the light decrease. Absolute light intensity is presently estimated through the integration of the photographic signals from individual channel segments, in order to simulate the calibrated, all-sky photoelectric data of Guo and Krider (1982).

Jordan, D. M.↗

Origin, evolution and present thermal state of the moon.

The relative absence of lunar volcanism in the last 3,000,000,000 years and the Apollo 15 heat flow measurement suggest that present-day temperatures in the moon are approximately steady-state to depths of about 100 km. An exponential distribution of heat sources with depth may then be scaled by equating the surface heat flow to the integrated heat production of this exterior shell. Presumed present-day interior temperatures, as well as the present-day surface heat flow of about 30 erg/sq cm-sec, may be obtained with an initial temperature roughly corresponding to the Apollo 11 basalt solidus and exponential scaling of heat sources. The concentrations of U for an originally homogeneous moon are estimated to be about 0.009 micrograms per gram, close to that measured for eucrites and inferred for primitive inclusions of the Allende meteorite. The estimated homogeneous concentrations of U, the chemistry of the lunar surface material and inferences to modest depth, and the short accretion time of the moon necessary to provide large-scale differentiation at 4.6 AE suggest that the moon had its origin in the rapid accretion of compounds first condensing from the protoplanetary nebula.

Hanks, T. C.↗

Development of OPTIMOD - The new Space Shuttle Orbiter/Payload Thermal Integration Model

The development of the Orbiter/Payload Thermal Integration Model (OPTIMOD), a simplified Space Shuttle Orbiter thermal math model to be used for payload thermal integration analysis, is described. The Orbiter lower fuselage math model has been reduced to the minimum level of detail consistent with accuracy requirements. Instead, an 'environment radiation sensor', a simple closed-end cylinder oriented in a fixed position with respect to the Orbiter, is introduced. Equations are presented which relate the cylinder surface heating rates to external fuselage steady-state temperatures. These equations were derived from a regression correlation analysis of steady-state Orbiter thermal solutions for 2372 points on orbit, representing a complete range of vehicle-earth-sun orientations. An exponential equation is developed to relate the steady-state data to time-variant temperatures for the particular transient attitude timeline being simulated. The accuracy of OPTIMOD is shown to be excellent in the payload bay and to compare favorably with previous models in Orbiter elements outside the bay.

Russell, David J.↗

A Diagnostic Analysis of the Kennedy Space Center LDAR Network: Data Characteristics - 1

An analytical framework is developed in which to analyze climatological VHF (66 MHz) radiation measurements taken by the Kennedy Space Center Lightning Detection and Ranging (LDAR) network. A 19-month noise-filtered sample of LDAR observations is examined using this framework. It is found that the climatological VHF source density as observed by LDAR falls off approximately 10 dB every 71 km of ground range away from the network centroid (a 31 km e-folding scale). From this framework it is inferred that the underlying source distribution is likely inverse exponential in amplitude, with a log-slope of 126-147 mW(sup -1/2). The underlying vertical distribution of VHF sources is approximately normally distributed with a mean altitude of 9 km and a standard deviation of 2.7 km; this implies that the loss of below-horizon sources has a negligible effect on column-integrated source densities within 200 km ground range. At medium to far ranges, location errors are primarily radial and have a slightly asymmetric distribution whose first moment increases as range squared. Error moments estimated from observed lightning are significantly higher than those from aircraft-based signal generator or analytic solution estimates, suggesting that timing errors arising from poor signal identification and discrimination may dominate over timing errors arising from nominal sensor resolution. The VHF source properties of individual LDAR-observed flashes are computed and an analytic expression for flash detection efficiency vs. range is derived. This reveals nearly constant flash detection efficiency to 80-94 km range from the network centroid.

Boccippio, D. J.↗

In-Space Transportation Sensitivity to Roundtrip Mission Duration and Mars Vicinity Stay Time

The National Aeronautics and Space Administration’s Exploration Systems Development Mission Directorate has been developing architecture concepts for human missions to Mars in alignment with the agency’s Moon-to-Mars Strategy & Objectives. One of the key components of a human Mars mission is the in-space transportation system that delivers crew and cargo to Mars vicinity and returns the crew safely back to Earth. The Mars Architecture Team within ESDMD’s Strategy and Architecture Office has been evaluating multiple in-space transportation options to satisfy this functional need. The first step of this evaluation is to understand key architectural drivers for in-space transportation systems. Lunar and Mars missions have unique challenges, and systems designed for one may not be directly applicable to the other as they have different energy and mission needs. In addition, the time and distance for a Mars mission will be drastically different from the historical and current programs, requiring a new paradigm for mission architects. This paper investigates the impact of mission duration on a variety of different transportation options. As mission duration decreases, the total energy required to perform a roundtrip Mars mission increases exponentially, significantly increasing the required Earth departure mass. Sensitivity analyses of various in-space transportation system concepts were shown in this paper to begin mapping out the in-space transportation trade space. Each of the transportation systems has its own strengths and weaknesses, especially in the context of the continuous mission duration and vicinity stay time trade space. Although mission duration and Mars vicinity stay time may be key architecture drivers from an in-space transportation perspective, the decision of acceptable mission duration cannot be made solely from this perspective. Mission duration also has tangible impacts to crew health and performance, as well as system reliability and other key architecture considerations. The integrated nature of the architecture decision road map will be the primary focus of the Mars Architecture Team in the coming analysis cycle. This in-space transportation analysis will serve as the first guidepost of this critical activity that will guide the Moon-to-Mars Architecture into the humans to Mars segment.

Mars↗

In-Space Transportation Sensitivity to Roundtrip Mission Duration and Mars Vicinity Stay Time

The National Aeronautics and Space Administration’s Exploration Systems Development Mission Directorate has been developing architecture concepts for human missions to Mars in alignment with the agency’s Moon-to-Mars Strategy & Objectives. One of the key components of a human Mars mission is the in-space transportation system that delivers crew and cargo to Mars vicinity and returns the crew safely back to Earth. The Mars Architecture Team within ESDMD’s Strategy and Architecture Office has been evaluating multiple in-space transportation options to satisfy this functional need. The first step of this evaluation is to understand key architectural drivers for in-space transportation systems. Lunar and Mars missions have unique challenges, and systems designed for one may not be directly applicable to the other as they have different energy and mission needs. In addition, the time and distance for a Mars mission will be drastically different from the historical and current programs, requiring a new paradigm for mission architects. This paper investigates the impact of mission duration on a variety of different transportation options. As mission duration decreases, the total energy required to perform a roundtrip Mars mission increases exponentially, significantly increasing the required Earth departure mass. Sensitivity analyses of various in-space transportation system concepts were shown in this paper to begin mapping out the in-space transportation trade space. Each of the transportation systems has its own strengths and weaknesses, especially in the context of the continuous mission duration and vicinity stay time trade space. Although mission duration and Mars vicinity stay time may be key architecture drivers from an in-space transportation perspective, the decision of acceptable mission duration cannot be made solely from this perspective. Mission duration also has tangible impacts to crew health and performance, as well as system reliability and other key architecture considerations. The integrated nature of the architecture decision road map will be the primary focus of the Mars Architecture Team in the coming analysis cycle. This in-space transportation analysis will serve as the first guidepost of this critical activity that will guide the Moon-to-Mars Architecture into the humans to Mars segment.

Mars↗

Measurement of Initial Conditions at Nozzle Exit of High Speed Jets

The time averaged and unsteady density fields close to the nozzle exit (0.1 less than or = x/D less than or = 2, x: downstream distance, D: jet diameter) of unheated free jets at Mach numbers of 0.95, 1.4, and 1.8 were measured using a molecular Rayleigh scattering based technique. The initial thickness of shear layer and its linear growth rate were determined from time-averaged density survey and a modeling process, which utilized the Crocco-Busemann equation to relate density profiles to velocity profiles. The model also corrected for the smearing effect caused by a relatively long probe length in the measured density data. The calculated shear layer thickness was further verified from a limited hot-wire measurement. Density fluctuations spectra, measured using a two-Photomultiplier-tube technique, were used to determine evolution of turbulent fluctuations in various Strouhal frequency bands. For this purpose spectra were obtained from a large number of points inside the flow; and at every axial station spectral data from all radial positions were integrated. The radially-integrated fluctuation data show an exponential growth with downstream distance and an eventual saturation in all Strouhal frequency bands. The initial level of density fluctuations was calculated by extrapolation to nozzle exit.

Panda, J.↗

4U 1909+07: A Hidden Pearl

We present a detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is a persistent accreting X-ray pulsar with a period of approximately 605 s. The period changes erratically consistent with a random walk expected for a wind accreting system. INTEGRAL detects the source with an average of 2.4 cps (corresponding to 15 mCrab), but sometimes exhibits flaring activity up to 50 cps (i.e. 300 mCrab). The strongly energy dependent pulse profile shows a double peaked structure at low energies and only a single narrow peak at energies above 20 keV. The phase averaged spectrum is well described by a powerlaw modified at higher energies by an exponential cutoff and photoelectric absorption at low energies. In addition at 6.4 keV a strong iron fluorescence line and at lower energies a black body component are present. We performed phase resolved spectroscopy to study the pulse phase dependence of the spectral parameters: while most spectral parameters are constant within uncertainties, the blackbody normalization and the cutoff folding energy vary strongly with phase.

Kreykenbohm, Ingo↗

Prediction of Turbulence-Generated Noise in Unheated Jets: JeNo Technical Manual (Version 1.0) - Part 1

The model-based approach, used by the JeNo code to predict jet noise spectral directivity, is described. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. A Reynolds-averaged Navier-Stokes (RANS) solution yields the required mean flow for the solution of the propagation Green s function in a locally parallel flow. The RANS solution also produces time- and length-scales needed to model the non-compact source, the turbulent velocity correlation tensor, with exponential temporal and spatial functions. It is shown that while an exact non-causal Green s function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at low to moderate Mach numbers, the polar directivity of radiated sound is not entirely captured by this Green s function at high subsonic and supersonic acoustic Mach numbers. Results presented for unheated jets in the Mach number range of 0.51 to 1.8 suggest that near the peak radiation angle of high-speed jets, a different source/Green s function convolution integral may be required in order to capture the peak observed directivity of jet noise. A sample Mach 0.90 heated jet is also discussed that highlights the requirements for a comprehensive jet noise prediction model.

Khavaran, Abbas↗

Modeling of Turbulence Generated Noise in Jets

A numerically calculated Green's function is used to predict jet noise spectrum and its far-field directivity. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. In this paper, contributions from the so-called self- and shear-noise source terms will be discussed. A Reynolds-averaged Navier-Stokes solution yields the required mean flow as well as time- and length scales of a noise-generating turbulent eddy. A non-compact source, with exponential temporal and spatial functions, is used to describe the turbulence velocity correlation tensors. It is shown that while an exact non-causal Green's function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at moderate Mach numbers, at high subsonic and supersonic acoustic Mach numbers the polar directivity of radiated sound is not entirely captured by this Green's function. Results presented for Mach 0.5 and 0.9 isothermal jets, as well as a Mach 0.8 hot jet conclude that near the peak radiation angle a different source/Green's function convolution integral may be required in order to capture the peak observed directivity of jet noise.

Khavaran, Abbas↗