Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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 1,333 records · Page 74

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗

ECCE sensitivity studies for single hadron transverse single spin asymmetry measurements

Here, we performed feasibility studies for various single transverse spin measurements that are related to the Sivers effect, transversity and the tensor charge, and the Collins fragmentation function. The processes studied include semi-inclusive deep inelastic scattering (SIDIS) where single hadrons (pions and kaons) were detected in addition to the scattered DIS lepton. The data were obtained in pythia6 and geant4 simulated e+p collisions at 18 GeV on 275 GeV, 18 on 100, 10 on 100, and 5 on 41 that use the ECCE detector configuration. Typical DIS kinematics were selected, most notably Q 2 >1 GeV 2 , and cover the x range from 10 -4 to 1. The single spin asymmetries were extracted as a function of x and Q 2 , as well as the semi-inclusive variables z, which corresponds to the momentum fraction the detected hadron carries relative to the struck parton, and P T , which corresponds to the transverse momentum of the detected hadron relative to the virtual photon. They are obtained in azimuthal moments in combinations of the azimuthal angles of the hadron transverse momentum and transverse spin of the nucleon relative to the lepton scattering plane. In order to extract asymmetries, the initially unpolarized MonteCarlo was re-weighted in the true kinematic variables, hadron types and parton flavors based on global fits of fixed target SIDIS experiments and e + e – annihilation data. The expected statistical precision of such measurements is extrapolated to 10 fb –1 and potential systematic uncertainties are approximated given the deviations between true and reconstructed yields. Similar neutron information is obtained by comparing the ECCE e+p pseudo-data with the same from the EIC Yellow Report and scaling the corresponding Yellow Report e+ 3 He pseudo-data uncertainties accordingly. The impact on the knowledge of the Sivers functions, transversity and tensor charges, and the Collins function has then been evaluated in the same phenomenological extractions as in the Yellow Report. Finally, the impact is found to be comparable to that obtained with the parametrized Yellow Report detector and shows that the ECCE detector configuration can fulfill the physics goals on these quantities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗

Numerical Analysis of High-Order Modes in SRF Resonators for Particle Accelerators

Over the past decades, superconducting technology has rapidly evolved towards high accelerating gradients and low surface resistance, making it possible to operate particle accelerators with high average beam currents and large duty factors. However, RF losses due to coherent excitation of the HOM become the limiting factor for these regimes. Unlike the cavity operating mode, which is tuned separately, the HOM parameters can significantly vary from one cavity to another due to finite mechanical tolerances during the manufacturing process. Thus, it is of utmost importance to know the HOM parameter spread in advance in order to predict unexpected cryogenic losses, overheating of beam line components and maintain stable beam dynamics. In this paper, we present a method for generating cavity geometry with an arbitrary spread of mechanical imperfections and numerically evaluating HOM statistics. Knowing the spread of HOM parameters, we calculated the probability of resonant HOM losses in SRF accelerating cavities used in CW beam current machines such as the PIP-II and LCLS-II linacs, as well as for the SRF crab-cavity for the ILC project. Finally, we present experimental results of HOM spectra measurements in hundreds of 1.3 GHz cavities installed in LCLS-II cryomodules. Studying the effects of HOM excitation results in specifications of the SRF cavity and cryomodule and can significantly impact the efficiency and reliability of the machine operation.

Lunin, Andrei [Fermilab] (ORCID:0000000290096792)↗

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗

Mars gravitational field estimation error

The error covariance matrices associated with a weighted least-squares differential correction process have been analyzed for accuracy in determining the gravitational coefficients through degree and order five in the Mars gravitational potential junction. The results are presented in terms of standard deviations for the assumed estimated parameters. The covariance matrices were calculated by assuming Doppler tracking data from a Mars orbiter, a priori statistics for the estimated parameters, and model error uncertainties for tracking-station locations, the Mars ephemeris, the astronomical unit, the Mars gravitational constant (G sub M), and the gravitational coefficients of degrees six and seven. Model errors were treated by using the concept of consider parameters.

Compton, H. R.↗

Parameter uncertainties in control system design.

Development of a design method for including the effects of parameter uncertainties in the design of linear control systems. The approach taken to this problem may be classified as a special case of the stochastic control problem. Thus the formulation is based on the minimization of the expected value of a quadratic performance index defined in terms of the system state vector. The uncertainty in the value of the performance index is the result of the statistical nature of the system parameters rather than a random input signal. It is shown that the expected value of the performance index may be written as a sum of two terms under the assumption of first-order variations of the system state. The first of these terms expresses the nominal performance of the system when the system parameters assume their mean values. The second term represents the effect of the uncertainties on the expected value of the performance index, and is interpreted as an index of system sensitivity.

Palsson, T.↗

Cratering on Mars. I - Cratering and obliteration history. II Implications for future cratering studies from Mariner 4 reanalysis

It is pointed out that Mars is especially well adapted to statistical studies of crater morphologies for deciphering its geological history. A framework for understanding planetary geomorphological histories from the diameter-frequency relations of different morphological classes of craters described by Chapmam et al. (1970) is extended in order to understand Martian cratering, erosional, and depositional history. The cratering-obliteration history derived is compared with global interpretations considered by Hartman (1973) and Soderblom et al. (1974). An idealized dust-filling model is employed.

Chapman, C. R.↗

Cloud fluctuation statistics

A space-time statistical analysis of total outgoing infrared radiation is used to determine the gross features of day-to-day cloudiness fluctuations over the Pacific Ocean in summer and winter. IR fluctuations arise from the passage of cloudiness systems through a grid box as well as the creation and destruction of cloudiness in the box. Which process dominates depends upon the size of the box relative to the size, speed and persistence time of a typical cloudiness system. In most regions the statistical analysis yields advection speeds characteristic of 700 mb mean flow with spatial dependence resembling the 300 mb mean flow. Spatial scales less than 2000 km predominate, smaller scales having less persistence. Characteristic time scales are on the order of one or two days, even for a grid box spanning the entire North Pacific storm track. This result is remarkable in view of the much longer time scales commonly associated with atmospheric disturbances. Apparently many cloudiness systems are created and destroyed during the lifetime of a single disturbance.

Cahalan, R. F.↗

Gravity waves

Gravity waves contributed to the establishment of the thermal structure, small scale (80 to 100 km) fluctuations in velocity (50 to 80 m/sec) and density (20 to 30%, 0 to peak). Dominant gravity wave spectrum in the middle atmosphere: x-scale, less than 100 km; z-scale, greater than 10 km; t-scale, less than 2 hr. Theorists are beginning to understand middle atmosphere motions. There are two classes: Planetary waves and equatorial motions, gravity waves and tidal motions. The former give rise to variability at large scales, which may alter apparent mean structure. Effects include density and velocity fluctuations, induced mean motions, and stratospheric warmings which lead to the breakup of the polar vortex and cooling of the mesosphere. On this scale are also equatorial quasi-biennial and semi-annual oscillations. Gravity wave and tidal motions produce large rms fluctuations in density and velocity. The magnitude of the density fluctuations compared to the mean density is of the order of the vertical wavelength, which grows with height. Relative density fluctuations are less than, or of the order of 30% below the mesopause. Such motions may cause significant and variable convection, and wind shear. There is a strong seasonal variation in gravity wave amplitude. Additional observations are needed to address and quantify mean and fluctuation statistics of both density and mean velocity, variability of the mean and fluctuations, and to identify dominant gravity wave scales and sources as well as causes of variability, both temporal and geographic.

Fritts, David↗

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

Satellite Image Mosaic Engine

A computer program automatically builds large, full-resolution mosaics of multispectral images of Earth landmasses from images acquired by Landsat 7, complete with matching of colors and blending between adjacent scenes. While the code has been used extensively for Landsat, it could also be used for other data sources. A single mosaic of as many as 8,000 scenes, represented by more than 5 terabytes of data and the largest set produced in this work, demonstrated what the code could do to provide global coverage. The program first statistically analyzes input images to determine areas of coverage and data-value distributions. It then transforms the input images from their original universal transverse Mercator coordinates to other geographical coordinates, with scaling. It applies a first-order polynomial brightness correction to each band in each scene. It uses a data-mask image for selecting data and blending of input scenes. Under control by a user, the program can be made to operate on small parts of the output image space, with check-point and restart capabilities. The program runs on SGI IRIX computers. It is capable of parallel processing using shared-memory code, large memories, and tens of central processing units. It can retrieve input data and store output data at locations remote from the processors on which it is executed.

Plesea, Lucian↗

Determination of the Limiting Magnitude

The limiting magnitude of an optical camera system is an important property to understand since it is used to find the completeness limit of observations. Limiting magnitude depends on the hardware and software of the system, current weather conditions, and the angular speed of the objects observed. If an object exhibits a substantial angular rate during the exposure, its light spreads out over more pixels than the stationary stars. This spreading causes the limiting magnitude to be brighter when compared to the stellar limiting magnitude. The effect, which begins to become important when the object moves a full width at half max during a single exposure or video frame. For targets with high angular speeds or camera systems with narrow field of view or long exposures, this correction can be significant, up to several magnitudes. The stars in an image are often used to measure the limiting magnitude since they are stationary, have known brightness, and are present in large numbers, making the determination of the limiting magnitude fairly simple. In order to transform stellar limiting magnitude to object limiting magnitude, a correction must be applied accounting for the angular velocity. This technique is adopted in meteor and other fast-moving object observations, as the lack of a statistically significant sample of targets makes it virtually impossible to determine the limiting magnitude before the weather conditions change. While the weather is the dominant factor in observing satellites, the limiting magnitude for meteors also changes throughout the night due to the motion of a meteor shower or sporadic source radiant across the sky. This paper presents methods for determining the limiting stellar magnitude and the conversion to the target limiting magnitude.

Kingery, Aaron↗

Landslide Hazard is Projected to Increase Across High Mountain Asia

High Mountain Asia has long been known as a hotspot for landslide risk, and studies have suggested that landslide hazard is likely to increase in this region over the coming decades. Extreme precipitation may become more frequent, with a nonlinear response relative to increasing global temperatures. However, these changes are geographically varied. This article maps probable changes to landslide hazard, as shown by a landslide hazard indicator (LHI) derived from downscaled precipitation and temperature. In order to capture the nonlinear response of slopes to extreme precipitation, a simple machine-learning model was trained on a database of landslides across High Mountain Asia to develop a regional LHI. This model was applied to statistically downscaled data from the 30 members of the Seamless System for Prediction and Earth System Research large ensembles to produce a range of possible outcomes under the Shared Socioeconomic Pathways 2-4.5 and 5-8.5. The LHI reveals that landslide hazard will increase in most parts of High Mountain Asia. Absolute increases will be highest in already hazardous areas such as the Central Himalaya, but relative change is greatest on the Tibetan Plateau. Even in regions where landslide hazard declines by year 2100, it will increase prior to the mid-century mark. However, the seasonal cycle of landslide occurrence will not change greatly across High Mountain Asia. Although substantial uncertainty remains in these projections, the overall direction of change seems reliable. These findings highlight the importance of continued analysis to inform disaster risk reduction strategies for stakeholders across High Mountain Asia.

Thomas A Stanley↗

Probing me Reverse Shock in an Oxygen-Rich Supernova Remnant

The aim of this project is to examine the O VI emission at three positions around the X- ray bright ring of the remnant in order to investigate the relation between the O VI emission, the X-ray O VII and O VIII emission, and the optical [OIII} emission, and how these vary around the rim of the remnant. All three pointings and the background pointing have now been observed; the archive notification for the most recent dataset was Oct 30, 2003. After reprocessing and screening, the net exposure time for the SE exposure is only 54 percent of the approved time (15 kilosec). for the SE exposure, the available statistics are not good enough for analysis. A request for reobservation to make up for the lost time in the SE pointing has been approved. Broad O VI 1032 and O VI 1038 emission is detected with velocity width of at least 800 km/s, and possibly exceeding 1000 km/s. The Flanagan et al. analysis of the Chandra grating data show bulk velocities in the X-ray gas of order +/- 1000 km/s. In the region of the FUSE E0lO2-SE pointing, the Chandra data indicate both blue-shifted and red-shifted emission. Analysis of the velocity structure of the O VI emission will provide additional constraints on the kinematics of the gas: is it emission from a tilted expanding barrel, or a more symmetric expansion? The O VI fluxes are also needed to assess whether the O VI is radiation from recombining O VII or instead fiom cooler gas ionizing toward O VII. The emission is faint, however, which complicates the analysis. Because the lines are so broad, absorption by intervening H2, CII, and foreground OVI must be considered. A number of stars in the SMC have been observed which provide information on foreground OVI absorption. The initial analyses have concentrated on the Li1F channel since the guidance is based on that channel. The Li2F channel is being examined exposure by exposure to see if any of the data can be used to improve the signal to noise in the Li1F data.

Gaetz, Terrance↗

Organometallic Polymeric Conductors

For aerospace applications, the use of polymers can result in tremendous weight savings over metals. Suitable polymeric materials for some applications like EMI shielding, spacecraft grounding, and charge dissipation must combine high electrical conductivity with long-term environmental stability, good processability, and good mechanical properties. Recently, other investigators have reported hybrid films made from an electrically conductive polymer combined with insulating polymers. In all of these instances, the films were prepared by infiltrating an insulating polymer with a precursor for a conductive polymer (either polypyrrole or polythiophene), and oxidatively polymerizing the precursor in situ. The resulting composite films have good electrical conductivity, while overcoming the brittleness inherent in most conductive polymers. Many aerospace applications require a combination of properties. Thus, hybrid films made from polyimides or other engineering resins are of primary interest, but only if conductivities on the same order as those obtained with a polystyrene base could be obtained. Hence, a series of experiments was performed to optimize the conductivity of polyimide-based composite films. The polyimide base chosen for this study was Kapton. 3-MethylThiophene (3MT) was used for the conductive phase. Three processing variables were identified for producing these composite films, namely time, temperature, and oxidant concentration for the in situ oxidation. Statistically designed experiments were used to examine the effects of these variables and synergistic/interactive effects among variables on the electrical conductivity and mechanical strength of the films. Multiple linear regression analysis of the tensile data revealed that temperature and time have the greatest effect on maximum stress. The response surface of maximum stress vs. temperature and time (for oxidant concentration at 1.2 M) is shown. Conductivity of the composite films was measured for over 150 days in air at ambient temperature. The conductivity of the films dropped only half an order of magnitude in that time. Films aged under vacuum at ambient temperature diminished slightly in conductivity in the first day, but did not change thereafter. An experimental design approach will be applied to maximize the efficiency of the laboratory effort. The material properties (initial and long term) will also be monitored and assessed. The experimental results will add to the existing database for electrically conductive polymer materials. Attachments: 1) Synthesis Crystal Structure, and Polymerization of 1,2:5,6:9,10-Tribenzo-3,7,11,13-tetradehydro(14) annulene. 2) Reinvestigation of the Photocyclization of 1,4-Phenylene Bis(phenylmaleic anhydride): Preparation and Structure of (5)Helicene 5,6:9,10-Dianhydride. 3) Preparation and Structure Charecterization of a Platinum Catecholate Complex Containing Two 3-Ethynyltheophone Groups. and 4) Rigid-Rod Polymers Based on Noncoplanar 4,4'-Biphenyldiamines: A Review of Polymer Properties vs Configuration of Diamines.

Youngs, Wiley J.↗

A statistical mechanical calculation of the thermodynamic properties of interstitial solid solutions involving second nearest neighbor interactions.

A previous calculation of the thermodynamic properties of interstitial solid solutions based on the technique of Kirkwood expansions has been extended to include the effects of second nearest neighbor solute atom mutual interactions. The error inherent in the first order (or quasi-chemical) counting of the degeneracy of the solution crystal is avoided. It is shown that, at high temperatures, even strong second nearest neighbor solute mutual interactions have a negligible effect on the entropy of the solution and a small, temperature-dependent effect on the solute partial enthalpy.

Alex, K.↗