Engineering PapersSearch

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 109 records · Page 6

Effects of Space Environment on Flow and Concentration During Directional Solidification

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we numerically investigate the effects of g-jitter on directional solidification. A background gravity of 1 micro-g has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.

Effect of g-jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter: the small, rapid fluctuations in gravitational acceleration which may be experienced in an orbiting space vehicle. A background gravity of 1 micro-g has been assumed, and new results for the effects of orientation angle of the periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Santiviriyapanich, P.

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.

Evidence from Galactic Cosmic Rays That the Sun Has Likely Entered A Secular Minimum in Solar Activity

Since the beginning of the space age, the Sun has been in a multi-cycle period of elevated activity (secular maximum). This secular maximum is the longest in the last 9300 years. Since the end of solar cycle 21 (SC21), however, the Sun has shown a decline in over-all activity, which has remarkably increased the fluxes of galactic cosmic rays (GCRs). Here, we investigate the correlation between the modulation of GCRs, the heliospheric magnetic field, and solar wind speed for the last 24 solar cycles to find trends that can potentially be used to predict future solar activity. Specifically, we develop a tool for predicting future magnetic field intensity, based on the hysteresis in the GCR variation, during the last phases of the current cycle. This method estimates that SC25 will be as weak as or weaker than SC24. This would mean that the Sun has likely entered a secular minimum, which, according to historical records, should last for another two cycles (SC2527and SC26).

Solar cycles

Calibration of High Frequency MEMS Microphones

Understanding and controlling aircraft noise is one of the major research topics of the NASA Fundamental Aeronautics Program. One of the measurement technologies used to acquire noise data is the microphone directional array (DA). Traditional direction array hardware, consisting of commercially available condenser microphones and preamplifiers can be too expensive and their installation in hard-walled wind tunnel test sections too complicated. An emerging micro-machining technology coupled with the latest cutting edge technologies for smaller and faster systems have opened the way for development of MEMS microphones. The MEMS microphone devices are available in the market but suffer from certain important shortcomings. Based on early experiments with array prototypes, it has been found that both the bandwidth and the sound pressure level dynamic range of the microphones should be increased significantly to improve the performance and flexibility of the overall array. Thus, in collaboration with an outside MEMS design vendor, NASA Langley modified commercially available MEMS microphone as shown in Figure 1 to meet the new requirements. Coupled with the design of the enhanced MEMS microphones was the development of a new calibration method for simultaneously obtaining the sensitivity and phase response of the devices over their entire broadband frequency range. Over the years, several methods have been used for microphone calibration. Some of the common methods of microphone calibration are Coupler (Reciprocity, Substitution, and Simultaneous), Pistonphone, Electrostatic actuator, and Free-field calibration (Reciprocity, Substitution, and Simultaneous). Traditionally, electrostatic actuators (EA) have been used to characterize air-condenser microphones for wideband frequency ranges; however, MEMS microphones are not adaptable to the EA method due to their construction and very small diaphragm size. Hence a substitution-based, free-field method was developed to calibrate these microphones at frequencies up to 80 kHz. The technique relied on the use of a random, ultrasonic broadband centrifugal sound source located in a small anechoic chamber. Phase calibrations of the MEMS microphones were derived from cross spectral phase comparisons between the reference and test substitution microphones and an adjacent and invariant grazing-incidence 1/8-inch standard microphone.

Shams, Qamar A.

Electronic heterodyne moire deflectometry: A method for transient and three dimensional density fields measurements

Effects of diffraction and nonlinear photographic emulsion characteristics on the performance of deferred electronic heterodyne moire deflectometry are investigated. The deferred deflectometry is used for measurements of nonsteady phase objects where it is difficult to complete the analysis of the field in real time. The sensitivity, accuracy and resolution of the system are calculated and it is shown that they are weakly affected by diffraction and by nonlinear recording. The feactures of the system are significantly improved compared with the conventional deferred intensity moire technique, and are comparable with the online heterodyne moire. The system was evaluated experimentally by deferred measurements of the refractive index gradients of a weak phase object consisting of a large KD*P crystal. This was done by photographing the phase object through a Ronchi grating and analyzing the tranparency with the electronic heterodyne readout system. The results are compared with the measurements performed on the same phase object with online heterodyne moire deflectometry and with heterodyne holographic interferometry methods. Some practical considerations for system improvement are discussed.

Stricker, Josef

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.

Neutral line chaos and phase space structure

Phase space structure and chaos near a neutral line are studied with numerical surface-of-section (SOS) techniques and analytic methods. Results are presented for a linear neutral line model with zero crosstail electric field. It was found that particle motion can be divided into three regimes dependening on the value of the conserved canonical momentum, Py, and the conserved Hamiltonian, h. The phase space structure, using Poincare SOS plots, is highly sensitive to bn = Bn/B0 variations, but not to h variations. It is verified that the slow motion preserves the action, Jz, as evaluated by Sonnerup (1971), when the period of the fast motion is smaller than the time scale of the slow motion. Results show that the phase space structure and particle chaos depend sensitively upon Py and bn, but are independent of h.

Burkhart, Grant R.

Development of N/P AlGaAs free-standing top solar cells for tandem applications

The combination of a free standing AlGaAs top solar cell and an existing bottom solar cell is the highest performance, lowest risk approach to implementing the tandem cell concept. The solar cell consists of an AlGaAs substrate layer, an AlGaAs base layer, an AlGaAs emitter, and an ultra-thin AlGaAs window layer. The window layer is compositionally graded which minimizes reflection at the window layer/emitter interface and creates a built-in electric field to improve quantum response in the blue region of the spectrum. Liquid phase epitaxy (LPE) is the only viable method to produce this free standing top solar cell. Small (0.125 sq cm), transparent p/n AlGaAs top solar cells were demonstrated with optimum bandgap for combination with a silicon bottom solar cell. The efficiency of an AlGaAs/Si stack using the free standing AlGaAs device upon an existing silicon bottom solar cell is 24 pct. (1X, Air Mass Zero (AM0). The n/p AlGaAs top solar cell is being developed in order to facilitate the wiring configuration. The two terminal tandem stack will retain fit, form, and function of existing silicon solar cells. Progress in the development of large area (8 and 16 sq cm), free standing AlGaAs top solar cells is discussed.

Negley, Gerald H.

Heat transfer in space systems; Proceedings of the Symposium, AIAA/ASME Thermophysics and Heat Transfer Conference, Seattle, WA, June 18-20, 1990

Theoretical and experimental studies of heat-tranfer in a space environment are discussed in reviews and reports. Topics addressed include a small-scale two-phase thermosiphon to cool high-power electronics, a low-pressure-drop heat exchanger with integral heat pipe, an analysis of the thermal performance of heat-pipe radiators, measurements of temperature and concentration fields in a rectangular heat pipe, and a simplified aerothermal heating method for axisymmetric blunt bodies. Consideration is given to entropy production in a shock wave, bubble-slug transition in a two-phase liquid-gas flow under microgravity, plasma arc welding under normal and zero gravity, the Microgravity Thaw Experiment, the flow of a thin film on stationary and rotating disks, an advanced ceramic fabric body-mounted radiator for Space Station Freedom phase 0 design, and lunar radiators with specular reflectors.

Chan, S. H.

Directionally Solidified Ceramics Produced

Produced Multiphase, interpenetrating structures are an alternative route to obtaining structural ceramic materials with adequate strength, toughness, and stability for high-temperature aerospace applications. The eutectic architecture, a continuous-reinforcing phase within a higher volume phase or matrix, can be described as a naturally occurring, in situ composite. The phases of a eutectic are thermodynamically compatible at high homologous temperatures. Strong and stable materials have been produced. Toughness, however, remains a technical obstacle. The potential for producing materials with enhanced toughness along with adequate strength and stability was demonstrated using the laser-heated float zone (LHFZ) growth method at the NASA Glenn Research Center at Lewis Field. LHFZ growth at Glenn provides a means to efficiently produce and record the underlying growth phenomena associated with two-phase structures. To initiate directional solidification, a seed of single-crystal sapphire (<0001> direction) was lowered onto the molten liquid until wetting occurred and then withdrawn at a constant rate. Neither the crystal nor the source rod was rotated. The materials produced were tested mechanically in tension, and the resulting microstructure was examined with a scanning electron microscope. Both the inherent properties of the constituent phases and the properties of the interface between them affect the mechanical behavior and the fracture surfaces. The following scanning electron micrographs show the microstructures of two different materials that were tested to failure in tension. In the left micrograph, the flat fracture surface is typical of a material that is strong but has low toughness. In the right micrograph, the crack is effectively deflected at the interface between the two phases, achieving higher toughness at moderately lower strength levels. Conducting mechanical tests to determine the high temperature properties of these materials is the next step in determining their eventual suitability.

Farmer, Serene C.

Linear-phase approximation in the triangular facet near-field physical optics computer program

Analyses of reflector antenna surfaces use a computer program based on a discrete approximation of the radiation integral. The calculation replaces the actual surface with a triangular facet representation; the physical optics current is assumed to be constant over each facet. Described here is a method of calculation using linear-phase approximation of the surface currents of parabolas, ellipses, and shaped subreflectors and compares results with a previous program that used a constant-phase approximation of the triangular facets. The results show that the linear-phase approximation is a significant improvement over the constant-phase approximation, and enables computation of 100 to 1,000 lambda reflectors within a reasonable time on a Cray computer.

Imbriale, W. A.

Electronic Compensation For Distortion Of Antenna Reflector

Proposed method of obtaining approximately desired radiation or reception pattern from antenna that includes reflector based on concept of superposition of electromagnetic fields generated by multiple feedhorns or feed antenna elements arrayed at various positions near reflector and excited at electronically adjustable magnitudes and phases. In intended application, reflector nominally paraboloidal, feed elements N feedhorns in hexagonal array, and method used to compensate for deviations of real reflector surface from nominal paraboloidal shape. Method and concept also applicable to electronic beam steering and electronic antenna compensation in other situations.

Bailey, M. C.

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.

Electron-beam flow visualization - Applications in the definition of configuration aerothermal characteristics.

Comparisons between flow visualization systems using electron-beam fluorescence, schlieren, and shadowgraph techniques illustrate the advantages associated with the electron beam. Specific applications of this method as an aid in defining the origin of erosion on a heat-transfer model are cited. Results of combined electron-beam oil-flow studies on configurations illustrate that the simultaneous definition of the external flow field and its surface flow can be obtained. Comparisons between the electron-beam oil-flow visualization method and phase-change coating heat-transfer tests on a shuttle ascent configuration indicate the complementary nature of these two testing techniques. Potential methods for improving the electron-beam technique are included.

Woods, W. C.

The surface brightness of reflection nebulae

Hubble's equation relating the maximum apparent angular extent of a reflection nebula to the apparent magnitude of the illuminating star has been reconsidered under a set of less restrictive assumptions. A computational technique is developed which permits the use of fits to observed m, log a values to determine the albedo of the particles composing reflection nebulae, providing only that one assumes a particular phase function. Despite the fact that all orders of scattering, anisotropic phase functions, and illumination by the general stellar field are considered, the albedo which is determined for reflection nebulae by this method appears larger than that for interstellar particles in general. The possibility that the higher surface brightness might be due to a continuous fluorescence mechanism is considered both theoretically and observationally.

Rush, W. F.

Acoustic levitation and manipulation for space applications

A wide spectrum of experiments to be performed in space in a microgravity environment require levitation and manipulation of liquid or molten samples. A novel acoustic method has been developed at JPL for controlling liquid samples without physical contacts. This method utilizes the static pressure generated by three orthogonal acoustic standing waves excited within an enclosure. Furthermore, this method will allow the sample to be rotated and/or oscillated by modifying the phase angles and/or the amplitude of the acoustic field. This technique has been proven both in our laboratory and in a microgravity environment provided by KC-135 flights. Samples placed within our chamber driven at (1,0,0), (0,1,0), and (0,0,1), modes were indeed levitated, rotated, and oscillated.

Wang, T. G.

Scan-independent slot arrays with parasitic wire arrays in stratified medium

The effect of a parasitic wire array on the scan admittance of a slot array has been investigated. Structures considered can consist of an infinite slot array and an arbitrary number of parasitic infinte arrays of piecewise linear wires, all arrays being embedded in a stratified medium. These include, as particular cases, phased arrays of Clavin elements. Expressing the fields from the arrays as plane waves, a procedure similar to the periodic moment method for infinite periodic structures is set up to obtain the the scan admittance of the slot array. Scan admittances are presented for a slot array with monopole arrays in free space, and a slot array with a tilted dipole array in a stratified medium. Blind spots at which the incident energy is mainly reflected rather than transmitted were found. Results obtained indicate the possibility of using parasitic wire arrays for scan compensation of active slot arrays.

Ng, Kwong T.