Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse waves”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Perpendicular blade vortex interaction and its implications for helicopter noise prediction: Wave-number frequency spectra in a trailing vortex for BWI noise prediction

Perpendicular blade vortex interactions are a common occurrence in helicopter rotor flows. Under certain conditions they produce a substantial proportion of the acoustic noise. However, the mechanism of noise generation is not well understood. Specifically, turbulence associated with the trailing vortices shed from the blade tips appears insufficient to account for the noise generated. The hypothesis that the first perpendicular interaction experienced by a trailing vortex alters its turbulence structure in such a way as to increase the acoustic noise generated by subsequent interactions is examined. To investigate this hypothesis a two-part investigation was carried out. In the first part, experiments were performed to examine the behavior of a streamwise vortex as it passed over and downstream of a spanwise blade in incompressible flow. Blade vortex separations between +/- one eighth chord were studied for at a chord Reynolds number of 200,000. Three-component velocity and turbulence measurements were made in the flow from 4 chord lengths upstream to 15 chordlengths downstream of the blade using miniature 4-sensor hot wire probes. These measurements show that the interaction of the vortex with the blade and its wake causes the vortex core to loose circulation and diffuse much more rapidly than it otherwise would. Core radius increases and peak tangential velocity decreases with distance downstream of the blade. True turbulence levels within the core are much larger downstream than upstream of the blade. The net result is a much larger and more intense region of turbulent flow than that presented by the original vortex and thus, by implication, a greater potential for generating acoustic noise. In the second part, the turbulence measurements described above were used to derive the necessary inputs to a Blade Wake Interaction (BWI) noise prediction scheme. This resulted in significantly improved agreement between measurements and calculations of the BWI noise spectrum especially for the spectral peak at low frequencies, which previously was poorly predicted.

Devenport, William J.↗
NASA NTRS1985

Electrical coupling effects on the temporal evolution of F layer plasma structure

A time dependent model of F region structure decay by 'classical' cross field diffusion and electrical coupling along magnetic field lines to the E region is examined. The temporal behavior of the ion concentration fluctuations is determined by the electric field in the coupled system as well as by the initial perturbation spectra and the E region recombination rate. The formation of image structure in the E region ion concentration affects the lifetime of F layer structure in a scale size dependent way. Once an image is formed, the image amplitude and the driving F region structure amplitude decay at the same rate. At large scale sizes of at least lambda(2pi/k), this rate is proportional to k2 and the ratio of the temperatures in each region. At small scale sizes it depends on the E region recombination rate and the temperatures of the two regions but is only very weakly dependent on k. The background E region concentration determines the wave number beyond which the structure amplitude decay rate is almost independent of its scale size.

Heelis, R. A.↗
NASA NTRS2012

Investigation of Jupiter's Equatorial Hotspots and Plumes Using Cassini ISS Observations

We present updated analysis of Jupiter's equatorial meteorology from Cassini observations. For two months preceding the spacecraft's closest approach, the ISS onboard regularly imaged the atmosphere. We created time-lapse movies from this period in order to analyze the dynamics of equatorial 5-micron hot spots and their interactions with adjacent latitudes. Hot spots are quasi-stable, rectangular dark areas on visible-wavelength images, with defined eastern edges that sharply contrast with surrounding clouds, but a diffuse western edge serving as a nebulous boundary with adjacent equatorial plumes. Hot spots exhibit significant variations in size and shape over timescales of days and weeks. Some of these changes correspond with passing vortex systems from adjacent latitudes interacting with hot spots. Strong anticyclonic gyres present to the south and southeast of the dark areas appear to circulate into hot spots. Impressive, bright white plumes occupy spaces in between hot spots. Compact cirrus-iike 'scooter' clouds flow rapidly through the plumes before disappearing within the dark areas. This raises the possibility that the plumes and fast-moving clouds are at higher altitudes, because their speed does not match previously published zonal wind profiles. Most profiles represent the drift speed of the hot spots at their latitude from pattern matching of the entire longitudinal image strip. If a downward branch of an equatorially-trapped Rossby waves controls the overall appearance of hot spots, however, the westward phase velocity of the wave leads to underestimates of the true jet stream speed. Instead, our expanded data set demonstrating the rapid flow of these scooter clouds may be more illustrative of the actual jet stream speed at these latitudes. This research was supported by a NASA JDAP grant and the NASA Postdoctoral Program.

Choi, David S.↗
DOE OSTI2023

Effects of elevated-temperature deposition on the atomic structure of amorphous Ta 2 O 5 films

Brownian thermal noise as a result of mechanical loss in optical coatings will become the dominant source of noise at the most sensitive frequencies of ground-based gravitational-wave detectors. Experiments found, however, that a candidate material, amorphous Ta 2 O 5 , is unable to form an ultrastable glass and, consequently, to yield a film with significantly reduced mechanical loss through elevated-temperature deposition alone. X-ray scattering PDF measurements are carried out on films deposited and subsequently annealed at various temperatures. Inverse atomic modeling is used to analyze the short and medium range features in the atomic structure of these films. Furthermore, in silico deposition simulations of Ta 2 O 5 are carried out at various substrate temperatures and an atomic level analysis of the growth at high temperatures is presented. It is observed that upon elevated-temperature deposition, short range features remain identical, whereas medium range order increases. After annealing, however, both the short and medium range orders of films deposited at different substrate temperatures are nearly identical. A discussion on the surface diffusion and glass transition temperatures indicates that future pursuits of ultrastable low-mechanical-loss films through elevated temperature deposition should focus on materials with a high surface mobility, and/or lower glass transition temperatures in the range of achievable deposition temperatures.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Acceleration of low-energy protons and alpha particles at interplanetary shock waves

The low-energy protons and alpha particles in the energy range 30 keV/charge to 150 keV/charge associated with three different interplanetary shock waves in the immediate preshock and postshock region are studied using data obtained by the ISEE 3. The spatial distributions in the preshock and postshock medium are presented, and the dependence of the phase space density at different energies on the distance from the shock and on the form of the distribution function of both species immediately at the shock is examined. It is found that in the preshock region the particles are flowing in the solar wind frame of reference away from the shock and in the postshock medium the distribution is more or less isotropic in this frame of reference. The distribution function in the postshock region can be represented by a power law in energy which has the same spectral exponent for both protons and alpha particles. It is concluded that the first-order Fermi acceleration process can consistently explain the data, although the spectra of diffuse bow shock associated particles are different from the spectra of the interplanetary shock-associated particles in the immediate vicinity of the shock. In addition, the mean free path of the low energy ions in the preshock medium is found to be considerably smaller than the mean free path determined by the turbulence of the background interplanetary medium.

Scholer, M.↗
NASA NTRS2017

Shock Radiation Tests for Saturn and Uranus Entry Probes

This paper describes a test series in the Electric Arc Shock Tube at NASA Ames Research Center with the objective of quantifying shock-layer radiative heating magnitudes for future probe entries into Saturn and Uranus atmospheres. Normal shock waves are measured in Hydrogen-Helium mixtures (89:11 by volume) at freestream pressures between 13-66 Pa (0.1-0.5 Torr) and velocities from 20-30 kms. No shock layer radiation is detected within measurement limits below 25 kms, a finding consistent with predictions for Uranus entries. Between 25-30 kms, radiance is quantified from the Vacuum Ultraviolet through Near Infrared, with focus on the Lyman-a and Balmer series lines of Hydrogen. Shock profiles are analyzed for electron number density and electronic state distribution. The shocks do not equilibrate over several cm, and in many cases the state distributions are non-Boltzmann. Radiation data are compared to simulations of Decadal Survey entries for Saturn and shown to be as much as 8x lower than predicted with the Boltzmann radiation model. Radiance is observed in front of the shock layer, the characteristics of which match the expected diffusion length.

Hydrogen↗
DOE OSTI2022

Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. As a result, we measure little heating across the shock (<10 % of the ion kinetic energy) which is consistent with an absence of viscous dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. Here, in the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Optical Flow-Field Techniques Used for Measurements in High-Speed Centrifugal Compressors

The overall performance of a centrifugal compressor depends on the performance of the impeller and diffuser as well as on the interactions occurring between these components. Accurate measurements of the flow fields in each component are needed to develop computational models that can be used in compressor design codes. These measurements must be made simultaneously over an area that covers both components so that researchers can understand the interactions occurring between the two components. Optical measurement techniques are being used at the NASA Lewis Research Center to measure the velocity fields present in both the impeller and diffuser of a 4:1 pressure ratio centrifugal compressor operating at several conditions ranging from design flow to surge. Laser Doppler Velocimetry (LDV) was used to measure the intrablade flows present in the impeller, and the results were compared with analyses obtained from two three-dimensional viscous codes. The development of a region of low throughflow velocity fluid within this high-speed impeller was examined and compared with a similar region first observed in a large low-speed centrifugal impeller at Lewis. Particle Image Velocimetry (PIV) is a relatively new technique that has been applied to measuring the diffuser flow fields. PIV can collect data rapidly in the diffuser while avoiding the light-reflection problems that are often encountered when LDV is used. The Particle Image Velocimeter employs a sheet of pulsed laser light that is introduced into the diffuser in a quasi-radial direction through an optical probe inserted near the diffuser discharge. The light sheet is positioned such that its centerline is parallel to the hub and shroud surfaces and such that it is parallel to the diffuser vane, thereby avoiding reflections from the solid surfaces. Seed particles small enough to follow the diffuser flow are introduced into the compressor at an upstream location. A high-speed charge-coupled discharge (CCD) camera is synchronized to the laser pulse rate; this allows it to capture images of seed particle position that are separated by a small increment in time. A crosscorrelation of a particle's position in two consecutive images provides an estimate of flow velocity and direction. Multiple image pairs obtained in rapid succession at a particular flow condition provide enough measurements for statistical significance. PIV provides simultaneous velocity measurements over the entire plane that is illuminated by the light sheet instead of at a single point, as is the case when LDV is used. PIV has a further advantage in that the laser light pulse can be triggered by an external source such as a high-response pressure transducer. This feature will allow PIV to synchronize flow imaging to physical phenomena such as rotating stall or stall precursor waves. We hope that this technique can be used to obtain images of the flow field during and just prior to stall.

Skoch, Gary J.↗
NASA NTRS2000

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗
DOE OSTI2023

Fast Multigrid Reduction-in-Time for Advection via Modified Semi-Lagrangian Coarse-Grid Operators

Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs) but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the application of the multigrid reduction-in-time (MGRIT) algorithm to linear advection PDEs. Here, the key to efficient time integration with this method is using a coarse-grid operator that provides a sufficiently accurate approximation to the so-called ideal coarse-grid operator. For certain classes of semi-Lagrangian discretizations, we present a novel semi-Lagrangian-based coarse-grid operator that leads to fast and scalable multilevel time integration of linear advection PDEs. The coarse-grid operator is composed of a semi-Lagrangian discretization followed by a correction term, with the correction designed so that the leading-order truncation error of the composite operator is approximately equal to that of the ideal coarse-grid operator. Parallel results show substantial speed-ups over sequential time integration for variable-wave-speed advection problems in one and two spatial dimensions, and using high-order discretizations up to order five. The proposed approach establishes the first practical method that provides small and scalable MGRIT iteration counts for advection problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Optimal absorber thickness in long-wave multiple-stage detector

The detectivity characteristics of interband cascade infrared type-II superlattice detectors for long-wave infrared detection (λ{sub cut-off} = 8 μm at room temperature) are discussed. We present comparison of two superlattices: InAs/GaSb and InAs/InAsSb, assuming the characteristic parameters—absorption coefficients α and carrier lifetimes τ published in literature. Dependence of the Johnson-noise limited detectivity on the absorber thickness for a different number of stages is reported. Higher detectivity D{sup *} value can be achieved by increasing the carrier lifetime. However, for detectors based on type-II InAs/GaSb superlattice increasing the carrier lifetime up to 25 ns leads to a situation in which one stage is preferred, i.e. for detector with a single absorber, we obtain the highest detectivity value. In the case of InAs/InAsSb material, the situation is similar for τ ≥ 80 ns. We have shown that the optimal absorber thickness at which the highest detectivity values are obtained depends not only on the absorption coefficient α and the number of stages N{sub S}, but also on the carrier diffusion length L. According to a calculations, cascade detectors based on Ga-free material should have much higher optimal absorber thicknesses than materials based on InAs/GaSb.

42 ENGINEERING↗
DOE OSTI2021

Transient heat transfer of superfluid He 4 in nonhomogeneous geometries: Second sound, rarefaction, and thermal layer

Transient heat transfer in superfluid 4 He (He II) is a complex process that involves the interplay of the unique counterflow heat-transfer mode, the emission of second-sound waves, and the creation of quantized vortices. Many past researches focused on homogeneous heat transfer of He II in a uniform channel driven by a planar heater. Here, we report our systematic study of He II transient heat transfer in nonhomogeneous geometries that are pertinent to emergent applications. By solving the He II two-fluid equations of motion coupled with Vinen's equation for vortex-line density, we examine and compare the characteristics of transient heat transfer from planar, cylindrical, and spherical heaters in He II. Our results show that as the heater turns on, an outgoing second-sound pulse emerges, within which the vortex-line density grows rapidly. These vortices attenuate the second sound and result in a heated He II layer in front of the heater, i.e., the thermal layer. In the planar case where the vortices are created throughout the space, the second-sound pulse is continuously attenuated, leading to a thick thermal layer that diffusely spreads following the heat pulse. On the contrary, in the cylindrical and the spherical heater cases, vortices are created mainly in a thin thermal layer near the heater surface. As the heat pulse ends, a rarefaction tail develops following the second-sound pulse, in which the temperature drops. This rarefaction tail can promptly suppress the thermal layer and take away the deposited thermal energy. The effects of the heater size, heat flux, pulse duration, and temperature on the thermal-layer dynamics are discussed. We also show how the peak heat flux for the onset of boiling in He II can be studied in our model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

The Filtered Abel Transform and Its Application in Combustion Diagnostics

Many non-intrusive combustion diagnosis methods generate line-of-sight projections of a flame field. To reconstruct the spatial field of the measured properties, these projections need to be deconvoluted. When the spatial field is axisymmetric, commonly used deconvolution method include the Abel transforms, the onion peeling method and the two-dimensional Fourier transform method and its derivatives such as the filtered back projection methods. This paper proposes a new approach for performing the Abel transform method is developed, which possesses the exactness of the Abel transform and the flexibility of incorporating various filters in the reconstruction process. The Abel transform is an exact method and the simplest among these commonly used methods. It is evinced in this paper that all the exact reconstruction methods for axisymmetric distributions must be equivalent to the Abel transform because of its uniqueness and exactness. Detailed proof is presented to show that the two dimensional Fourier methods when applied to axisymmetric cases is identical to the Abel transform. Discrepancies among various reconstruction method stem from the different approximations made to perform numerical calculations. An equation relating the spectrum of a set of projection date to that of the corresponding spatial distribution is obtained, which shows that the spectrum of the projection is equal to the Abel transform of the spectrum of the corresponding spatial distribution. From the equation, if either the projection or the distribution is bandwidth limited, the other is also bandwidth limited, and both have the same bandwidth. If the two are not bandwidth limited, the Abel transform has a bias against low wave number components in most practical cases. This explains why the Abel transform and all exact deconvolution methods are sensitive to high wave number noises. The filtered Abel transform is based on the fact that the Abel transform of filtered projection data is equal to an integral transform of the original projection data with the kernel function being the Abel transform of the filtering function. The kernel function is independent of the projection data and can be obtained separately when the filtering function is selected. Users can select the best filtering function for a particular set of experimental data. When the kernal function is obtained, it can be used repeatedly to a number of projection data sets (rovs) from the same experiment. When an entire flame image that contains a large number of projection lines needs to be processed, the new approach significantly reduces computational effort in comparison with the conventional approach in which each projection data set is deconvoluted separately. Computer codes have been developed to perform the filter Abel transform for an entire flame field. Measured soot volume fraction data of a jet diffusion flame are processed as an example.

Simons, Stephen N.↗
NASA NTRS2013

Meteorology of Jupiter's Equatorial Hot Spots and Plumes from Cassini

We present an updated analysis of Jupiter's equatorial meteorology from Cassini observations. For two months preceding the spacecraft's closest approach, the Imaging Science Subsystem (ISS) onboard regularly imaged the atmosphere. We created time-lapse movies from this period in order to analyze the dynamics of equatorial hot spots and their interactions with adjacent latitudes. Hot spots are relatively cloud-free regions that emit strongly at 5 lm; improved knowledge of these features is crucial for fully understanding Galileo probe measurements taken during its descent through one. Hot spots are quasistable, rectangular dark areas on visible-wavelength images, with defined eastern edges that sharply contrast with surrounding clouds, but diffuse western edges serving as nebulous boundaries with adjacent equatorial plumes. Hot spots exhibit significant variations in size and shape over timescales of days and weeks. Some of these changes correspond with passing vortex systems from adjacent latitudes interacting with hot spots. Strong anticyclonic gyres present to the south and southeast of the dark areas appear to circulate into hot spots. Impressive, bright white plumes occupy spaces in between hot spots. Compact cirrus-like 'scooter' clouds flow rapidly through the plumes before disappearing within the dark areas. These clouds travel at 150-200 m/s, much faster than the 100 m/s hot spot and plume drift speed. This raises the possibility that the scooter clouds may be more illustrative of the actual jet stream speed at these latitudes. Most previously published zonal wind profiles represent the drift speed of the hot spots at their latitude from pattern matching of the entire longitudinal image strip. If a downward branch of an equatorially-trapped Rossby wave controls the overall appearance of hot spots, however, the westward phase velocity of the wave leads to underestimates of the true jet stream speed.

Jupitor's Equatorial Zone↗
DOE OSTI2018

Structure and Dynamics of Three-Dimensional Magnetotail Reconnection

In this study, we have used a global magnetohydrodynamic simulation and embedded particle-in-cell (PIC) simulation to analyze magnetotail reconnection in and near the electron diffusion region (EDR). Results from the magnetohydrodynamic simulation were used to set the initial and boundary conditions for the large-scale implicit PIC simulation. We examined proxies for the EDR (the nongyrotropy of the electron distribution function, slippage, and the nonideal terms in Ohm's law and the work done [J·E'] in the electron frame). The reconnection was well organized by the B x = 0 surface. All of the proxies gave false positive values but together, along with the magnetic field B z component and calculations of magnetic field lines, we were able to locate the EDR. Three of the proxies (the slippage, the nonideal terms in Ohm's law, and agyrotropy) give consistent results in the EDR. The EDR is structured and time dependent. Multiple EDRs aligned roughly parallel to the X axis extended from X ~ -31R E to X ~ -38R E with structure in the Y direction. There are regions with both J·E' < 0 and J·' > 0. Wave-like behavior with a scale of ~0.5R E or 2–3d i develops in the reconnected plasma sheet. The structure is closely related to strong electron flows in the YZ plane and to the density gradient at the outer edge of the plasma sheet. Finally, these results are consistent with expectations for the lower hybrid drift instability coupled with a shear-type instability such as the Kelvin-Helmholtz instability. Similar results were found in a Harris sheet PIC simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Comparative study on the formation of Cr and Ti ohmic contacts to (001) β -Ga 2 O 3

Here, a comparative study of Cr/Au and Ti/Au ohmic contacts on (001) β-Ga 2 O 3 was conducted. The electrical behavior from current-voltage measurements and the interfacial composition and microstructure as determined from high-resolution transmission electron microscopy (TEM) with energy dispersive x-ray analysis were compared for the different contacts at selected points in an annealing series (300–700 °C, 1 min. anneals in N 2 ). Cr/Au contacts became ohmic at temperatures (300–350 °C) approximately 50–100 °C lower than Ti/Au contacts (400–450 °C). Cr/Au and Ti/Au contacts demonstrated optimal ohmic behavior (lowest resistance) when annealed to 450–500 and 500–600 °C, respectively, with Ti/Au contacts yielding a lower total resistance than Cr/Au. Cross-sectional TEM images of Cr/Au contacts annealed at 450 °C revealed the presence of Au nanoclusters at the Ga 2 O 3 interface and CrO x layers at both the top of the contact and the Ga 2 O 3 interface. Whereas TiO x also formed at the top and bottom interfaces in 450 °C-annealed Ti/Au contacts, the TiO x surface layer appeared to be variable in thickness and/or discontinuous, unlike the CrO x surface layer. Au nanoclusters were not detected at the interface in the Ti/Au contacts. The interdiffusion and oxidation observed in both contact metallizations point to the need for diffusion barriers that may allow these contacts to be used in future Ga 2 O 3 -based devices that operate at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Surprising robustness of particle-hole symmetry for composite-fermion liquids

Here, we report on fixed phase diffusion Monte Carlo calculations that show that, even for a large amount of Landau level mixing, the energies of the Pfaffian and anti-Pfaffian phases remain very nearly the same, as also do the excitation gaps at 1/3 and 2/3. These results, combined with previous theoretical and experimental investigations, indicate that particle hole (PH) symmetry for composite fermion states is much more robust than a priori expected, emerging even in models that explicitly break PH symmetry. We provide insight into this fact by showing that the low energy physics of a generic repulsive 3-body interaction is captured, to a large extent and over a range of filling factors, by a mean field approximation that maps it into a PH symmetric 2-body interaction. This explains why Landau level mixing, which effectively generates such a generic 3-body interaction, is inefficient in breaking PH symmetry. As a byproduct, our results provide a systematic construction of a 2-body interaction which produces, to a good approximation, the Pfaffian wave function as its ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗