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 55 records · Page 3

The photochemistry of ammonia in the Jovian atmosphere

A reinvestigation of the ammonia photochemistry in the Jovian lower atmosphere above the cloud tops is conducted in connection with certain questions concerning the model proposed by Strobel (1973). Chemical reactions and their rate coefficients are considered along with the model atmosphere, solar fluxes, the continuity equations, and numerical details. Differences between the newmodel and Strobel's model are related to more moderate ammonia depletion relations. The height profiles of H, NH2, N2H3, N2H4, and NH3 concentrations are presented in a graph.

Prasad, S. S.↗

A laboratory study of the subjective response to helicopter blade-slap noise

The test stimuli recorded during a recent field study consisted of 16 sounds, each presented at 4 peak noise levels. Two helicopters and a fixed-wing aircraft were used. The impulsive characteristics of one helicopter were varied by operating at different rotor speeds, whereas the other helicopter, the noise of which was dominated by the tail rotor, displayed little variation in blade-slap noise. Thirty-two subjects made noisiness judgments on a continuous, 11 point, numerical scale. Preliminary results indicate that proposed impulsiveness corrections provide no significant improvement in the noisiness predictive ability of Effective Perceived Noise Levels (EPNL). For equal EPNL, the two categories of helicopter stimuli, one of which was far more impulsive than the other, showed no difference in judged noisiness. Examination of the physical characteristics of the sounds presented in the laboratory highlighted the difficulty of reproducing acoustical signals with high-crest factors.

Shepherd, K. P.↗

An analytical transformation technique for generating uniformly spaced computational mesh

An analytical transformation method which can map arbitrary physical coordinate grid distribution into desired computational coordinate with uniform grid distribution is derived. The transformation function and its higher derivatives are differentiable. Salient features include; (1) precise control of grid sizes; (2) more than one location of clustered grids; (3) exact positioning of particular computational nodes in the physical plane; and (4) ensuring several patches of uniformly spaced grids in the physical plane for the higher accuracies (such as at the boundaries), etc., while keeping the variation of grid spacing continuous to avoid numerical instability.

Oh, Y. H.↗

Basaltic Pressure Ridges: Formation and Relationship to Flow Emplacement

Pressure ridges and pressure plateau are common but not ubiquitous surface features on terrestrial basaltic lava flows and may reflect a specific flow emplacement mechanism. Many of the ridges are large enough to be detected on high-resolution orbital images; therefore, understanding these features could provide a means for interpreting volcanic flows on the mechanism of pressure ridge formation and how their formation relates to flow emplacement. Results from detailed field studies indicate that pressure ridges and plateaus: (1) are emplaced as individual flow lobes, (2) can be composed of primary or secondary material, (3) are dependent on duration and volume of activity within the flow unit, (4) are penecontemporaneous in formation, and (5) are indicative of slowly advancing flows with numerous flow lobes continuously forming and overriding each other.

Theilig, E.↗

Planet-Crossing Asteroid Search (PCAS)

511 independent fields were photographed with the 0.46m Schmidt during 15 observing runs along with a limited number of 1.2m Schmidt exposures. From these fields 122 new asteroids were reported. Special sets of field using the 1.2m Schmidt are continuing to yield numerous new faint asteroids. The combined result from both telescopes include 4 Mars-crossers, 3 Hungarias, 2 Phocaeas, a unique high inclination object (1985XB), 1 Apollo (1985PA) and Aten (3362) Khufu. Halley was regularly monitored (astrometry and large-scale phenomena) IHW participation. Scientifically value photographs were obtained from both 0.46m and 1.2m Schmidt Telescopes. Other accessible comets were observed including G-Z. Recently a new comet (1986d) was discovered.

Helin, E. F.↗

Thermochemical erosion

A quasi-steady multidimensional Stefan problem is treated by means of similarity, asymptotic methods, and numerical analysis. The continuous wearing away of a solid 'erodee' is effected by gradual consumption of a solid 'eroder'. In the process under study, the erodee melts and the (liquid) melt diffuses to the eroder, where an exothermic, indefinitely rapid, diffusion-controlled surface reaction occurs. The products of reaction convey some of the chemically-generated heat to the surface of the erodee, so that further melting may occur, and the process may proceed. The analysis seeks to estimate the quasi-steady rates of consumption of the erodee and eroder, for the (convenient) self-similar case of confocal parabolas defining the boundaries of the melt layer. The analysis is simplified by observing that the flow is of lubrication-theory type: creeping, primarily unidirectional motion owing to the thinness of the melt layer between the erodee and eroder, with diffusion of heat and mass predominantly perpendicular to the principal motion. The analysis is completed by consideration of the force balance on the eroder.

Butler, Gerald↗

Precision orbit determination at the NASA Goddard Space Flight Center

This paper describes the GEODYN computer program developed by the Geodynamics Branch at the NASA Goddard Space Flight Center and outlines the procedure for accurate satellite orbit and tracking-data analyses. The capabilities of the program allow the development of gravity fields as large as 90 by 90, and a complete modeling of tidal parameters. It is also feasible to numerically integrate a continuous orbit of a satellite such as Lageos for up to 12 years. The evolution of the orbit can be studied, and, by comparison with locally determined orbits, force model improvements can be made. The GEODYN flow diagrams are presented.

Putney, B.↗

Analytical study of the liquid phase transient behavior of a high temperature heat pipe

The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.

Roche, Gregory Lawrence↗

Program For Simulation Of Trajectories And Events

Universal Simulation Executive (USE) program accelerates and eases generation of application programs for numerical simulation of continuous trajectories interrupted by or containing discrete events. Developed for simulation of multiple spacecraft trajectories with events as one spacecraft crossing the equator, two spacecraft meeting or parting, or firing rocket engine. USE also simulates operation of chemical batch processing factory. Written in Ada.

Gottlieb, Robert G.↗

Minimizing Segregation During the Controlled Directional Solidification of Dendritic Alloys Publication: Metallurgical and Materials Transactions

Gravity-driven thermosolutal convection that arises during controlled directional solidification (DS) of dendritic alloys promotes detrimental macro-segregation (e.g. freckles and steepling) in products such as turbine blades. Considerable time and effort has been spent to experimentally and theoretically investigate this phenomena; although our knowledge has advanced to the point where convection can be modeled and accurately compared to experimental results, little has been done to minimize its onset and deleterious effects. The experimental work demonstrates that segregation can be. minimized and microstructural uniformity promoted when a slow axial rotation is applied to the sample crucible during controlled directional solidification processing. Numerical modeling utilizing continuation and bifurcation methods have been employed to develop accurate physical and mathematical models with the intent of identifying and optimizing processing parameters.

Grugel, R. N.↗

Fiber-Optic Continuous Liquid Sensor for Cryogenic Propellant Gauging

An innovative fiber-optic sensor has been developed for low-thrust-level settled mass gauging with measurement uncertainty <0.5 percent over cryogenic propellant tank fill levels from 2 to 98 percent. The proposed sensor uses a single optical fiber to measure liquid level and liquid distribution of cryogenic propellants. Every point of the sensing fiber is a point sensor that not only distinguishes liquid and vapor, but also measures temperature. This sensor is able to determine the physical location of each point sensor with 1-mm spatial resolution. Acting as a continuous array of numerous liquid/vapor point sensors, the truly distributed optical sensing fiber can be installed in a propellant tank in the same manner as silicon diode point sensor stripes using only a single feedthrough to connect to an optical signal interrogation unit outside the tank. Either water or liquid nitrogen levels can be measured within 1-mm spatial resolution up to a distance of 70 meters from the optical interrogation unit. This liquid-level sensing technique was also compared to the pressure gauge measurement technique in water and liquid nitrogen contained in a vertical copper pipe with a reasonable degree of accuracy. It has been demonstrated that the sensor can measure liquid levels in multiple containers containing water or liquid nitrogen with one signal interrogation unit. The liquid levels measured by the multiple fiber sensors were consistent with those virtually measured by a ruler. The sensing performance of various optical fibers has been measured, and has demonstrated that they can survive after immersion at cryogenic temperatures. The fiber strength in liquid nitrogen has also been measured. Multiple water level tests were also conducted under various actual and theoretical vibration conditions, and demonstrated that the signal-to-noise ratio under these vibration conditions, insofar as it affects measurement accuracy, is manageable and robust enough for a wide variety of spacecraft applications. A simple solution has been developed to absorb optical energy at the termination of the optical sensor, thereby avoiding any feedback to the optical interrogation unit

Xu. Wei↗

Model of Mixing Layer With Multicomponent Evaporating Drops

A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The study is motivated by the fact that typical real petroleum fuels contain hundreds of chemical species. Previously, for the sake of computational efficiency, spray studies were performed using either models based on a single representative species or models based on surrogate fuels of at most 15 species. The present multicomponent model makes it possible to perform more realistic simulations by accounting for hundreds of chemical species in a computationally efficient manner. The model is used to perform Direct Numerical Simulations in continuing studies directed toward understanding the behavior of liquid petroleum fuel sprays. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and the drops, respectively. This representation is consistent with the expected volumetrically small loading of the drops in gas (of the order of 10 3), although the mass loading can be substantial because of the high ratio (of the order of 103) between the densities of liquid and gas. The drops are treated as point sources of mass, momentum, and energy; this representation is consistent with the drop size being smaller than the Kolmogorov scale. Unsteady drag, added-mass effects, Basset history forces, and collisions between the drops are neglected, and the gas is assumed calorically perfect. The model incorporates the concept of continuous thermodynamics, according to which the chemical composition of a fuel is described probabilistically, by use of a distribution function. Distribution functions generally depend on many parameters. However, for mixtures of homologous species, the distribution can be approximated with acceptable accuracy as a sole function of the molecular weight. The mixing layer is initially laden with drops in its lower stream, and the drops are colder than the gas. Drop evaporation leads to a change in the gas-phase composition, which, like the composition of the drops, is described in a probabilistic manner

Bellan, Josette↗

Crack Growth Behavior in the Threshold Region for High Cyclic Loading

The present studies show that fanning in the threshold regime is likely caused by other factors than a plastic wake developed during load shedding. The cause of fanning at low R-values is a result of localized roughness, mainly formation of a faceted crack surface morphology , plus crack bifurcations which alters the crack closure at low R-values. The crack growth behavior in the threshold regime involves both crack closure theory and the dislocation theory of metals. Research will continue in studying numerous other metal alloys and performing more extensive analysis, such as the variation in dislocation properties (e.g., stacking fault energy) and its effects in different materials.

Forman, R.↗

An Analysis and Review of Measures and Relationships in Space Transportation Affordability

The affordability of transportation to or from space is of continued interest across numerous and diverse stakeholders in our aerospace industry. Such an important metric as affordability deserves a clear understanding among stakeholders about what is meant by affordability, costs, and related terms, as otherwise it's difficult to see where specific improvements are needed or where to target specific investments. As captured in the famous words of Lewis Carroll, "If you don't know where you are going, any road will get you there". As important as understanding a metric may be, with terms such as costs, prices, specific costs, average costs, marginal costs, etc., it is equally important to understand the relationship among these measures. In turn, these measures intermingle with caveats and factors that introduce more measures in need of a common understanding among stakeholders. These factors include flight rates, capability, and payload. This paper seeks to review the costs of space transportation systems and the relationships among the many factors involved in costs from the points of view of diverse decision makers. A decision maker may have an interest in acquiring a single launch considering the best price (along with other factors in their business case), or an interest in many launches over time. Alternately, a decision maker may have a specific interest in developing a space transportation system that will offer certain prices, or flight rate capability, or both, at a certain up-front cost. The question arises for the later, to reuse or to expend? As it is necessary in thinking about the future to clearly understand the past and the present, this paper will present data and graphics to assist stakeholders in visualizing trends and the current state of affairs in the launch industry. At all times, raw data will be referenced (or made available separately) alongside detailed explanations about the data, so as to avoid the confusion or misleading conclusions that occur more often than not with complex graphs or statements when such context is lacking.

affordability↗

QWIPs, SLS, Landsat and the International Space Station

In 1988 DARPA provided funding to NASA’s Goddard Space Flight Center to support the development of GaAs Quantum Well Infrared Photodetectors (QWIP). The goal was to make a single element photodetector that might be expandable to a two-dimensional array format. Ultimately, this led to the development of a 128 x 128 element array in collaboration with AT&T Bell Labs and Rockwell Science Center in 1990. We continued to develop numerous generations of QWIP arrays most recently resulting in the multi-QWIP focal plane for the NASA-US Geological Survey (USGS) Landsat 8 mission launched in 2013 and a similar instrument on the Landsat 9 mission to be launched in 2020. Toward the end of the Landsat 8 QWIP-based Thermal Infrared Sensor (TIRS) instrument the potential of the newly developed Strained Layer Superlattice (SLS) detector array technology became of great interest to NASA for three primary reasons: 1) higher operating temperature; 2) broad spectral response and; 3) higher sensitivity. We have collaborated extensively with QmagiQ, LLC and Northwestern University to further pursue and advance the SLS technology ever since we started back in 2012. In December of 2018 we launched the first SLS-based IR camera system to the International Space Station on board the Robotic Refueling Mission #3 (RRM3). This paper will describe the evolution of QWIP technology leading to the current development of SLS-based imaging systems at the Goddard Space Flight Center over the past 30 years.

Jhabvala, M.↗

A Computational Analysis of Boundary Layer Instability over the BOLT Configuration

The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle-of-attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane. At a selected condition from the ascent part of the anticipated flight trajectory, either type of departure from the design condition is shown to have a considerable impact on the structure of the laminar basic state and, hence, also on the amplification characteristics of the instability waves. A broader set of computations, including N-factor predictions for the relevant instabilities, will be presented in the final paper.

Boundary layer transition↗

A Computational Analysis of Boundary Layer Instability over the BOLT Configuration

The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle of attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane, specifically near both ends of the azimuthal region of a thick boundary layer in the middle, where this region rapidly changes to a thinner boundary layer on either side. At the t = 28.8767 s condition from the ascent part of the anticipated flight trajectory with a flight Mach number of M ∞ = 5.53 and unit Reynolds number of 4.25 x 10 6 /m, either type of departure from the design condition is shown to have a considerable impact on the structure of the basic state rollup within the region of interest and, hence, also on the amplification characteristics of instability waves within the resulting streaks. Yet, for a yaw angle of β = 4 degrees, the computations indicate only a slight reduction with respect to the peak N-factor of nearly 11 at the design condition of zero degrees yaw and zero degrees angle of attack. In contrast, an angle of attack equal to α = 4 degrees, the peak N-factor decreases to nearly 6 on the leeward side and increases above 16 on the windward side, making the onset of transition highly likely on the windward side. Computations also highlight the role of streak instabilities that originate as Mack mode disturbances and also demonstrate the potential pitfalls in using surface pressure sensors alone to gauge the magnitude of instability amplification.

Boundary layer transition↗

A Computational Analysis of Boundary Layer Instability over the BOLT Configuration (Effect of Nonzero AoA and Yaw)

The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle of attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane, specifically near both ends of the azimuthal region of a thick boundary layer in the middle, where this region rapidly changes to a thinner boundary layer on either side. At the t = 28.8767 s condition from the ascent part of the anticipated flight trajectory with a flight Mach number of M ∞ = 5.53 and unit Reynolds number of 4.25 x 10 6 /m, either type of departure from the design condition is shown to have a considerable impact on the structure of the basic state rollup within the region of interest and, hence, also on the amplification characteristics of instability waves within the resulting streaks. Yet, for a yaw angle of β = 4 degrees, the computations indicate only a slight reduction with respect to the peak N-factor of nearly 11 at the design condition of zero degrees yaw and zero degrees angle of attack. In contrast, an angle of attack equal to α = 4 degrees, the peak N-factor decreases to nearly 6 on the leeward side and increases above 16 on the windward side, making the onset of transition highly likely on the windward side. Computations also highlight the role of streak instabilities that originate as Mack mode disturbances and also demonstrate the potential pitfalls in using surface pressure sensors alone to gauge the magnitude of instability amplification.

Boundary layer transition↗