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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1999

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.↗

NASA NTRS1987

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗

NASA NTRS1989

Scaling and modeling of turbulent suspension flows

Scaling factors determining various aspects of particle-fluid interactions and the development of physical models to predict gas-solid turbulent suspension flow fields are discussed based on two-fluid, continua formulation. The modes of particle-fluid interactions are discussed based on the length and time scale ratio, which depends on the properties of the particles and the characteristics of the flow turbulence. For particle size smaller than or comparable with the Kolmogorov length scale and concentration low enough for neglecting direct particle-particle interaction, scaling rules can be established in various parameter ranges. The various particle-fluid interactions give rise to additional mechanisms which affect the fluid mechanics of the conveying gas phase. These extra mechanisms are incorporated into a turbulence modeling method based on the scaling rules. A multiple-scale two-phase turbulence model is developed, which gives reasonable predictions for dilute suspension flow. Much work still needs to be done to account for the poly-dispersed effects and the extension to dense suspension flows.

Chen, C. P.↗

NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗

NASA NTRS1998

Coulomb Crystallization of Charged Microspheres Levitated in a Gas Discharge Plasma

The technical topic of the project was the experimental observation of Coulomb crystallization of charged microspheres levitated in a gas discharge plasma. This suspension, sometimes termed a dusty plasma, is closely analogous to a colloidal suspension, except that it has a much faster time response, is more optically thin, and has no buoyancy forces to suspend the particles. The particles are levitated by electric fields. Through their collective Coulomb repulsions, the particles arrange themselves in a lattice with a crystalline symmetry, which undergoes an order-disorder phase transition analogous to melting when the effective temperature of the system is increased. Due to gravitational sedimentation, the particles form a thin layer in the laboratory, so that the experimental system is nearly 2D, whereas in future microgravity experiments they are expected to fill a larger volume and behave like a 3D solid or liquid. The particles are imaged using a video camera by illuminating them with a sheet of laser light. Because the suspension is optically thin, this imaging method will work as well in a 3D microgravity experiment as it does in a 2D laboratory system.

Goree, John↗

NASA NTRS1994

Visualizing Flow Over Parametric Surfaces Using Line Integral Convolution

Line Integral Convolution (LIC) is a powerful technique for imaging and animating vector fields. We extend the LIC paradigm in three ways: (1) The existing technique is limited to vector fields over a regular Cartesian grid. We extend it to vector fields over parametric surfaces, such as those found in curvilinear grids, used in computational fluid dynamics simulations; (2) Periodic motion filters can be used to animate the flow visualization. When the flow lies on a parametric surface, however, the motion appears misleading. We explain why this problem arises and show how to adjust the LIC algorithm to handle it; (3) We introduce a technique to visualize vector magnitudes as well as vector direction. Cabral and Leedom have suggested a method for variable-speed animation, which is based on varying the frequency of the filter function. We develop a different technique based on kernel phase shifts which we have found to show substantially better results. Our implementation of these algorithms utilizes texture-mapping hardware to run in real time, which allows them to be included in interactive applications.

Forssell, Lisa↗

NASA NTRS2003

Planar Particle Imaging and Doppler Velocimetry System and Method

A planar velocity measurement system (100) is operative to measure all three velocity components of a flowing fluid (106) across an illuminated plane (108) using only a single line of sight. The fluid flow is seeded with small particles which accurately follow the flow field fluctuations. The seeded flow field is illuminated with pulsed laser light source (102) and the positions of the particles in the flow are recorded on CCD cameras (122,124). The in-plane velocities are measured by determining the in-plane particle displacements. The out-of-plane velocity component is determined by measuring the Doppler shift of the light scattered by the particles. Both gas and liquid velocities can be measured, as well as two-phase flows.

Wernet, Mark P.↗

NASA NTRS1985

Middle atmosphere (60-110 km) tidal oscillations at Saskatoon, Canada (52 deg N, 107 deg W) during 1983/84

Since mid-1978, the medium frequency (2.2 MHz; M.F.) radar at Saskatoon has run continuously, measuring winds by the spaced antenna method in the upper middle atmosphere. As there is so much dynamics activity during the years of MAP it seems very important to establish the features of the tidal field at Saskatoon during 1983 to 1984, and compare these with the climatology from the four years 1978 to 1979 and 1981 to 1982. Wind profiles with 60 to 110 km (80 to 110 km at night) with 3 km resolution are obtained every 5 minutes by means of a full-correlation analysis. Hourly means are formed and harmonic analysis applied to sets of data to obtain the mean winds and tidal oscillation amplitudes and phases. The zonal wind field at Saskatoon is structurally similar to CIRA-72; however, the wind eastward cell is weaker and much more structured, and the summer eastward cell above 85 km is also weaker. For the semidiurnal tide, a strong seasonal change exists, with larger amplitudes and shorter wavelengths, in winter-like months and the reverse in summer-like months. A summary of the tidal parameters at 90 km is given, where mean values from 1978 to 1982 are also shown. This seasonal variation of tides is a midlatitude climatology, as comparisons with Monpazier, France and Christchurch, New Zealand show very similar behavior, as does Durham; and requires antisymmetric tidal modes such as 2,3 and 2,5. Considering the diurnal tide, for most months centered on summer there are long wavelength modes, or superposition of several models. However, some winter months show a short wavelength mode. A summary of the 90 km characteristics is presented.

Manson, A. H.↗

NASA NTRS1989

A combined Eulerian-Lagrangian two-phase analysis of the SSME HPOTP nozzle plug trajectories

As a result of high cycle fatigue, hydrogen embrittlement, and extended engine use, it was observed in testing that the trailing edge on the first stage nozzle plug in the High Pressure Oxygen Turbopump (HPOTP) could detach. The objective was to predict the trajectories followed by particles exiting the turbine. Experiments had shown that the heat exchanger soils, which lie downstream of the turbine, would be ruptured by particles traveling in the order of 360 ft/sec. An axisymmetric solution of the flow was obtained from the work of Lin et. al., who used INS3D to obtain the solution. The particle trajectories were obtained using the method of de Jong et. al., which employs Lagrangian tracking of the particle through the Eulerian flow field. The collision parameters were obtained from experiments conducted by Rocketdyne using problem specific alloys, speeds, and projectile geometries. A complete 3-D analysis using the most likely collision parameters shows maximum particle velocities of 200 ft/sec. in the heat exchanger region. Subsequent to this analysis, an engine level test was conducted in which seven particles passed through the turbine but no damage was observed on the heat exchanger coils.

Garcia, Robert↗

NASA NTRS2005

Field demonstration of DNAPL dehalogenation using emulsified zero-valent iron

This paper describes the results of the first field-scale demonstration conducted to evaluate the performance of nanoscale emulsified zero-valent iron (EZVI) injected into the saturated zone to enhance in situ dehalogenation of dense, nonaqueous phase liquids (DNAPLs) containing trichloroethene (TCE). EZVI is an innovative and emerging remediation technology. EZVI is a surfactant-stabilized, biodegradable emulsion that forms emulsion droplets consisting of an oil-liquid membrane surrounding zero-valent iron (ZVI) particles in water. EZVI was injected over a five day period into eight wells in a demonstration test area within a larger DNAPL source area at NASA's Launch Complex 34 (LC34) using a pressure pulse injection method. Soil and groundwater samples were collected before and after treatment and analyzed for volatile organic compounds (VOCs) to evaluate the changes in VOC mass, concentration and mass flux. Significant reductions in TCE soil concentrations (>80%) were observed at four of the six soil sampling locations within 90 days of EZVI injection. Somewhat lower reductions were observed at the other two soil sampling locations where visual observations suggest that most of the EZVI migrated up above the target treatment depth. Significant reductions in TCE groundwater concentrations (57 to 100%) were observed at all depths targeted with EZVI. Groundwater samples from the treatment area also showed significant increases in the concentrations of cis-1,2-dichloroethene (cDCE), vinyl chloride (VC) and ethene. The decrease in concentrations of TCE in soil and groundwater samples following treatment with EZVI is believed to be due to abiotic degradation associated with the ZVI as well as biodegradation enhanced by the presence of the oil and surfactant in the EZVI emulsion.

Water Pollutants/isolation & purification↗

NASA NTRS2006

Coherent Detection of High-Rate Optical PPM Signals

A method of coherent detection of high-rate pulse-position modulation (PPM) on a received laser beam has been conceived as a means of reducing the deleterious effects of noise and atmospheric turbulence in free-space optical communication using focal-plane detector array technologies. In comparison with a receiver based on direct detection of the intensity modulation of a PPM signal, a receiver based on the present method of coherent detection performs well at much higher background levels. In principle, the coherent-detection receiver can exhibit quantum-limited performance despite atmospheric turbulence. The key components of such a receiver include standard receiver optics, a laser that serves as a local oscillator, a focal-plane array of photodetectors, and a signal-processing and data-acquisition assembly needed to sample the focal-plane fields and reconstruct the pulsed signal prior to detection. The received PPM-modulated laser beam and the local-oscillator beam are focused onto the photodetector array, where they are mixed in the detection process. The two lasers are of the same or nearly the same frequency. If the two lasers are of different frequencies, then the coherent detection process is characterized as heterodyne and, using traditional heterodyne-detection terminology, the difference between the two laser frequencies is denoted the intermediate frequency (IF). If the two laser beams are of the same frequency and remain aligned in phase, then the coherent detection process is characterized as homodyne (essentially, heterodyne detection at zero IF). As a result of the inherent squaring operation of each photodetector, the output current includes an IF component that contains the signal modulation. The amplitude of the IF component is proportional to the product of the local-oscillator signal amplitude and the PPM signal amplitude. Hence, by using a sufficiently strong local-oscillator signal, one can make the PPM-modulated IF signal strong enough to overcome thermal noise in the receiver circuits: this is what makes it possible to achieve near-quantum-limited detection in the presence of strong background. Following quantum-limited coherent detection, the outputs of the individual photodetectors are automatically aligned in phase by use of one or more adaptive array compensation algorithms [e.g., the least-mean-square (LMS) algorithm]. Then the outputs are combined and the resulting signal is processed to extract the high-rate information, as though the PPM signal were received by a single photodetector. In a continuing series of experiments to test this method (see Fig. 1), the local oscillator has a wavelength of 1,064 nm, and another laser is used as a signal transmitter at a slightly different wavelength to establish an IF of about 6 MHz. There are 16 photodetectors in a 4 4 focal-plane array; the detector outputs are digitized at a sampling rate of 25 MHz, and the signals in digital form are combined by use of the LMS algorithm. Convergence of the adaptive combining algorithm in the presence of simulated atmospheric turbulence for optical PPM signals has already been demonstrated in the laboratory; the combined output is shown in Fig. 2(a), and Fig. 2(b) shows the behavior of the phase of the combining weights as a function of time (or samples). We observe that the phase of the weights has a sawtooth shape due to the continuously changing phase in the down-converted output, which is not exactly at zero frequency. Detailed performance analysis of this coherent free-space optical communication system in the presence of simulated atmospheric turbulence is currently under way.

Vilnrotter, Victor↗

NASA NTRS1972

A critical discussion of theories of flame spread across solid and liquid fuels.

Mathematical descriptions of flames spreading over liquid and solid fuels are obtained, using basic assumptions derived from observations or physical reasoning. A review of existing theories shows that they are incomplete in that they either treat an uncoupled problem of the condensed phase where the spreading rate and heat flux at the surface are given, or they merely determine the spreading rate in terms of a new vaguely defined eigenvalue. An important difference between the liquid and solid cases, due to convection, is pointed out, and it is shown that solid-fuel flame-spread theories which claim to apply to the liquid case, in reality do not apply to it. A mathematical formulation and a method of solution are presented for the phenomenon of flame spread over solid fuels with forward heat conduction in both the solid and the gas. The method uses an energy integral over the field to determine the spreading rate in terms of the basic properties of the fuel and air.

Sirignano, W. A.↗

NASA NTRS1993

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗

NASA NTRS2001

Noise Sources in a Low-Reynolds-Number Turbulent Jet at Mach 0.9

The mechanisms of sound generation in a Mach 0.9, Reynolds number 3600 turbulent jet are investigated by direct numerical simulation. Details of the numerical method are briefly outlined and results are validated against an experiment at the same flow conditions. Lighthill's theory is used to define a nominal acoustic source in the jet, and a numerical solution of Lighthill's equation is compared to the simulation to verify the computational procedures. The acoustic source is Fourier transformed in the axial coordinate and time and then filtered in order to identify and separate components capable of radiating to the far field. This procedure indicates that the peak radiating component of the source is coincident with neither the peak of the full unfiltered source nor that of the turbulent kinetic energy. The phase velocities of significant components range from approximately 5% to 50% of the ambient sound speed which calls into question the commonly made assumption that the noise sources convect at a single velocity. Space-time correlations demonstrate that the sources are not acoustically compact in the streamwise direction and that the portion of the source that radiates at angles greater than 45 deg. is stationary. Filtering non-radiating wavenumber components of the source at single frequencies reveals that a simple modulated wave forms for the source, as might be predicted by linear stability analysis. At small angles from the jet axis the noise from these modes is highly directional, better described by an exponential than a standard Doppler factor.

Freund, Jonathan B.↗

NASA NTRS1983

Microwave holographic surface measurement of the Tidbinbilla 64-m antenna

The use of a holographic technique for accurately measuring the surface figure of the Tidbinbilla (Canberra, Australia) 64-m antenna is described. Described by Bennet et al. (1976) and Scott and Ryle (1977), the technique consists of measuring the complex far-field response of the antenna at a single frequency using a terrestrial, satellite-borne or celestial radiation source of small angular diameter. This two-dimensional pattern is then Fourier transformed to yield the complex illumination function across the antenna aperture; antenna surface deviations are manifested as phase fluctuations in this function. The speed and low cost of using the technique are real advantages where observing schedules with the telescopes are heavily subscribed. For the 64-m antenna, the unweighted rms surface deviation is found to be about 1.8 mm; when weighted by the radial amplitude taper it is about 1.1 mm. Two earlier estimates by other methods had yielded tolerances of about 1.2 mm rms and 1.1 mm rms. It is planned to use the technique to survey the 22-m antennas of the Australian Telescope accurate to approximately 0.1 mm rms near 10 Ghz.

Rahmat-Samii, Y.↗

NASA NTRS1988

Multi-point measurements of ULF wave phases using a multi-channel energetic ion detector

The oscillation of differential fluxes of energetic ions modulated by a ULF wave often shows a phase shift between measurements in back-to-back detectors with look directions perpendicular to the ambient magnetic field. In a plasma of a single ion species, the phase difference is caused by displacement of the effective measurement positions by one ion gyroradius to each side of the detector. As the wave phase is periodic, the observed phase shift can correspond to a family of possible wavelengths. Simultaneous measurements of the flux modulations in different energy channels, which are equivalent to measurements of the wave phase at different positions, may make it possible to single out a unique wavelength consistent with all the measurements. Using the medium-energy-particle experiments ISEE-1 and 2, each of which may serve as a back-to-back detector, the above method was applied to a compressional Pc 5 wave observed near the equatorial plane at L between about 7 and 11. The transverse propagation properties of the wave were determined unambiguously.

Lin, N.↗

NASA NTRS1986

Coherent lidar signal fluctuation reduction by means of frequency diversity technique

The atmospheric return measured by a coherent lidar is typically characterized by rapid and deep fluctuations in signal strength. These fluctuations result from the interference of the fields backscattered to the lidar from randomly located aerosol particles which move relative to the lidar pulse. In many applications, it is necessary to determine the average value of the lidar signal intensity at some range. A new method utilizes frequency diversity initially suggested by Goldstein and subsequently studied in the microwave radar domain by others. It is expected that the application of the frequency diversity method in the coherent lidar domain will eventually provide greater efficiency and speed in the return signal averaging needed to obtain accurate intensity estimates. The frequency diversity method recognizes that the transmitted lidar pulse is very long compared to a wavelength and consequently a given phase, theta sub i, is repeated many times within the pulse. In order to test this concept, a fairly simple laboratory experiment was designed which simulates scattering of a lidar pulse from atmospheric aerosol. The testing of the frequency diversity method is discussed.

Schotland, R. M.↗

NASA NTRS1990

Investigation of mechanical and thermal properties of microwave-sintered lunar simulant materials using 2.45 GHz radiation

The mechanical and thermal properties of lunar simulant material were investigated. An alternative method of examining thermal shock in microwave-sintered lunar samples was researched. A computer code was developed that models how the fracture toughness of a thermally shocked lunar simulant sample is related to the sample hardness as measured by a micro-hardness indentor apparatus. This technique enables much data to be gathered from a few samples. Several samples were sintered at different temperatures and for different times at the temperatures. The melting and recrystallization characteristics of a well-studied binary system were also investigated to see if the thermodynamic barrier for the nucleation of a crystalline phase may be affected by the presence of a microwave field. The system chosen was the albite (sodium alumino silicate) anorthite system (calcium alumino silicate). The results of these investigations are presented.

Meek, T. T.↗