Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Decay Curve”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Ginga observations of the dipping low-mass X-ray binaries XB 1916-053 and EXO 1748-676

Ginga observations of XB 1916-053 and EXO 0748-676, low-mass X-ray binaries with pronounced dips of variable depth and duration in their X-ray light curves, are reported. Periodicity studies of XB-1916-053 show results consistent with those from previous Ginga observations, specifically, that the derived X-ray period is shorter than the optical period by about 1 percent. This problem is discussed, and a more refined model for the system is described. Three bursts were detected during the observation of XB 1916-053, the second of which occurred during a dip. Spectral analysis suggests that the second burst caused an almost instantaneous ionization of the absorbing medium in the line of sight. This material then returned to its equilibrium state on the same time scale as the burst decay. Loose limits are placed upon the possible clumpiness of the material responsible for the accretion disk structure in XB 1916-053 and EXO 0748-676.

Smale, Alan P.↗
NASA NTRS2021

TOI-2109: An Ultrahot Gas Giant on a 16 hr Orbit

We report the discovery of an ultrahot Jupiter with an extremely short orbital period of 0.67247414 ± 0.00000028 days (∼16 hr). The 1.347 ± 0.047 RJup planet, initially identified by the Transiting Exoplanet Survey Satellite (TESS) mission, orbits TOI-2109 (TIC 392476080)—a T(eff) ∼ 6500 K F-type star with a mass of 1.447 ± 0.077 Mꙩ, a radius of 1.698 ± 0.060 Rꙩ, and a rotational velocity of v sin i =81.9 ± 1.7 km/s. The planetary nature of TOI-2109b was confirmed through radial-velocity measurements, which yielded a planet mass of 5.02 ± 0.75 M(Jup). Analysis of the Doppler shadow in spectroscopic transit observations indicates a well-aligned system, with a sky-projected obliquity of λ = 1.°7± 1.°7. From the TESS full-orbit light curve, we measured a secondary eclipse depth of 731 ± 46 ppm, as well as phase-curve variations from the planet's longitudinal brightness modulation and ellipsoidal distortion of the host star. Combining the TESS-band occultation measurement with a K(s)-band secondary eclipse depth (2012 ± 80 ppm) derived from ground-based observations, we find that the dayside emission of TOI-2109b is consistent with a brightness temperature of 3631 ± 69 K, making it the second hottest exoplanet hitherto discovered. By virtue of its extreme irradiation and strong planet–star gravitational interaction, TOI-2109b is an exceptionally promising target for intensive follow-up studies using current and near-future telescope facilities to probe for orbital decay, detect tidally driven atmospheric escape, and assess the impacts of H2 dissociation and recombination on the global heat transport.

Ian Wong↗
NASA NTRS2003

GRB 011121: A Collimated Outflow into Wind-Blown Surroundings

We report optical and near-infrared follow-up observations of GRB 011121 collected predominantly at ESO telescopes in Chile. We discover a break in the afterglow light curve after 1.3 days, which implies an initial jet opening angle of about 9 deg. The jet origin of this break is supported by the fact that the spectral energy distribution is achromatic during the first four days. During later phases, GRB 011121 shows significant excess emission above the flux predicted by a power law, which we interpret as additional light from an underlying supernova. In particular, the spectral energy distribution of the optical transient approximately 2 weeks after the burst is clearly not of power-law type, but can be presented by a black body with a temperature of approx. 6000 K. The deduced parameters for the decay slope as well as the spectral index favor a wind scenario, i.e. an outflow into a circum-burst environment shaped by the stellar wind of a massive GRB progenitor. Due to its low redshift of z=0.36, GRB 011121 has been the best example for the GRB-supernova connection until GRB 030329, and provides compelling evidence for a circum-burster wind region expected to exist if the progenitor was a massive star.

Greiner, J.↗
NASA NTRS1994

The X-ray selected cataclysmic variables H0459 + 246 and H0857 - 242

We report the discovery of two cataclysmic variables located with the assistance of X-ray positions from the HEAO 1 Modulation Collimator and the Large-Area Sky Survey. Each case is distinguished by evidence of two periodic modulations that appear to represent the orbital period and the spin period of the white dwarf, respectively. The first case, H0459 + 246, has been observed optically during faint states (V approximately 16), in which there are spectral features of both an accretion disk and a K star. The light curves in the V and I bands are consistent with ellipsoidal variations in the secondary stars with a binary period of 9.952 hr. X-ray observations with EXOSAT reveal a strong pulsation with a period approximately 62 minutes. This result confirms an 'intermediate polar' classification for H0459 + 246. The pulsation is observed at 63.2 minutes in the optical I band. The long orbital period opens the possibility that H0459 + 246 is a relatively young intermediate polar that might evolve into a polar. The second case, H0857--242, shows radial velocity modulations at 1.78 hr along with photometric variations at a period approximately 49 minutes. The latter are observed only during the decay phase of dwarf nova outbursts (13 less than V less than 17), which apparently recur frequently. Given the lack of X-ray monitoring observations and the absence of proof that the 49 minute periodicity is coherent over long time-scales, we regard H0857--242 as a candidate intermediate polar. Photographic records from th e Harvard Observatory Plate Library further reveal superoutbursts for H0857--242 (V approximately 11). A bright X-ray source that is both an intermediate polar and a continually cycling dwarf nova may provide an effective means of measuring the time delay for the arrival of accreting matter at the white dwarf surface, relative to the onset of optical brightening.

Remillard, R. A.↗
NASA NTRS1989

Size distribution of dust grains: A problem of self-similarity

Distribution functions describing the results of natural processes frequently show the shape of power laws, e.g., mass functions of stars and molecular clouds, velocity spectrum of turbulence, size distributions of asteroids, micrometeorites and also interstellar dust grains. It is an open question whether this behavior is a result simply coming about by the chosen mathematical representation of the observational data or reflects a deep-seated principle of nature. The authors suppose the latter being the case. Using a dust model consisting of silicate and graphite grains Mathis et al. (1977) showed that the interstellar extinction curve can be represented by taking a grain radii distribution of power law type n(a) varies as a(exp -p) with 3.3 less than or equal to p less than or equal to 3.6 (example 1) as a basis. A different approach to understanding power laws like that in example 1 becomes possible by the theory of self-similar processes (scale invariance). The beta model of turbulence (Frisch et al., 1978) leads in an elementary way to the concept of the self-similarity dimension D, a special case of Mandelbrot's (1977) fractal dimension. In the frame of this beta model, it is supposed that on each stage of a cascade the system decays to N clumps and that only the portion beta N remains active further on. An important feature of this model is that the active eddies become less and less space-filling. In the following, the authors assume that grain-grain collisions are such a scale-invarient process and that the remaining grains are the inactive (frozen) clumps of the cascade. In this way, a size distribution n(a) da varies as a(exp -(D+1))da (example 2) results. It seems to be highly probable that the power law character of the size distribution of interstellar dust grains is the result of a self-similarity process. We can, however, not exclude that the process leading to the interstellar grain size distribution is not fragmentation at all. It could be, e.g., diffusion-limited growth discussed by Sander (1986), who applied the theory of fractal geometry to the classification of non-equilibrium growth processes. He received D=2.4 for diffusion-limited aggregation in 3d-space.

Henning, TH.↗
NASA NTRS2019

Bowles-Tatnall Wake Vortex Encounter Hazard Metric

Wake vortex spacing standards constrict the terminal area throughput and impose severe constraints on the overall capacity and efficiency of the National Airspace System. For more than two decades starting in the early 1990s, the National Aeronautics and Space Administration conducted extensive research on characterizing the formation and evolution of aircraft wakes. This multidisciplinary work included comprehensive field experiments (Pruis et al. 2016), flight tests (Vicroy et al. 1998), and wind tunnel tests (Rossow 1994; Chow et al. 1997). Parametric studies using large eddy simulations (Proctor 1998; Proctor et al. 2006) were conducted in order to develop fast-time models for the prediction of wake transport and decay (Ahmad et al. 2016). Substantial effort was spent on the formulation of acceptable vortex hazard metrics (Tatnall 1995; Hinton and Tatnall 1997). Several wake encounter severity metrics have been suggested in the past, which include the wake circulation strength, vortex-induced rolling moment coefficient (Clv), bank angle, and the roll control ratio (Tatnall 1995; Hinton and Tatnall 1997; Van der Geest 2012). The vortex-induced rolling moment coefficient introduced by Bowles and Tatnall (Tatnall 1995; Gloudemans et al. 2016) has been used extensively for risk and safety analysis of newly proposed air traffic management concepts and procedures. The original method of Bowles and Tatnall assumed a constant wing loading (the wing lift-curve slope, CL is constant), which resulted in an overestimation of the vortexinduced rolling moment coefficient. Bowles (2014) suggested a correction to the original method that provides more accurate values of Clv and which is also consistent with the underlying physics of the problem. The overestimation of Clv in the original method can be corrected by assuming an elliptical lift distribution. Figure 1.1 illustrates the correction in Clv achieved by the modified method.

Joel Malissa↗
NASA NTRS1997

The Decay of Optical Emission from the gamma-Ray Burst GRB970228

The origin of gamma-ray bursts has been one of the great unsolved mysteries in high-energy astrophysics for almost 30 years. The recent discovery of fading sources at X-ray and optical wavelengths coincident with the location of the gamma-ray burst GRB970228 therefore provides an unprecedented opportunity to probe the nature of these high-energy events. The optical counterpart appears to be a transient point source embedded in a region of extended nebulosity, the latter having been tentatively identified as a high-redshift galaxy. This would seem to favour models that place gamma-ray bursts at cosmological distances, although a range of mechanisms for producing the bursts is still allowed. A crucial piece of information for distinguishing between such models is how the brightness of the optical counterpart evolves with time. Here we re-evaluate the existing photometry of the optical counterpart of GRB970228 to construct an optical light curve for the transient event. We find that between 21 hours and six days after the burst, the R-band brightness decreased by a factor of approximately 40, with any subsequent decrease in brightness occurring at a much slower rate. As the point source faded, it also became redder. The initial behaviour of the source appears to be consistent with the 'fireball' model, but the subsequent decrease in the rate of fading may prove harder to explain.

Galama, T.↗
NASA NTRS1997

The Decay of Optical Emission Form the Gamma-Ray Burst GRB 970228

The origin of gamma-ray bursts has been one of the great unsolved mysteries in high-energy astrophysics for almost 30 years. The recent discovery of fading sources at X-ray and optical wave-lengths coincident with the location of the gamma-ray burst GRB970228 therefore provides an unprecedented opportunity to probe the nature of these high-energy events. The optical counterpart appears to be a transient point source embedded in a region of extended nebulosity, the latter having been tentatively identified as a high-redshift galaxy. This would seem to favour models that place gamma-ray bursts at cosmological distances, although a range of mechanisms for producing the bursts is still allowed. A crucial piece of information for distinguishing between such models is how the brightness of the optical counterpart evolves with time. Here we re-evaluate the existing photometry of the optical counterpart of GRB970228 to construct an optical light curve for the transient event. We find that between 21 hours and six days after the burst, the R-band brightness decreased by a factor of approx. 50, with any subsequent decrease in brightness occurring at a much slower rate. As the point source faded, it also became redder. The initial behaviour of the source appears to be consistent with the 'fireball' model, in which the burst results from the merger of two neutron stars, but the subsequent decrease in the rate of fading may prove harder to explain. The gamma-ray burst of 28 February 1997, detected with the Gamma-Ray Burst Monitor on board the BeppoSAX satellite, and located with an approx. 3 feet radius position with the Wide Field Camera on the same satellite, was the first for which a fading X-ray and optical counterpart were discovered. The optical Counterpart was discovered from a comparison of V- and I-band images taken with the William Herschel Telescope (WHT) on February 28.99 UT, and the Isaac Newton Telescope (INT; V band) and the WHT (I band) on March 8.8 uT.

Galama, T.↗
NASA NTRS1988

Novae as distance indicators

Nova shells are characteristically prolate with equatorial bands and polar caps. Failure to account for the geometry can lead to large errors in expansion parallaxes for individual novae. When simple prescriptions are used for deriving expansion parallaxes from an ensemble of randomly oriented prolate spheroids, the average distance will be too small by factors of 10 to 15 percent. The absolute magnitudes of the novae will be underestimated and the resulting distance scale will be too small by the same factors. If observations of partially resolved nova shells select for large inclinations, the systematic error in the resulting distance scale could easily be 20 to 30 percent. Extinction by dust in the bulge of M31 may broaden and shift the intrinsic distribution of maximum nova magnitudes versus decay rates. We investigated this possibility by projecting Arp's and Rosino's novae onto a composite B - 6200A color map of M31's bulge. Thirty two of the 86 novae projected onto a smooth background with no underlying structure due to the presence of a dust cloud along the line of sight. The distribution of maximum magnitudes versus fade rates for these unreddened novae is indistinguishable from the distribution for the entire set of novae. It is concluded that novae suffer very little extinction from the filamentary and patchy distribution of dust seen in the bulge of M31. Time average B and H alpha nova luminosity functions are potentially powerful new ways to use novae as standard candles. Modern CCD observations and the photographic light curves of M31 novae found during the last 60 years were analyzed to show that these functions are power laws. Consequently, unless the eruption times for novae are known, the data cannot be used to obtain distances.

Ford, Holland C.↗
NASA NTRS2003

Giant Molecular Cloud Structure and Evolution

Bodenheimer and Burkert extended earlier calculations of cloud core models to study collapse and fragmentation. The initial condition for an SPH collapse calculation is the density distribution of a Bonnor-Ebert sphere, with near balance between turbulent plus thermal energy and gravitational energy. The main parameter is the turbulent Mach number. For each Mach number several runs are made, each with a different random realization of the initial turbulent velocity field. The turbulence decays on a dynamical time scale, leading the cloud into collapse. The collapse proceeds isothermally until the density has increased to about 10(exp 13) g cm(exp -3). Then heating is included in the dense regions. The nature of the fragmentation is investigated. About 15 different runs have been performed with Mach numbers ranging from 0.3 to 3.5 (the typical value observed in molecular cloud cores is 0.7). The results show a definite trend of increasing multiplicity with increasing Mach number (M), with the number of fragments approximately proportional to (1 + M). In general, this result agrees with that of Fisher, Klein, and McKee who published three cases with an AMR grid code. However our results show that there is a large spread about this curve. For example, for M=0.3 one case resulted in no fragmentation while a second produced three fragments. Thus it is not only the value of M but also the details of the superposition of the various velocity modes that play a critical role in the formation of binaries. Also, the simulations produce a wide range of separations (10-1000 AU) for the multiple systems, in rough agreement with observations. These results are discussed in two conference proceedings.

Hollenbach, David↗
NASA NTRS1975

Photomultiplier window materials under electron irradiation - Fluorescence and phosphorescence

The fluorescence and phosphorescence of photomultiplier window materials under electron irradiation have been investigated using a Sr-90/Y-90 beta emitter as the electron source. Spectral emission curves of UV-grade, optical-grade, and electron-irradiated samples of MgF2 and LiF, and of CaF2, BaF2, sapphire, fused silica, and UV-transmitting glasses were obtained over the 200-650-nm spectral range. Fluorescence yields were determined on these materials utilizing photomultiplier tubes with cesium telluride, bialkali, and trialkali (S-20) photocathodes, respectively. Optical-grade MgF2 and LiF, as well as electron-irradiated UV-grade samples of these two materials, show enhanced fluorescence due to color-center formation and associated emission bands in the blue and red wavelength regions. Large variations in fluorescence intensities were found in UV-grade sapphire samples of different origins, particularly in the red end of the spectrum, presumably due to various amounts of chromium-ion content. Phosphorescence decay with time is best described by a sum of exponential terms, with time constants ranging from a few minutes to several days.

Viehmann, W.↗
NASA NTRS2013

Ultraviolet and Extreme-Ultraviolet Emissions at the Flare Footpoints Observed by Atmosphere Imaging Assembly

A solar flare is composed of impulsive energy release events by magnetic reconnection, which forms and heats flare loops. Recent studies have revealed a two-phase evolution pattern of UV 1600 A emission at the feet of these loops: a rapid pulse lasting for a few seconds to a few minutes, followed by a gradual decay on timescales of a few tens of minutes. Multiple band EUV observations by the Atmosphere Imaging Assembly further reveal very similar signatures. These two phases represent different but related signatures of an impulsive energy release in the corona. The rapid pulse is an immediate response of the lower atmosphere to an intense thermal conduction flux resulting from the sudden heating of the corona to high temperatures (we rule out energetic particles due to a lack of significant hard X-ray emission). The gradual phase is associated with the cooling of hot plasma that has been evaporated into the corona. The observed footpoint emission is again powered by thermal conduction (and enthalpy), but now during a period when approximate steady-state conditions are established in the loop. UV and EUV light curves of individual pixels may therefore be separated into contributions from two distinct physical mechanisms to shed light on the nature of energy transport in a flare.We demonstrate this technique using coordinated, spatially resolved observations of UV and EUV emissions from the footpoints of a C3.2 thermal flare.

Flares↗
NASA NTRS2006

Switching Characteristics of a 4H-SiC Based Bipolar Junction Transistor to 200 C

Static curves and resistive load switching characteristics of a 600 V, 4 A rated, SiC-based NPN bipolar power transistor (BJT) were observed at selected temperatures from room to 200 C. All testing was done in a pulse mode at low duty cycle (approx.0.1 percent). Turn-on was driven by an adjustable base current pulse and turn-off was accelerated by a negative base voltage pulse of 7 V. These base drive signals were implemented by 850 V, gated power pulsers, having rise-times of roughly 10 ns, or less. Base charge sweep-out with a 7 V negative pulse did not produce the large reverse base current pulse seen in a comparably rated Si-based BJT. This may be due to a very low charge storage time. The decay of the collector current was more linear than its exponential-like rise. Switching observations were done at base drive currents (I(sub B)) up to 400 mA and collector currents (I(sub C)) up to 4 A, using a 100 Omega non-inductive load. At I(sub B) = 400 mA and I(sub C) = 4 A, turn-on times typically varied from 80 to 94 ns, over temperatures from 23 to 200 C. As expected, lowering the base drive greatly extended the turn-on time. Similarly, decreasing the load current to I(sub C) = 1 A with I(sub B) = 400 mA produced turn-on times as short as 34 ns. Over the 23 to 200 C range, with I(sub B) = 400 mA and I(sub C) = 4 A, turn-off times were in the range of 72 to 84 ns with the 7 V sweep-out.

Niedra, Janis M.↗
NASA NTRS1990

The effects of streamline curvature and swirl on turbulent flows in curved ducts

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.

Cheng, Chih-Hsiung↗
NASA NTRS1994

Disk instability and the spectral evolution of the 1992 outburst of the intermediate polar GK Persei

The disk instability model can explain the previous history of dwarf-nova-like outbursts in the intermediate polar GK Per, which occur about once every three years. Disk models that reproduce the recurrence time and outburst light curves suggest that GK Per has a large effective inner disk radius (approx. 30-40 white dwarf radii) truncated by a strong magnetic field (10(exp 7) G). In this context, the effective radius is that of the portion of the disk that participates in the disk thermal instability. The radius derived is larger than the corotation radius, which must be an upper limit on the true dynamical inner radius of the disk. Disk instability models with this large effective inner radius predict that the ultraviolet continuum should be rather flat. Here we compare the predictions of the disk instability model to IUE observations of the 1981 outburst and to IUE and ROSAT observations of the recent 1992 outburst of GK Per. The model disk continuum spectral evolution is consistent with the observed UV and optical spectra, especially at maximum and in the early decay phase of the outburst. The consistency of the model with the observed UV spectra suggests that the effective inner radius of the disk is almost constant, independent of mass accretion rate, and that whatever structure lies between the effective inner radius and the corotation radius neither participates in the disk instability nor radiates substantially in the UV. The related physics of the inner disk region will be briefly discussed.

Kim, S.-W.↗
NASA NTRS1990

Observations of gamma-ray line profiles from SN 1987A

Fully resolved gamma-ray line observations from the decay of Co-56 in SN 1987A are presented. Data are from the first two balloon flights of the Gamma-Ray Imaging Spectrometer. On day 433 (after the initial optical sighting) the 847 and 1238 keV lines from Co-56 were observed at 2.3 and 4.3 sigma significance. On day 613, the lines at 847, 1238, and 2599 keV were observed at 4.6, 3.4, and 1.9 sigma, respectively. The combined significance for the three-line complex in both flights is 7.8 sigma. Gaussian profiles yield acceptable least-squares fits to the lines. The line profiles are centered on the red side of the rest energy with typical velocity dispersions of about 3500 km/s FWHM, consistent with an optically thin source, but the line intensities are less than about 30 percent of those produced by the 0.075 solar mass of Co-56 determined from the bolometric light curve.

Tueller, J.↗
NASA NTRS2002

RHESSI and Trace Observations of the 21 April 2002 X1.5 Flare

Observations of the X1.5 flare on 21 April 2002 are reviewed using the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Transition Region and Coronal Explorer (TRACE). The major findings are as follows: 1. The 3-25 keV X-rays started 54 mins before the EUV (195 A) emission suggesting that the initial energy release heated plasma directly to 220 MK, well above the 1.6 MK needed to produce the Fe XII (195 A) line. 2. Using coaligned 12-25 keV RHESSI and TRACE images, further evidence is found for the existence of hot (15-20 MK) plasma in the 195 A passband. This hot, diffuse emission is attributed to the presence of the Fe XXIV (192 A) line within the TRACE 195 A passband. 3. The 12-25 keV source centroid moves away from the limb with an apparent velocity of approx. 9.9 km/ s, slowing to approx. 1.7 km/ s after 3 hours, its find altitude being approx. 120 Mm after approx. 12 hours. This suggests that the energy release site moves to higher altitudes in agreement with classical flare models. 4. The 50-100 keV emission correlates well with EUV flare ribbons, suggesting thick-target interactions at the footpoints of the magnetic arcade. The 50-100 keV time profile matches the time derivative of the GOES light curve (Neupert effect), which suggests that the same electrons that produced the thick-target hard X-ray emission also heat the plasma seen in soft X-rays. 5. X-ray footpoint emission has an E(sup -3) spectrum down to approx. 10 keV suggesting a lower electron cutoff energy than previously thought. 6. The hard X-ray (25-200 keV) peaks have FWHM durations of approx. 1 min suggesting a more gradual energy release process than expected. 7. The TRACE images reveal a bright symmetric front propagating away from the main flare site at speeds of greater than or = 120 km/ s. This may be associated with fast CME observed several minutes later by LASCO. 8. Dark sinuous lanes are observed in the TRACE images that extend almost radially from the post-flare loop system. This "fan of spines" becomes visible well into the decay phase of the flare and shows evidence for both lateral and downward motions.

Gallagher, Peter T.↗
NASA NTRS1991

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗