Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Decoherence and Exponential Law: A Solvable Model

We analyze a modified version of the 'AgBr' Hamiltonian, solve exactly the equations of motion in terms of SU(2) coherent states, and study the weak-coupling, macroscopic limit of the model, obtaining an exponential behavior at all times. The asymptotic dominance of the exponential behavior is representative of a purely stochastic evolution and can be derived quantum mechanically in the so-called van Hove's limit (which is a weak-coupling, macroscopic limit). At the same time, a temporal behavior of the exponential type, yielding a 'probability dissipation' is closely related to dephasing ('decoherence') effects and one can expect a close connection with a dissipative and irreversible behavior. We stress the central relevance of the problem of dissipation to the quantum measurement theory and to the general topic of decoherence.

Pascazio, Saverio↗

Optical amplifiers for coherent lidar

We examine application of optical amplification to coherent lidar for the case of a weak return signal (a number of quanta of the return optical field close to unity). We consider the option that has been explored to date, namely, incorporation of an optical amplifier operated in a linear manner located after reception of the signal and immediately prior to heterodyning and photodetection. We also consider alternative strategies where the coherent interaction, the nonlinear processes, and the amplification are not necessarily constrained to occur in the manner investigated to date. We include the complications that occur because of mechanisms that occur at the level of a few, or one, quantum excitation. Two factors combine in the work to date that limit the value of the approach. These are: (1) the weak signal tends to require operation of the amplifier in the linear regime where the important advantages of nonlinear optical processing are not accessed, (2) the linear optical amplifier has a -3dB noise figure (SN(out)/SN(in)) that necessarily degrades the signal. Some improvement is gained because the gain provided by the optical amplifier can be used to overcome losses in the heterodyned process and photodetection. The result, however, is that introduction of an optical amplifier in a well optimized coherent lidar system results in, at best, a modest improvement in signal to noise. Some improvement may also be realized on incorporating more optical components in a coherent lidar system for purely practical reasons. For example, more compact, lighter weight, components, more robust alignment, or more rapid processing may be gained. We further find that there remain a number of potentially valuable, but unexplored options offered both by the rapidly expanding base of optical technology and the recent investigation of novel nonlinear coherent interference phenomena occurring at the single quantum excitation level. Key findings are: (1) insertion of linear optical amplifiers in well optimized conventional lidar systems offers modest improvements, at best, (2) the practical advantages of optical amplifiers, especially fiber amplifiers, such as ease of alignment, compactness, efficiency, lightweight, etc., warrant further investigation for coherent lidar, (3) the possibility of more fully optical lidar systems should be explored, (4) advantages gained by use of coherent interference of optical fields at the level of one, or a few, signal quanta should be explored, (5) amplification without inversion, population trapping, and use of electromagnetic induced transparency warrant investigation in connection with coherent lidar, (6) these new findings are probably more applicable to earth related NASA work, although applications to deep space should not be excluded, and (7) our own work in the Ultrafast Laboratory at UAH along some of the above lines of investigation, may be useful.

Fork, Richard↗

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗

Effects of Spatial Gradients on Electron Runaway Acceleration

The runaway process is known to accelerate electrons in many laboratory plasmas and has been suggested as an acceleration mechanism in some astrophysical plasmas, including solar flares. Current calculations of the electron velocity distributions resulting from the runaway process are greatly restricted because they impose spatial homogeneity on the distribution. We have computed runaway distributions which include consistent development of spatial gradients in the energetic tail. Our solution for the electron velocity distribution is presented as a function of distance along a finite length acceleration region, and is compared with the equivalent distribution for the infinitely long homogenous system (i.e., no spatial gradients), as considered in the existing literature. All these results are for the weak field regime. We also discuss the severe restrictiveness of this weak field assumption.

MacNeice, Peter↗

Performance of Optimized Actuator and Sensor Arrays in an Active Noise Control System

Experiments have been conducted in NASA Langley's Acoustics and Dynamics Laboratory to determine the effectiveness of optimized actuator/sensor architectures and controller algorithms for active control of harmonic interior noise. Tests were conducted in a large scale fuselage model - a composite cylinder which simulates a commuter class aircraft fuselage with three sections of trim panel and a floor. Using an optimization technique based on the component transfer functions, combinations of 4 out of 8 piezoceramic actuators and 8 out of 462 microphone locations were evaluated against predicted performance. A combinatorial optimization technique called tabu search was employed to select the optimum transducer arrays. Three test frequencies represent the cases of a strong acoustic and strong structural response, a weak acoustic and strong structural response and a strong acoustic and weak structural response. Noise reduction was obtained using a Time Averaged/Gradient Descent (TAGD) controller. Results indicate that the optimization technique successfully predicted best and worst case performance. An enhancement of the TAGD control algorithm was also evaluated. The principal components of the actuator/sensor transfer functions were used in the PC-TAGD controller. The principal components are shown to be independent of each other while providing control as effective as the standard TAGD.

Palumbo, D. L.↗

Turbulent Fluid Motion 6: Turbulence, Nonlinear Dynamics, and Deterministic Chaos

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The turbulence is sustained by a nonrandom time-independent external force. The solutions, on the average, separate exponentially with time, having a positive Liapunov exponent. Thus, the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number (or strength of the forcing) is reduced. Several qualitatively different flows are noted. These are, respectively, fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, we classify only the fully chaotic flows as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows, although having positive Liapunov exponents, have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, as generally understood, is both time-dependent and aperiodic.

Deissler, Robert G.↗

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(exp -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗

The Mid-Infrared Transmission Spectra of Antarctic Ureilites

The mid-infrared (4000-450 1/cm; 2.5-22.2 micron) transmission spectra of seven Antarctic ureilites and 10 Antarctic H-5 ordinary chondrites are presented. The ureilite spectra show a number of absorption bands, the strongest of which is a wide, complex feature centered near 1000 1/cm (10 micron) due to Si-O stretching vibrations in silicates. The profiles and positions of the substructure in this feature indicate that Mg-rich olivines and pyroxenes are the main silicates responsible. The relative abundances of these two minerals, as inferred from the spectra, show substantial variation from meteorite to meteorite, but generally indicate olivine is the most abundant (olivine:pyroxene = 60:40 to 95:5). Both the predominance of olivine and the variable olivine-to-pyroxene ratio are consistent with the known composition and heterogeneity of ureilites. The H-5 ordinary chondrites spanned a range of weathering classes and were used to provide a means of addressing the extent to which the ureilite spectra may have been altered by weathering processes. It was found that, while weathering of these meteorites produces some weak bands due to the formation of small amounts of carbonates and hydrates, the profile of the main silicate feature has been little affected by Antarctic exposure in the meteorites studied here. The mid-infrared ureilite spectra provide an additional means of testing potential asteroidal parent bodies for the ureilites. At present, the best candidates include the subset of S-type asteroids having low albedos and weak absorption features in the near infrared.

Sandford, Scott A.↗

Structure and Dynamics of Coronal Plasmas

During the past year this grant has funded research in the interaction between magnetic fields and the hot plasma in the solar outer atmosphere. The following is a brief summary of the published papers, abstracts and talks which have been supported. The paper 'Coronal Structures Observed in X-rays and H-alpha Structures' was published in the Kofu Symposium proceedings. The study analyzes cool and hot behavior of two x-ray events, a small flare and a surge. We find that a large H-alpha surge appears in x-rays as a very weak event, while a weak H-alpha feature corresponds to the brightest x-ray emission on the disk at the time of the observation. Calculations of the heating necessary to produce these signatures, and implications for the driving and heating mechanisms of flares vs. surges are presented. A copy of the paper is appended to this report. The paper 'Differential Magnetic Field Shear in an Active Region' has been published in The Astrophysical Journal. We have compared the 3D extrapolation of magnetic fields with the observed coronal structure in an active region. Based on the fit between observed coronal structure throughout the volume of the region and the calculated magnetic field configurations, we propose a differential magnetic field shear model for this active region. The decreasing field shear in the outer portions of the AR may indicate a continual relaxation of the magnetic field with time, corresponding to a net transport of helicity outward. The paper 'Difficulties in Observing Coronal Structure' has been accepted for publication in the journal Solar Physics. In this paper we discuss the evidence that the temperature and density structure of the corona are far more complicated than had previously been thought. The discussion is based on five studies carried out by our group on coronal plasma properties, showing that any one x-ray instrument does see all of the plasma present in the corona, that hot and cool material may appear to be co-spatial at a given location in the corona, and that simple magnetic field extrapolations provide only a poor fit to the observed structure. A copy of the paper is appended to this report.

Golub, Leon↗

Neutral-Line Magnetic Shear and Enhanced Coronal Heating in Solar Active Regions

By examining the magnetic structure at sites in the bright coronal interiors of active regions that are not flaring but exhibit persistent strong coronal heating, we establish some new characteristics of the magnetic origins of this heating. We have examined the magnetic structure of these sites in five active regions, each of which was well observed by both the Yohkoh SXT and the Marshall Space Flight Center Vector Magnetograph and showed strong shear in its magnetic field along part of at least one neutral line (polarity inversion). Thus, we can assess whether this form of nonpotential field structure in active regions is a characteristic of the enhanced coronal heating and vice versa. From 27 orbits' worth of Yohkoh SXT images of the five active regions, we have obtained a sample of 94 persistently bright coronal features (bright in all images from a given orbit), 40 long (greater than or approximately equals 20,000 km) neutral-line segments having strong magnetic shear throughout (shear angle greater than 45 deg), and 39 long neutral-line segments having weak magnetic shear throughout (shear angle less than 45 deg). From this sample, we find that: (1) all of our persistently bright coronal features are rooted in magnetic fields that are stronger than 150 G; (2) nearly all (95%) of these enhanced coronal features are rooted near neutral lines (closer than 10,000 km); (3) a great majority (80%) of the bright features are rooted near strong-shear portions of neutral lines; (4) a great majority (85%) of long strong-shear segments of neutral lines have persistently bright coronal features rooted near them; (5) a large minority (40%) of long weak-shear segments of neutral lines have persistently bright coronal features rooted near them; and (6) the brightness of a persistently bright Coronal feature often changes greatly over a few hours. From these results, we conclude that most persistent enhanced heating of coronal loops in active regions: (1) requires the presence of a polarity inversion in the magnetic field near at least one of the loop footpoints; (2) is greatly aided by the presence of strong shear in the core magnetic field along that neutral line; and (3) is controlled by some variable process that acts in this magnetic environment. We infer that this variable process is low-lying reconnection accompanying flux cancellation.

Falconer, D. A.↗

Enhancing Adhesion: Relative Merits of Different Approaches

Adhesive performance is improved mainly by manipulation of the bimaterials interface zone, which is only a few molecules thick. There are three approaches to enhancement of interfacial adhesion at the molecular level. They are 1) changing the nonchemically bonded interactions across the interface from weak ones to strong ones, 2) making the true interfacial area much larger than the simple geometric area, and 3) inducing chemical bonding between the two materials forming the interface. Our goal this summer was to question some of the built-in assumptions contained within these approaches and to determine the most promising approach, both theoretically and practically, for enhancing adhesion in NASA structures. Our computations revealed that all three of these approaches have, in theory, the potential to enhance molecular adhesion approximately ten-fold. Experiments, however, revealed that this excellent level of enhancement is not likely to be reached in practice. Each approach was found to be severely limited by practical problems. In addition, some of the built-in assumptions associated with these approaches were found to be insufficient or inadequate. The first approach, changing the nonchemically bonded interactions from weak to strong, Is an example of one containing inadequate assumptions. The extensive literature on intermolecular interactions, based on solution studies, shows that certain functional group pairs interact much more strongly than others. It has always been assumed that these data can be reliably extended to systems where only one member of the pair is in solution and the other Is contained in a solid surface. Our experiments this summer demonstrated that solution data do not adequately predict the strength of functional group interaction at the solid-liquid interface. Furthermore, the strong solvents needed to dissolve the monomers or polymers to which the functional groups of interest are attached compete successfully with the solid surface for the functional group. As a result, functional groups in solution cannot pair with the complementary groups in the solid surface, and the expected enhancement of nonchemically bonded interactions is not realized. The second approach, increasing the true interfacial area, is an example of one containing inadequate assumptions and suffering from numerous practical problems. First, practitioners have assumed that material removal, such as bead blasting or etching, increases true surface area (and therefore interfacial area) in a meaningful way. Our geometric analysis demonstrated that removal methods increase area by a factor of two at most. To increase interfacial area by an order of magnitude or more. a thin layer of high porosity must be added to the substrate surface prior to application of the adhesive phase. Consistent with this finding, we attempted to create a thin layer of rigid, highly porous glass on the surface of our smooth glass substrate by means of sol-gel technology. We were unable to surmount a wide variety of practical problems and obtained only collapsed, nonporous layers. Thus this approach, appealing in principle, would require long term development and is not promising in the near term. The third approach, inducing chemical bonding at the interface, is an example of one having neither inadequate assumptions nor insurmountable practical problems. When silicate glass is the substrate, there are only a few chemical reactions that can be successfully conducted to create these chemical bonds, and these reactions usually involve silicon-containing reagents. We compared the silazane reagents to the silane reagents and found through experiment that the silazanes react with the glass surface much more readily, and under milder conditions, than the silanes. The functional groups attached to the glass surface by silazane reactions were not able to be removed by solvent extraction, elevated temperature exposure, or mechanical action. This clearly indicates that the formation of chemical bonds at the interface is the most effective approach for enhancing molecular adhesion.

Penn, L. S.↗

The source and distribution of Galactic positrons

The oriented scintillation spectrometer experiment (OSSE) observations of the Galactic plane and the Galactic center region were combined with observations acquired with other instruments in order to produce a map of the Galactic 511 keV annihilation radiation. Two mapping techniques were applied to the data: the maximum entropy method, and the basis pursuit inversion method. The resulting maps are qualitatively similar and show evidence for a central bulge and a weak galactic disk component. The weak disk is consistent with that expected from positrons produced by the decay of radioactive Al-26 in the interstellar medium. Both maps suggest an enhanced region of emission near l = -4 deg, b = 7 deg, with a flux of approximately 50 percent of that of the bulge. The existence of this emission appears significant, although the location is not well determined. The source of this enhanced emission is presently unknown.

Purcell, W. R.↗

Hard X-ray spectra of neutron stars and black hole candidates

The hard X-ray behavior of several X-ray binary systems containing a neutron star or a black hole candidate is analyzed in an attempt to determine the specific signature of these categories of compact objects. Limiting the consideration to two subclasses of neutron stars, Atoll sources and non-pulsating Z sources, it appears that only the Atoll sources have a spectral behavior similar to black holes. It is proposed that Atoll sources are weakly magnetized neutron stars, whereas Z sources are small radius moderate magnetized neutron stars. Large magnetic fields funnel the accreting matter, thus preventing spherical accretion and free fall if the neutron star radius is smaller than the last stable accreting orbit. Weak magnetic fields do not have this effect, and blackbody soft photons from the stellar surface are upscattered on the relativistic infalling matter, leading to excess hard X-rays. This excess is visible in two of the observed Atoll sources and in the spectrum of a black hole candidate. In the case of a Z source, a lack of photons was remarked, providing a possible signature to distinguish between these classes of objects.

Durouchoux, P.↗

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(sup -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗

Numerical Computation of Sensitivities and the Adjoint Approach

We discuss the numerical computation of sensitivities via the adjoint approach in optimization problems governed by differential equations. We focus on the adjoint problem in its weak form. We show how one can avoid some of the problems with the adjoint approach, such as deriving suitable boundary conditions for the adjoint equation. We discuss the convergence of numerical approximations of the costate computed via the weak form of the adjoint problem and show the significance for the discrete adjoint problem.

Lewis, Robert Michael↗

Advanced High Brilliance X-Ray Source

The possibility to dramatically increase the efficiency of laboratory based protein structure measurements through the use of polycapillary X-ray optics was investigated. This project initiated April 1, 1993 and concluded December 31, 1996 (including a no cost extension from June 31, 1996). This is a final report of the project. The basis for the project is the ability to collect X-rays from divergent electron bombardment laboratory X-ray sources and redirect them into quasiparallel or convergent (focused) beams. For example, a 0.1 radian (approx. 6 deg) portion of a divergent beam collected by a polycapillary collimator and transformed into a quasiparallel beam of 3 millradian (0.2 deg) could give a gain of 6(exp 2)/0.2(exp 2) x T for the intensity of a diffracted beam from a crystal with a 0.2 deg diffraction width. T is the transmission efficiency of the polycapillary diffraction optic, and for T=0.5, the gain would be 36/0.04 x O.5=45. In practice, the effective collection angle will depend on the source spot size, the input focal length of the optic (usually limited by the source spot-to-window distance on the x-ray tube) and the size of the crystal relative to the output diameter of the optic. The transmission efficiency, T, depends on the characteristics (fractional open area, surface roughness, shape and channel diameter) of the polycapillary optic and is typically in the range 0.2-0.4. These effects could substantially reduce the expected efficiency gain. During the course of this study, the possibility to use a weakly focused beam (0.5 deg convergence) was suggested which could give an additional 10-20 X efficiency gain for small samples . Weakly focused beams from double focusing mirrors are frequently used for macromolecular crystallography studies. Furthermore the crystals are typically oscillated by as much as 2 deg during each X-ray exposure in order to increase the reciprocal space (number of crystal planes) sampled and use of a slightly convergent beam could, in principle, provide a similar sampling benefit without oscillation. Although more problematic, because of complications in analyzing the diffraction patterns, it was also suggested that even more extreme beam convergence might be used to give another order of magnitude intensity gain and even smaller focused spot size which could make it possible to study smaller protein crystals than can be studied using standard laboratory based X-ray diffraction systems. This project represents the first systematic investigation of these possibilities. As initially proposed, the contract included requirements for design, purchase, evaluation and delivery of three polycapillary lenses to the Laboratory for Structural Biology at MSFC and demonstration of such optics at MSFC for selected protein crystal diffraction applications.

Gibson, Walter M.↗

Some Results Relevant to Statistical Closures for Compressible Turbulence

For weakly compressible turbulent fluctuations there exists a small parameter, the square of the fluctuating Mach number, that allows an investigation using a perturbative treatment. The consequences of such a perturbative analysis in three different subject areas are described: 1) initial conditions in direct numerical simulations, 2) an explanation for the oscillations seen in the compressible pressure in the direct numerical simulations of homogeneous shear, and 3) for turbulence closures accounting for the compressibility of velocity fluctuations. Initial conditions consistent with small turbulent Mach number asymptotics are constructed. The importance of consistent initial conditions in the direct numerical simulation of compressible turbulence is dramatically illustrated: spurious oscillations associated with inconsistent initial conditions are avoided, and the fluctuating dilatational field is some two orders of magnitude smaller for a compressible isotropic turbulence. For the isotropic decay it is shown that the choice of initial conditions can change the scaling law for the compressible dissipation. A two-time expansion of the Navier-Stokes equations is used to distinguish compressible acoustic and compressible advective modes. A simple conceptual model for weakly compressible turbulence - a forced linear oscillator is described. It is shown that the evolution equations for the compressible portions of turbulence can be understood as a forced wave equation with refraction. Acoustic modes of the flow can be amplified by refraction and are able to manifest themselves in large fluctuations of the compressible pressure.

Ristorcelli, J. R.↗