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.
Two apparently disparate lines of inquiry in kinetic theory are shown to be equivalent: (1) Brownian motion as treated by the (stochastic) Langevin equation and Fokker-Planck equation; and (2) Boltzmann's equation. The method is to derive the kinetic equation for Brownian motion from the Boltzmann equation for a two-component neutral gas by a simultaneous expansion in the density and mass ratios.
Solution of the kinetic dispersion equation shows that heavy pickup ions with large perpendicular energies excite a new hydromagnetic wave mode with maximum growth, high (magnetic field and mass density) compression ratios, and almost linear polarization coexisting slightly away 6-15 deg from parallel propagation. The mode arises from the coupling of the heavy ion beam to the modified left-hand wave dispersion caused by the multiion nature of the medium. The ion beam can also bring about nonoscillatory, purely growing structures distinct from the mirror wave through the coupling of a counter-streaming left-hand mode to a costreaming right-hand wave. The free energy source parameters are varied to determine their relative influence and define the domain of existence of the instability. Its characteristics may prove helpful in the interpretation of recent cometary (Halley and Giacobini-Zinner) and artificial release (AMPTE) observations.
A numerical definition of the thermal evolution of spherically symmetric models of the nucleus in the orbit of Comet P/Tempel-1 is presently used to ascertain the properties of the outer layers of comets under consideration for the future Comet Rendezvous and Asteroid Flyby and the Comet Nucleus Sample Return missions. Evolutionary sequences are computed for different values of density, dust/ice mass ratio, and the dust fraction not lost with ice sublimation. It is found that inner and outer surface temperatures of the dust mantle are comparatively insensitive to parameter changes, and that the total thickness of the crystalline ice layer between the dust mantle and the amorphous ice core will make it difficult for the comet-mission probes to sample pristine ice.
The distribution of Si II emission in the Galactic center has been measured and found to extend beyond 4 pc from the center. The observed forbidden Si II 34.8 micron line, which is thermally excited by H I and H2 collisions, originates from the photodissociation region formed in the neutral gas ring surrounding the ionized core. Two peaks in the Si II emission are found, one lying 20 arcsec northeast of Sgr A-asterisk along the Galactic plane and the other located about 80 arcsec north of Sgr A-asterisk. The latter peak is located beyond the northern arm seen in the ionized gas. The ratio of Si II to dust mass density varies in the ring, indicating changes in the gas-to-dust ratio. The mass density of Si II and dust are anticorrelated, indicating that Si II tracks the destruction of dust.
Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.
We introduce a novel nonparametric method to find solar neighborhood analogs (SNAs) in extragalactic integral field spectroscopic surveys. The main ansatz is that the physical properties of the solar neighborhood (SN) should be encoded in its optical stellar spectrum. We assume that our best estimate of such a spectrum is the one extracted from the analysis performed by the Code for Stellar properties Heuristic Assignment (CoSHA) from the MaStar stellar library. It follows that finding SNAs in other galaxies consist in matching, in a χ 2 sense, the SN reference spectrum across the optical extent of the observed galaxies. We applied this procedure to a selection of CALIFA galaxies, by requiring a close to face-on projection, relative isolation, and non-active galactic nucleus. We explore how the local and global properties of the SNAs (stellar age, metallicity, dust extinction, mass-to-light ratio, stellar surface mass density, star-formation density, and galactocentric distance) and their corresponding host galaxies (morphological type, total stellar mass, star-formation rate, and effective radius) compared with those of the SN and the Milky Way (MW). We find that SNAs are located preferentially in S(B)a–S(B)c galaxies, in a ring-like structure, which radii seem to scale with the galaxy size. Despite the known sources of systematics and errors, most properties present a considerable agreement with the literature on the SN. We conclude that the solar neighborhood is relatively common in our sample of SNAs. Our results warrant a systematic exploration of correlations among the physical properties of the SNAs and their host galaxies. We reckon that our method should inform current models of the galactic habitable zone in our MW and other galaxies.
ABSTRACT Dark matter as scalar particles consisting of multiple species is well motivated in string theory where axion fields are ubiquitous. A two-field fuzzy dark matter (FDM) model features two species of ultralight axion particles with different masses, m1 ≠ m2, which is extended from the standard one-field model with $m_a \sim 10^{-22} \, {\rm eV}$. Here we perform numerical simulations to explore the properties of two-field FDM haloes. We find that the central soliton has a nested structure when m2 ≫ m1, which is distinguishable from the generic flat-core soliton in one-field haloes. However, the formation of this nested soliton is subject to many factors, including the density fraction and mass ratio of the two fields. Finally, we study non-linear structure formation in two-field cosmological simulations with self-consistent initial conditions and find that the small-scale structure in two-field cosmology is also distinct from the one-field model in terms of DM halo counts and soliton formation time.
The extended method of equivalent inclusion developed is applied to study the specific wave problems of the transmission of elastic waves in an infinite medium containing a layer of inhomogeneity, and of the scattering of elastic waves in an infinite medium containing a perfect spherical inhomogeneity. The eigenstrains are expanded as a geometric series and the method of integration for the inhomogeneous Helmholtz operator given by Fu and Mura is adopted. The results obtained by using a limited number of terms in the eigenstrain expansion are compared with exact solutions for the layer problem and for a perfect sphere. Two parameters are singled out for this comparison: the ratio of elastic moduli, and the ratio of the mass densities. General trends for three different situations are shown.
A production 3-D elliptic flow code has been used to calculate non-reacting and reacting flow fields in an experimental mixing section relevant to a rich burn/quick mix/lean burn (RQL) combustion system. A number of test cases have been run to assess the effects of the variation in the number of orifices, mass flow ratio, and rich-zone equivalence ratio on the flow field and mixing rates. The calculated normalized temperature profiles for the non-reacting flow field agree qualitatively well with the normalized conserved variable isopleths for the reacting flow field indicating that non-reacting mixing experiments are appropriate for screening and ranking potential rapid mixing concepts. For a given set of jet momentum-flux ratio, mass flow ratio, and density ratio (J, MR, and DR), the reacting flow calculations show a reduced level of mixing compared to the non-reacting cases. In addition, the rich-zone equivalence ratio has noticeable effect on the mixing flow characteristics for reacting flows.
Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A three-dimensional computational fluid dynamics (CFD) code has been used to predict the mixing flow field characteristics and NOx emission in a quench section of a rich-burn/quick-mix/lean-burn (RQL) combustor. Sixty configurations have been analyzed in both circular and annular geometries in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying the number of orifices per row and orifice shape. Other parameters such as J (momentum-flux ratio), MR (mass flowrate ratio), DR (density ratio), and mixer sector orifice ACd (effective orifice area) were maintained constant throughout the entire study. The results indicate that the mixing flow field can be correlated with the NOx production if they are referenced with the stoichiometric equivalence ratio value and not the equilibrium value. The mixing flowfields in both circular and annular mixers are different. The penetration of equal jets in both annular and circular geometries is vastly different which significantly affects the performance of the mixing section. In the computational results with the circular mixer, most of the NOx formation occurred behind the orifice starting at the orifice wake region. General trends have been observed in the NOx production as the number of orifices is changed and this appears to be common for all hole configurations and mixer types (circular or annular). The performance of any orifice shape (in producing minimum NOx) appears to be acceptable if the number of orifices can be freely varied in order to attain the optimum jet penetration.
ABSTRACT We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new multi-Gaussian expansion photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a = 1.294 ± 0.039, b = 0.912 ± 0.025, and zero-point c = 7.067 ± 0.078. We leverage the high signal-to-noise ratio of SAMI integral field spectroscopy, to determine how structural and stellar population observables affect the scatter about the FP. The FP residuals correlate most strongly (8σ significance) with luminosity-weighted simple stellar population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, Sérsic index, and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio Υ⋆ ) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti)correlate with the residuals of the relation between surface density and Υ⋆ . This correlation implies that part of the FP scatter is due to the broad age and Υ⋆ distribution at any given surface mass density. Using virial mass and Υ⋆, we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75 per cent) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane.
The comparison of perigee density data near 230 km for satellite Aeros-A has been extended to the complete mission time. The average ratio between orbital drag derived density and mass spectrometer measurements is very near to 1 with a large scatter reflected by a standard deviation of the order of 20%. A method of comparison was developed and tested which uses the observed rate of change of orbital period in comparison with the orbital decay computed from the actual mass spectrometer data measured at discrete positions along the orbit. This method proves that deviations from the average ratio of 1 between perigee densities from drag and mass spectrometers are due to the smoothing and poor resolution of the orbital drag technique.
Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.
Prediction of uncontrolled satellite reentry time and trajectory are essential to avoid damages to human being and properties over land. Reentry time and trajectory are largely controlled by aerodynamic drag and mass neutral density below 200 km. In this paper, using a numerical model of the whole atmosphere, WACCM‐X, and a Precision Orbital Determination software, we demonstrate the major sudden stratospheric warming (SSW) in January 2009, during the deep solar minimum of 2009, may cause 0.2–0.3 day of variation in reentry time. The SSW can alter the mass density and composition mixing ratio in the lower thermosphere below 200 km through changed wave forcing and general circulation in the lower thermosphere. This is comparable to the effect of a 2‐SFU (Solar Flux Unit) solar radiance variation on the neutral density. 0.1 day of reentry time uncertainty corresponds to a global shift in impact location, that is, thousands of kilometers distance. For comparison, a major SSW during solar maximum contributes little to reentry time variation, which is predominantly controlled by the solar and geomagnetic forcing. We suggest a few proxies such as O/N 2 in the lower thermosphere or mean meridional wind to account for the lower atmosphere impacts in an operational model. This needs to be investigated in rigorous future work.
Photoelectric U, B, and V, photovoltaic K, and four-filter CCD surface photometry of the small, M81-linked Irr II galaxy NGC 3077 are reported. The observations are compared with data from previous studies, and the mass, density, and optical depth ratios are determined for a prominent dust cloud near the center of the galaxy. These characteristics are found to be similar to those of galactic dust clouds associated with star-formation sites. The observed intensity distribution at various wavelengths indicates that NGC 3077 was an elliptical galaxy before its interaction with M81. During the interaction NGC 3077 acquired the interstellar material which is currently forming stars in its central regions.
We present a study of the high-energy diffuse emission observed toward Orion by the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray Observatory. The total exposure by EGRET in this region has increased by more than a factor of two since a previous study. A simple model for the diffuse emission adequately fits the data; no significant point sources are detected in the region studied (1 = 195 deg to 220 deg and b = -25 deg to -10 deg) in either the composite dataset or in two separate groups of EGRET viewing periods considered. The gamma-ray emissivity in Orion is found to be (1.65 +/- 0.11) x 10(exp -26)/s.sr for E > 100 MeV, and the differential emissivity is well-described as a combination of contributions from cosmic-ray electrons and protons with approximately the local density. The molecular mass calibrating ratio is N(H2)/W(sub CO) = (1.35 +/- 0.15) x 10(exp 20)/sq cm.(K.km/s).
We use the NASA/Ames grid Mars general circulation model with an updated radiation code and cloud nlicrophysics scheme. To speed up the simulations, we run the model at fairly coarse resolution (7.5" latitude x 22.5" longitude). Radiation Code Fluxes and heating rates are calculated from a radiation code based on the two-stream solution to radiative transfer that fully accounts for multiple scattering in the presence of gaseous absorp- tion. The model has 12 spectral intervals. Dust and water ice scattering properties are included. For dust, we use the Ockert-Bell [9] values in the visible, and Forget [IO] values in the infrared. For ice, we can either compute them online as the cloud evolves, or we can specify them. Gaseous opacities for water vapor and CO2 are calculated from correlated k-distributions taken from full line-by-line models. Cloud Microphysics Our cloud scheme is based on momentlorder scheme in which the mass mixing ratio and number density of the cloud ensemble are the advected species. From these we obtain a mean particle size and an estimate or the particle size distribution assuming a variance) which we then divide into 8 bins. Cloud microphysics is performed in each of these bins and includes nucleation, condensation, and gravitational settling. Dust is treated as a tracer and serves as condensation nuclei. The altered size distribution is then converted back into a mean size, a mixing ratio, and a particle number density.
Bi-component protection-systems need to withstand a variety of extreme loading conditions such as high-velocity projectile impact to shock and blast mitigation. The large number of parameters at each stage of bi-component protection system design process results in a near-infinite number of possible variables for optimization. The aim of this project is to identify the scaling laws linking material and geometric scales of the problem to the final impact performance, further including the effects of wave propagation/interaction and inter-component bonding. Specifically, the impact dynamics for different optimal projectile/target density ratios and relative component mass fractions are quantified for bonded ceramic/titanium alloy systems. The impact behavior such as wave mechanics and damage propagation within the target system is characterized and quantified via digital image correlation or via x-ray phase-contrast imaging.