Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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 739 records · Page 41

Sensitivity Analysis of Particle-In-Cell Modeling Parameters in Settling Bed, Bubbling Fluidized Bed and Circulating Fluidized Bed

The objective of the work presented is to perform a preliminary sensitivity analysis of particle-in-cell (PIC) model parameters when applied to settling bed, bubbling fluidized bed, and circulating fluidized bed simulations. These examples correspond to widely different flow conditions commonly seen in chemical engineering applications. Simulations were performed using the PIC method in the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). As part of the non-intrusive uncertainty quantification (UQ) analysis, simulation campaigns were generated using Nodeworks. Sampling locations or settings for PIC model parameters were determined using the Latin Hypercube method. Response surfaces were created using radial basis functions (RBF), and Sobol’ indices were estimated to quantify the influence of model parameters on the quantities of interest (QoI). This study marks a first step towards systematically determining optimal ranges for model parameters used in MFiX-PIC. Based on limited experience, it is expected that these values would depend strongly on flow conditions. Given the complexity of the multiphase flow systems under analysis, a non-intrusive UQ based approach is used to identify the most influential parameters in each case. This prior knowledge will help in proposing an effective design of experiments (DoE) and determine optimal parameters through techniques such as deterministic or Bayesian calibration, which will be pursued in the future.

42 ENGINEERING↗

Improving the Advancement of Women in Computer and Computational Science Research with the CRA-W Career Mentoring Workshops (Final Report)

The mission of the Computing Research Association’s Committee on Widening Participation in Computing Research (CRA-WP) is to widen the participation and improve the access, opportunities, and positive experiences of individuals from groups underrepresented in computing research and education. CRA-WP programs serve this overarching goal at all career stages; in addition, CRA-WP, through the formation of the CRA Center for Evaluating the Research Pipeline (CERP), has developed a methodology for thoroughly evaluating the success of its programs by comparing a nationwide sample of students, researchers, and faculty (non-participants) to program participants. Achieving these objectives requires that an increasing number of individuals from populations underrepresented in computing start and progress to the next stage while understanding and supporting the myriad computing pathways. CRA-WP offers programs for participants from undergraduate to senior professional levels. Different career stages need different types of interventions, and the goal of all Alliance program activities can be described within the unifying framework of Social Cognitive Career Theory, which finds that interest in and choice of a particular career path will be increased by interventions that improve one or more of the following: (1) outcome expectations (understanding and valuing the rewards of a particular outcome), (2) self-efficacy (a belief that one can successfully achieve an outcome), and (3) social supports that help one persist and overcome obstacles.

97 MATHEMATICS AND COMPUTING↗

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Meteorological data from the Billy Bar site at the East River, Colorado

This dataset includes Quality Assessed and Quality Controlled (QA/QC) meteorological data from the Billy Bar field site in the East River Watershed, Colorado in order to inform watershed hydrobiogeochemical processes. The data includes 1-hour aligned time series of Solar Radiation, Wind Speed, Wind Direction, Air Temperature, Relative Humidity, Barometric Pressure, and Precipitation. For each parameter, the QA flags are given: 1 is for extreme (potential abnormal values), and 0 is for values in the expected range. The Rmarkdown document is a stand-alone file/notebook that preserves the text, code, and code results, as well as formatting contained in the original R script. The Data profiling report is a summary and graphical presentation of the data exploration process for data analysis and model building, so that users could focus on understanding data and extracting insights. The report provides a summary of each variable and does data profiling.

54 ENVIRONMENTAL SCIENCES↗

Reinterpretation of the Pioneer 6 bow shock crossing.

This paper re-examines and reinterprets the interesting bow shock crossing of Pioneer 6 on December 16, 1965, by combining the high-resolution data for the magnetic field measurements with all the available plasma data. It is shown that, although a small correction of the magnetic field data improves the validity test of the MHD Rankine-Hugoniot (R-H) conditions, it is the plasma density measurements that are the principal source of the remaining inconsistencies. With the exception of the densities, the best agreement between the measurements and the R-H relations is obtained when the specific heat ratio is slightly less than the classical value of 5/3 expected for a magnetized plasma in the magnetosheath. A more accurate shock normal determination has been made, and it is found that the shock itself is possibly in motion with respect to the earth and that locally it may be approximated as a normal shock.

Mariani, F.↗

F-15 flight flutter test program

The modes to be observed during the F-15 flight flutter test program were selected on the basis of the results of analytical studies, wind tunnel tests, and ground vibration tests. The modes (both symmetrical and antisymmetrical) tracked on this basis were: fin first bending, fin torsion, fin tip roll, stabilator bending, stabilator pitch, boom lateral bending, boom torsion, boom vertical bending, wing first bending, wing second bending, wing first torsion, outer wing torsion, and aileron rotation. Data obtained for these various modes were evaluated in terms of damping versus airspeed at 1525 m (5000 ft), damping versus altitude at the cross-section Mach numbers (to extrapolate to the damping value to be expected at sea level), and flutter boundaries on the basis of flutter margin of various modal pairs representing potential flutter mechanisms. Results of these evaluations are summarized in terms of minimum predicted flutter margin for the various mechanisms.

Katz, H.↗

Shear dampers for deflections of micro-inch magnitudes

The rotation of the TV camera on the Mariner Jupiter/Saturn 1977 (MJS77) spacecraft will excite the long magnetometer and heavy RTG booms which are slightly damped perturbances that have long settling times. To ensure clear TV pictures additional damping was required over micro-inch displacements and at low frequencies. The solution incorporated elastomeric shear dampers at the base of these booms which possessed high loss tangent values in the expected frequency and temperature ranges. Presented are the shear damper design, results of full scale horizontal boom tests and relative effects of damper sizing and temperature. Boom instrumentation and excitation used for the full scale tests are also discussed.

Otth, D.↗

Positron annihilation spectroscopy with magnetically analyzed beams

Lifetime measurements with magnetically analyzed positron beams were made in condensed media with uniform and non-uniform properties. As expected, the lifetime values with magnetically analyzed positron beams in uniform targets are similar to those obtained with conventional positron sources. The lifetime values with magnetically analyzed beams in targets which have non-uniform properties vary with positron energy and are different from the conventional positron source derived lifetime values in these targets.

Singh, J. J.↗

Discrimination of coherent features in turbulent boundary layers by the entropy method

Entropy in information theory is defined as the expected or mean value of the measure of the amount of self-information contained in the ith point of a distribution series x sub i, based on its probability of occurrence p(x sub i). If p(x sub i) is the probability of the ith state of the system in probability space, then the entropy, E(X) = - sigma p(x sub i) logp (x sub i), is a measure of the disorder in the system. Based on this concept, a method was devised which sought to minimize the entropy in a time series in order to construct the signature of the most coherent motions. The constrained minimization was performed using a Lagrange multiplier approach which resulted in the solution of a simultaneous set of non-linear coupled equations to obtain the coherent time series. The application of the method to space-time data taken by a rake of sensors in the near-wall region of a turbulent boundary layer was presented. The results yielded coherent velocity motions made up of locally decelerated or accelerated fluid having a streamwise scale of approximately 100 nu/u(tau), which is in qualitative agreement with the results from other less objective discrimination methods.

Corke, T. C.↗

Impact of GRM: New evidence from the Soviet Union

Gravity information released by the Soviet Union allows the quantitative assessment of how the geopotential research mission (GRM) mission might effect the ability to use global gravity data for continental tectonic interpretation. The information is of an isostatic response spectra for eight individual tectonic units in the USSR. The regions examined include the Caroathians, Caucasus, Urals, Pamirs, Tien-Shan, Altal, Chersky Ridge, and East Siberian Platform. The 1 deg x 1 deg gravity data are used to calculate the admittances are used in two different sorts of tectonic studies of mountain belts in the USSR: (1) interpretation of isostatic responses in terms of plate models of compensation for mountainous terrain. Using geologic information concerning time of the orogeny, lithospheric plates involved, and polarity of subduction in collision zones, they convert the best-fitting flexural rigidity to an elastic plate thickness for the lithospheric plate inferred to underlie the mountains; the isostatic admittance functions is an attempt to directly model gravity and topography data for a few select regions in the Soviet Union. By knowing the value of the expected correlation between topography and gravity from the admittances, the Artemjev's map in mountainous areas can be calibrated, and the maps are converted back to Bouguer gravity. This procedure is applied to the Caucasus and southern Urals.

Mcnutt, M.↗

Predictive monitoring research: Summary of the PREMON system

Traditional approaches to monitoring are proving inadequate in the face of two important issues: the dynamic adjustment of expectations about sensor values when the behavior of the device is too complex to enumerate beforehand, and the selective but effective interpretation of sensor readings when the number of sensors becomes overwhelming. This system addresses these issues by building an explicit model of a device and applying common-sense theories of physics to model causality in the device. The resulting causal simulation of the device supports planning decisions about how to efficiently yet reliably utilize a limited number of sensors to verify correct operation of the device.

Doyle, Richard J.↗

He I 10830-A observations of Seyfert 2 galaxies

Measurements of the 10830-A transition of neutral helium are presented for 12 Seyfert 2 galaxies and the emission-line radio galaxy NGC 1275. The observed 10830-A to He I 5876-A ratios for these objects range from 4 to 24, showing that the objects are reddened. For a number of these galaxies, the observed H-alpha/H-beta values exceed that expected with uniform mixing of dust and gas. The results suggest that one or more mechanisms are destroying the 10830-A photons or are depopulating the 2 3S level. Three previously proposed mechanisms to explain this process are evaluated and it is found that all of the proposed mechanisms are lacking.

Rudy, Richard J.↗

Approximation algorithms for planning and control

A control system operating in a complex environment will encounter a variety of different situations, with varying amounts of time available to respond to critical events. Ideally, such a control system will do the best possible with the time available. In other words, its responses should approximate those that would result from having unlimited time for computation, where the degree of the approximation depends on the amount of time it actually has. There exist approximation algorithms for a wide variety of problems. Unfortunately, the solution to any reasonably complex control problem will require solving several computationally intensive problems. Algorithms for successive approximation are a subclass of the class of anytime algorithms, algorithms that return answers for any amount of computation time, where the answers improve as more time is allotted. An architecture is described for allocating computation time to a set of anytime algorithms, based on expectations regarding the value of the answers they return. The architecture described is quite general, producing optimal schedules for a set of algorithms under widely varying conditions.

Boddy, Mark↗

A model for the Space Shuttle Main Engine High Pressure Oxidizer Turbopump shaft seal system

A model of the High Pressure Oxidizer Turbopump (HPOTP) shaft seal system on the Space Shuttle Main Engine (SSME) is described. The model predicts the fluid properties and flow rates throughout this system for a number of conditions simulating failed seals. The results agree well with qualitative expectations and redline values but cannot be verified with actual data due to the lack thereof. The results indicate that each failure mode results in a unique distribution of properties throughout the seal system and can therefore be individually identified given the proper instrumentation. Furthermore, the detection process can be built on the principle of qualitative reasoning without the use of exact fluid property values. A simplified implementation of the model which does not include the slinger/labyrinth seal combination has been developed and will be useful for inclusion in a real-time diagnostic system.

Paxson, Daniel E.↗

The best-fit universe

Inflation provides very strong motivation for a flat Universe, Harrison-Zel'dovich (constant-curvature) perturbations, and cold dark matter. However, there are a number of cosmological observations that conflict with the predictions of the simplest such model: one with zero cosmological constant. They include the age of the Universe, dynamical determinations of Omega, galaxy-number counts, and the apparent abundance of large-scale structure in the Universe. While the discrepancies are not yet serious enough to rule out the simplest and most well motivated model, the current data point to a best-fit model with the following parameters: Omega(sub B) approximately equal to 0.03, Omega(sub CDM) approximately equal to 0.17, Omega(sub Lambda) approximately equal to 0.8, and H(sub 0) approximately equal to 70 km/(sec x Mpc) which improves significantly the concordance with observations. While there is no good reason to expect such a value for the cosmological constant, there is no physical principle that would rule out such.

Turner, Michael S.↗

The abundance of boron in three halo stars

B abundances for three halo stars: HD 140283, HD 19445, and HD 201891 are presented. Using recent determinations of the Be abundance in HD 140283, B/Be of 10 +5/-4 is found for this star, and similar ratios are inferred for HD 19445 and HD 201891. This ratio is equal to the minimum value of 10 expected from a synthesis of B and Be by high-energy cosmic-ray spallation reactions in the interstellar medium. It is shown that the accompanying synthesis of Li by alpha on alpha fusion reactions is probably a minor contributor to the observed 'primordial' Li of halo stars. The observed constant ratios of B/O and Be/O are expected if the principal channel of synthesis involves cosmic-ray CNO nuclei from the supernovae colliding with interstellar protons.

Duncan, Douglas K.↗

Theoretical determination of the alkali-metal superoxide bond energies

The bond dissociation energies for the alkali-metal superoxides have been computed using extensive Gaussian basis sets and treating electron correlation at the modified coupled-pair functional level. Our computed D0 values are 61.4, 37.2, 40.6, and 38.4 kcal/mol for LiO2, NaO2, KO2, and RbO2, respectively. These values, which are expected to be lower bounds and accurate to 2 kcal/mol, agree well with some of the older flame data, but rule out several recent experimental measurements.

Partridge, Harry↗

On the decay of outbursts in dwarf novae nad X-ray novae

We perform computations using a time-dependent model for the accretion disk limit-cycle mechanism to examine the decay of the optical light following the peak of a dwarf nova outburst. We present the results of a parameter study of the physical input variables which affect the decay rate. In the model, the decay is brought about by a cooling transition front which begins at large radii in the disk and moves inward. The nature of the decay is strongly influenced by the radial dependence of the accretion disk viscosity parameter alpha. To obtain exponential decays for typical dwarf nova parameters, we require alpha proportional to r(exp epsilon(sub 0)), where epsilon(sub 0) approximately = 0.3-0.4. The exact value of epsilon(sub 0) which produces exponential decays depends on factors such as the mass of the accreting star and the inner radius of the accretion disk. Therefore, the observed ubiquity of exponential decays in two different types of systems (dwarf novae and X-ray novae) leads us to believe that alpha is an unnatural scaling for the viscosity. The physics of the cooling transition front must be self-regulating in that the timescale (-parital derivative of lnSigma(r)/partial derivative +)(exp -1) (where Sigma is the surface density) for mass extraction across the front remains constant. This may be consistent with a scaling alpha proportional to (h/r)(exp n), where h is the local disk semi-thickness and n approximately 1-2. As regards the speed of the cooling front, we find v(sub F)(r) proportional to r(exp p), where p approximately 3 at large radii, with an abrupt transition to p approximately 0 at some smaller radius. The r(exp 3) dependence is much steeper than has been found by previous workers and appears to result from the strong variation of specific heat within the cooling front when the front resides at a large radius in the disk. The outflow of disk material across the cooling front causes a significant departure of dln T(sub dff0/dln r from the standard value of -0.75 (expected from steady state accretion) within about 0.2 dex in radius of the break associated with the cooling front -- T(sub eff) aproximately 10(exp 3.9) K (r/10(exp 10 cm)) (exp -0.1). These effects should be observable with eclipse mapping. Finally, it appears that the relatively slow decay rate for the optical flux in the 1975 outburst of A0620-00 can be accounted for if the primary is a approximately 10 Solar mass black hole.

Cannizzo, John K.↗