SEARCH · Search NASA
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.
The method for X-ray wavelength tuning using an adaptive element (a quartz monolithic resonator modulated by a standing longitudinal low-frequency ultrasonic wave) has been proposed and implemented. This method is characterized by high tuning accuracy and possibility of scanning X-ray beam parameters with a time resolution of up to 3 μs.
This project was a continuation of the project entitled "Resolution Earth Surface Features from Repeat Moderate Resolution Satellite Imagery". In the previous study, a Bayesian Maximum Posterior Estimate (BMPE) algorithm was used to obtain a composite series of repeat imagery from the Advanced Very High Resolution Radiometer (AVHRR). The spatial resolution of the resulting composite was significantly greater than the 1 km resolution of the individual AVHRR images. The BMPE algorithm utilized a simple, no-atmosphere geometrical model for the short-wave radiation budget at the Earth's surface. A necessary assumption of the algorithm is that all non geometrical parameters remain static over the compositing period. This assumption is of course violated by temporal variations in both the surface albedo and the atmospheric medium. The effect of the albedo variations is expected to be minimal since the variations are on a fairly long time scale compared to the compositing period, however, the atmospheric variability occurs on a relatively short time scale and can be expected to cause significant errors in the surface reconstruction. The current project proposed to incorporate an atmospheric correction into the BMPE algorithm for the purpose of investigating the effects of a variable atmosphere on the surface reconstructions. Once the atmospheric effects were determined, the investigation could be extended to include corrections various cloud effects, including short wave radiation through thin cirrus clouds. The original proposal was written for a three year project, funded one year at a time. The first year of the project focused on developing an understanding of atmospheric corrections and choosing an appropriate correction model. Several models were considered and the list was narrowed to the two best suited. These were the 5S and 6S shortwave radiation models developed at NASA/GODDARD and tested extensively with data from the AVHRR instrument. Although the 6S model was a successor to the 5S and slightly more advanced, the 5S was selected because outputs from the individual components comprising the short-wave radiation budget were more easily separated. The separation was necessary since both the 5S and 6S did not include geometrical corrections for terrain, a fundamental constituent of the BMPE algorithm. The 5S correction code was incorporated into the BMPE algorithm and many sensitivity studies were performed.
The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW–LISM interaction. By performing three-dimensional, MHD plasma/kinetic neutral atom simulations, we determine the width of the outer heliosheath—the LISM plasma region affected by the presence of the heliosphere—and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by Interstellar Boundary Explorer .
Planar arrays of waveguide-fed slots have been employed in many radar and remote sensing applications. Such arrays are designed in the standing wave configuration because of high efficiency. Traveling wave arrays can produce greater bandwidth at the expense of efficiency due to power loss in the load or loads. Traveling wave planar slot arrays may be designed with a long feed waveguide consisting of centered-inclined coupling slots. The feed waveguide is terminated in a matched load, and the element spacing in the feed waveguide is chosen to produce a beam squinted from the broadside. The traveling wave planar slot array consists of a long feed waveguide containing resonant-centered inclined coupling slots in the broad wall, coupling power into an array of stacked radiating waveguides orthogonal to it. The radiating waveguides consist of longitudinal offset radiating slots in a standing wave configuration. For the traveling wave feed of a planar slot array, one has to design the tilt angle and length of each coupling slot such that the amplitude and phase of excitation of each radiating waveguide are close to the desired values. The coupling slot spacing is chosen for an appropriate beam squint. Scattering matrix parameters of resonant coupling slots are used in the design process to produce appropriate excitations of radiating waveguides with constraints placed only on amplitudes. Since the radiating slots in each radiating waveguide are designed to produce a certain total admittance, the scattering (S) matrix of each coupling slot is reduced to a 2x2 matrix. Elements of each 2x2 S-matrix and the amount of coupling into the corresponding radiating waveguide are expressed in terms of the element S11. S matrices are converted into transmission (T) matrices, and the T matrices are multiplied to cascade the coupling slots and waveguide sections, starting from the load end and proceeding towards the source. While the use of non-resonant coupling slots may provide an additional degree of freedom in the design, resonant coupling slots simplify the design process. The amplitude of the wave going to the load is set at unity. The S11 parameter, r of the coupling slot closest to the load, is assigned an arbitrary value. A larger value of r will reduce the power dissipated in the load while increasing the reflection coefficient at the input port. It is now possible to obtain the excitation of the radiating waveguide closest to the load and the coefficients of the wave incident and reflected at the input port of this coupling slot. The next coupling slot parameter, r , is chosen to realize the excitation of that radiating waveguide. One continues this process moving towards the source, until all the coupling slot parameters r and hence the S11 parameter of the 4-port coupler, r, are known for each coupling slot. The goal is to produce the desired array aperture distribution in the feed direction. From an interpolation of the computed moment method data for the slot parameters, all the coupling slot tilt angles and lengths are obtained. From the excitations of the radiating waveguides computed from the coupling values, radiating slot parameters may be obtained so as to attain the desired total normalized slot admittances. This process yields the radiating slot parameters, offsets, and lengths. The design is repeated by choosing different values of r for the last coupling slot until the percentage of power dissipated in the load and the input reflection coefficient values are satisfactory. Numerical results computed for the radiation pattern, the tilt angles and lengths of coupling slots, and excitation phases of the radiating waveguides, are presented for an array with uniform amplitude excitation. The design process has been validated using computer simulations. This design procedure is valid for non-uniform amplitude excitations as well.
The electronic band structure and related band parameters such as, lattice parameters, energy band gaps, density of states and magnetic moment of Zn{sub 1−x}Mn{sub x}O ternary system in the hypothetical hexagonal wurtzite structure have been investigated. The calculations are performed using a full-potential linearized augmented plane wave method based on density functional theory in the generalized gradient approximation of Perdew-Burke-Ernzerhof. The evolution of the features of interest with increasing the Mn doping concentration from 0 to 1 in the ZnO structure has been examined and discussed. Our findings show that the behaviour of all studied band parameters with respect to the alloy composition x is consistent with those reported experimentally in the literature. The present study may be useful for spintronics applications.
An extensive aircraft underflight program was conducted along the Skylab groundpath for the purpose of documenting wind, wave, and atmospheric conditions affecting the amplitude of the active and passive microwave signatures. The S-193 microwave system senses a roughness parameter at the ocean surface that is proportional to the surface windspeed. The aircraft program is described, and some comparisons of satellite and aircraft results are presented. The principal result of the comparison of active radar is that direct inferences of the surface windspeed are possible, but subject to considerable scatter, and that this scatter appears to be due to interaction between long gravity and short Bragg waves and backscatter due to rain as well as errors in correcting for azimuth dependence. An unforeseen opportunity to observe a Pacific hurricane by both Skylab and NOAA aircraft has contributed to the development of a simplified wave forecasting scheme applicable to hurricanes, and more general conditions, which combines the better qualities of both spectral and height/period forecasting techniques. Horizontal polarization data obtained by the aircraft in Hurricane Ava, and in other experiments, which led to this development are presented.
To better understand the relationship between P-wave velocities and ice content in saturated, unconsolidated saline permafrost, we constructed an effective-medium model based upon ultrasonic P-wave data that were obtained from earlier laboratory studies. The model uses a two-end-member mixing approach in which an ice-filled, fully frozen end member and a water-filled, fully unfrozen end member are mixed together to form the effective medium of partially frozen sediments. This mixing approach has two key advantages: (1) It does not require parameter tuning of the mixing ratios, and (2) it inherently assumes mixed pore-scale distributions of ice that consist of frame-strengthening (i.e., cementing and/or loadbearing) ice and pore-filling ice. The model-predicted P-wave velocities agree well with our laboratory data, demonstrating the effectiveness of the model for quantitatively inferring ice content from P-wave velocities. The modeling workflow is simple and is largely free of calibration parameters -attributes that ease its application in interpreting field data sets.
The mean curvature of the influence surface of the space-time point (x, t) appears in linear supersonic propeller noise theory and in the Kirchhoff formula for a supersonic surface. Both these problems are governed by the linear wave equation with sources on a moving surface. The influence surface is also called the Sigma - surface in the aeroacoustic literature. This surface is the locus, in a frame fixed to the quiescent medium, of all the points of a radiating surface f(x, t) = 0 whose acoustic signals arrive simultaneously to an observer at position x and at the time t. Mathematically, the Sigma- surface is produced by the intersection of the characteristic conoid of the space-time point (x, t) and the moving surface. In this paper, we derive the expression for the local mean curvature of the Sigma - space of the space-time point for a moving rigid or deformable surface f(x, t) = 0. This expression is a complicated function of the geometric and kinematic parameters of the surface f(x, t) = 0. Using the results of this paper, the solution of the governing wave equation of high speed propeller noise radiation as well as the Kirchhoff formula for a supersonic surface can be written as very compact analytic expression.
Spectroscopic parameters are accurately determined for the alkali and alkaline-earth monofluorides and monochlorides by means of ab initio self-consistent field and correlated wave function calculations. Numerical Hartree-Fock calculations are performed on selected systems to ensure that the extended Slater basis sets employed are near the Hartree-Fock limit. Since the bonding is predominantly electrostatic in origin, a strong correlation exists between the dissociation energy (to ions) and the spectroscopic parameter r(e). By dissociating to the ionic limits, most of the differential correlation effects can be embedded in the accurate experimental electron affinities and ionization potentials.
We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.
Various systems can be modeled as a point-like explosion of ionized debris into a magnetized, collisionless background plasma-including astrophysical examples, active experiments in space, and laser-driven laboratory experiments. Debris streaming from the explosion parallel to the magnetic field may drive multiple resonant and non-resonant ion-ion beam instabilities, some of which can efficiently couple the debris energy to the background and may even support the formation of shocks. We present a large-scale hybrid (kinetic ions + fluid electrons) particle-in-cell simulation, extending hundreds of ion inertial lengths from a 3D explosion, that resolves these instabilities. We show that the character of these instabilities differs notably from the 1D equivalent by the presence of unique transverse structure. Additional 2D simulations explore how the debris beam length, width, density, and speed affect debris-background coupling, with implications for the generation of quasi-parallel shocks.
The steady-state amplitude of the output of an ultrasonic through-transmission measurement is analyzed and the result is given in closed form. Provided that the product of the input and output transduction ratios, the specimen-transducer reflection coefficient, the specimen-transducer phase-shift parameter, and the material phase velocity are known, this analysis gives a means for determining the through-thickness attenuation of an individual thin sample. Multiple stress-wave reflections are taken into account, and so signal echoes do not represent a difficulty. An example is presented for a graphite fiber epoxy composite (Hercules AS/3501-6). Thus, the technique provides a direct method for continuous or intermittent monitoring of through-thickness attenuation of plate structures which may be subject to service structural degradation.
Recently, a braking index for the pulsar PSR J1640-4631 has been measured. With a braking index of n = 3.15 ± 0.03, this pulsar has the highest braking index ever measured. As it is well known, a pure magnetic dipole brake yields n = 3, whereas a pure gravitational wave (GW) brake yields n = 5. Therefore, each of these mechanisms alone can not account for the braking index found for PSR J1640-4631. Here we consider in detail that such a braking index could be accounted for if the spindown model combines magnetic dipole and GW brakes. Then, we briefly discuss the detectability of this pulsar by aLIGO and the planned Einstein Telescope. In particular, we show that the amplitude of the GW that comes from our model is around a factor four lower than the amplitude modeled exclusively by GW energy loss. Another interesting outcome of our modeling is that it is possible to obtain the ellipticity from the braking index and other pulsar parameters.
Here, we consider various models realizing baryogenesis during the electroweak phase transition (EWBG). Our focus is their possible detection in future collider experiments and possible observation of gravitational waves emitted during the phase transition. We also discuss the possibility of a non-standard cosmological history which can facilitate EWBG. We show how acceptable parameter space can be extended due to such a modification and conclude that next generation precision experiments such as the ILC will be able to confirm or falsify many models realizing EWBG. We also show that, in general, collider searches are a more powerful probe than gravitational wave searches. However, observation of a deviation from the SM without any hints of gravitational waves can point to models with modified cosmological history that generically enable EWBG with weaker phase transition and thus, smaller GW signals.
Solar energetic particle acceleration by a shock wave accompanying a coronal mass ejection (CME) is studied. The description of the accelerated particle spectrum evolution is based on the numerical calculation of the diffusive transport equation with a set of realistic parameters. The relation between the CME and shock speeds, which depend on the initial CME radius, is determined. Depending on the initial CME radius, its speed, and the magnetic energy of the scattering Alfvén waves, the accelerated particle spectrum is established 10–60 minutes from the beginning of CME motion. The maximum energies of particles reach 0.1–10 GeV. The CME radii of 3–5 R{sub ⊙} and the shock radii of 5–10 R{sub ⊙} agree with observations. The calculated particle spectra agree with the observed ones in events registered by ground-based detectors if the turbulence spectrum in the solar corona significantly differs from the Kolmogorov one.
Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.
The (low frequency) zonal flow-ion temperature gradient (ITG) wave system, constructed on Braginskii's fluid model in tokamak, is shown here to be a reaction-diffusion-advection system; it is derived by making use of a multiple spatiotemporal scale technique and two-dimensional (2D) ballooning theory. For real regular group velocities of ITG waves, two distinct temporal processes, sharing a very similar meso-scale radial structure, are identified in the nonlinear self-organized stage. The stationary and quasi-stationary structures reflect a particular feature of the poloidal group velocity. The equation set posed to be an initial value problem is numerically solved for JET low mode parameters; the results are presented in several figures and two movies that show the spatiotemporal evolutions as well as the spectrum analysis—frequency-wave number spectrum, auto power spectrum, and Lissajous diagram. This approach reveals that the zonal flow in tokamak is a local traveling wave. For the quasi-stationary process, the cycle of ITG wave energy is composed of two consecutive phases in distinct spatiotemporal structures: a pair of Cavitons growing and breathing slowly without long range propagation, followed by a sudden decay into many Instantons that carry negative wave energy rapidly into infinity. A spotlight onto the motion of Instantons for a given radial position reproduces a Blob-Hole temporal structure; the occurrence as well as the rapid decay of Caviton into Instantons is triggered by zero-crossing of radial group velocity. A sample of the radial profile of zonal flow contributed from 31 nonlinearly coupled rational surfaces near plasma edge is found to be very similar to that observed in the JET Ohmic phase [J. C. Hillesheim et al., Phys. Rev. Lett. 116, 165002 (2016)]. The theory predicts an interior asymmetric dipole structure associated with the zonal flow that is driven by the gradients of ITG turbulence intensity.
A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfvén wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at β ≳ 1 and by magnetosonic/whistler fluctuations at lower β . The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma β , and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions.