Short-term frequency stability - characterization, theory, and measurement.
Short-term frequency stability of sinusoidal oscillation sources and loop-controlled oscillations discussing characterization, theory and measurement
SEARCH · Engineering Papers
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.
Short-term frequency stability of sinusoidal oscillation sources and loop-controlled oscillations discussing characterization, theory and measurement
Seismic refraction tomography provides images of the elastic properties of subsurface materials in landslide settings. Seismic velocities are sensitive to changes in moisture content, which is a triggering factor in the initiation of many landslides. However, the application of the method to long-term monitoring of landslides is rarely used, given the challenges in undertaking repeat surveys and in handling and minimizing the errors arising from processing time-lapse surveys. This work presents a simple method and workflow for producing a reliable time-series of inverted seismic velocity models. This method is tested using data acquired during a recent, novel, long-term seismic refraction monitoring campaign at an active landslide in the UK. Potential sources of error include those arising from inaccurate and inconsistent determination of first-arrival times, inaccurate receiver positioning, and selection of inappropriate inversion starting models. At our site, a comparative analysis of variations in seismic velocity to real-world variations in topography over time shows that topographic error alone can account for changes in seismic velocity of greater than ±10% in a significant proportion (23%) of the data acquired. The seismic velocity variations arising from real material property changes at the near-surface of the landslide, linked to other sources of environmental data, are demonstrated to be of a similar magnitude. Over the monitoring period we observe subtle variations in the bulk seismic velocity of the sliding layer that are demonstrably related to variations in moisture content. This highlights the need to incorporate accurate topographic information for each time-step in the monitoring time-series. The goal of the proposed workflow is to minimize the sources of potential errors, and to preserve the changes observed by real variations in the subsurface. Following the workflow produces spatially comparable, time-lapse velocity cross-sections formulated from disparate, discretely-acquired datasets. These practical steps aim to aid the use of the seismic refraction tomography method for the long-term monitoring of landslides prone to hydrological destabilization.
An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.
Linear wave equation models which have been used over the years at NASA Langley for describing noise emissions from high speed rotating blades are summarized. The noise sources are assumed to lie on a moving surface, and analysis of the situation has been based on the Ffowcs Williams-Hawkings (FW-H) equation. Although the equation accounts for two surface and one volume source, the NASA analyses have considered only the surface terms. Several variations on the FW-H model are delineated for various types of applications, noting the computational benefits of removing the frequency dependence of the calculations. Formulations are also provided for compact and noncompact sources, and features of Long's subsonic integral equation and Farassat's high speed integral equation are discussed. The selection of subsonic or high speed models is dependent on the Mach number of the blade surface where the source is located.
A new analytic result in acoustics called 'Formulation 1B,' proposed by Farassat, is used to compute broadband trailing edge noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term, and has been shown in previous research to provide time domain predictions of broadband noise that are in excellent agreement with experiment. Furthermore, this formulation lends itself readily to rotating reference frames and statistical analysis of broadband trailing edge noise. Formulation 1B is used to calculate the far field noise radiated from the trailing edge of a NACA 0012 airfoil in low Mach number flows, using both analytical and experimental data on the airfoil surface. The results are compared to analytical results and experimental measurements that are available in the literature. Good agreement between predictions and measurements is obtained.
Predictions of community noise exposure from a NASA urban air mobility (UAM) concept vehicle have been conducted for representative operations using the FAA Aviation Environmental Design Tool (AEDT) in order to demonstrate modeling tool interoperability and assess applicability, capabilities and limitations of integrated noise modeling tools. To both quantify limitations and highlight other capabilities, a comparative analysis is performed using a time-domain simulation method, in particular, using the Volpe Advanced Acoustic Model (AAM). Using the same source noise model, the 3D directivity of a UAM concept vehicle is predicted in terms of aeroacoustic pressure time histories at a sphere of observer positions near the vehicle. In addition to distilling those data to a set of noise-power-distance data for input to AEDT, the data are processed preserving spectral directivity, into one-third and one-twelfth octave bands for input to AAM. Results from AEDT and AAM modeling are provided for a variety of metrics to demonstrate the effect that source noise directivity and propagation modeling fidelity have on predicted results at receptors over a study area.
Predictions of community noise exposure from a NASA urban air mobility (UAM) concept vehicle have been conducted for representative operations using the FAA Aviation Environmental Design Tool (AEDT) in order to demonstrate modeling tool interoperability and assess applicability, capabilities and limitations of integrated noise modeling tools. To both quantify limitations and highlight other capabilities, a comparative analysis is performed using a time-domain simulation method, in particular, using the Volpe Advanced Acoustic Model (AAM). Using the same source noise model, the 3D directivity of a UAM concept vehicle is predicted in terms of aeroacoustic pressure time histories at a sphere of observer positions near the vehicle. In addition to distilling those data to a set of noise-power-distance data for input to AEDT, the data are processed preserving spectral directivity, into one-third and one-twelfth octave bands for input to AAM. Results from AEDT and AAM modeling are provided for a variety of metrics to demonstrate the effect that source noise directivity and propagation modeling fidelity have on predicted results at receptors over a study area.
The linearized unsteady aerodynamic response of a cascade of airfoils subjected to entropic, vortical, and acoustic gusts is analyzed. Field equations for the first-order unsteady perturbation flow are obtained by linearizing the full time-dependent mass, momentum, and energy conservation equations about a nonlinear, isentropic, and irrotational mean or steady flow. A splitting technique is then used to decompose the unsteady velocity field into irrotational and rotational parts leading to field equations for the unsteady entropy, rotational velocity, and irrotational velocity fluctuations that are coupled only sequentially. The entropic and rotational velocity fluctuations can be described in terms of the mean-flow drift and stream functions which can be computed numerically. The irrotational unsteady velocity is described by an inhomogeneous linearized potential equation which contains a source term that depends on the rotational velocity field. This equation is solved via a finite difference technique. Results are presented to indicate the status of the numerical solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described leads to very efficient predictions of cascade unsteady aerodynamic phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.
The linearized unsteady aerodynamic response of a cascade of airfoils subjected to entropic, vortical, and acoustic gusts is analyzed. Field equations for the first-order unsteady perturbation flow are obtained by linearizing the full time-dependent mass, momentum, and energy conservation equations about a nonlinear, isentropic, and irrotational mean or steady flow. A splitting technique is then used to decompose the unsteady velocity field into irrotational and rotational parts leading to field equations for the unsteady entropy, rotational velocity, and irrotational velocity fluctuations that are coupled only sequentially. The entropic and rotational velocity fluctuations can be described in terms of the mean-flow drift and stream functions which can be computed numerically. The irrotational unsteady velocity is described by an inhomogeneous linearized potential equation which contains a source term that depends on the rotational velocity field. This equation is solved via a finite difference technique. Results are presented to indicate the status of the numerical solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described leads to very efficient predictions of cascade unsteady aerodynamics phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.
Many efforts have been made to illuminate the nature of past hydroclimates in semi-arid and arid regions, where current and future shifts in water availability have enormous consequences on human subsistence. Deep desert aquifers, where groundwater is stored for prolonged periods, might serve as a direct record of major paleo-recharge events. To date, groundwater-based paleoclimate reconstructions have mainly focused on a relatively narrow timescale (up to ∼40 kyr), limited by the relatively short half-life of the widely used radiocarbon (5.73 kyr). Here we demonstrate the usage of deep regional aquifers in the arid southeastern Mediterranean as a hydroclimate archive for earlier Mid-to-Late Pleistocene epochs. State-of-the-art dating tools, primarily the 81 Kr radioisotope (t 1/2 = 229 kyr), were combined with other atmosphere-derived tracers to illuminate the impact of four distinguishable wetter episodes over the past 400 kyr, with differences in climatic conditions and paleo-recharge locations. Variations in stable water isotope composition suggest moisture transport from more proximal (Mediterranean) and distal (Atlantic) sources to different parts of the region at distinct times. Large variability in the computed noble gas-based recharge temperature (NGT), ranging ~15–30 °C, cannot be explained by climate variations solely, and points to different recharge pathways, including geothermal heating in the deep unsaturated zone and recharge from high-elevation (colder) regions. The obtained groundwater record complements and enhances the interpretation of other terrestrial archives in the arid region, including a contribution of valuable information regarding the moisture source origin as reflected in the deuterium-excess values, which is unattainable from the common practice analysis of calcitic cave deposits. We conclude that similar applications in other deep (hundred-m-order) regional groundwater systems (e.g., the Sahara desert aquifers) can significantly advance our understanding of long-term (up to 1 Myr) paleo-hydroclimate in arid regions, including places where no terrestrial remnants, such as cave, lake, and spring sediments, are available.
Data Analytics applications have made successful strides in the business world where co-analyzing extremely large sets of independent variables have proven profitable. Today, most data analytics tools and techniques, sometimes applicable to Earth science, have targeted the business industry. In fact, the literature is nearly absent of discussion about Earth science data analytics. Earth science data analytics (ESDA) is the process of examining large amounts of data from a variety of sources to uncover hidden patterns, unknown correlations, and other useful information. ESDA is most often applied to data preparation, data reduction, and data analysis. Co-analysis of increasing number and volume of Earth science data has become more prevalent ushered by the plethora of Earth science data sources generated by US programs, international programs, field experiments, ground stations, and citizen scientists.Through work associated with the Earth Science Information Partners (ESIP) Federation, ESDA types have been defined in terms of data analytics end goals. Goals of which are very different than those in business, requiring different tools and techniques. A sampling of use cases have been collected and analyzed in terms of data analytics end goal types, volume, specialized processing, and other attributes. The goal of collecting these use cases is to be able to better understand and specify requirements for data analytics tools and techniques yet to be implemented. This presentation will describe the attributes and preliminary findings of ESDA use cases, as well as provide early analysis of data analytics toolstechniques requirements that would support specific ESDA type goals. Representative existing data analytics toolstechniques relevant to ESDA will also be addressed.
Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.
A formulation based on a generalized reciprocity theorem is developed for analyzing the external high frequency EM scattering by a complex obstacle inside a relatively arbitrary open-ended waveguide cavity when it is illuminated by an external source. This formulation is also extended to include EM fields whose time dependence may be nonperiodic. A significant advantage of this formulation is that it allows one to break up the analysis into two independent parts; one deals with the waveguide cavity shape alone and the other with the obstacle alone. The external scattered field produced by the obstacle (in the presence of the waveguide cavity structure) is given in terms of a generalized reciprocity integral over a surface S(T) corresponding to the interior waveguide cavity cross section located conveniently but sufficiently close to the obstacle. Furthermore, the fields coupled into the cavity from the source in the exterior region generally need to propagate only one-way via the open front end (which is directly illuminated) to the interior surface S(T) in this approach, and not back, in order to find the external field scattered by the obstacle.
We present a method to determine the leading-order (LO) contact term contributing to the nn → ppe - e - amplitude through the exchange of light Majorana neutrinos. Our approach is based on the representation of the amplitude as the momentum integral of a known kernel (proportional to the neutrino propagator) times the generalized forward Compton scattering amplitude n ( p 1 ) n ( p 2 ) W + ( k ) → \( p\left({p}_1^{\prime}\right)p\left({p}_2^{\prime}\right){W}^{-}(k) \) , in analogy to the Cottingham formula for the electromagnetic contribution to hadron masses. We construct model-independent representations of the integrand in the low- and high-momentum regions, through chiral EFT and the operator product expansion, respectively. We then construct a model for the full amplitude by interpolating between these two regions, using appropriate nucleon factors for the weak currents and information on nucleon-nucleon ( NN ) scattering in the 1 S 0 channel away from threshold. By matching the amplitude obtained in this way to the LO chiral EFT amplitude we obtain the relevant LO contact term and discuss various sources of uncertainty. We validate the approach by computing the analog I = 2 NN contact term and by reproducing, within uncertainties, the charge-independence-breaking contribution to the 1 S 0 NN scattering lengths. While our analysis is performed in the \( \overline{\mathrm{MS}} \) scheme, we express our final result in terms of the scheme-independent renormalized amplitude \( {\mathcal{A}}_{\nu}\left(\left|\mathbf{p}\right|,\left|\mathbf{p}^{\prime}\right|\right) \) at a set of kinematic points near threshold. We illustrate for two cutoff schemes how, using our synthetic data for \( {\mathcal{A}}_{\nu } \) , one can determine the contact-term contribution in any regularization scheme, in particular the ones employed in nuclear-structure calculations for isotopes of experimental interest.
A framework is developed for the accurate calculation of Compton and synchrotron processes in a special class of spherically-symmetric sources. The models considered take a distribution in magnetic field values and high-energy electron density that varies smoothly from central values to values at the outer boundary of the source. Further, the assumption is made that the magnetic field has local disorder (isotropy) and the high-energy electron distribution has local isotropy. Then, in terms of parameters of the source, the synchrotron brightness distribution, total synchrotron flux, Compton-synchrotron flux, and synchrotron self-absorption turnover are all computed accurately. Results are given with a view toward application to the analysis of compact non-thermal sources.
The power output from given sources is usually ascertained via an energy flux integral over the normal directions to a remote (farfield) surface; an alternative procedure, which utilizes an integral that specifies the direct rate of working by the source on the resultant field, is described and illustrated for both point and continuous source distributions. A comparison between the respective procedures is made in the analysis of sound radiated from a periodic dipole source whose axis rotates in a plane, on a full or partial angular range, with prescribed frequency. Thus, adopting a conventional approach, Sretenskii (1956) characterizes the rotating dipole in terms of an infinite number of stationary ones along a pair of orthogonal directions in the plane and, through the farfield representation of the latter, arrives at a series development for the instantaneous radiated power, whereas the local manner of power calculation dispenses with the equivalent infinite aggregate of sources and yields a compact analytical result.
A high accuracy convection scheme using a sequential solution technique has been developed and applied to simulate the longitudinal combustion instability and its active control. The scheme has been devised in the spirit of the Total Variation Diminishing (TVD) concept with special source term treatment. Due to the substantial heat release effect, a clear delineation of the key elements employed by the scheme, i.e., the adjustable damping factor and the source term treatment has been made. By comparing with the first-order upwind scheme previously utilized, the present results exhibit less damping and are free from spurious oscillations, offering improved quantitative accuracy while confirming the spectral analysis reported earlier. A simple feedback type of active control has been found to be capable of enhancing or attenuating the magnitude of the combustion instability.
The goal of this one-year investigation was to determine whether N2O is formed in atmospherically significant quantities by the reaction of vibrationally excited levels of the O2((A3 Sigma(sub u)(sup +)) state with nitrogen. O2(A3 Sigma(sub u)(sup +)) is made throughout the upper stratosphere in considerable amounts by solar photoabsorption, and only a very small reactive yield is necessary for this mechanism to be a major N2O source. By long-term 245-252 nm irradiation of O2/N2 mixtures on- and off-resonance with absorption lines in the O2(A3 Sigma(sub u)(sup +) - X3 Sigma(sub g)(sup -)) transition, followed by N2O analysis by frequency-modulated diode laser absorption spectroscopy, we determined an upper limit for the N2O yield of the candidate reaction. This limit, 3 x 10(exp -5), eliminates O2(A3 Sigma(sub u)(sup +)) + N2 as a significant channel for the generation of stratospheric N2O. In further measurements, we established that N2O is stable under our photolysis conditions, showing that the small amounts of ozone generated from the reaction of O2(A) and O2 do not indirectly lead to destruction of N2O.