Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 397 records · Page 22

Large-scale electric fields resulting from magnetic reconnection in the corona

The method of Forbes and Priest (2-D model) is applied to the large two-ribbon flare of July 29, 1973, for which both detailed H observations and magnetic data are available. For this flare the ribbons were long, nearly straight, and parallel to each other, and the 2-D model for the coronal field geometry may be adequate. The temporal profile E(t) is calculated and indicates that reconnection sets in at the beginning of the decay phase. From this time the electric field grows rapidly to a maximum value of about 2 V/cm within just a few minutes. Thereafter it decreases monotonically with time.

Kopp, R. A.↗

Constructing the 'Best' Reliability Data for the Job - Developing Generic Reliability Data from Alternative Sources Early in a Product's Development Phase

Reliability practitioners advocate getting reliability involved early in a product development process. However, when assigned to estimate or assess the (potential) reliability of a product or system early in the design and development phase, they are faced with lack of reasonable models or methods for useful reliability estimation. Developing specific data is costly and time consuming. Instead, analysts rely on available data to assess reliability. Finding data relevant to the specific use and environment for any project is difficult, if not impossible. Instead, analysts attempt to develop the "best" or composite analog data to support the assessments. Industries, consortia and vendors across many areas have spent decades collecting, analyzing and tabulating fielded item and component reliability performance in terms of observed failures and operational use. This data resource provides a huge compendium of information for potential use, but can also be compartmented by industry, difficult to find out about, access, or manipulate. One method used incorporates processes for reviewing these existing data sources and identifying the available information based on similar equipment, then using that generic data to derive an analog composite. Dissimilarities in equipment descriptions, environment of intended use, quality and even failure modes impact the "best" data incorporated in an analog composite. Once developed, this composite analog data provides a "better" representation of the reliability of the equipment or component. It can be used to support early risk or reliability trade studies, or analytical models to establish the predicted reliability data points. It also establishes a baseline prior that may updated based on test data or observed operational constraints and failures, i.e., using Bayesian techniques. This tutorial presents a descriptive compilation of historical data sources across numerous industries and disciplines, along with examples of contents and data characteristics. It then presents methods for combining failure information from different sources and mathematical use of this data in early reliability estimation and analyses.

Kleinhammer, Roger K.↗

UAS-Based P-Band Signals of Opportunity for Remote Sensing of Snow and Root Zone Soil Moisture

We have developed the P-band Signals of Opportunity (SoOp) sensor based on the Unmanned Aircraft System (UAS) to remotely sense Snow Water Equivalent (SWE) and Root Zone Soil Moisture (RZSM). The P-band UAS SoOp sensor for Hydrology (UASHydro) would operate on the S2 aircraft developed by Black Swift Technologies for sensing of SWE and RZSM with a spatial resolution of about 10m. Root-zone soil moisture and snow water storage in land are critical parameters of the water cycle. The long-term goal of our development would be to use small UAS to perform regional high resolution observation of two key hydrological measurements to improve the estimation of terrestrial water storage for water management, crop production and forecasts of natural hazard. The UASHydro concept utilizes passive receivers to detect the reflection of strong existing P-band radio signals at the 360-380 MHz band from geostationary Mobile Use Objective System (MUOS) communication satellites launched by the US Navy. The SWE remote sensing measurement principle using the P-band SoOp is based on the propagation delay (or phase change) of radio signals through the snowpack. The time delay of the reflected signal due to the snowpack with respect to snow-free conditions is directly proportional to the snowpack SWE, while the soil moisture can be retrieved from the reflectivity at the P-band frequencies for MUOS. We have been conducting ground-based campaigns to test the instrumentation and data processing methods at the Fraser Experimental Forest in Colorado since February 2016. The field campaign data has provided support to the measurement concept. To install the SoOp technologies on the UAS, a lightweight antenna has been built and interfaces with the S2 built by Black Swift Technologies have been completed. A set of flights have been planned starting April 2018 through the end of 2018 in Colorado.

Elston, Jack↗

Turbofan forced mixer lobe flow modeling. 1: Experimental and analytical assessment

A joint analytical and experimental investigation of three-dimensional flowfield development within the lobe region of turbofan forced mixer nozzles is described. The objective was to develop a method for predicting the lobe exit flowfield. In the analytical approach, a linearized inviscid aerodynamical theory was used for representing the axial and secondary flows within the three-dimensional convoluted mixer lobes and three-dimensional boundary layer analysis was applied thereafter to account for viscous effects. The experimental phase of the program employed three planar mixer lobe models having different waveform shapes and lobe heights for which detailed measurements were made of the three-dimensional velocity field and total pressure field at the lobe exit plane. Velocity data was obtained using Laser Doppler Velocimetry (LDV) and total pressure probing and hot wire anemometry were employed to define exit plane total pressure and boundary layer development. Comparison of data and analysis was performed to assess analytical model prediction accuracy. As a result of this study a planar mixed geometry analysis was developed. A principal conclusion is that the global mixer lobe flowfield is inviscid and can be predicted from an inviscid analysis and Kutta condition.

Barber, T.↗

The kink instability in infinite cylindrical flux tubes - Eigenvalues for power-law twist profiles

Simple, accurate methods of calculating ideal MHD instability eigenvalues for infinitely long cylindrical tubes with twist function T(r) are developed. The results show that the most rapidly growing and energetic instabilities occur in the Gold-Hoyle v = 0 field, with the instability progressively weakening with increasing v. However, the maximum force eigenvalue is always small, so that even in the Gold-Hoyle case only a small proportion of the available magnetic energy can be released in the linear phase. The results also confirm that the linear pinch is remarkably weak yet relatively resistant to line-tying. It is shown that the weakness of the force eigenvalue implies that the influence of uniform gas pressure on stability is negligible. Implications for the energy-release mechanism in solar flares are discussed.

Craig, I. J. D.↗

Deep Learning Models for Planetary Seismicity Detection

Research in planetary seismology is fundamentally constrained by a lack of data. Seismo-logical science products of future missions can typically only be informed by theoretical signal/noise characteristics of the environment or likely Earth-analogues. Although objectives can be re-assessed after some initial data-collection upon lander arrival, transfer of high-resolution data back to Earth is costly on lander power usage. Over the last several years, development of GPU computing techniques and open-source high-level APIs have led to rapid advances in deep learning within the fields of computer vision, natural language processing, and collaborative filtering. These techniques are actively being adapted in seismology for a variety of tasks, including: earthquake detection, seismic phase discrimination, and ground-motion prediction. Until the recent detection of mars quakes during the Mars InSight mission, the only other measurements of seismicity recorded outside of Earth was on the Moon during the Apollo missions between 1969 to 1977. These unique data sets have been periodically revisited using new seismological methods, including ambient noise interferometry and Hidden Markov Models. Our objective is to develop a deep learning seismic detector and use it to catalog moonquakes from the Apollo 17 Lunar Seismic Profiling Experiment (LSPE) and compare the results with those obtained by other methods. Additionally, we will assess the accuracy tradeoff between using a training set of lunar data and one composed of Earth seismicity. In this document, we present preliminary results using a prototype classifier trained on a small set of earthquakes that was able to obtain detections for LSPE moonquakes with a greater accuracy than a recent study using Hidden Markov Models.

Civilini, F.↗

Computational technology for flight vehicles; Proceedings of the Symposium on Computational Technology on Flight Vehicles, Washington, DC, Nov. 5-7, 1990

The present conference on computational methods for aeronautics applications discusses topics in the fields of parallel computing, multidisciplinary computational methods, grid generation, visualization methods for CFD, probabilistic modeling, numerical simulations and methodologies for different flow regimes, computational strategies and adaptive methods in CFD, and computational strategies and dynamics and control. Attention is given to the MACH system-software kernel, a multidisciplinary approach to aeroelastic analysis, interactive grid generation with control points, interactive flow visualization using stream surfaces, numerical simulations of dynamic/aerodynamic interactions, implicit mathods for the Navier-Stokes equations, and the automatic phase-space analysis of dynamical systems.

Noor, Ahmed K.↗

Proposed dynamic phase difference method for the detection of tile debonding from the space shuttle orbiter

A noncontracting, semi-global, dynamic technique was developed for detecting loose tiles on the space shuttle orbiter. In laboratory tests on a single tile, the substrate was excited into lateral motion at a constant frequency and amplitude of 2g. The phase relationship between the motions of tile and substrate was examined by noncontacting probes in order to relate the dynamic properties of the tile SIP system to its fatigue history; by a visual technique using a stroboscope and split screen video monitor for practical application in the field. When the substrate is excited at an appropriate frequency (between 30 and 60 Hz), a good tile moves in phase and a loose tile out of phase with the substrate. The out of phase motion is readily observable in the form of a "beat" between the tile and a reference marker on the substrate.

Zuckerwar, A. J.↗

Multiple scattering from finite inhomogeneous media

Utilizing the characteristic information concerning the apparent phase constant difference between the electric and magnetic fields propagating in an inhomogeneous medium, a theoretical analysis of the multiple scattering of electromagnetic waves in finite inhomogeneous media is presented. The solution is obtained by first approximating the coefficients of a pair of exact coupled first-order differential equations and then solving the equations by first-order iteration. The present first-order approximate solution with multiple scattering considerations is shown to be more accurate than the WKB solution. Methods to improve the accuracy of the first-order solution further are discussed. Application of the solution to slowly varying finite media with periodic properties demonstrates the validity of the solution. The same approach can be extended to frequencies in the optical region by retaining additional terms in the coefficients of the coupled differential equations.

Tang, C. C. H.↗

Astrochemistry: Recent Advances in the Study of Carbon Molecules in Space

Carbon molecules and ions play an important role in space. Polycyclic Aromatic Hydrocarbons (PAHs) are the best-known candidates to account for the infrared emission bands (UIR bands) and PAH spectral features are now being used as probes of the interstellar medium in Galactic and extra-galactic environments. PAHs are also thought to be among the carriers of the diffuse interstellar absorption bands (DIBs). In the model dealing with the interstellar spectral features, PAHs are present as a mixture of radicals, ions and neutral species. PAH ionization states reflect the ionization balance of the medium while PAH size, composition, and structure reflect the energetic and chemical history of the medium. A major challenge for laboratory Astrochemistry is to reproduce (in a realistic way) the physical conditions that exist in the emission and absorption interstellar zones. An extensive laboratory program has been developed in various laboratories to characterize the physical and chemical properties of PAHs in astrophysical environments and to describe how they influence the radiation and energy balance in space and the interstellar chemistry. In particular, laboratory experiments provide measurements of the spectral characteristics of interstellar PAH analogs from the ultraviolet and visible range to the infrared range for comparison with astronomical data. The harsh physical conditions of the interstellar medium - characterized by a low temperature, an absence of collisions and strong ultraviolet radiation fields - are simulated in the laboratory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature (of the order of 100 K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase.

Salama, Farid↗

Progress in Turbulence Detection via GNSS Occultation Data

The increased availability of radio occultation (RO) data offers the ability to detect and study turbulence in the Earth's atmosphere. An analysis of how RO data can be used to determine the strength and location of turbulent regions is presented. This includes the derivation of a model for the power spectrum of the log-amplitude and phase fluctuations of the permittivity (or index of refraction) field. The bulk of the paper is then concerned with the estimation of the model parameters. Parameter estimators are introduced and some of their statistical properties are studied. These estimators are then applied to simulated log-amplitude RO signals. This includes the analysis of global statistics derived from a large number of realizations, as well as case studies that illustrate various specific aspects of the problem. Improvements to the basic estimation methods are discussed, and their beneficial properties are illustrated. The estimation techniques are then applied to real occultation data. Only two cases are presented, but they illustrate some of the salient features inherent in real data.

Cornman, L. B.↗

X-ray Mapping of Terrestrial and Extraterrestrial Materials Using the Electron Microprobe

Lunar samples returned from the Apollo program motivated development of the Bence-Albee algorithm for the rapid and accurate analysis of lunar materials, and established interlaboratory comparability through its common use. In the analysis of mineral and rock fragments it became necessary to combine micro- and macroscopic analysis by coupling electron-probe microanalysis (EPMA) with automated stage point counting. A coarse grid that included several thousand points was used, and initially wavelength-dispersive (WDS) and later energydispersive (EDS) data were acquired at discrete stage points using approx. 5 sec count times. A approx 50 micrometer beam diameter was used for WDS and up to 500 micrometer beam diameter for EDS analysis. Average analyses of discretely sampled phases were coupled with the point count data to calculate the bulk composition using matrix algebra. Use of a defocused beam resulted in a contribution from multiple phases to each analytical point, and the analytical data were deconvolved relative to end-member phase chemistry on the fly. Impressive agreement was obtained between WDS and EDS measurements as well as comparison with bulk chemistry obtained by other methods. In the 30 years since these methods were developed, significant improvements in EPMA automation and computer processing have taken place. Digital beam control allows routine collection of x-ray maps by EDS, and stage mapping for WDS is conducted continuously at slew speed and incrementally by sampling at discrete points. Digital pulse processing in EDS systems has significantly increased the throughput for EDS mapping, and the ongoing development of Si-drift detector systems promises mapping capabilities rivaling WDS systems. Spectrum imaging allows a data cube of EDS spectra to be acquired and sophisticated processing of the original data is possible using matrix algebra techniques. The study of lunar and meteoritic materials includes the need to conveniently: (1) Characterize the sample at microscopic and macroscopic scales with relatively high sensitivity, (2) Determine the modal abundance of minerals, and (3) Identify and relocate discrete features of interest in terms of size and chemistry. The coupled substitution of cations in minerals can result in significant variation in mineral chemistry, but at similar average Z, leading to poor backscattered-electron (BSE) contrast discrimination of mineralogy. It is necessary to discriminate phase chemistry at both the trace element level and the major element level. To date, the WDS of microprobe systems is preferred for mapping due to high throughput and the ability to obtain the necessary intensity to discriminate phases at both trace and major element concentrations. It is desirable to produce fully quantitative compositional maps of geological materials, which requires the acquisition of k-ratio maps that are background and dead-time corrected, and which have been corrected by phi(delta z> or an equivalent algorithm at each pixel. To date, turnkey systems do not allow the acquisition of k-ratio maps and the rigorous correction in this manner. X-ray maps of a chondrule from the Ourique meteorite, and a comb-layered xenolith from the San Francisco volcanic field, have been analyzed and processed to extract phase information. The Ourique meteorite presents a challenge due to relatively low BSE contrast, and has been studied using spectrum imaging. X-ray maps for Si, Mg, and FeK(alpha) were used to produce RGB images. The xenolith sample contains sector-zoned augite, olivine, plagioclase, and basaltic glass. X-ray maps were processed using Lispix and ImageJ software to produce mineral phase maps. The x-ray maps for Mg, Ca, and Ti were used with traceback to generate binary images that were converted to RGB images. These approaches are successful in discriminating phases, but it is desirable to achieve the methods that were used on lunar samples 30 years ago on current microprobe systems. Curnt research includes x-ray mapping analysis of the Dalgety Downs chondrite by micro x-ray fluorescence and spectrum imaging, in collaboration with Kenny Witherspoon of IXRF Systems and Dale Newbury of NIST.

Carpenter, P.↗

Waveform synthesis of surface waves in a laterally heterogeneous earth by the Gaussian beam method

The present investigation is concerned with an application of the Gaussian beam method to surface waves in the laterally heterogeneous earth. The employed method has been developed for ray tracing and synthesizing seismograms of surface waves in cases involving the laterally heterogeneous earth. The procedure is based on formulations derived by Yomogida (1985). Vertical structure of the wave field is represented by the eigenfunctions of normal mode theory, while lateral variation is expressed by the parabolic equation as in two-dimensional acoustic waves or elastic body waves. It is demonstrated that a large-amplitude change can result from a slight perturbation in the phase velocity model.

Yomogida, K.↗

The neutral atomic phases of the interstellar medium

We calculate the thermal equilibrium gas temperature of the diffuse interstellar medium. Our method incorporates a new photoelectric heating rate from small grains and polycyclic aromatic hydrocarbons (PAHs) that accounts for a size distribution of particles extending from 100 to 3 A radius. We also include a detailed treatment of the ionization rates and heating due to the soft X-ray background and due to cosmic rays. Phase diagrams (thermal pressure P versus hydrogen density n) are presented for gas that is illuminated by local interstellar far-ultraviolet (FUV) and X-ray radiation fields. A stable two-phase medium is produced with thermal pressure in the range P/k approximately = to 10(exp 3-4) K/cc. We demonstrate that photoelectric heating from PAHs dominates in the warm neutral phase (WNM) and cold neutral phase (CNM). If the C II (158 micrometers cooling per hydrogen nucleus in the solar neighborhood represents an average value for the Galaxy, we predict L(sub CII) approximately = to 7 x 10(exp 7) solar luminosities from the CNM in the Galaxy, comparable to that observed by the Cosmic Background Explorer (COBE). We discuss the dependence of the results on absorbing column density, gas phase abundances, dust abundances and metallicity, FUV field, and the X-ray radiation field. These results will be useful in modeling the multiphase structure of high-velocity clouds in the halo, the interstellar matter (ISM) at other galactocentric radii, and the ISM in external galaxies and galactic nuclei.

Wolfire, M. G.↗

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower

Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.

Computational fluid dynamics↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

Feasibility of wake vortex monitoring systems for air terminals

Wake vortex monitoring systems, especially those using laser Doppler sensors, were investigated. The initial phases of the effort involved talking with potential users (air traffic controllers, pilots, etc.) of a wake vortex monitoring system to determine system requirements from the user's viewpoint. These discussions involved the volumes of airspace to be monitored for vortices, and potential methods of using the monitored vortex data once the data are available. A subsequent task led to determining a suitable mathematical model of the vortex phenomena and developing a mathematical model of the laser Doppler sensor for monitoring the vortex flow field. The mathematical models were used in combination to help evaluate the capability of laser Doppler instrumentation in monitoring vortex flow fields both in the near vicinity of the sensor (within 1 kilometer and at long ranges(10 kilometers).

Wilson, D. J.↗