Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOLE DISTRIBUTION”

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 595 records · Page 33

Long Term Variability of Cygnus X-1. VII. Orbital Variability of the Focused Wind in Cyg X-1/HDE 226868 System

Binary systems with an accreting compact object offer a unique opportunity to investigate the strong, clumpy, line-driven winds of early-type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1, using 4.77 Ms Rossi X-ray Timing Explorer (RXTE) observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott- Klein wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could either be due to lack of a focussed wind component in the model or to a more complicated clump structure.

binaries↗

Extragalactic scenarios for gamma-ray bursts

The cosmological distance scale to gamma-ray bursts is indicated by the observed distribution of the bursts in angle and intensity. The redshifts of the weakest and the strongest bursts detected by the BATSE experiment are estimated to be z = 1.5 and z = 0.2, respectively, while the strongest bursts detected by the PVO in its 11 years of operation were at z = 0.05. The bursts' rate is about 10 exp -6 per year in a galaxy like ours, and a typical energy release is about 10 exp 51 ergs, assuming isotropic emission. The rate of the neutron star and black hole collisions in close binaries is somewhat higher, and those collisions release 10 exp 53 - 10 exp 54 ergs in a burstlike form, making them the primary candidates for the observed gamma-ray events. Contrary to the common belief, no observations available to date are in conflict with the cosmological distance to the bursters. A detection of gravitational lensing of some bursts, the redshift of the weakest events, and the absorption of soft X-rays by the interstellar matter in the plane of our galaxy, are the specific predictions of the cosmological scenario.

Paczynski, Bohdan↗

Experimental studies on effervescent atomizers with wide spray angles

An experimental investigation was conducted to examine the operating and spray characteristics of two internal-mixing twin-fluid atomizers that were designed to produce effervescent atomization at low air/liquid mass ratios (ALR's). These two experimental atomizers ejected the two phase flow so as to produce a wide spray angle. One atomizer was a plain orifice design that used a four-hole exit orifice which divided and turned the two phase flow just prior to ejection. The other atomizer, called the conical sheet atomizer, ejected the two phase flow through an annular passage in such a way as to form a hollow cone spray. The atomizer operating parameters varied during this investigation were the air/liquid mass ratio, atomizer operating pressure, and, in the case of the conical sheet atomizer, the exit gap width. Studies of spray characteristics included measurements of the spray Sauter mean diameter (SMD), drop size distribution, and, for the conical sheet atomizer, circumferential distribution of the liquid mass within the spray. For both atomizers it was found that SMD decreases with an increase in either ALR or operating pressure. The effect of ALR on SMD diminishes as the value of ALR increases. For the conical sheet atomizer, when operating at low values of pressure and ALR, SMD increases with increase in gap width, but the influence of gap width on SMD diminishes with an increase in either pressure or ALR. At the highest operating pressure of the conical sheet atomizer (552 kPa), SMD is independent of gap width at all ALR's. For both atomizers, changes in operating pressure and ALR have little effect on the distribution of drop sizes in the spray.

Whitlow, J. D.↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗

Chapter 13. Atmospheric Dynamics and Meteorology

Titan, after Venus, is the second example in the solar system of an atmosphere with a global cyclostrophic circulation, but in this case a circulation that has a strong seasonal modulation in the middle atmosphere. Direct measurement of Titan's winds, particularly observations tracking the Huygens probe at 10 deg S, indicate that the zonal winds are mostly in the sense of the satellite's rotation. They generally increase with altitude and become cyclostrophic near 35 km above the surface. An exception to this is a sharp minimum centered near 75 km, where the wind velocity decreases to nearly zero. Zonal winds derived from temperatures retrieved from Cassini orbiter measurements, using the thermal wind equation, indicate a strong winter circumpolar vortex, with maximum winds of 190 m/s at mid northern latitudes near 300 km. Above this level, the vortex decays. Curiously, the stratospheric zonal winds and temperatures in both hemispheres are symmetric about a pole that is offset from the surface pole by about 4 deg. The cause of this is not well understood, but it may reflect the response of a cyclostrophic circulation to the onset between the equator, where the distance to the rotation axis is greatest, and the seasonally varying subsolar latitude. The mean meridional circulation can be inferred from the temperature field and the meridional distribution of organic molecules and condensates and hazes. Both the warm temperatures near 400 km and the enhanced concentration of several organic molecules suggest subsidence in the north polar region during winter and early spring. Stratospheric condensates are localized at high northern latitudes, with a sharp cut-off near 50 deg N. Titan's winter polar vortex appears to share many of the same characteristics of isolating high and low-latitude air masses as do the winter polar vortices on Earth that envelop the ozone holes. Global mapping of temperatures, winds, and composition in the troposphere, by contrast, is incomplete. The few suitable discrete clouds that have been found for tracking indicate smaller velocities than aloft, consistent: with the Huygens measurements, Along the descent trajectory, the Huygens measurements indicate eastward zonal winds down to 7 km, where they shift westward, and then eastward again below 1 km dawn to the surface. The low-latitude dune fields seen in Cassini RADAR images have been interpreted as longitudinal dunes occurring in a mean eastward zonal wind. This is not like Earth, where the low-latitude winds are westward above the surface. Because the net zonal-mean time-averaged torque exerted by the surface on the atmosphere should vanish, there must be westward flow over part of the surface; the question is where and when. The meridional contrast in tropospheric temperatures deduced from radio occultations at low, mid, and high latitudes. is small, approximately 5 K at the tropopause and approximately 3 K at the surface. This implies efficient heat transport, probably by axisymmetric meridional circulations. The effect of the methane "hydrological" cycle on the atmospheric circulation is not well constrained by existing measurements, Understanding the mature of the surface-atmosphere coupling will be critical to elucidating the atmospheric transports of momentum, heat, and volatiles.

Flasar, F. M.↗

Theoretical study of stress concentrations at circular holes and inclusions in strain hardening materials.

Nonlinear boundary value problems of an infinite elastic-plastic plate with a circular hole subjected to pure tension and pure shear at infinity are solved by a method involving Fourier series and finite difference. On the basis of these solutions, the validity of Neuber's relationship between the stress and strain concentration factors for the plane stress problems is examined and a generalized Stowell formula for the stress concentration factor is proposed for problems in which the applied loading may be pure shear as well as pure tension and, furthermore, other stress states. By the same method of solution, the stress distributions around a rigid circular cylindrical inclusion embedded in an infinite rigid-plastic matrix subjected to uniform transverse pure shear and tension are obtained.

Huang, W. C.↗

The NANOGrav Nine-Year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr(exp -1) for a GWB from supermassive black hole binaries of A(sub gw) less than 1.5 x 10(exp -15). We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, OMEGA(sub gw) (f) h squared less than 4.2 x 10(exp -10). Our limit on the cosmic string GWB, OMEGA(sub gw) (f) h squared less than 2.2 x 10(exp -10), translates to a conservative limit on the cosmic string tension with G mu less than 3.3 x 10(exp -8), a factor of four better than the joint Planck and high-l‚ cosmic microwave background data from other experiments.

Arzoumanian, Z.↗

Full-coverage film cooling. I - Three-dimensional measurements of turbulence structure. II - Prediction of the recovery-region hydrodynamics

Hydrodynamic measurements of turbulence structure were performed with a triaxial hot wire in the full coverage and the recovery regions following an array of injection holes under isothermal conditions at ambient temperature and pressure for blowing ratios of 0.9 and 0.4. High levels of turbulence kinetic energy (TKE) were determined for low blowing, and low TKE levels were found for the high blowing levels; in the recovery region, the flow can be represented by a model with an outer boundary layer and a 2-dimensional inner boundary layer. Recovery region hydrodynamics can be modelled by considering that a new boundary layer started to grow immediately after the end of blowing; the Prandtl mixing length distributions calculated from the values of mean velocity and turbulent shear stresses were consistent with the presence of a dual boundary layer structure in the recovery region. The program used here contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length; this 2-dimensional, finite difference program can predict the mean velocity and turbulence kinetic energy profiles based on initial values, boundary conditions, and a closure condition.

Yavuzkurt, S.↗

Compact Instrument for Measuring Profile of a Light Beam

The beamviewer is an optical device designed to be attached to a charge-coupled-device (CCD) image detector for measuring the spatial distribution of intensity of a beam of light (the beam profile ) at a designated plane intersecting the beam. The beamviewer-and-CCD combination is particularly well suited for measuring the radiant- power profile (for a steady beam) or the radiant-energy profile (for a pulsed beam) impinging on the input face or emerging from the output face of a bundle of optical fibers. The beamviewer and-CCD combination could also be used as a general laboratory instrument for profiling light beams, including beams emerging through small holes and laser beams in free space.

Papanyan, Valeri↗

Polar Coronal Hole Ephemeral Regions, the Fast Solar Wind and the Global Magnetic Dynamo

The X-Ray Telescope aboard Hinode has been regularly observing both the north and south solar polar coronal holes from November 2006 through March 2009. We use the observations of emerged flux regions within the coronal hole as evidenced by small x-ray bright points to study the physical properties of these regions. The width of the emerged flux region loop footpoints, the duration of the x-ray emission lifetime for the emerged flux region, the latitude of formation and whether an x-ray or EUV jet was observed were all recorded. In the present work we detail these observations and show a dependence on the width of the emerged flux region (bright point) to the number of x-ray jets observed. The distribution of base width is then related to a power law for number of emerged flux regions as a function of base width.

Cirtain, Jonathan W.↗

On the physics of waves in the solar atmosphere: Wave heating and wind acceleration

In the area of solar physics, new calculations of the acoustic wave energy fluxes generated in the solar convective zone was performed. The original theory developed was corrected by including a new frequency factor describing temporal variations of the turbulent energy spectrum. We have modified the original Stein code by including this new frequency factor, and tested the code extensively. Another possible source of the mechanical energy generated in the solar convective zone is the excitation of magnetic flux tube waves which can carry energy along the tubes far away from the region. The problem as to how efficiently those waves are generated in the Sun was recently solved. The propagation of nonlinear magnetic tube waves in the solar atmosphere was calculated, and mode coupling, shock formation, and heating of the local medium was studied. The wave trapping problems and evaluation of critical frequencies for wave reflection in the solar atmosphere was studied. It was shown that the role played by Alfven waves in the wind accelerations and the coronal hole heating is dominant. Presently, we are performing calculations of wave energy fluxes generated in late-type dwarf stars and studying physical processes responsible for the heating of stellar chromospheres and coronae. In the area of physics of waves, a new analytical approach for studying linear Alfven waves in smoothly nonuniform media was recently developed. This approach is presently being extended to study the propagation of linear and nonlinear magnetohydrodynamic (MHD) waves in stratified, nonisothermal and solar atmosphere. The Lighthill theory of sound generation to nonisothermal media (with a special temperature distribution) was extended. Energy cascade by nonlinear MHD waves and possible chaos driven by these waves are presently considered.

Musielak, Z. E.↗

Star-jet Interactions and Gamma-ray Outbursts from 3C454.3

We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (a red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelope lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long that varies on timescales of hours, GeV gamma-ray flare from 3C454.3, observed during 2010 November on top of a plateau lasting weeks. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the stellar atmosphere due to nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: Blackhole Mass is approx. equal to 10(exp 9) Solar Mass, the total jet power: L(j) is approx. equal to 10(exp 48) erg s(exp −1), and the Doppler factor of the gamma-ray emitting clouds: Delta is approx. equal to 20. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model.We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directlyaccelerated electrons. An injected proton distribution varies as E(exp −1) or harder below the relevant energies would be favored to alleviate the tight energetic constraints and to avoid the violation of the observational low-energy constraints.

radiation mechanisms↗

A theoretical analysis of steady-state charge collection in simple diodes under high-injection conditions

A previous rigorous mathematical analysis of drift-diffusion equations was used to investigate collected charge in a simple reverse-biased p-n junction diode exposed to an ionization source that liberates carriers (electron-hole pairs) in a quasi-neutral region within the diode. Each of two simple models was found to agree with the more rigorous analysis when carrier liberation is sufficiently intense. One is the sensitive volume (SV) model, and the other was called “ambipolar diffusion with a cutoff” (ADC). The earlier rigorous analysis was worked out in detail only for a localized source, i.e., a point source of carrier liberation, so it was able to validate the applicability of each simple model only for that case. The present paper treats an arbitrary spatial distribution of carrier generation and concludes that the ADC model remains valid for this more general case, but the SV model does

Edmonds, Larry D.↗

EUV energy distributions of accretion disks in active galactic nuclei

We have calculated UV/EUV (300 A less than or = lambda less than or = 1500 A) continuous energy distributions of accretion disks in the centers of active galactic nuclei (AGNs) for disk luminosities in the range 0.1 L(sub Edd) and central masses ranging from 10(exp 8) solar mass to 10(exp 10) solar mass. The vertical gas pressure structure of the disk is obtained analytically, the temperature stratification and the resulting continuum radiation fields are calculated numerically. We show that weak Lyman edges are an intrinsic feature of such disks. The strength of the H I Lyman edge decreases for increasing accretion rate and fixed mass of the central black hole. It increases for increasing central mass and fixed luminosity in terms of the Eddington luminosity.

Storzer, H.↗

Mass and Size Frequency Distribution of the Impact Debris from Disruption of Chondritic Meteorites

Since direct observation of the collision of asteroids in space is not always convenient for earthbound observers, we have undertaken simulations of these collisions using the NASA Ames Vertical Gun Range (AVGR). To simulate the collision of asteroids in space, and aluminum projectiles with velocities ranging from approx.1 to approx.6 km/sec were fired at 70g to approx.200 g fragments of chondritic meteorites. The target meteorite was placed in an evacuated chamber at the AVGR. Detectors, usually four, were set up around the target meteorite. These detectors consisted of aerogel and aluminum foil of varying thickness. The aerogel's purpose was to catch debris after the collision, and the aluminum foil.s purpose was to show the size of the debris particles through the size of the holes in the aluminum foil. Outside the chamber, a camera was set up to record high-speed film of the collision. This camera recorded at either 500 frames per second or 1000 frames per second. Three different types of targets were used for these tests. The first were actual meteorites, which varied in mineralogical composition, density, and porosity. The second type of target was a Hawaiian basalt, consisting of olivine phenocrysts in a porous matrix, which we thought might be similar to the chondritic meteorites, thus providing data for comparison. The final type was made out of Styrofoam. The Styrofoam was thought to simulate very low-density asteroids and comets.

VanVeghten, T. W.↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

Non-Geiger-Mode Single-Photon Avalanche Detector with Low Excess Noise

This design constitutes a self-resetting (gain quenching), room-temperature operational semiconductor single-photon-sensitive detector that is sensitive to telecommunications optical wavelengths and is scalable to large areas (millimeter diameter) with high bandwidth and efficiencies. The device can detect single photons at a 1,550-nm wavelength at a gain of 1 x 10(exp 6). Unlike conventional single photon avalanche detectors (SPADs), where gain is an extremely sensitive function to the bias voltage, the multiplication gain of this device is stable at 1 x 10(exp 6) over a wide range of bias from 30.2 to 30.9 V. Here, the multiplication gain is defined as the total number of charge carriers contained in one output pulse that is triggered by the absorption of a single photon. The statistics of magnitude of output signals also shows that the device has a very narrow pulse height distribution, which demonstrates a greatly suppressed gain fluctuation. From the histograms of both pulse height and pulse charge, the equivalent gain variance (excess noise) is between 1.001 and 1.007 at a gain of 1 x 10(exp 6). With these advantages, the device holds promise to function as a PMT-like photon counter at a 1,550- nm wavelength. The epitaxial layer structure of the device allows photons to be absorbed in the InGaAs layer, generating electron/hole (e-h) pairs. Driven by an electrical field in InGaAs, electrons are collected at the anode while holes reach the multiplication region (InAlAs p-i-n structure) and trigger the avalanche process. As a result, a large number of e-h pairs are created, and the holes move toward the cathode. Holes created by the avalanche process gain large kinetic energy through the electric field, and are considered hot. These hot holes are cooled as they travel across a p -InAlAs low field region, and are eventually blocked by energy barriers formed by the InGaAsP/ InAlAs heterojunctions. The composition of the InGaAsP alloy was chosen to have an 80 meV valance band offset with InAlAs, which is high enough to hinder the transport of the already cooled holes. Being stopped by the energy barrier, holes are accumulated at the junctions to shield the electric field, resulting in a decrease of the electric field in the multiplication region. Because the impact ionization rate is extremely sensitive to the magnitude of the electric field, the field-screening effect drastically reduces the impact ionization rate and quenches the output signals. After the avalanche pulse signal is self-quenched, the accumulated holes at the InGaAsP/ InAlAs interface escape the energy barrier through thermal excitation and tunneling and finally leave the device. The device is thus reset and ready for subsequent photon detection. This recovery time is controlled by the height of the energy barrier and the hole-cooling rate.

Zhao, Kai↗