Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time dependent fluctuations”

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 325 records · Page 18

The structure, energetics and evolution of the dominant frequency-dependent three-dimensional atmospheric modes

The three-dimensional structure of the dominant modes of wintertime variability are examined, and the mode structures are determined by an EOF analysis of the coupled vertical mean and shear streamfunction fields for a wide range of time-scale classes. The time mean energy conversions are determined from a two-layer quasigeostrophic model. Lag-cross correlations between the dominant low-frequency modes tend to be small and/or symmetric about zero, the high frequency modes are found to show highly asymmetric sinusoidal cross correlation functions similar to traveling waves, and some intermediate and short-time scale modes exhibit correlations indicating an association with the decaying phases of blocking in the Pacific. It is suggested that the mean flow barotropic (baroclinic) instabilities are an important source of variance for some of the lowest (highest) fluctuations considered.

Schubert, S. D.↗

Propagation velocity of perturbations in turbulent channel flow

A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the streamwise component of the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) less than 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity Uc for velocity and vorticity perturbations and 75 percent of U sub c for pressure perturbations. This is equal to U at y(+) = 15 for velocity and vorticity perturbations, and equal to U at y(+) = 20 for pressure perturbations, indicating that the dynamics of the nearwall turbulence is controlled by turbulence structures present near y(+) about 15-20. Scale dependence of V is also examined by analyzing the bandpass-filtered flow fields. This paper contains comprehensive documentation on the propagation velocities, which should prove useful in the evaluation of Taylor's hypothesis.

Kim, John↗

Sporadic E Morphology from GPS-CHAMP Radio Occultation

The scintillations of phase and amplitude in terms of signal-to-noise ratio (SNR) of the GPS radio occultation signal are caused by thin ionization layers. These thin irregular electron density layers in the E region ionosphere are often called sporadic E (Es). For a monthly retrieval of Es morphology we use the variances of the phase and SNR fluctuations of worldwide ~6000 GPS/CHAMP occultations in the E region. The Es climatology is studied globally with the SNR and phase variances in terms of monthly zonal means, seasonal maps, and diurnal and long-term variations. The zonal mean variances reveal strong, extended Es activities at summertime midlatitudes but weak, confined activities in wintertime high latitudes, peaking at ~105 km. Global maps at 105-km altitude show clear dependence of Es activities on the geomagnetic dip angle, where the summertime midlatitude Es occurs mostly at dip angles of 30 deg. - 60 deg. and the wintertime high-latitude enhancement occurs mostly at dip angles greater than 80 deg. The midlatitude Es variances exhibit a strong semidiurnal variation with peak hours near 0800 1000 and 2000 local solar time, respectively. The peak hours are delayed slightly with decreasing height, suggesting influences from the semidiurnal tide. To provide more insights on the observed SNR and phase variances, we model radio wave propagation for the CHAMP observing geometry under several perturbed cases in the E region ionosphere. The model simulations indicate that the SNR variance has the maximum response to Es perturbations at vertical wavelengths of 1.2 km, whereas the phase response maximizes at ~2 km (for the 1-s variance analysis). The characteristic scale depends little on the truncation time used in the SNR variance analysis, but it increases with the truncation time for the phase variances. Initial studies show that reasonable global Es morphology can be produced on a monthly and seasonal basis with the CHAMP one-antenna occultations. Better results from other existing and upcoming GPS occultation missions are anticipated in future studies, and they will significantly improve our understanding of this important phenomenon.

signal to noise ration (SNR)↗

Dominance of ENSO-Like Variability in Controlling Tropical Ocean Surface Energy Fluxes in the Satellite Era

Ocean surface turbulent and radiative fluxes are critical links in the climate system since they mediate energy exchange between the two fluid systems (ocean and atmosphere) whose combined heat transport determines the basic character of Earth's climate. Moreover, interannual to decadal climate variability depends crucially on the nature of these exchange processes. For example, addressing the question of the degree to which the global hydrologic cycle is changing depends on our ability to observe and model these fluxes accurately. In this work we investigate the interannual to decadal variation of fluxes over the global tropics, especially the tropical oceans. Recent versions of satellite-derived fresh water flux estimates as well as some reanalyses (e.g. products from Remote Sensing Systems, the Woods Hole Oceanographic Institute, and Global Precipitation Climatology Project) suggest that increases in evaporation and precipitation over the past 20 years exceed those expected on the basis of climate model projected responses to greenhouse gas forcing. At the same time, it is well known that E1 Nino / Southern Oscillation behavior in the Pacific exhibits significant variability at scales longer than interannual. We examine here the degree to which surface fluxes attending these interannual to decadal fluctuations are related to ENSO. We examine consistency between these data sets and explore relationships between SST variations, flux changes and modulation of tropical Walker and Hadley circulations.

Robertson, F. R.↗

Turbulent Transport in the Gray Zone: A Large Eddy Model Intercomparison Study of the CONSTRAIN Cold Air Outbreak Case

To quantify the turbulent transport at gray zone length scales between 1 and 10 km the Lagrangian evolution of the CONSTRAIN cold air outbreak (CAO) case was simulated with seven large eddy models. The case is characterized by rather large latent and sensible heat fluxes, and a rapid deepening rate of the boundary layer. In some models the entrainment velocity exceeds 4 cm/s. A significant fraction of this growth is attributed to a strong longwave radiative cooling of the inversion layer. The evolution and the timing of the breakup of the stratocumulus cloud deck differ significantly among the models. Sensitivity experiments demonstrate that a decrease in the prescribed cloud droplet number concentration, and the inclusion of ice microphysics, both act to speed up the thinning of the stratocumulus by enhancing the production of precipitation. In all models the formation of mesoscale fluctuations is clearly evident in the cloud fields but also in the horizontal wind velocity. Resolved vertical fluxes remain important for scales up to 10 km. The simulation results show that the resolved vertical velocity variance gradually diminishes with a coarsening of the horizontal mesh, but the total vertical fluxes of heat, moisture, and momentum are only weakly affected. This is a promising result as it demonstrates the potential use of a mesh size dependent turbulent length scale for convective boundary layers at gray zone model resolutions.

gray zone↗

Penrose pair production as a power source of quasars and active galactic nuclei

Penrose pair production in massive canonical Kerr black holes (those with a/M equal to 0.998) is proposed as a way to explain the nature of the vast fluctuating energy production associated with active galactic nuclei and quasars. It is assumed that a Kerr black hole with a mass of the order of 100 million solar masses lies at the center of an active nucleus and that an accretion disk is formed. Penrose pair production in the inner ergosphere of such a massive canonical Kerr black hole is analyzed. The results indicate that: (1) particle pairs are ejected within a 40 deg angle relative to the equator; (2) the particle energy is of the order of 1 GeV per pair; (3) the pressure of the electron-positron relativistic gas is proportional to the electron-positron number density; (4) pair production may occur in bursts; and (5) the overall lifetime of an active nucleus would depend on the time required to exhaust the disk of its matter content. A test of the theory is suggested which involves observation of the 0.5-MeV pair-annihilation gamma rays that would be generated by annihilating particle pairs.

Kafatos, M.↗

Comparative Ecology of H2 Cycling in Organotrophic and Phototrophic Ecosystems

The simple biochemistry of H2 is critical to a large number of microbial processes, affecting the interaction of organisms with each other and with the environment. The sensitivity of these many processes to H2 can be described quantitatively, at a basic thermodynamic level. This shared dependence on H2 may provide a means for interpreting the ecology and system-level biogeochemistry of widely variant microbial ecosystems on a common (and quantitative) level. Understanding the factors that control H2 itself is a critical prerequisite. Here, we examine two ecosystems that vary widely with respect to H2 cycling. In anoxic, 'organotrophic' sediments from Cape Lookout Bight (North Carolina, USA), H2 partial pressures are strictly maintained at low, steady-state levels by H2-consuming organisms, in a fashion that can be quantitatively predicted by simple thermodynamic calculations. In phototrophic microbial mats from Baja, Mexico, H2 partial pressures are instead controlled by the activity of light-sensitive H2-producing organisms. In consequence, H2 partial pressures within the system fluctuate by orders of magnitude on hour-long time scales. The differences in H2 cycling subsequently impact H2-sensitive microbial processes, such as methanogenesis. For example, the presence of sulfate in the organotrophic system always yielded low levels of H2 that were inhibitory to methanogenesis; however, the elevated levels of H2 in the phototrophic system favored methane production at significant levels, even in the presence of high sulfate concentrations. The myriad of other H2-sensitive microbial processes are expected to exhibit similar behavior.

Hoehler, Tori M.↗

Use of cross-correlation measurements to investigate noise generating regions of a real jet engine and a model jet

Cross-correlations are reported of the jet static pressure fluctuations (as measured with a B and K microphone fitted with a nose cone), with the far-field radiated sound pressure. These measurements were made for various probe positions and a large number of far-field positions (at various angles). In addition, the tests were run for a number of different jet exit velocities. The measured, normalized cross-correlation functions vary between 0.004 and 0.155. These values depend upon the angular position of the far-field microphone, the jet exit Mach number, and the position of the probe. In addition, the cross-correlation technique was employed to study the symmetry of the far-field radiated sound about the jet axis. Third-octave analyses of both the probe signal and the far-field radiated sound were made. This is the first time correlation measurements have been made on a jet engine. In addition, a report is given on an extensive noise survey of a model jet. The correlations are related to sound source functions and jet source regions are discussed.

Meecham, W. C.↗

Shear Thinning Near the Critical Point of Xenon

We measured shear thinning, a viscosity decrease ordinarily associated with complex liquids, near the critical point of xenon. The data span a wide range of reduced shear rate: 10(exp -3) < gamma-dot tau < 700, where gamma-dot tau is the shear rate scaled by the relaxation time tau of critical fluctuations. The measurements had a temperature resolution of 0.01 mK and were conducted in microgravity aboard the Space Shuttle Columbia to avoid the density stratification caused by Earth's gravity. The viscometer measured the drag on a delicate nickel screen as it oscillated in the xenon at amplitudes 3 mu,m < chi (sub 0) >430 mu, and frequencies 1 Hz < omega/2 pi < 5 Hz. To separate shear thinning from other nonlinearities, we computed the ratio of the viscous force on the screen at gamma-dot tau to the force at gamma-dot tau approximates 0: C(sub gamma) is identical with F(chi(sub 0), omega tau, gamma-dot tau )/F)(chi(sub 0, omega tau, 0). At low frequencies, (omega tau)(exp 2) < gamma-dot tau, C(sub gamma) depends only on gamma-dot tau, as predicted by dynamic critical scaling. At high frequencies, (omega tau)(exp 2) > gamma-dot tau, C(sub gamma) depends also on both x(sub 0) and omega. The data were compared with numerical calculations based on the Carreau-Yasuda relation for complex fluids: eta(gamma-dot)/eta(0)=[1+A(sub gamma)|gamma-dot tau|](exp - chi(sub eta)/3+chi(sub eta)), where chi(sub eta) =0.069 is the critical exponent for viscosity and mode-coupling theory predicts A(sub gamma) =0.121. For xenon we find A(sub gamma) =0.137 +/- 0.029, in agreement with the mode coupling value. Remarkably, the xenon data close to the critical temperature T(sub c) were independent of the cooling rate (both above and below T(sub c) and these data were symmetric about T(sub c) to within a temperature scale factor. The scale factors for the magnitude of the oscillator s response differed from those for the oscillator's phase; this suggests that the surface tension of the two-phase domains affected the drag on the screen below T(sub c).

Zimmerli, Gregory A.↗

Permeation of halide anions through phospholipid bilayers occurs by the solubility-diffusion mechanism

Two alternative mechanisms are frequently used to describe ionic permeation of lipid bilayers. In the first, ions partition into the hydrophobic phase and then diffuse across (the solubility-diffusion mechanism). The second mechanism assumes that ions traverse the bilayer through transient hydrophilic defects caused by thermal fluctuations (the pore mechanism). The theoretical predictions made by both models were tested for halide anions by measuring the permeability coefficients for chloride, bromide, and iodide as a function of bilayer thickness, ionic radius, and sign of charge. To vary the bilayer thickness systematically, liposomes were prepared from monounsaturated phosphatidylcholines (PC) with chain lengths between 16 and 24 carbon atoms. The fluorescent dye MQAE (N-(ethoxycarbonylmethyl)-6-methoxyquinolinium bromide) served as an indicator for halide concentration inside the liposomes and was used to follow the kinetics of halide flux across the bilayer membranes. The observed permeability coefficients ranged from 10(-9) to 10(-7) cm/s and increased as the bilayer thickness was reduced. Bromide was found to permeate approximately six times faster than chloride through bilayers of identical thickness, and iodide permeated three to four times faster than bromide. The dependence of the halide permeability coefficients on bilayer thickness and on ionic size were consistent with permeation of hydrated ions by a solubility-diffusion mechanism rather than through transient pores. Halide permeation therefore differs from that of a monovalent cation such as potassium, which has been accounted for by a combination of the two mechanisms depending on bilayer thickness.

NASA Discipline Exobiology↗

Terrestrial Waters and Sea Level Variations on Interannual Time Scale

On decadal to multi-decadal time scales, thermal expansion of sea waters and land ice loss are the main contributors to sea level variations. However, modification of the terrestrial water cycle due to climate variability and direct anthropogenic forcing may also affect sea level. For the past decades, variations in land water storage and corresponding effects on sea level cannot be directly estimated from observations because these are almost non-existent at global continental scale. However, global hydrological models developed for atmospheric and climatic studies can be used for estimating total water storage. For the recent years (since mid-2002), terrestrial water storage change can be directly estimated from observations of the GRACE space gravimetry mission. In this study, we analyse the interannual variability of total land water storage, and investigate its contribution to mean sea level variability at interannual time scale. We consider three different periods that, each, depend on data availability: (1) GRACE era (2003-2009), (2) 1993-2003 and (3) 1955-1995. For the GRACE era (period 1), change in land water storage is estimated using different GRACE products over the 33 largest river basins worldwide. For periods 2 and 3, we use outputs from the ISBA-TRIP (Interactions between Soil, Biosphere, and Atmosphere-Total Runoff Integrating Pathways) global hydrological model. For each time span, we compare change in land water storage (expressed in sea level equivalent) to observed mean sea level, either from satellite altimetry (periods 1 and 2) or tide gauge records (period 3). For each data set and each time span, a trend has been removed as we focus on the interannual variability. We show that whatever the period considered, interannual variability of the mean sea level is essentially explained by interannual fluctuations in land water storage, with the largest contributions arising from tropical river basins.

Llovel, W.↗

The Indirect Effect of UV: Some Good News for Microbes?

Ultraviolet (UV) radiation is of great concern because its biological effects are predominantly harmful. UV damage may be direct or indirect, the latter mediated through the photochemical production of reactive oxygen species such as hydrogen peroxide. We measured the effect of H2O2 on various microbes both in the lab and in nature. At our study site in Yellowstone National Park, there is a UV-induced diurnal fluctuation of H2O2 extending up to one micron. Levels of DNA synthesis resulting from exposure to H2O2 were measured in several algal mats. Within naturally-occurring concentrations of H2O2, DNA synthesis increased. Laboratory studies showed that similar concentrations of H2O2 induce mitosis. We hypothesize that the low levels of H2O2 encountered in nature are inducing mitotic division. At higher levels of H2O2 a second peak in DNA synthesis was found which we interpret to represent DNA damage repair. These experiments suggest that in nature, the low levels of H2O2 produced may have a mitogenic rather than damaging effect. Assuming early levels of UV radiation were substantially higher at the time protists evolved, differential mitogenic effects could have influenced protistan evolution. With H2O2 likely to be present on such bodies as Mars and Europa, the response of organisms will be concentration-dependent and not linear.

Rothschild, Lynn J.↗

Analysis of shadowing effects on MIR photovoltaic and solar dynamic power systems

The NASA Lewis Research Center is currently working with RSC-Energia, the Russian Space Agency, and Allied Signal in developing a flight demonstration solar dynamic power system. This type of power system is dependent upon solar flux that is reflected and concentrated into a thermal storage system to provide the thermal energy input to a closed-cycle Brayton heat engine. The solar dynamic unit will be flown on the Russian Mir space station in anticipation of use on the International Space Station Alpha. By the time the power system is launched, the Mir will be a spatially complex configuration which will have, in addition to the three-gimbaled solar dynamic unit, eleven solar array wings that are either fixed or track the Sun along one axis and a variety or repositionable habitation and experiment modules. The proximity of arrays to modules creates a situation which makes it highly probable that there will be varying solar flux due to shadowing on the solar dynamic unit and some of the arrays throughout the orbit. Shadowing causes fluctuations in the power output from the arrays and the solar dynamic power system, thus reducing the energy capabilities of the spacecraft. An assessment of the capabilities of the power system under these conditions is an important part in influencing the design and operations of the spacecraft and predicting its energy performance. This paper describes the results obtained from using the Orbiting Spacecraft Shadowing Analysis Station program that was integrated into the Station Power Analysis for Capability Evaluation (SPACE) electrical power system computer program. OSSA allows one to consider the numerous complex factors for analyzing the shadowing effects on the electrical power system including the variety of spacecraft hardware geometric configurations, yearly and daily orbital variations in the vehicle attitude and orbital maneuvers (for communications coverage, payload pointing requirements and rendezvous/docking with other vehicles). The geometric models of the MIR with a solar dynamic power unit that were used in performing shadowing analyses are described. Also presented in this paper are results for individual orbits for several flight attitude cases which include assessments of the shadowing impacts upon the solar dynamic unit and the solar arrays. These cases depict typical MIR flight attitudes likely to have shadowing impact. Because of the time varying nature of the Mir orientation with respect to the Sun and the lack of knowledge of the precise timing of the attitude changes, strategies must be devised to assess and depict the shadowing impacts on power generation throughout the year. To address this, the best, nominal and worst impacts of shadowing considering a wide possible range of parameter changes for typical mission operation period are shown.

Fincannon, James↗

Improved Estimates of Clear Sky Longwave Flux and Application to the Tropical Greenhouse Effect

The first objective of this investigation is to eliminate the clear-sky offset introduced by the scene-identification procedures developed for the Earth Radiation Budget Experiment (ERBE). Estimates of this systematic bias range from 10 to as high as 30 W/sq m. The initial version of the ScaRaB data is being processed with the original ERBE algorithm. Since the ERBE procedure for scene identification is based upon zonal flux averages, clear scenes with longwave emission well below the zonal mean value are mistakenly classified as cloudy. The erroneous classification is more frequent in regions with deep convection and enhanced mid- and upper-tropospheric humidity. We will develop scene identification parameters with zonal and/or time dependence to reduce or eliminate the bias in the clear- sky data. The modified scene identification procedure could be used for the ScaRaB-specific version of the Earth-radiation products. The second objective is to investigate changes in the clear-sky Outgoing Longwave Radiation (OLR) associated with decadal variations in the tropical and subtropical climate. There is considerable evidence for a shift in the climate state starting in approximately 1977. The shift is accompanied by higher SSTs in the equatorial Pacific, increased tropical convection, and higher values of atmospheric humidity. Other evidence indicates that the humidity in the tropical troposphere has been steadily increasing over the last 30 years. It is not known whether the atmospheric greenhouse effect has increased during this period in response to these changes in SST and precipitable water. We will investigate the decadal-scale fluctuations in the greenhouse effect using Nimbus-7, ERBE, and ScaRaB measurements spaning 1979 to the present. The data from the different satellites will be intercalibrated by comparison with model calculations based upon ship radiosonde observations. The fluxes calculated from the radiation model will also be used for validation of the ScaRaB fluxes.

Collins, W. D.↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

Solar Forced Dansgaard/Oeschger Events?

Climate records for the last ice age (which ended 11,500 years ago) show enormous climate fluctuations in the North Atlantic region - the so-called Dansgaard/Oeschger events. During these events air temperatures in Greenland changed on the order of 10 degrees Celsius within a few decades. These changes were attributed to shifts in ocean circulation which influences the warm water supply from lower latitudes to the North Atlantic region. Interestingly, the rapid warmings tend to recur approximately every 1500 years or multiples thereof. This has led researchers to speculate about an external cause for these changes with the variable Sun being one possible candidate. Support for this hypothesis came from climate reconstructions, which suggested that the Sun influenced the climate in the North Atlantic region on these time scales during the last approximately 12,000 years of relatively stable Holocene climate. However, Be-10 measurements in ice cores do not indicate that the Sun caused or triggered the Dansgaard/Oeschger events. Depending on the solar magnetic shielding more or less Be-10 is produced in the Earth's atmosphere. Therefore, 10Be can be used as a proxy for solar activity changes. Since Be-10 can be measured in ice cores, it is possible to compare the variable solar forcing directly with the climate record from the same ice core. This removes any uncertainties in the relative dating, and the solar-climate link can be reliably studied. Notwithstanding that some Dansgaard/Oeschger warmings could be related to increased solar activity, there is no indication that this is the case for all of the Dansgaard/Oeschger events. Therefore, during the last ice age the Be-10 and ice core climate data do not indicate a persistent solar influence on North Atlantic climate.

Muscheler, R.↗

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗