Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

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 559 records · Page 31

Strong X-Ray Absorption in a Broad Absorption Line Quasar: PHL 5200

We present ASCA observations of the z=1.98 prototype broad absorption line quasar (BALQSO): PHL 5200. The source detected in both SIS and GIS. A power-law spectrum (alpha(sub epsilon)=0.6(sup 0.9)(sub -0.6)) with large intrinsic absorption (N(sub H)=1.3(sup 2.3)(sub -1.1)x10(exp 23)cm(exp -2)) best describes the spectrum. Excess column density over the local Galactic value is required at the 99% confidence level. This detection suggests that although BALQSOs are X-ray-quiet, it is strong absorption in the BAL region that makes them appear faint to low-energy X-ray experiments. The required intrinsic absorbing column density is 2-3 orders of magnitude larger than earlier estimates of column densities in BALQSOs. This implies that the BAL systems are much more highly ionized than was previously thought.

Mathur, Smita↗

Plasmapause diffusion

The Bohm diffusion coefficient and observed electrostatic wave scattering are used as the bases of estimates of the smoothing effect that diffusion may have on steep plasmapause density gradients. The estimate for diffusion resulting from scattering by observed electrostatic waves is found to be much lower than that of the perpendicular Bohm diffusion coefficient for characteristic plasma temperatures and magnetic fields. This diffusion rate estimate may be too small, however, if the wave amplitudes are significantly higher for steep plasmapauses. The effects are therefore negligible for most considerations of macroscopic plasmapause dynamics, but may be significant in limiting drift wave instabilities and similar phenomena driven by the steepness of the plasmapause density gradient.

Horwitz, J. L.↗

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced forLarge Synoptic Survey Telescope The Rubin Observatory Legacy Survey of Space and Time (lsst) Dark Energy Science Collaboration (desc). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. Additionally, we report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

79 ASTRONOMY AND ASTROPHYSICS↗

The distribution of particulate material on Mars

The surface materials on Mars were extensively studied using a variety of spacecraft and Earth-based remote sensing observations. These measurements include: (1) diurnal thermal measurements, used to determine average particle size, rock abundance, and the presence of crusts; (2) radar observations, used to estimate the surface slope distributions, wavelength scale roughness, and density; (3) radio emission observations, used to estimate subsurface density; (4) broadband albedo measurements, used to study the time variation of surface brightness and dust deposition and removal; and (5) color observations, used to infer composition, mixing, and the presence of crusts. Remote sensing observations generally require some degree of modeling to interpret, making them more difficult to interpret than direct observations from the surface. They do, however, provide a means for examining the surface properties over the entire planet and a means of sampling varying depths within the regolith. Albedo and color observations only indicate the properties of the upper-most few microns, but are very sensitive to thin, sometimes emphemeral dust coatings. Thermal observations sample the upper skin depth, generally 2 to 10 cm. Rock abundance measurements give an indirect indication of surface mantling, where the absence of rocks suggests mantles of several meters. Finally, radar and radio emission data can penetrate several meters into the surface, providing an estimate of subsurface density and roughness.

Christensen, Philip R.↗

Numerical simulation of the hot-tail runaway electron production mechanism using CQL3D and comparison with Smith–Verwichte analytical model

Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.

Physics↗

Identified charged hadron production in Au+Au collisions at $\sqrt{s_{NN}}$ = 54.4 GeV with the STAR detector

Here, we present results on the production of 𝜋±, 𝐾±, 𝑝, and $\bar{𝑝}$ in Au+Au collisions at $\sqrt{𝑠_{N⁢N}}$ = 54.4~GeV using the STAR detector at RHIC, at mid-rapidity (|𝑦|< 0.1). Invariant yields of these particles as a function of transverse momentum are shown. We determine bulk properties such as integrated particle yields (𝑑 ⁢𝑁 /𝑑 ⁢𝑦 ), mean transverse momentum (⟨𝑝 𝑇 ⟩), particle ratios, which provide insight into the particle production mechanisms. Additionally, the kinetic freeze-out parameters (𝑇 kin and ⟨𝛽 𝑇 ⟩), which provide information about the dynamics of the system at the time of freeze-out, are obtained. The Bjorken energy density (𝜖 𝐵⁢𝐽 ), which gives an estimate of the energy density in the central rapidity region of the collision zone at the formation time 𝜏, is calculated and presented as a function of multiplicity for various energies. The results are compared with those from the models such as A Multi-Phase Transport (AMPT) and Heavy Ion Jet INteraction Generator (HIJING) for further insights.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗

Martian cratering V - Toward an empirical Martian chronology, and its implications

This paper estimates ages of Martian features by dividing least-square-fitted crater densities by an estimated Martian crater production rate. This method is more rigorous than methods based purely on assumptions about correspondence of lunar and Martian cratering rates. Results are interpreted as supporting the conclusion of Burns et al (in press) that Mars shifted climatic states due to a change in obliguity caused by buildup of massive volcanics. Prior to a few gy. ago, conditions favored erosion, deposition, and fluvial channel formation. Recent volcanics postdate this era, and the surfaces of the major shields are fairly well constrained in age to a few hundred my, in agreement with Masursky et al.

Hartmann, W. K.↗

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗

Estimation and Modeling of Enceladus Plume Density Using Attitude Control Data Collected by the Cassini Spacecraft During Low-Altitude Enceladus Flybys

The Cassini spacecraft was launched on 15 October 1997. After an interplanetary cruise of almost seven years, it arrived at Saturn on June 30, 2004. Major science objectives of the Cassini mission include investigations of the configuration and dynamics of Saturn's magnetosphere, the structure and composition of the rings, the characterization of several of Saturn's icy satellites, and Titan's atmosphere constituent abundance

Enceladus↗

Predicting Space Weather Effects on Close Approach Events

The NASA Robotic Conjunction Assessment Risk Analysis (CARA) team sends ephemeris data to the Joint Space Operations Center (JSpOC) for conjunction assessment screening against the JSpOC high accuracy catalog and then assesses risk posed to protected assets from predicted close approaches. Since most spacecraft supported by the CARA team are located in LEO orbits, atmospheric drag is the primary source of state estimate uncertainty. Drag magnitude and uncertainty is directly governed by atmospheric density and thus space weather. At present the actual effect of space weather on atmospheric density cannot be accurately predicted because most atmospheric density models are empirical in nature, which do not perform well in prediction. The Jacchia-Bowman-HASDM 2009 (JBH09) atmospheric density model used at the JSpOC employs a solar storm active compensation feature that predicts storm sizes and arrival times and thus the resulting neutral density alterations. With this feature, estimation errors can occur in either direction (i.e., over- or under-estimation of density and thus drag). Although the exact effect of a solar storm on atmospheric drag cannot be determined, one can explore the effects of JBH09 model error on conjuncting objects' trajectories to determine if a conjunction is likely to become riskier, less risky, or pass unaffected. The CARA team has constructed a Space Weather Trade-Space tool that systematically alters the drag situation for the conjuncting objects and recalculates the probability of collision for each case to determine the range of possible effects on the collision risk. In addition to a review of the theory and the particulars of the tool, the different types of observed output will be explained, along with statistics of their frequency.

Space Weather↗

Mixture densities, maximum likelihood, and the EM algorithm

The problem of estimating the parameters which determine a mixture density is reviewed as well as maximum likelihood estimation for it. A particular iterative procedure for numerically approximating maximum likelihood estimates for mixture density problems is considered. This EM algorithm, is a specialization to the mixture density context of a general algorithm of the same name used to approximate maximum likelihood estimates for incomplete data problems. The formulation and theoretical and practical properties of the EM algorithm for mixture densities are discussed focussing in particular on mixtures of densities from exponential families.

Redner, R. A.↗

Satellite Geopotential Anomaly Constraints for the Crust of the Greenland-Iceland Region

Satellite magnetometer observations of the Greenland-Iceland region compare quite well with lower altitude data. The satellite magnetic data suggest magnetically enhanced crust was emplaced by the Iceland Plume. Crustal thicknesses, which may be more than 30 km for the Greenland-Scotland Ridge, were obtained from inversion of the compensating terrain gravity effects that were estimated by spectral correlation analysis of the free-air gravity anomalies and terrain gravity effects. Regional magnetic anomaly maxima overlie possible thickened crust from eastern Iceland to the Greenland Coast. The Iceland-Faroe Ridge may involve thinner crust than the Greenland-Iceland portion of the Greenland-Scotland Ridge. The gravity derived crustal model exceeds a 0.7 correlation with available seismic estimates. In thermally active areas our gravity Moho estimates are systematically deeper than the seismic estimates suggesting local density reductions of the underlying lower crust/upper mantle. In south central Greenland, on the other hand, the gravity Moho estimates are shallower than seismic estimates to suggest a local enhancement of the lower crust/upper mantle density. The dichotomous crust of the Greenland-Iceland and Iceland-Faroe Ridges suggests unequal crustal development by the Iceland Plume and the Mid-Atlantic Ridge, where more crustal material may have been contributed to the North Atlantic Plate than the Eurasian Plate. A new thermal modeling scheme based on Poisson's relation between point pole gravity and thermal potentials allows estimation of magnetic crustal thicknesses. Subsequent magnetic anomaly inversion for susceptibility contrasts infers crustal development of the Greenland-Scotland Ridge by temporally variable pulses in plume strength.

vonFrese, R. R.↗

Constrained Kalman Filtering Via Density Function Truncation for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops an analytic method of incorporating state variable inequality constraints in the Kalman filter. The resultant filter truncates the PDF (probability density function) of the Kalman filter estimate at the known constraints and then computes the constrained filter estimate as the mean of the truncated PDF. The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is demonstrated via simulation results obtained from a turbofan engine model. The turbofan engine model contains 3 state variables, 11 measurements, and 10 component health parameters. It is also shown that the truncated Kalman filter may be a more accurate way of incorporating inequality constraints than other constrained filters (e.g., the projection approach to constrained filtering).

Simon, Dan↗

Small grains and IRAS colors

The paper studies how infrared colors of dust emission from the interstellar medium vary with the energy density of the radiation field on the basis of IRAS observation of the California Nebula. The data suggest that color variations result from a combinatin of equilibrium emission from large grains, and nonequilibrium emission from small grains, with destruction of the small grains emitting at 12 microns at high energy density; it is estimated that 80 percent of these small particles are destroyed for an energy density in ultraviolet photons larger than 50 times that of the average interstellar radiation field in the solar neighborhood. In a color-color diagram, I(v)(60 microns)/I(v)(100 microns) versus I(v)(12 microns)/I(v)(25 microns), the California Nebula measurements at various distances to the ionizing star Zeta Per follow a sequence similar to that of galaxies. This result shows that the position of a galaxy along this sequence is a measure of the intensity of the radiation field in the regions responsible for the infrared emission.

Boulanger, F.↗

Quantifying the Storm Time Thermospheric Neutral Density Variations Using Model and Observations

Accurate determination of thermospheric neutral density holds crucial importance for satellited rag calculations. The problem is twofold and involves the correct estimation of the quiet time climatology and storm time variations. In this work, neutral density estimations from two empirical and threephysics based models of the ionosphere thermosphere are compared with the neutral densities along the Challenging MicroSatellite Payload satellite track for six geomagnetic storms. Storm time variations are extracted from neutral density by (1) subtracting the mean difference between model and observation (bias),(2) setting climatological variations to zero, and (3) multiplying model data with the quiet time ratio between the model and observation. Several metrics are employed to evaluate the model performances. We find that the removal of bias or climatology reveals actual performance of the model in simulating the storm time variations. When bias is removed, depending on event and model, storm time errors in neutral density can decrease by an amount of 113% or can increase by an amount of 12% with respect to error in models with quiet time bias. It is shown that using only average and maximum values of neutral density to determine the model performances can be misleading since a model can estimate the averages fairly well but may not capture the maximum value or vice versa. Since each of the metrics used for determining model performances provides different aspects of the error, among these, we suggest employing mean absolute error, prediction efficiency, and normalized root mean square error together as a standard set ofmetrics for the neutral density.

Eyiguler, E. Ceren Kalafatoglu↗

Efficient star formation in the spiral arms of M51

The molecular, neutral, and ionized hydrogen distributions in the Sbc galaxy M51 (NGC 5194) are compared. To estimate H2 surface densities observations of the CO (J = 1 - 0) transition were made in 60 positions out to a radius of 155 arcsec. Extinction-corrected H-alpha intensities were used to compute the detailed massive star formation rates (MSFRs) in the disk. Estimates of the gas surface density, the MSFR, and the ratio of these quantities, MSFR/sigma(p), were then examined. The spiral arms were found to exhibit an excess gas density, measuring between 1.4 and 1.6 times the interarm values at 45 arcsec resolution. The total (arm and interarm) gas content and massive star formation rates in concentric annuli in the disk of M51 were computed. The two quantities fall off together with radius, yielding a relatively constant MSFR/sigma(p) with radius. This behavior is not explained by current models of star formation in galactic disks.

Lord, Steven D.↗