Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic processes.”

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 721 records · Page 40

FORECAST FOR THE PLANCK PRECISION ON THE TENSOR-TO-SCALAR RATIO AND OTHER COSMOLOGICAL PARAMETERS

The Planck satellite, successfully launched on 2009 May 14 to measure with unprecedented accuracy the primary cosmic microwave background (CMB) anisotropies, is operating as expected. The Standard Model of the Universe ('concordance' model) provides the current realistic context to analyze the CMB and other cosmological/astrophysical data, inflation in the early universe being part of it. The Planck performance for the crucial primordial parameter r, the tensor-to-scalar ratio related to primordial B-mode polarization, will depend on the quality of data analysis and interpretation. The Ginzburg-Landau (G-L) approach to inflation allows us to take high benefit of the CMB data. The fourth-degree double-well inflaton potential gives an excellent fit to the current CMB+LSS data. We evaluate the Planck precision to the recovery of cosmological parameters, taking into account a reasonable toy model for residuals of systematic effects of instrumental and astrophysical origin based on publicly available information. We use and test two relevant models: the {Lambda}CDMr model, i.e., the standard {Lambda}CDM model augmented by r, and the {Lambda}CDMrT model, where the scalar spectral index, n{sub s} , and r are related through the theoretical 'banana-shaped' curve r = r(n{sub s}) coming from the G-L theory with a double-well inflaton potential. In the latter case, the analytical expressions for n{sub s} and r are imposed as a hard constraint in a Monte Carlo Markov Chain (MCMC) data analysis. We consider two C{sub l}-likelihoods (with and without B modes) and take into account the white noise sensitivity of Planck (LFI and HFI) in the 70, 100, and 143 GHz channels as well as the residuals from systematic errors and foregrounds. We also consider a cumulative channel of the three mentioned. We produce the sky (mock data) for the CMB multipoles C{sup TT}{sub l} , C{sup TE}{sub l} , C{sup EE}{sub l} , and C{sup BB}{sub l} from the {Lambda}CDMr and {Lambda}CDMrT models and obtain the cosmological parameter marginalized likelihood distributions for the two models. Foreground residuals affect only the cosmological parameters sensitive to the B modes. As expected, the likelihood r distribution is more clearly peaked near the fiducial value (r = 0.0427) in the {Lambda}CDMrT model than in the {Lambda}CDMr model. The best value for r in the presence of residuals turns out to be about r {approx_equal} 0.04 for both the {Lambda}CDMr and the {Lambda}CDMrT models. The {Lambda}CDMrT model is very stable; its distributions do not change by including residuals and the B modes. For r we find 0.028 < r < 0.116 at a 95% confidence level (CL) with the best value r = 0.04. We also compute the B mode detection probability by the most sensitive HFI-143 channel. At the level of foreground residual equal to 30% of our toy model, only a 68% CL (1{sigma}) detection is very likely. For a 95% CL detection (2{sigma}), the level of foreground residual should be reduced to 10% or lower of the adopted toy model. The lower bounds (and most probable value) we infer for r support the searching of CMB B-mode polarization in the current data as well as the planned CMB missions oriented toward B polarization.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics and quantum Zeno effect for a qubit in either a low- or high-frequency bath beyond the rotating-wave approximation

We use a non-Markovian approach to study the decoherence dynamics of a qubit in either a low- or high-frequency bath modeling the qubit environment. This is done for two separate cases: either with measurements or without them. This approach is based on a unitary transformation and does not require the rotating-wave approximation. In the case without measurement, we show that, for low-frequency noise, the bath shifts the qubit energy toward higher energies (blue shift), while the ordinary high-frequency cutoff Ohmic bath shifts the qubit energy toward lower energies (red shift). In order to preserve the coherence of the qubit, we also investigate the dynamics of the qubit subject to measurements (quantum Zeno regime) in two cases: low- and high-frequency baths. For very frequent projective measurements, the low-frequency bath gives rise to the quantum anti-Zeno effect on the qubit. The quantum Zeno effect only occurs in the high-frequency-cutoff Ohmic bath, after counterrotating terms are considered. In the condition that the decay rate due to the two kinds of baths are equal under the Wigner-Weisskopf approximation, we find that without the approximation, for a high-frequency environment, the decay rate should be faster (without measurements) or slower (with frequent measurements, in the Zeno regime), compared to the low-frequency bath case. The experimental implementation of our results here could distinguish the type of bath (either a low- or high-frequency one) and protect the coherence of the qubit by modulating the dominant frequency of its environment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deterministic creation and stabilization of entanglement in circuit QED by homodyne-mediated feedback control

In a solid-state circuit QED system, we demonstrate that a homodyne-current-based feedback can create and stabilize highly entangled two-qubit states in the presence of a moderate noisy environment. Particularly, we present an extended analysis for the current-based Markovian feedback, which leads to an improved feedback scheme. We show that this is essential to achieve a desirable control effect by the use of dispersive measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cosmic discordance: Detection of a modulation in the primordial fluctuation spectrum

As a test of the standard inflationary cosmology, which generically predicts a nearly scale-invariant spectrum of primordial curvature fluctuations, we perform a Markov-Chain Monte-Carlo analysis to search for possible modulations in the power spectrum and determine its shape together with the cosmological parameters using cosmic microwave background radiation data. By incorporating various three-parameter features on the simple power-law spectrum, we find an oscillatory modulation localized around the comoving wave number k{approx_equal}0.009 Mpc{sup -1} at 99.995% confidence level which improves the log-likelihood as much as -{Delta}2lnL{identical_to}{Delta}{chi}{sub eff}{sup 2}=-22. This feature can be detected even if we use only the cross correlation between the temperature and the E-mode polarization anisotropies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Current observational constraints to the holographic dark energy model with a new infrared cutoff via the Markov chain Monte Carlo method

In this paper, the holographic dark energy models with the new infrared cutoff for both the flat case and the nonflat case are confronted with the combined constraints of current cosmological observations: type Ia supernovae, baryon acoustic oscillations, current cosmic microwave background, and the observational Hubble data. By utilizing the Markov chain Monte Carlo method, we obtain the best fit values of the parameters with 1{sigma}, 2{sigma} errors in the flat model: {Omega}{sub b}h{sup 2}=0.0230{sub -0.0010-0.0014}{sup +0.0008+0.0012}, {alpha}=0.9788{sub -0.0927-0.1249}{sup +0.1297+0.1354}, {beta}=0.4739{sub -0.0723-0.0984}{sup +0.0793+0.1055}, {Omega}{sub de0}=0.7869{sub -0.0304-0.0455}{sup +0.0291+0.0370}, {Omega}{sub m0}=0.2131{sub -0.0291-0.0370}{sup +0.0304+0.0455}, H{sub 0}=70.46{sub -2.97-4.02}{sup +2.82+3.89}. In the nonflat model, the constraint results are found in 1{sigma}, 2{sigma} regions: {Omega}{sub b}h{sup 2}=0.0229{sub -0.0010-0.0014}{sup +0.0010+0.0013}, {Omega}{sub k}=0.0014{sub -0.0597-0.0743}{sup +0.0604+0.0604}, {alpha}=0.9637{sub -0.2840-0.3333}{sup +0.2291+0.2894}, {beta}=0.4712{sub -0.0756-0.0961}{sup +0.1412+0.1703}, {Omega}{sub de0}=0.7829{sub -0.2130-0.2386}{sup +0.1588+0.1901}, {Omega}{sub m0}=0.2157{sub -0.1012-0.1241}{sup +0.1562+0.1800}, H{sub 0}=70.64{sub -3.32-4.45}{sup +2.74+3.52}. In the best fit holographic dark energy models, the equation of state of dark energy and the deceleration parameter at present are characterized by w{sub de0}=-0.9278{+-}0.0626, q{sub 0}=-0.5951{+-}0.0586 (flat case) and w{sub de0}=-0.9501{+-}0.1442, q{sub 0}=-0.6164{+-}0.0805 (nonflat case). Compared to the {Lambda}CDM model, it is found the current combined data sets do not favor the holographic dark energy model over the {Lambda}CDM model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Testing general relativity with current cosmological data

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large-scale structure and the deflection of light by that structure. We clarify the relations between several different model-independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. A Markov chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MDLoader: A Hybrid Model-Driven Data Loader for Distributed Graph Neural Network Training

Scalable data management is essential for processing large scientific dataset on HPC platforms for distributed deep learning. In-memory distributed storage is preferred for its speed, enabling rapid, random, and frequent data access required by stochastic optimizers. Processes use one-sided or collective communication to fetch remote data, with optimal performance depending on (i) dataset characteristics, (ii) training scale, and (iii) interconnection network. Empirical analysis shows collective communication excels with larger mini-batch sizes and/or fewer processes, whereas one-sided communication outperforms at larger scales. We propose MDLoader, a hybrid in-memory data loader for distributed graph neural network training. MDLoader features a model-driven performance estimator that dynamically selects between one-sided and collective communication at the beginning of training using Tree of Parzen Estimators (TPE). Evaluations on NERSC Perlmutter and OLCF Summit show MDLoader outperforms single-backend loaders by up to 2.83 × and predicts the suitable communication method with 96.3% (Perlmutter) and 94.3% (Summit) success rate.

Bae, Jonghyun↗

PERIOD ESTIMATION FOR SPARSELY SAMPLED QUASI-PERIODIC LIGHT CURVES APPLIED TO MIRAS

We develop a nonlinear semi-parametric Gaussian process model to estimate periods of Miras with sparsely sampled light curves. The model uses a sinusoidal basis for the periodic variation and a Gaussian process for the stochastic changes. We use maximum likelihood to estimate the period and the parameters of the Gaussian process, while integrating out the effects of other nuisance parameters in the model with respect to a suitable prior distribution obtained from earlier studies. Since the likelihood is highly multimodal for period, we implement a hybrid method that applies the quasi-Newton algorithm for Gaussian process parameters and search the period/frequency parameter space over a dense grid. A large-scale, high-fidelity simulation is conducted to mimic the sampling quality of Mira light curves obtained by the M33 Synoptic Stellar Survey. The simulated data set is publicly available and can serve as a testbed for future evaluation of different period estimation methods. The semi-parametric model outperforms an existing algorithm on this simulated test data set as measured by period recovery rate and quality of the resulting period–luminosity relations.

79 ASTRONOMY AND ASTROPHYSICS↗

BBGKY kinetic approach for an e{sup -}e{sup +}{gamma} plasma created from the vacuum in a strong laser-generated electric field: The one-photon annihilation channel

In the present work a closed system of kinetic equations is obtained from the truncation of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy for the description of the vacuum creation of an electron-positron plasma and secondary photons due to a strong laser field. This truncation is performed in the Markovian approximation for the one-photon annihilation channel which is accessible due to the presence of the strong external field. Estimates of the photon production rate are obtained for different domains of laser field parameters (frequency {nu} and field strength E). A huge quantity of optical photons of the quasiclassical laser field is necessary to satisfy the conservation laws of the energy and momentum of the constituents (e{sup -}, e{sup +} and {gamma}) in this channel. Since the number of these optical photons corresponds to the order of perturbation theory, a vanishingly small photon production rate results for the optical region and strongly subcritical fields E<<E{sub c}. In the {gamma}-ray region {nu} < or approx. m the required number of laser photons is small and the production rate of photons from the one-photon annihilation process becomes accessible to observations for subcritical fields E < or approx. E{sub c}. In the infrared region the photon distribution has a 1/k spectrum typical for flicker noise.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Algorithms and analyses for stochastic optimization for turbofan noise reduction using parallel reduced-order modeling

Simulation-based optimization of acoustic liner design in a turbofan engine nacelle for noise reduction purposes can dramatically reduce the cost and time needed for experimental designs. Because uncertainties are inevitable in the design process, a stochastic optimization algorithm is posed based on the conditional value-at-risk measure so that an ideal acoustic liner impedance is determined that is robust in the presence of uncertainties. A parallel reduced-order modeling framework is developed that dramatically improves the computational efficiency of the stochastic optimization solver for a realistic nacelle geometry. The reduced stochastic optimization solver takes less than 500 s to execute. Additionally, well-posedness and finite element error analyses of the state system and optimization problem are provided.

97 MATHEMATICS AND COMPUTING↗

GPU-BASED MONTE CARLO DUST RADIATIVE TRANSFER SCHEME APPLIED TO ACTIVE GALACTIC NUCLEI

A three-dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing-time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons. Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one- and two-dimensional benchmarks are reproduced at high precision. The parallelization capability of various MC algorithms is analyzed and included in our treatment. We utilize the Lucy algorithm for the optical thin case where the Poisson noise is high, the iteration-free Bjorkman and Wood method to reduce the calculation time, and the Fleck and Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGNs) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 {mu}m silicate features in absorption or emission is discussed. The SED of the radio-loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far-infrared emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Programmable entanglement oscillations in a non-Markovian channel

We suggest and demonstrate an all-optical experimental setup to observe and engineer entanglement oscillations of a pair of polarization qubits in an effective non-Markovian channel. We generate entangled photon pairs by spontaneous parametric down-conversion (SPDC), and then insert a programmable spatial light modulator in order to impose a polarization-dependent phase shift on the spatial domain of the SPDC output. This creates an effective programmable non-Markovian environment where modulation of the environment spectrum is obtained by inserting a spatial grating on the signal arm. In our experiment, programmable oscillations of entanglement are achieved, where the entangled state obtained at the maximum of the revival after sudden death violates Bell's inequality by 17 standard deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum algorithm for simulating the dynamics of an open quantum system

In the study of open quantum systems, one typically obtains the decoherence dynamics by solving a master equation. The master equation is derived using knowledge of some basic properties of the system, the environment, and their interaction: One basically needs to know the operators through which the system couples to the environment and the spectral density of the environment. For a large system, it could become prohibitively difficult to even write down the appropriate master equation, let alone solve it on a classical computer. In this paper, we present a quantum algorithm for simulating the dynamics of an open quantum system. On a quantum computer, the environment can be simulated using ancilla qubits with properly chosen single-qubit frequencies and with properly designed coupling to the system qubits. The parameters used in the simulation are easily derived from the parameters of the system + environment Hamiltonian. The algorithm is designed to simulate Markovian dynamics, but it can also be used to simulate non-Markovian dynamics provided that this dynamics can be obtained by embedding the system of interest into a larger system that obeys Markovian dynamics. We estimate the resource requirements for the algorithm. In particular, we show that for sufficiently slow decoherence a single ancilla qubit could be sufficient to represent the entire environment, in principle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cosmic microwave background with Brans-Dicke gravity. II. Constraints with the WMAP and SDSS data

Using the covariant formalism developed in a companion paper [F.-Q. Wu, L. E. Qiang, X. Wang, and X. Chen, preceding Article, Phys. Rev. D 82, 083002 (2010)] (paper I), we derive observational constraints on the Brans-Dicke model in a flat Friedmann-Lemaitre-Robertson-Walker universe with a cosmological constant and cold dark matter. The CMB observations we use include the Wilkinson Microwave Anisotropy Probe 5 yr data, the Arcminute Cosmology Bolometer Array Receiver 2007 data, the Cosmic Background Imager polarization data, and the Balloon Observations of Millimetric Extragalactic Radiation and Geophysics 2003 flight data. For the large scale structure we use the matter power spectrum data measured with the luminous red galaxy survey of the Sloan Digital Sky Survey Data Release 4. We parametrize the Brans-Dicke parameter {omega} with a new parameter {zeta}=ln(1/{omega}+1), and use the Markov-Chain Monte Carlo method to explore the parameter space. We find that using CMB data alone, one could place some constraints on positive {zeta} or {omega}, but negative {zeta} or {omega} could not be constrained effectively. However, with additional large scale structure data, one could break the degeneracy at {zeta}<0. The 2{sigma} (95.5%) bound on {zeta} is -0.008 37<{zeta}<0.010 18 (corresponding to {omega}<-120.0 or {omega}>97.8). We also obtained constraints on G/G, the rate of change of G at present, as -1.75x10{sup -12} yr{sup -1}<G/G<1.05x10{sup -12} yr{sup -1}, and {delta}G/G, the total variation of G since the epoch of recombination, as -0.083<{delta}G/G<0.095 at the 2{sigma} confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the Use of the Beta Distribution in Probabilistic Resource Assessments

The triangular distribution is a popular choice when it comes to modeling bounded continuous random variables. Its wide acceptance derives mostly from its simple analytic properties and the ease with which modelers can specify its three parameters through the extremes and the mode. On the negative side, hardly any real process follows a triangular distribution, which from the outset puts at a disadvantage any model employing triangular distributions. At a time when numerical techniques such as the Monte Carlo method are displacing analytic approaches in stochastic resource assessments, easy specification remains the most attractive characteristic of the triangular distribution. The beta distribution is another continuous distribution defined within a finite interval offering wider flexibility in style of variation, thus allowing consideration of models in which the random variables closely follow the observed or expected styles of variation. Despite its more complex definition, generation of values following a beta distribution is as straightforward as generating values following a triangular distribution, leaving the selection of parameters as the main impediment to practically considering beta distributions. This contribution intends to promote the acceptance of the beta distribution by explaining its properties and offering several suggestions to facilitate the specification of its two shape parameters. In general, given the same distributional parameters, use of the beta distributions in stochastic modeling may yield significantly different results, yet better estimates, than the triangular distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A stepped leader model for lightning including charge distribution in branched channels

The stepped leader process in negative cloud-to-ground lightning plays a vital role in lightning protection analysis. As lightning discharge usually presents significant branched or tortuous channels, the charge distribution along the branched channels and the stochastic feature of stepped leader propagation were investigated in this paper. The charge density along the leader channel and the charge in the leader tip for each lightning branch were approximated by introducing branch correlation coefficients. In combination with geometric characteristics of natural lightning discharge, a stochastic stepped leader propagation model was presented based on the fractal theory. By comparing simulation results with the statistics of natural lightning discharges, it was found that the fractal dimension of lightning trajectory in simulation was in the range of that observed in nature and the calculation results of electric field at ground level were in good agreement with the measurements of a negative flash, which shows the validity of this proposed model. Furthermore, a new equation to estimate the lightning striking distance to flat ground was suggested based on the present model. The striking distance obtained by this new equation is smaller than the value estimated by previous equations, which indicates that the traditional equations may somewhat overestimate the attractive effect of the ground.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Application of the MCMC Method for the Calibration of DSMC Parameters

A Markov Chain Monte Carlo (MCMC) algorithm was employed to obtain a calibrated distribution for the hard sphere diameter, the VHS reference diameter, and the temperature viscosity exponent of argon, for use in the Direct Simulation Monte Carlo (DSMC) method. Shock-tube experiments from Alsmeyer [1] were used to provide the necessary calibration data for use in the MCMC method. The DSMC method used in this work employs the algorithm of Bird [2], with modifications to allow for the efficient simulation of a 1D shock. When calibrating for the hard sphere diameter (the temperature viscosity exponent is set to 0.5 for the hard sphere method), the results of the MCMC method agree with a simple brute-force method, and a single value for the hard-sphere diameter is obtained. For the VHS method, however, when simultaneously calibrating the VHS reference diameter and the temperature viscosity exponent, we find that normalized density data alone does not provide sufficient information to obtain a single solution for both parameters. Instead we find a band in parameter space where acceptable solutions are obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modelling and analysis of turbulent datasets using Auto Regressive Moving Average processes

We introduce a novel way to extract information from turbulent datasets by applying an Auto Regressive Moving Average (ARMA) statistical analysis. Such analysis goes well beyond the analysis of the mean flow and of the fluctuations and links the behavior of the recorded time series to a discrete version of a stochastic differential equation which is able to describe the correlation structure in the dataset. We introduce a new index Υ that measures the difference between the resulting analysis and the Obukhov model of turbulence, the simplest stochastic model reproducing both Richardson law and the Kolmogorov spectrum. We test the method on datasets measured in a von Kármán swirling flow experiment. We found that the ARMA analysis is well correlated with spatial structures of the flow, and can discriminate between two different flows with comparable mean velocities, obtained by changing the forcing. Moreover, we show that the Υ is highest in regions where shear layer vortices are present, thereby establishing a link between deviations from the Kolmogorov model and coherent structures. These deviations are consistent with the ones observed by computing the Hurst exponents for the same time series. We show that some salient features of the analysis are preserved when considering global instead of local observables. Finally, we analyze flow configurations with multistability features where the ARMA technique is efficient in discriminating different stability branches of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗