Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete phase model”

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 217 records · Page 12

SpF: Enabling Petascale Performance for Pseudospectral Dynamo Models

Pseudospectral (PS) methods possess a number of characteristics (e.g., efficiency, accuracy, natural boundary conditions) that are extremely desirable for dynamo models. Unfortunately, dynamo models based upon PS methods face a number of daunting challenges, which include exposing additional parallelism, leveraging hardware accelerators, exploiting hybrid parallelism, and improving the scalability of global memory transposes. Although these issues are a concern for most models, solutions for PS methods tend to require far more pervasive changes to underlying data and control structures. Further, improvements in performance in one model are difficult to transfer to other models, resulting in significant duplication of effort across the research community.We have developed an extensible software framework for pseudospectral methods called SpF that is intended to enable extreme scalability and optimal performance. High-level abstractions provided by SpF unburden applications of the responsibility of managing domain decomposition and load balance while reducing the changes in code required to adapt to new computing architectures. The key design concept in SpF is that each phase of the numerical calculation is partitioned into disjoint numerical kernels that can be performed entirely in-processor. The granularity of domain-decomposition provided by SpF is only constrained by the data-locality requirements of these kernels. SpF builds on top of optimized vendor libraries for common numerical operations such as transforms, matrix solvers, etc., but can also be configured to use open source alternatives for portability. SpF includes several alternative schemes for global data redistribution and is expected to serve as an ideal testbed for further research into optimal approaches for different network architectures.In this presentation, we will describe the basic architecture of SpF as well as preliminary performance data and experience with adapting legacy dynamo codes. We will conclude with a discussion of planned extensions to SpF that will provide pseudospectral applications with additional flexibility with regard to time integration, linear solvers, and discretization in the radial direction.

Pseudospectral (PS)↗

EFQPSK Versus CERN: A Comparative Study

This report presents a comparative study on Enhanced Feher's Quadrature Phase Shift Keying (EFQPSK) and Constrained Envelope Root Nyquist (CERN) techniques. These two techniques have been developed in recent times to provide high spectral and power efficiencies under nonlinear amplifier environment. The purpose of this study is to gain insights into these techniques and to help system planners and designers with an appropriate set of guidelines for using these techniques. The comparative study presented in this report relies on effective simulation models and procedures. Therefore, a significant part of this report is devoted to understanding the mathematical and simulation models of the techniques and their set-up procedures. In particular, mathematical models of EFQPSK and CERN, effects of the sampling rate in discrete time signal representation, and modeling of nonlinear amplifiers and predistorters have been considered in detail. The results of this study show that both EFQPSK and CERN signals provide spectrally efficient communications compared to filtered conventional linear modulation techniques when a nonlinear power amplifier is used. However, there are important differences. The spectral efficiency of CERN signals, with a small amount of input backoff, is significantly better than that of EFQPSK signals if the nonlinear amplifier is an ideal clipper. However, to achieve such spectral efficiencies with a practical nonlinear amplifier, CERN processing requires a predistorter which effectively translates the amplifier's characteristics close to those of an ideal clipper. Thus, the spectral performance of CERN signals strongly depends on the predistorter. EFQPSK signals, on the other hand, do not need such predistorters since their spectra are almost unaffected by the nonlinear amplifier, Ibis report discusses several receiver structures for EFQPSK signals. It is observed that optimal receiver structures can be realized for both coded and uncoded EFQPSK signals with not too much increase in computational complexity. When a nonlinear amplifier is used, the bit error rate (BER) performance of the CERN signals with a matched filter receiver is found to be more than one decibel (dB) worse compared to the bit error performance of EFQPSK signals. Although channel coding is found to provide BER performance improvement for both EFQPSK and CERN signals, the performance of EFQPSK signals remains better than that of CERN. Optimal receiver structures for CERN signals with nonlinear equalization is left as a possible future work. Based on the numerical results, it is concluded that, in nonlinear channels, CERN processing leads towards better bandwidth efficiency with a compromise in power efficiency. Hence for bandwidth efficient communications needs, CERN is a good solution provided effective adaptive predistorters can be realized. On the other hand, EFQPSK signals provide a good power efficient solution with a compromise in band width efficiency.

Borah, Deva K.↗

Comet Halley near-nucleus phenomena in 1986

Computer enhanced high resolution groundbased images taken during Halley's most active phases in 1985-86 exhibit all of the coma structures identified in 1910 photographs. The structures result mostly from discrete dust jet activity producing spirals, arcs, and jagged-edged envelopes. Multiple nucleus and antisunward jet phenomena are also observed. The curvature of the jets is consistent with the location of the Earth with respect to the nucleus equatorial plane. The variation in jet activity during the Vega and Giotto flybys is easily seen on the groundbased images. Images of CN jets show that they usually do not correlate with dust jets. Wide field, medium resolution images of faint outer envelopes are required to augment the smaller field, high resolution images of near-nucleus jets in modeling the jets.

Larson, S.↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar↗

Plasma Pressure in the Topside Ionosphere

A previous three year NASA-funded project resulted in the first 2-D maps of magnetotail pressure, density and temperature. A proposal to continue the work was declined, but modest funding was provided for one year to ramp down of the work. During the phase-out year, we used a time when 5 DMSP satellites were simultaneously active to produce the first instantaneous partial image of the magnetotail. The results have been submitted to the proceedings of the 1998 Huntsville Meeting on "The New Millennium Magnetosphere: Integrating Imaging, Discrete Observations and Global Simulations". A method of inferring central plasma sheet (CPS) temperature, density, and pressure from ionospheric observations was developed under a previous 3-year grant. These particles properties are calculated from data taken by particle instruments on DMSP satellites. Ion spectra occurring in conjunction with electron acceleration events are excluded. Because of the variability of magnetotail stretching, mapping to the plasma sheet was done using a modified Tsyganenko 1989 magnetic field model adjusted to agree with the actual magnetotail stretch. On May 25, 1997, five DMSP satellites (F10-F14) passed through the southern hemisphere nightside oval within a 19 minute period. Attached is the first magnetotail image, which results from applying our technique to that data set.

Newell, Patrick T.↗

A performance comparison of integration algorithms in simulating flexible structures

Asymptotic formulas for the characteristic root errors as well as transfer function gain and phase errors are presented for a number of traditional and new integration methods. Normalized stability regions in the lambda h plane are compared for the various methods. In particular, it is shown that a modified form of Euler integration with root matching is an especially efficient method for simulating lightly-damped structural modes. The method has been used successfully for structural bending modes in the real-time simulation of missiles. Performance of this algorithm is compared with other special algorithms, including the state-transition method. A predictor-corrector version of the modified Euler algorithm permits it to be extended to the simulation of nonlinear models of the type likely to be obtained when using the discretized structure approach. Performance of the different integration methods is also compared for integration step sizes larger than those for which the asymptotic formulas are valid. It is concluded that many traditional integration methods, such as RD-4, are not competitive in the simulation of lightly damped structures.

Howe, R. M.↗

CFD Lagrangian Modeling of Water Droplet Transport for ISS Hygiene Activity Application

The goal of this study was to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC) installed in Node 3. Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow enable identifying the paths of water transport. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 5-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain.

Son, Chang H.↗

Coherent Backscattering and Opposition Effects Observed in Some Atmosphereless Bodies of the Solar System

The results of photometric and polarimetric observations carried out for some bright atmosphere-less bodies of the Solar system near the zero phase angle reveal the simultaneous existence of two spectacular optical phenomena, the so-called brightness and polarization opposition effects. In a number of studies, these phenomena were explained by the influence of coherent backscattering. However, in general, the interference concept of coherent backscattering can be used only in the case where the particles are in the far-field zones of each other, i.e., when the scattering medium is rather rarefied. Because of this, it is important to prove rigorously and to demonstrate that the coherent backscattering effect may also exist in densely packed scattering media like regolith surface layers of celestial bodies. From the results of the computer modeling performed with the use of numerically exact solutions of the macroscopic Maxwell equations for discrete random media with different packing densities of particles, we studied the origin and evolution of all the opposition phenomena predicted by the coherent backscattering theory for low-packing-density media. It has been shown that the predictions of this theory remain valid for rather high-packing densities of particles that are typical, in particular, of regolith surfaces of the Solar system bodies. The results allow us to conclude that both opposition effects observed simultaneously in some high-albedo atmosphereless bodies of the Solar system are caused precisely by coherent backscattering of solar light in the regolith layers composed of microscopic particles.

Regolith↗

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗

The general linear inverse problem - Implication of surface waves and free oscillations for earth structure.

The discrete general linear inverse problem reduces to a set of m equations in n unknowns. There is generally no unique solution, but we can find k linear combinations of parameters for which restraints are determined. The parameter combinations are given by the eigenvectors of the coefficient matrix. The number k is determined by the ratio of the standard deviations of the observations to the allowable standard deviations in the resulting solution. Various linear combinations of the eigenvectors can be used to determine parameter resolution and information distribution among the observations. Thus we can determine where information comes from among the observations and exactly how it constraints the set of possible models. The application of such analyses to surface-wave and free-oscillation observations indicates that (1) phase, group, and amplitude observations for any particular mode provide basically the same type of information about the model; (2) observations of overtones can enhance the resolution considerably; and (3) the degree of resolution has generally been overestimated for many model determinations made from surface waves.

Wiggins, R. A.↗

Mars aqueous chemistry experiment

Mars Aqueous Chemistry Experiment (MACE) is designed to conduct a variety of measurements on regolith samples, encompassing mineral phase analyses, chemical interactions with H2O, and physical properties determinations. From these data, much can be learned or inferred regarding the past weathering environment, the contemporaneous soil micro-environments, and the general chemical and physical state of the Martian regolith. By analyzing both soil and duricrust samples, the nature of the latter may become more apparent. Sites may be characterized for comparative purposes and criteria could be set for selection of high priority materials on future sample return missions. The second year of the MACE project has shown significant progress in two major areas. MACE Instrument concept definition is a baseline design that has been generated for the complete MACE instrument, including definition of analysis modes, mass estimates and thermal model. The design includes multiple reagent reservoirs, 10 discrete analysis cells, sample manipulation capability, and thermal control. The MACE Measurement subsystems development progress is reported regarding measurement capabilities for aqueous ion sensing, evolved gas sensing, solution conductivity measurement, reagent addition (titration) capabilities, and optical sensing of suspended particles.

Clark, Benton C.↗

Effect of Boundary-Layer Bleed Hole Inclination Angle and Scaling on Flow Coefficient Behavior

Phase II data results of the Fundamental Inlet Bleed Experiments study at NASA Glenn Research Center are presented which include flow coefficient behavior for 21 bleed hole configurations. The bleed configurations are all round holes with hole diameters ranging from 0.795 to 6.35 mm, hole inclination angles from 20deg to 90deg, and thickness-to-diameter ratios from 0.25 to 2.0. All configurations were tested at a unit Reynolds number of 2.46 10(exp 7)/m and at discrete local Mach numbers of 1.33, 1.62, 1.98, 2.46, and 2.92. Interactions between the design parameters of hole diameter, hole inclination angle, and thickness-to-diameter as well as the interactions between the flow parameters of pressure ratio and Mach number upon the flow coefficient are examined, and a preliminary statistical model is proposed. An existing correlation is also examined with respect to this data.

Eichorn, Michael B.↗

CFD Model of Water Droplet Transport for ISS Hygiene Activity

The goal of the study is to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC). Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow made possible to identify the paths of water transport. The Node 3 airflow was computed for several ventilation scenarios. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 2-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain. The probability of the droplet transport to the adjacent rack surface with electronic equipment was predicted.

Son, Chang H.↗

Pulse Phase Thermography for Defect Detection and Visualization

Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Marinetti, Sergio↗

A Study of the 20 Day Superorbital Modulation in the High-Mass X-Ray Binary IGR J16493-4348

We report on Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory(Swift) X-ray Telescope (XRT), and Swift Burst Alert Telescope (BAT) observations of IGR J16493-4348, a wind-fed supergiant X-ray binary showing significant superorbital variability. From a discrete Fourier transform of the BAT light curve,we refine its superorbital period to be 20.058±0.007 days. The BAT dynamic power spectrum and a fractional root mean square analysis both show strong variations in the amplitude of the superorbital modulation, but no observed changes in the period are found. The superorbital modulation is significantly weaker between MJD 55,700 and MJD 56,300. The joint NuSTAR and XRT observations, which were performed near the minimum and maximum of one cycle of the 20 day superorbital modulation, show that the flux increases by more than a factor of two between superorbital minimum and maximum. We find no significant changes in the 3-50 keV pulse profiles between superorbital minimum and maximum, which suggests a similar accretion regime. Modeling the pulse-phase-averaged spectra we find a possible Fe Kα emission line at 6.4 keV at superorbital maximum. This feature is not significant at superorbital minimum. While we do not observe any significant differences between the pulse-phase-averaged spectral continua apart from the overall flux change, we find that the hardness ratio near the broad main peak of the pulse profile increases from superorbital minimum to maximum. This suggests the spectral shape hardens with increasing luminosity. We discuss different mechanisms that might drive the observed superorbital modulation.

Coley, Joel B.↗

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗

Design and simulation of a descent controller for strategic four-dimensional aircraft navigation

A time-controlled navigation system applicable to the descent phase of flight for airline transport aircraft was developed and simulated. The design incorporates the linear discrete-time sampled-data version of the linearized continuous-time system describing the aircraft's aerodynamics. Using optimal linear quadratic control techniques, an optimal deterministic control regulator which is implementable on an airborne computer is designed. The navigation controller assists the pilot in complying with assigned times of arrival along a four-dimensional flight path in the presence of wind disturbances. The strategic air traffic control concept is also described, followed by the design of a strategic control descent path. A strategy for determining possible times of arrival at specified waypoints along the descent path and for generating the corresponding route-time profiles that are within the performance capabilities of the aircraft is presented. Using a mathematical model of the Boeing 707-320B aircraft along with a Boeing 707 cockpit simulator interfaced with an Adage AGT-30 digital computer, a real-time simulation of the complete aircraft aerodynamics was achieved. The strategic four-dimensional navigation controller for longitudinal dynamics was tested on the nonlinear aircraft model in the presence of 15, 30, and 45 knot head-winds. The results indicate that the controller preserved the desired accuracy and precision of a time-controlled aircraft navigation system.

Lax, F. M.↗

Three-dimensional aspects of radiative transfer in remote sensing of precipitation: Application to the 1986 COHMEX storm

Several multifrequency techniques for passive microwave estimation of precipitation based on the absorption and scattering properties of hydrometers have been proposed in the literature. In the present study, plane-parallel limitations are overcome by using a model based on the discrete-ordinates method to solve the radiative transfer equation in three-dimensional rectangular domains. This effectively accounts for the complexity and variety of radiation problems encountered in the atmosphere. This investigation presents result for plane-parallel and three-dimensional radiative transfer for a precipitating system, discusses differences between these results, and suggests possible explanations for these differences. Microphysical properties were obtained from the Colorado State University Regional Atmospehric Modeling System and represent a hailstorm observed during the 1986 Cooperative Huntsville Meteorological Experiment. These properties are used as input to a three-dimensional radiative transfer model in order to simulate satellite observation of the storm. The model output consists of upwelling brightness temperatures at several of the frequencies on the Special Sensor Microwave/Imager. The radiative transfer model accounts for scattering and emission of atmospheric gases and hydrometers in liquid and ice phases. Brightness temperatures obtained from the three-dimensional model of this investigation indicate that horizontal inhomogeneities give rise to brightness temperature fields that can be quite different from fields obtained using plane-parallel radiative transfer theory. These differences are examined for various resolutions of the satellite sensor field of view. In adddition, the issue of boundary conditions for three-dimensional atmospheric radiative transfer is addressed.

Haferman, J. L.↗