Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 163 records · Page 9

Influence of transient baroclinic eddies on planetary-scale waves

Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.

Branscome, L. E.↗

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

Babin, A.↗

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

Atmospheric Diabatic Heating in Different Weather States and the General Circulation

Analysis of multiple global satellite products identifies distinctive weather states of the atmosphere from the mesoscale pattern of cloud properties and quantifies the associated diabatic heating/cooling by radiative flux divergence, precipitation, and surface sensible heat flux. The results show that the forcing for the atmospheric general circulation is a very dynamic process, varying strongly at weather space-time scales, comprising relatively infrequent, strong heating events by ''stormy'' weather and more nearly continuous, weak cooling by ''fair'' weather. Such behavior undercuts the value of analyses of time-averaged energy exchanges in observations or numerical models. It is proposed that an analysis of the joint time-related variations of the global weather states and the general circulation on weather space-time scales might be used to establish useful ''feedback like'' relationships between cloud processes and the large-scale circulation.

Cooling↗

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗

Recurrent episodic mass loss in a B2e star - 66 Ophiuchi, 1982-1985

The results of three intensive ultraviolet monitoring campaigns for 66 Oph in September-October 1982, September 1983, and September 1984, supplemented by the available IUE archival data from April, 1982 to June, 1985 are reported. The stellar wind went through at least two and probably three episodes of weak-wind absorption. At wind minimum, the resonance profiles and equivalent widths of C IV, Si IV, Si III, and Al III are similar to those observed in normal B stars of comparable luminosity, effective temperature, and v sin i. The strong-wind absorption profiles show one or more shortward-shifted discrete components superposed upon the weak-wind absorption profiles. Significant profile variations are seen on time scales spanning over two orders of magnitude, ranging from small-amplitude changes detected on time scales of hours to the transitions from wind minimum to maximum, which take months.

Grady, C. A.↗

Cosmic gamma-ray bursts

A review of the cosmic gamma-ray burst phenomenon is presented. Both the light curves and the energy spectra of these short transient events display a great diversity. However, rapid rise times and periodicities sometimes observed in the light curves suggest a compact object origin. Similarly, absorption and emission features in the energy spectra argue strongly in favor of this interpretation. Counterparts to gamma-bursters in other energy ranges, such as optical and sort x-ray, have still not been identified, however, leading to a large uncertainty in the distances to bursters. Although gamma-ray burst sources have not yet been observed to repeat, numerous bursts from three objects which may be related to the gamma-bursters, called Soft Gamma Repeaters, have been recorded; there is weak evidence that they may be relatively distant on a galactic scale. Future missions, particularly those emphasizing high energy, time, and/or spatial resolution, as well as a multiwavelength approach, are likely to advance our understanding of this enigmatic phenomenon.

Hurley, K.↗

Weakly nonlinear models for turbulent mixing in a plane mixing layer

New closure models for turbulent free shear flows are presented in this paper. They are based on a weakly nonlinear theory with a description of the dominant large-scale structures as instability waves. Two models are presented that describe the evolution of the free shear flows in terms of the time-averaged mean flow and the dominant large-scale turbulent structure. The local characteristics of the large-scale motions are described using linear theory. Their amplitude is determined from an energy integral analysis. The models have been applied to the study of an incompressible mixing layer. For both models, predictions of the mean flow developed are made. In the second model, predictions of the time-dependent motion of the large-scale structures in the mixing layer are made. The predictions show good agreement with experimental observations.

Liou, William W.↗

Transport Coefficients in weakly compressible turbulence

A theory of transport coefficients in weakly compressible turbulence is derived by applying Yoshizawa's two-scale direct interaction approximation to the compressible equations of motion linearized about a state of incompressible turbulence. The result is a generalization of the eddy viscosity representation of incompressible turbulence. In addition to the usual incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluctuations also generate an effective turbulent mean pressure and corrections to the speed of sound. Finally, a prediction unique to Yoshizawa's two-scale approximation is that terms containing gradients of incompressible turbulence quantities also appear in the mean flow equations. The form these terms take is described.

Rubinstein, Robert↗

Hydrodynamic Force on a Cylinder Oscillating at Low Frequency

The hydrodynamic force on a cylinder oscillating transversely to its axis is a nonlinear function of the displacement amplitude x0. We report measurements and numerical calculations of the force at frequencies low enough that delta > R, where delta is the viscous penetration length and R is the cylinder radius. For small amplitudes, the numerically calculated Fourier transform of the force per unit length, F(sub small), agrees with Stokes' analytical calculation. For larger amplitudes, the force per unit length found by both calculation and measurement is F = F(sub small)C (x(sub 0)/delta,R/delta). The complex function C depends only weakly on R/delta, indicating that x0/delta is more appropriate as a scaling variable than the Keulegan-Carpenter number KC = pi*x(sub 0)/R. The measurements used a torsion oscillator driven at frequencies from 1 to 12 Hz while immersed in dense xenon. The oscillator comprised cylinders with an effective radius of R = 13.4 micron and oscillation amplitudes as large as x(sub 0)/delta = 4 (corresponding to KC as large as 71). The calculations used similar conditions except that the amplitudes were as large as x0/delta = 28.

Berg, Robert F.↗

Surface Height and Sea Ice Freeboard of the Arctic Ocean from ICESat-2: Characteristics and Early Results

We present the first winter season of surface height and sea ice freeboards of the Arctic Ocean from the new ICESat-2 (IS-2) mission. The Advanced Topographic Laser Altimeter System (ATLAS) onboard has six photon-counting beams for surface profiling with a 10-kHz pulse-rate (inter-pulse distance ~0.7m) and footprints of ~17m. Geolocated heights assigned to individual photons scattered from the surface allow significant flexibility in the construction of height distributions used in surface finding. For IS-2 sea ice products, a fixed 150-photon aggregate is used to control height precision and obtain better along-track resolution over high reflectance surfaces. Quasi-specular returns in openings as narrow as ~27m, crucial for freeboard calculations, are resolved. The fixed photon aggregate results in unique variable along-track resolutions and non-uniform sampling (17m x 27m to 17m x 200m for the strong beams) of the surface. The six profiling beams – three pairs separated by 3.3 km with a strong and weak beam in each pair – provide correlated statistics at regional length scales for assessment of beam-to-beam retrieval consistency and accuracy. Analysis shows along-track height precisions of ~2cm, and agreement in the monthly freeboard distributions across the strong beams to 1-2cm. In this paper, we describe briefly the approaches used in surface height and freeboard retrievals from ATLAS photon clouds, and detail the key features of these along-track sea ice products, focusing on the first release of data collected over the Arctic Ocean, which spans the period between October 14, 2018 – the start of data collection – and the end of March 2019.

R Kwok↗

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

Analytical estimates of turbulent MHD transport coefficients

Turbulent transfer rates from small-scale MHD excitations to large-scale Fourier modes are calculated algebraically, using the method of Biskamp and Welter. Three cases are considered: two-dimensional Navier-Stokes flows, two-dimensional incompressible MHD, and the weakly three-dimensional Strauss equations. In all cases, an initially large spectral gap between the small-scale and large-scale excitations is assumed, and attention focusses on the initial values of the back-transfer rates. The sign of the transfer is determined by the sign of an analytically calculable eddy viscosity and/or anomalous resistivity. We are able to confirm the results of Biskamp and Welter for the case of two-dimensional MHD, but find some differences for the case of the Strauss equations. It is argued that the Strauss equations may not exhibit an inverse cascade phenomenon for the spatially periodic case unless their initial spectra are such that the behavior is essentially that of two-dimensional MHD.

Montgomery, D.↗

Structure and dating errors in the geologic time scale and periodicity in mass extinctions

Structure in the geologic time scale reflects a partly paleontological origin. As a result, ages of Cenozoic and Mesozoic stage boundaries exhibit a weak 28-Myr periodicity that is similar to the strong 26-Myr periodicity detected in mass extinctions of marine life by Raup and Sepkoski. Radiometric dating errors in the geologic time scale, to which the mass extinctions are stratigraphically tied, do not necessarily lessen the likelihood of a significant periodicity in mass extinctions, but do spread the acceptable values of the period over the range 25-27 Myr for the Harland et al. time scale or 25-30 Myr for the DNAG time scale. If the Odin time scale is adopted, acceptable periods fall between 24 and 33 Myr, but are not robust against dating errors. Some indirect evidence from independently-dated flood-basalt volcanic horizons tends to favor the Odin time scale.

Stothers, Richard B.↗

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

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

Heelis, R. A.↗

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗

Small-scale Interstellar Structure Toward the Open Cluster Chi Persei; Program ID: BO29 (Cycle 2)

The primary objective of this program was to obtain FUSE observations of the interstellar H2 absorption toward six early-type stars in the core of the open star cluster Chi Persei. High resolution optical observations of the interstellar Na I absorption toward these stars and others in the core of Chi Per have revealed a rich variety of small-scale diffuse ISM structure in both the distant Perseus and more local Orion spiral arms. At the 1 arc minute angular scales typically probed by the separations of these stars, this structure corresponds to respective physical length scales of 0.6 and 0.15 pc in the Perseus and Orion arm gas. The principal scientific goal of the FUSE observations was to determine if this small-scale Na I structure toward the core of Chi Per was reflected in the abundance and/or excitation of interstellar H2. In September 2001, we obtained FUSE observations of the Chi Per stars BD +56 563, BD +56 571, BD +56 573, BD +56 574, BD +56 575, and BD +56 578 for this program. The data reveal that the Perseus arm gas exhibits significant small-scale variations in the strength of its molecular hydrogen absorption while the Orion arm H2 absorption is relatively uniform. We presented these results at the January 2003 American Astronomical Society meeting in Seattle, WA (points, Meyer, Lauroesch, Nguyen, and Friedman 2002, Bulletin of the American Astronomical Society, 34,1229). Since Na I is not a dominant ion in H I clouds, the variations detected in the Perseus arm Na I line profiles could be due to spatial differences in either the physical structure or the physical conditions in the gas. Although many of the H2 absorption profiles are too saturated or too weak to be definitive in this regard, it is clear that small-scale variations are present in H2 gas with excitation as high as the J=4 rotational level. In contrast, lower spatial resolution H I 21 cm studies of the Perseus arm gas across the face of the h and Chi Per double cluster have shown smoother, continuous variations on the larger scale. In order to investigate if the global velocity structure of H2 toward h and Chi Per more closely follows the chaotic structure seen in interstellar Na I absorption or the smoother, continuous H I structure, we have expanded our FUSE study with a Cycle 4 program (D902, PI: Points) to survey the H2 absorption toward a number of other stars across the face of the double cluster. We have obtained FUSE spectra of 14 additional h and Chi Per stars through this program and are now combining all of the data into a comprehensive study that will be submitted for publication next year.

Meyer, David M.↗

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗