Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local density approximations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Explosions of infalling comets in Jupiter's atmosphere

In view of the expected collision of comet Shoemaker-Levy 9 (1993e) with Jupiter in 1994 July, we calculate basic properties of the initial interaction for a simplified Jovian atmosphere. The comet is expected to impact Jupiter at 60 km/sec and at an angle of 45 deg to the zenith. The shock wave generated by the bolide should be optically thick once it has penetrated to an atmospheric density approximately 10(exp -6) gr/cm(exp 3), and we calculate the post-shock conditions assuming local thermodynamic equilibrium (LTE) for shock velocities v(sub sh) in the range 10 to 60 km/sec and preshock densities rho(sub a) = 10(exp -6) to 10(exp -2) gr/cm(exp 3). Our shock calculations include molecular hydrogen, atomic hydrogen, ionized hydrogen, neutral helium, and singly ionized helium. Even at the highest shock velocity, the gas is only partially ionized and the postshock temperature rises with preshock density in order to maintain the ionization. The value of the effective shock adiabatic index gamma(sub sh) varies from 1.17 (at low v(sub sh) and rho(sub a) to 1.40 (at high v(sub sh) and rho(sub a). The ablation rate is limited by the radiative flux that reaches the bolide surface. We argue that the ablated gas does not efficiently transfer its kinetic energy to the atmosphere, and it ultimately slows in a similar fashion to the comet material. As the bolide initially falls through the atmosphere, the character of the shock emission changes. At rho(sub a) approximately 10(exp -8) gr/cm (exp 3), the gas is optically thin and we expect line emission; in the optical spectrum, Balmer emission is expected from the shocked atmosphere and low-ionization metal lines from ablated cometary material. At rho(sub a) approximately 10(exp -6) gr/m(exp 3), the shocked gas is optically thick and the shock front near the bolide produces a blackbody spectrum. The temperature is favorable for ultraviolet (1000 to 3000 A) emission and the luminosity may be approximately 5 x 10(exp 23) ergs/sec for approximately 0.6 sec for a bolide 1 km in radius. At rho(sub a) approximately 10(exp -4) gr/cm(exp 3), the bolide has passed below the ultraviolet photosphere. The shock front emits considerable ionizing radiation, but it is absorbed in a narrow preshock region. The bolometric correction for the optical luminosity is large and we expect a 3000 to 8000 A luminosity of approximately 3 x 10(exp 23) ergs/sec for approximately 1 sec. The optical emission is strongly peaked in the vicinity of the bolide. The bolide does have a somewhat less luminous, optically thick trail extending greater than or equal to 10 km, but the radiation is characterized by a temperature of 4000 to 5000 K. From the fragmentation model of Chyba, Thomas, & Zahnle (1993), the bolide deposits most of its kinetic energy at rho(sub a) approximately 10(exp -3) gr/cm(exp 3) and this is the effective explosion site. The shock wave from such an explosion can move up about one density scale height. We examine the breakout of the shock front from the Jovian atmosphere and find that the shock acceleration in the decreasing density region is slow, so that the energy flux in the shock front is small. Higher velocities might be generated by shock acceleration along the channel left by the bolide if the shock motion can occur before the channel closes off as a result of radiative cooling. Hot gas created by the explosion ultimately rises due to buoyancy on a timescale of a minute. The luminosity is highest when the bubble first rises into the optically thin part of the atmosphere and may be approximately 1 x 10(exp 25) ergs/sec in the near-infrared. Roughly 1% of the initial bolide energy may be radiated in this way; the rest of the energy is lost to sound waves from the initial explosion and to work done by the bubble on the surrounding atmosphere.

Chevalier, Roger A.↗
NASA NTRS1976

Ultraviolet observations of local gas

From satellite measurements of ultraviolet spectra of stars, an average density of approximately 1.1/cu cm for hydrogen atoms, in both atomic and molecular form, is estimated for regions of space along the galactic plane within about 1 kpc of the sun. About 20% of the atoms are bound in molecular form although this figure is uncertain since the ultraviolet measurements avoid the very dense interstellar clouds. Discrete values for this percentage are observed to vary markedly; regions with less than average density seem to have fractional abundances of H2 several orders of magnitude lower than average. A ratio of CO/H is observed for regions in front of stars observed by the Copernicus satellite.

Jenkins, E. B.↗
NASA NTRS1977

Ultraviolet observations of local gas

From satellite measurements of ultraviolet spectra of stars, an average density of approximately 1.1 cu m/sec for hydrogen atoms, in both atomic and molecular form, is estimated for regions of space along the galactic plane within about 1 kpc of the sun. About 20 percent of the atoms are bound in molecular form, although this figure is uncertain because the ultraviolet measurements avoid the very dense interstellar clouds. Discrete values for this percentage are observed to vary markedly; regions with less than average density appear to have fractional abundances of H2 several orders of magnitude lower than average. A ratio of CO/H which ranges from 10 to the minus 9th to 10 to the minus 6th power is observed for regions in front of stars observed by the Copernicus satellite.

Jenkins, E. B.↗
NASA NTRS1982

Far-infrared mapping of the double-lobed H II region S106

Photometry of S106 at 50 and 100 microns, as well as its mapping with 35 arcsec resolution and a search for a nearby H2O maser source at 2 and 10 microns, have yielded data that are consistent with a 900 pc distance and an O9 (ZAMS) spectral type for a single, central exciting star. The correlation of dust temperature with optical emission, and anticorrelation of dust optical depth with optical emission, support models of S106 which envision an exciting star surrounded by a dusty disk having very low optical extinction along its approximately north-south axis. A local maximum in the dust column density was found to the west of the optical nebula, near the position of an H2O maser, although no local near- or far-IR flux maximum was detected.

Harvey, P. M.↗
NASA NTRS1995

Are Unidentified Extreme-Ultraviolet Sources the Closest Neutron Stars?

Unidentified extreme-ultraviolet (EUV) sources, detected in the EUVE and ROSAT WFC all-sky surveys, could be isolated old neutron stars, accreting material from the interstellar medium (ISM). The closest neutron stars, which are located in the local ISM bubble of unusually low density, are faint and cool (L approximately 10(exp 27) erg/s, T approximately < 6 eV). The extreme-UV spectrum of these sources is very sensitive to the H I column density, since a large fraction of the energy is emitted just below the hydrogen Lyman edge. The EUVE sources with large count rates in the long-wavelength bandpass (600 A) seem to be the most promising candidates. These sources should have low H I column density (N(sub H I) approximately < 10(exp 18) sq cm), constraining their distances to a few tens of parsecs. Otherwise, their spectra would be significantly modified by ISM absorption, and inevitably they would appear stronger in the short-wavelength (100 A, 200 A) EUVE bandpasses. If these unidentified objects are familiar EUV sources rather than neutron stars, i.e. white dwarfs, late-type stars or cataclysmic variables, they are expected to be identifiable, and generally brighter than V approximately 14.

Shemi, A.↗
NASA NTRS1994

The local distribution of Na I interstellar gas

We present high-resolution absorption measurements (lambda/Delta lambda approximately 75,000) of the interstellar Na I D lines at 5890 A toward 80 southern hemisphere early-type stars located in the local interstellar medium (LISM). Combining these results with other sodium measurements taken from the literature, we produce galactic maps of the distribution of neutral sodium column density for a total of 293 stars generally lying within approximately 250 pc of the Sun. These maps reveal the approximate shape of the mid-plane contours of the rarefied region of interstellar space termed the Local Bubble. Its shape is seen as highly asymmetric, with a radius ranging from 30 to 300 pc, and with an average radius of 60 pc. Similar plots of the Galactic mid-plane distribution of sources emitting extreme ultraviolet radiation show that they also trace out similar contours of the Local Bubble derived from Na I absorption measurements. We conclude that the Local Bubble absorption interface can be represented by a hydrogen column density, Nu(sub ETA) = 2 x 10(exp 19) cm(exp -2), which explains both the local distribution of Na I absorption and the observed galactic distribution of extreme ultraviolet sources. The derived mid-plane contours of the Bubble generally reproduce the large-scale features carved out in the interstellar medium by several nearby galactic shell structures.

Welsh, B. Y.↗
NASA NTRS2013

Concurrent Supermassive Black Hole and Galazy Growth: Linking Environment and Nuclear Activity in Zeta Equals 2.23 H Alpha Emitters

We present results from an approximately equal 100 ks Chandra observation of the 2QZ Cluster 1004+00 structure at z = 2.23 (hereafter 2QZ Clus). 2QZ Clus was originally identified as an overdensity of four optically-selected QSOs at z = 2.23 within a 15 × 15 arcmin square region. Narrow-band imaging in the near-IR (within the K band) revealed that the structure contains an additional overdensity of 22 z = 2.23 H alpha-emitting galaxies (HAEs), resulting in 23 unique z = 2.23 HAEs/QSOs (22 within the Chandra field of view). Our Chandra observations reveal that three HAEs in addition to the four QSOs harbor powerfully accreting supermassive black holes (SMBHs), with 2-10 keV luminosities of approximately equal (8-60) × 10(exp 43) erg s(exp−1) and X-ray spectral slopes consistent with unobscured active galactic nucleus (AGN). Using a large comparison sample of 210 z = 2.23 HAEs in the Chandra-COSMOS field (C-COSMOS), we find suggestive evidence that the AGN fraction increases with local HAE galaxy density. The 2QZ Clus HAEs reside in a moderately overdense environment (a factor of approximately equal 2 times over the field), and after excluding optically-selected QSOs, we find that the AGN fraction is a factor of approximately equal 3.5(+3.8/ −2.2) times higher than C-COSMOS HAEs in similar environments. Using stacking analyses of the Chandra data and Herschel SPIRE observations at 250micrometers, we respectively estimate mean SMBH accretion rates ( M(BH)) and star formation rates (SFRs) for the 2QZ Clus and C-COSMOS samples. We find that the mean 2QZ Clus HAE stacked X-ray luminosity is QSO-like (L(2−10 keV) approximately equal [6-10] × 10(exp 43) erg s(exp −1)), and the implied M(BH)/SFR approximately equal (1.6-3.2) × 10(exp −3) is broadly consistent with the local M(BH)/Stellar Mass relation and z approximately equal 2 X-ray selected AGN. In contrast, the C-COSMOS HAEs are on average an order of magnitude less X-ray luminous and have M(BH)/SFR approximately equal (0.2-0.4) × 10(exp −3), somewhat lower than the local MBH/M relation, but comparable to that found for z approximately equal 1-2 star-forming galaxies with similar mean X-ray luminosities. We estimate that a periodic QSO phase with duty cycle approximately 2%-8% would be sufficient to bring star-forming galaxies onto the local M(BH)/Stellar Mass relation. This duty cycle is broadly consistent with the observed C-COSMOS HAE AGN fraction (Approximately equal 0.4%-2.3%) for powerful AGN with LX approximately greater than 10(exp 44) erg s(exp −1). Future observations of 2QZ Clus will be needed to identify key factors responsible for driving the mutual growth of the SMBHs and galaxies.

Lehmer, B. D.↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Mesospheric and lower thermospheric winds, temperatures densities, and volume emission rates

This memo has been written to report on the progress made under grant NAG1-1315 and to request a continuation of funding for this work. Our proposal involved a plan to utilize an existing ground-based chain of optical and radar facilities to assemble a comprehensive, long-term, multi-station base of upper-mesospheric and lower thermospheric measurements of neutral winds, temperatures, and volume emission rates that can be used to make comparisons with data from the HRDI and WINDII instruments that are flying aboard the UARS spacecraft. The ground-based, optical data were to be obtained on a routine basis from five geographically separated observatories at: Thule, Greenland; Sondrestrom, Greenland; Watson Lake, Yukon (replacing Calgary, Alberta); Ann Arbor, Michigan; and Maynooth, Eire. Several different optical instruments are present at these sites: the total array of instruments is comprised of five Fabry-Perot interferometers, two half meter Ebert-Fastie spectrometers, one all-sky CCD imager, and a near infra-red Michelson Fourier transform spectrometer. In addition to these optical instruments, data were to be obtained from the incoherent scatter radar at Sondrestrom, Greenland. These radar measurements are comprised of neutral winds, temperatures, and densities from altitudes between approximately 70 - 120 km. The optical measurements are obtained locally from specific altitudes depending on the emission line studied. For example, red line optical data come from about 220 km. In this report we summarize the progress made in obtaining these data and relate it to the specific tasks outlined in the original grant application. These tasks are discussed in the next section. Progress towards their completion is discussed in section three, while future plans and summary are described in section four.

Killeen, Timothy L.↗
NASA NTRS1994

On the determination of the cosmic infrared background radiation from the high-energy spectrum of extragalactic gamma-ray sources

In a recent paper Stecker, De Jager, & Salamon have suggested using the observed approximately MeV to TeV spectra of extragalactic gamma-ray sources as probes of the local density of the cosmic infrared background radiation (CIBR) and have subsequently claimed a first possible measurement of the CIBR from the analysis of the gamma-ray spectrum of Mrk 421 (De Jager, Stecker, & Salamon). The CIBR from normal galaxies consists of two components: a stellar emission component (CIBRs), and a thermal dust emission component (CIBRd). Photons with energies in the approximately 0.1-2 TeV range interact primarily with the CIBRs, whereas interactions with CIBRd dominate the absorption of photons in the approximately 2-100 TeV energy range. SDS 92 and DSS94 considered only the interaction of the gamma-rays with the dust emission component of the CIBR. We present here an improved analysis of the absorption of extragalactic TeV gamma rays by the CIBR, taking the dual nature of its origin into account. Applying the analysis to the observed gamma-ray spectrum of Mrk 421, a BL Lac object at z = 0.031, we find agreement with DSS94 tentative evidence for absorption by the CINRs. Our analysis therefore limits the detection of the CIBR to the approximately 15-40 micron wavelength regime which, considering the uncertainties in the highest energy (greater than 4 TeV) data and ion the possibility of absorption inside the source, many turn out to be an upper limit on its energy density. At shorter wavelengths (lambda approximately = 1-15 microns), where the gamma-ray interactions are dominated by the CIBRs, our analysis definitely yields only an upper limit on the energy density of the CIBR. In contrast, DSS94 have claimed a possible first measurement of the CIBR over the entire 1-120 micron wavelength region. The upper limit on the CIBRs and tentative detection of the CIBRd are consistent with normal galaxies contributing most of the energy to the CIBR, and constrain the contribution of some exotic sources. With careful modeling of infrared foreground emissions, these constraints on the CIBR are above the values measurable by the DIRBE experiment on board the Cosmic Background Explorer (COBE) satellite.

Dwek, Eli↗
NASA NTRS2024

Variability and Distribution of Nighttime Equatorial to Mid Latitude Ionospheric Irregularities and Vertical Plasma Drift Observed by FORMOSAT-5 Advanced Ionospheric Probe In-Situ Measurements from 2017 – 2020

Irregularities in ionospheric plasma distribution can result in severe scintillation and disruption to the radio frequencies utilized for satellite communications and navigation. In the low and mid latitudes, these irregularities can include Equatorial Plasma Bubbles (EPBs) and Travelling Ionospheric Disturbances (TIDs). EPBs are irregularities manifesting in low latitude nighttime ionosphere plasma density that can extend along magnetic field lines with zonal scales on the order of 100 km or less, while TIDs are propagating wave disturbances. High frequency in-situ measurements of ionospheric plasma aboard spacecraft in Low Earth Orbit (LEO) are a direct measurement of irregularities in plasma density and are therefore valuable for resolving EPB and TID occurrences, variability, and relation to other ionospheric parameters that are believed to play a driving role in the formation of such irregularities. In this study, we utilize observations taken over a three-year period between 2017 and 2020 by the Advanced Ionospheric Probe (AIP) carried aboard the FORMOSAT-5 satellite to examine the spatial, seasonal, and interannual variability of equatorial to mid latitude ionospheric irregularities and vertical ion drift during this time. AIP provides in-situ measurements of ion density and vertical ion drift in the equatorial to mid latitude ionosphere at approximately 720 km altitude with local times between 22:00 – 23:00 local time. Our global scale results resolve distinct and inter-annually recurrent seasonal patterns in the distribution of nighttime ionospheric irregularities and vertical plasma drift during this time. Elevated occurrences of ion density irregularities are resolved along the Equatorial Ionization Anomaly (EIA) latitudes, while notable occurrences with variability consistent with EPBs also observed along the low and equatorial magnetic latitudes. Zonal variability of equatorial irregularities consistent with the signatures of nonmigrating atmospheric tides are observed. It is also notable that the occurrences and geographic distribution of ion density irregularities showed a considerable level of interannual variability, especially at mid latitudes over the South Atlantic and Southern African sectors, which showed much higher levels of irregularities in 2017–––2018, compared to 2019 and 2020. In comparison, the spatial and interannual variation of the co-located vertical ion drifts were much more consistent during the years examined, indicating that the driver for the observed interannual variability in ion density irregularities cannot be attributed to the vertical ion drift at the same time and location of the observations. This highlights the need for in-situ instruments distributed across multiple satellites in different local time zones.

ionosphere↗
DOE OSTI2024

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1977

Biaxial load effects on the crack border elastic strain energy and strain energy rate

The validity of the singular solution (first term of a series representation) is investigated for the crack tip stress and displacement field in an infinite sheet with a flat line crack with biaxial loads applied to the outer boundaries. It is shown that if one retains the second contribution to the series approximations for stress and displacement in the calculation of the local elastic strain energy density and elastic strain energy rate in the crack border region, both these quantities have significant biaxial load dependency. The value of the J-integral does not depend on the presence of the second term of the series expansion for stress and displacement. Thus J(I) is insensitive to the presence of loads applied parallel to the plane of the crack.

Eftis, J.↗
NASA NTRS1989

Probability density and exceedance rate functions of locally Gaussian turbulence

A locally Gaussian model of turbulence velocities is postulated which consists of the superposition of a slowly varying strictly Gaussian component representing slow temporal changes in the mean wind speed and a more rapidly varying locally Gaussian turbulence component possessing a temporally fluctuating local variance. Series expansions of the probability density and exceedance rate functions of the turbulence velocity model, based on Taylor's series, are derived. Comparisons of the resulting two-term approximations with measured probability density and exceedance rate functions of atmospheric turbulence velocity records show encouraging agreement, thereby confirming the consistency of the measured records with the locally Gaussian model. Explicit formulas are derived for computing all required expansion coefficients from measured turbulence records.

Mark, W. D.↗
NASA NTRS1981

Results of a new approach to determining the density function in the galactic plane - The local system

Star counts in the directions of 15 opaque CO clouds of known distance are used to determine the density function in the galactic plane of the local system of apparently enhanced stellar density. Counts were made on Palomar Sky Survey blue prints in the directions of the highest opacity portions of 15 molecular clouds. Average values of the density functions obtained are found to be significantly less than unity for clouds more distant than 800 pc regardless of galactic longitude, confirming the presence of the local mass density enhancement. Data suggest that the local system represents a density enhancement of a factor of 2 over a region of approximately 500 pc. An empirical relationship is also determined between the distance to a molecular cloud and the number of foreground stars projected on it, which may be useful in determining distances to clouds nearer than 3 kpc. The potential of the Space Telescope for complete stellar mass determinations out to a distance of 1 kpc is pointed out.

Herbst, W.↗
NASA NTRS1989

Problems in determining the surface density of the Galactic disk

A new method is presented for determining the local surface density of the Galactic disk from distance and velocity measurements of stars toward the Galactic poles. The procedure is fully three-dimensional, approximating the Galactic potential by a potential of Staeckel form and using the analytic third integral to treat the tilt and the change of shape of the velocity ellipsoid consistently. Applying the procedure to artificial data superficially resembling the K dwarf sample of Kuijken and Gilmore (1988, 1989), it is shown that the current best estimates of local disk surface density are uncertain by at least 30 percent. Of this, about 25 percent is due to the size of the velocity sample, about 15 percent comes from uncertainties in the rotation curve and the solar galactocentric distance, and about 10 percent from ignorance of the shape of the velocity distribution above z = 1 kpc, the errors adding in quadrature. Increasing the sample size by a factor of 3 will reduce the error to 20 percent. To achieve 10 percent accuracy, observations will be needed along other lines of sight to constrain the shape of the velocity ellipsoid.

Statler, Thomas S.↗