Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PROCESS”

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 793 records · Page 44

Zero Temperature Limit for Directed Polymers and Inviscid Limit for Stationary Solutions of Stochastic Burgers Equation

We consider a space-continuous and time-discrete polymer model for positive temperature and the associated zero temperature model of last passage percolation type. In our previous work, we constructed and studied infinite-volume polymer measures and one-sided infinite minimizers for the associated variational principle, and used these objects for the study of global stationary solutions of the Burgers equation with positive or zero viscosity and random kick forcing, on the entire real line. In this paper, we prove that in the zero temperature limit, the infinite-volume polymer measures concentrate on the one-sided minimizers and that the associated global solutions of the viscous Burgers equation with random kick forcing converge to the global solutions of the inviscid equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of the Flows of Admixtures in a Random Layered Strip with Probable Arrangement of Inclusions Near the Boundaries of the Body

We study the random flow of admixtures in a two-phase stochastically inhomogeneous strip with the most probable arrangement of inclusions in the vicinity of the surfaces of the body. A mathematical model is formulated for the function of diffusion flow with nonzero constant initial concentration. A random diffusion flow is represented in the form of a Neumann series. The procedure of averaging of the random mass flow over the ensemble of phase configurations with arcsine distribution function is performed. The influence of the characteristics of the medium on the distribution of mass flow is analyzed. It is shown that if the diffusion coefficient of admixtures in the inclusion is higher than for the matrix, then the increase in the characteristic thickness of the layers causes a decrease in the value of the diffusion flow, whereas the mass flow in the entire body increases with the volume fraction of the inclusions.

97 MATHEMATICS AND COMPUTING↗

Optimal guidance law in quantum mechanics

Following de Broglie’s idea of a pilot wave, this paper treats quantum mechanics as a problem of stochastic optimal guidance law design. The guidance scenario considered in the quantum world is that an electron is the flight vehicle to be guided and its accompanying pilot wave is the guidance law to be designed so as to guide the electron to a random target driven by the Wiener process, while minimizing a cost-to-go function. After solving the stochastic optimal guidance problem by differential dynamic programming, we point out that the optimal pilot wave guiding the particle’s motion is just the wavefunction Ψ(t,x), a solution to the Schrödinger equation; meanwhile, the closed-loop guidance system forms a complex state–space dynamics for Ψ(t,x), from which quantum operators emerge naturally. Quantum trajectories under the action of the optimal guidance law are solved and their statistical distribution is shown to coincide with the prediction of the probability density function Ψ{sup ∗}Ψ. -- Highlights: •Treating quantum mechanics as a pursuit-evasion game. •Reveal an interesting analogy between guided flight motion and guided quantum motion. •Solve optimal quantum guidance problem by dynamic programming. •Gives a formal proof of de Broglie–Bohm’s idea of a pilot wave. •The optimal pilot wave is shown to be a wavefunction solved from Schrödinger equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A semi-intrusive deterministic approach to uncertainty quantification in non-linear fluid flow problems

This paper deals with the formulation of a semi-intrusive (SI) method allowing the computation of statistics of linear and non linear PDEs solutions. This method shows to be very efficient to deal with probability density function of whatsoever form, long-term integration and discontinuities in stochastic space. Given a stochastic PDE where randomness is defined on Ω, starting from (i) a description of the solution in term of a space variables, (ii) a numerical scheme defined for any event ω∈Ω and (iii) a (family) of random variables that may be correlated, the solution is numerically described by its conditional expectancies of point values or cell averages and its evaluation constructed from the deterministic scheme. One of the tools is a tessellation of the random space as in finite volume methods for the space variables. Then, using these conditional expectancies and the geometrical description of the tessellation, a piecewise polynomial approximation in the random variables is computed using a reconstruction method that is standard for high order finite volume space, except that the measure is no longer the standard Lebesgue measure but the probability measure. This reconstruction is then used to formulate a scheme on the numerical approximation of the solution from the deterministic scheme. This new approach is said semi-intrusive because it requires only a limited amount of modification in a deterministic solver to quantify uncertainty on the state when the solver includes uncertain variables. The effectiveness of this method is illustrated for a modified version of Kraichnan–Orszag three-mode problem where a discontinuous pdf is associated to the stochastic variable, and for a nozzle flow with shocks. The results have been analyzed in terms of accuracy and probability measure flexibility. Finally, the importance of the probabilistic reconstruction in the stochastic space is shown up on an example where the exact solution is computable, the viscous Burgers equation.

97 MATHEMATICS AND COMPUTING↗

Reweighted ℓ{sub 1} minimization method for stochastic elliptic differential equations

We consider elliptic stochastic partial differential equations (SPDEs) with random coefficients and solve them by expanding the solution using generalized polynomial chaos (gPC). Under some mild conditions on the coefficients, the solution is “sparse” in the random space, i.e., only a small number of gPC basis makes considerable contribution to the solution. To exploit this sparsity, we employ reweighted l{sub 1} minimization to recover the coefficients of the gPC expansion. We also combine this method with random sampling points based on the Chebyshev probability measure to further increase the accuracy of the recovery of the gPC coefficients. We first present a one-dimensional test to demonstrate the main idea, and then we consider 14 and 40 dimensional elliptic SPDEs to demonstrate the significant improvement of this method over the standard l{sub 1} minimization method. For moderately high dimensional (∼10) problems, the combination of Chebyshev measure with reweighted l{sub 1} minimization performs well while for higher dimensional problems, reweighted l{sub 1} only is sufficient. The proposed approach is especially suitable for problems for which the deterministic solver is very expensive since it reuses the sampling results and exploits all the information available from limited sources.

97 MATHEMATICS AND COMPUTING↗

A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media

A new computational method is proposed for efficient uncertainty quantification of multiphase flow in porous media with stochastic permeability. For pressure estimation, it combines the dimension-adaptive stochastic collocation method on Smolyak sparse grids and the Kronrod–Patterson–Hermite nested quadrature formulas. For saturation estimation, an additional stage is developed, in which the pressure and velocity samples are first generated by the sparse grid interpolation and then substituted into the transport equation to solve for the saturation samples, to address the low regularity problem of the saturation. Numerical examples are presented for multiphase flow with stochastic permeability fields to demonstrate accuracy and efficiency of the proposed two-stage adaptive stochastic collocation method on nested sparse grids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Segmentation, statistical analysis, and modelling of the wall system in ceramic foams

Closed walls in otherwise open foam structures may have a great impact on macroscopic properties of the materials. In this paper, we present two algorithms for the segmentation of such closed walls from micro-computed tomography images of the foam structure. The techniques are compared on simulated data and applied to tomographic images of ceramic filters. This allows for a detailed statistical analysis of the normal directions and sizes of the walls. Finally, we explain how the information derived from the segmented wall system can be included in a stochastic microstructure model for the foam.

36 MATERIALS SCIENCE↗

T24 HRAS transformed NIH/3T3 mouse cells (GhrasT-NIH/3T3) in serial tumorigenic in vitro/in vivo passages give rise to increasingly aggressive tumorigenic cell lines T1-A and T2-A and metastatic cell lines T3-HA and T4-PA

Cancer cells often arise progressively from “normal” to “pre-cancer” to “transformed” to “local metastasis” to “metastatic disease” to “aggressive metastatic disease”. Recent whole genome sequencing (WGS) and spectral karyotyping (SKY) of cancer cells and tumorigenic models have shown this progression involves three major types of genome rearrangements: ordered small step-wise changes, more dramatic “punctuated evolution” (chromoplexy), and large catastrophic steps (chromothripsis) which all occur in random combinations to generate near infinite numbers of stochastically rearranged metastatic cancer cell genomes. This paper describes a series of mouse cell lines developed sequentially to mimic this type of progression. This starts with the new GhrasT-NIH/Swiss cell line that was produced from the NIH/3T3 cell line that had been transformed by transfection with HRAS oncogene DNA from the T24 human bladder carcinoma. These GhrasT-NIH/Swiss cells were injected s.c. into NIH/Swiss mice to produce primary tumors from which one was used to establish the T1-A cell line. T1-A cells injected i.v. into the tail vein of a NIH/Swiss mouse produced a local metastatic tumor near the base of the tail from which the T2-A cell line was established. T2-A cells injected i.v. into the tail vein of a nude NIH/Swiss mouse produced metastases in the liver and one lung from which the T3-HA (H=hepatic) and T3-PA (P=pulmonary) cell lines were developed, respectively. T3-HA cells injected i.v. into a nude mouse produced a metastasis in the lung from which the T4-PA cell line was established. PCR analysis indicated the human T24 HRAS oncogene was carried along with each in vitro/in vivo transfer step and found in the T2-A and T4-PA cell lines. Light photomicrographs indicate that all transformed cells are morphologically similar. GhrasT-NIH/Swiss cells injected s.c. produced tumors in 4% of NIH/Swiss mice in 6–10 weeks; T1-A cells injected s.c. produced tumors in 100% of NIH/Swiss mice in 7–10 days. T1-A, T-2A, T3-HA and T4-PA cells when injected i.v. into the tail produced local metastasis in non-nude or nude NIH/Swiss mice. T4-PA cells were more widely metastatic than T3-HA cells when injected i.v. into nude mice. Evaluation of the injected mice indicated a general increase in metastatic potential of each cell line in the progression as compared to the GhrasT-NIH/3T3 transformed cells. A new photomicrographic technique to follow growth rates within six preselected 2×2 mm{sup 2} grids per plate is described. Average doubling times of the transformed cells GhrasT-NIH/3T3 (17 h), T1A (17.5 h), T2A (15.5 h), T3-HA (17.5 h) and T4-PA (18.5 h) (average 17.2 h) were significantly faster (P=0.006) than NIH Swiss primary embryonic cells and NIH/3T3 cells (22 h each). This cell series is currently used in this lab for studies of cancer cell inhibitors, mitochondrial biogenesis and gene expression and is available for further study by other investigators for intra- and inter-laboratory comparisons of WGS, transcriptome sequencing, SKY and other analyses. The genome rearrangements in these cells together with their phenotypic properties may help provide more insights into how one tumorigenic progression occurred to produce the various cell lines that led to the highly metastatic T4-PA cell line. - Highlights: • A tumor progression model derived from a new T24 HRAS transformed NIH/3T3 cell line. • T1-A, T2-A, T3-HA, T4PA lines derived from GhrasT-NIH/3T3, T1-A, T2-A, T3-HA tumors. • PCR detected the human HRAS oncogene in T2-A and T4-PA metastatic tumors. • These five new cell lines showed a progressive increase in tumorigenic potential. • New method obtained doubling times using daily photomicrographs of gridded plates.

60 APPLIED LIFE SCIENCES↗

Doppler-resolved kinetics of saturation recovery

Frequency modulated laser transient absorption has been used to monitor the ground state rotational energy transfer rates of CN radicals in a double-resonance, depletion recovery experiment. When a pulsed laser is used to burn a hole in the equilibrium ground state population of one rotational state without velocity selection, the population recovery rate is found to depend strongly on the Doppler detuning of a narrow-band probe laser. Similar effects should be apparent for any relaxation rate process that competes effectively with velocity randomization. Alternative methods of extracting thermal rate constants in the presence of these non-thermal conditions are evaluated. Total recovery rate constants, analogous to total removal rate constants in an experiment preparing a single initial rotational level, are in good agreement with quantum scattering calculations, but are slower than previously reported experiments and show qualitatively different rotational state dependence between Ar and He collision partners. As a result, quasi-classical trajectory studies confirm that the differing rotational state dependence is primarily a kinematic effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The structural origin of the hard-sphere glass transition in granular packing

Glass transition is accompanied by a rapid growth of the structural relaxation time and a concomitant decrease of configurational entropy. It remains unclear whether the transition has a thermodynamic origin, and whether the dynamic arrest is associated with the growth of a certain static order. Using granular packing as a model hard-sphere glass, we show the glass transition as a thermodynamic phase transition with a ‘hidden’ polytetrahedral order. This polytetrahedral order is spatially correlated with the slow dynamics. It is geometrically frustrated and has a peculiar fractal dimension. Additionally, as the packing fraction increases, its growth follows an entropy-driven nucleation process, similar to that of the random first-order transition theory. In conclusion, our study essentially identifies a long-sought-after structural glass order in hard-sphere glasses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The spectral expansion of the elasticity random field

We consider a deformable body that occupies a region D in the plane. In our model, the body’s elasticity tensor H(x) is the restriction to D of a second-order mean-square continuous random field. Under translation, the expected value and the correlation tensor of the field H(x) do not change. Under action of an arbitrary element k of the orthogonal group O(2), they transform according to the reducible orthogonal representation k ⟼ S{sup 2}(S{sup 2}(k)) of the above group. We find the spectral expansion of the correlation tensor R(x) of the elasticity field as well as the expansion of the field itself in terms of stochastic integrals with respect to a family of orthogonal scattered random measures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sublinear scaling for time-dependent stochastic density functional theory

A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (≈16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From laser particle acceleration to the synthesis of extremely neutron rich isotopes via the novel fission-fusion mechanism

High-power, short pulse lasers have emerged in the last decade as attractive tools for accelerating charged particles (electrons, ions) to high energies over mm-scale acceleration lengths, thus promising to rival conventional acceleration techniques in the years ahead. In the first part of the article, the principles of laser-plasma interaction as well as the techniques and the current status of the acceleration of electron and ion beams will be briefly introduced. In particular with the upcoming next generation of multi-PW class laser systems, such as the one under construction for the ELI-Nuclear Physics project in Bucharest (ELI-NP), very efficient acceleration mechanisms for brilliant ion beams like radiation pressure acceleration (RPA) come into reach. Here, ultra-dense ion beams reaching solid-state density can be accelerated from thin target foils, exceeding the density of conventionally accelerated ion beams by about 14 orders of magnitude. This unique property of laser-accelerated ion beams can be exploited to explore the scenario of a new reaction mechanism called ‘fission-fusion’, which will be introduced in the second part of the article. Accelerating fissile species (e.g. {sup 232}Th) towards a second layer of the same material will lead to fission both of the beam-like and target-like particles. Due to the close to solid-state density of the accelerated ion bunches, fusion may occur between neutron-rich (light) fission products. This may open an access path towards extremely neutron-rich nuclides in the vicinity of the N=126 waiting point of the astrophysical r process. ‘Waiting points’ at closed nucleon shells play a crucial role in controlling the reaction rates. However, since most of the pathway of heavy-element formation via the rapid-neutron capture process (r-process) runs in ‘terra incognita’ of the nuclear landscape, in particular the waiting point at N=126 is yet unexplored and will remain largely inaccessible to conventional nuclear reaction schemes even at next-generation radioactive beam facilities, underlining the attractive perspectives offered, e.g., by ELI-NP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A multilevel stochastic collocation method for SPDEs

We present a multilevel stochastic collocation method that, as do multilevel Monte Carlo methods, uses a hierarchy of spatial approximations to reduce the overall computational complexity when solving partial differential equations with random inputs. For approximation in parameter space, a hierarchy of multi-dimensional interpolants of increasing fidelity are used. Rigorous convergence and computational cost estimates for the new multilevel stochastic collocation method are derived and used to demonstrate its advantages compared to standard single-level stochastic collocation approximations as well as multilevel Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical analysis of randomly forced glycolitic oscillations

Randomly forced glycolytic oscillations in Higgins model are studied both numerically and analytically. Numerical analysis is based on the direct simulation of the solutions of stochastic system. Non-uniformity of the stochastic bundle along the deterministic cycle is shown. For the analytical investigation of the randomly forced Higgins model, the stochastic sensitivity function technique and confidence domains method are applied. Results of the influence of additive noise on the cycle of this model are given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical analysis of thermally assisted spin-transfer torque magnetization reversal in synthetic ferrimagnetic free layers

The spin transfer torque magnetization reversal of synthetic ferrimagnetic free layers under pulsed temperature rise was numerically studied by solving the Landau–Lifshitz–Gilbert equation, taking into account the stochastic random fields, the temperature dependence of magnetic parameters, and the spin torque terms. The anti-parallel magnetization configuration was retained at the elevated temperature, due to interlayer dipole coupling. A significant thermal assistance effect, resulting in a 40% reduction in the switching current, was demonstrated during a nanosecond pulsed temperature rise up to 77% of the Curie temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theoretically informed Monte Carlo simulation of liquid crystals by sampling of alignment-tensor fields

A theoretically informed coarse-grained Monte Carlo method is proposed for studying liquid crystals. The free energy functional of the system is described in the framework of the Landau-de Gennes formalism. The alignment field and its gradients are approximated by finite differences, and the free energy is minimized through a stochastic sampling technique. The validity of the proposed method is established by comparing the results of the proposed approach to those of traditional free energy minimization techniques. Its usefulness is illustrated in the context of three systems, namely, a nematic liquid crystal confined in a slit channel, a nematic liquid crystal droplet, and a chiral liquid crystal in the bulk. It is found that for systems that exhibit multiple metastable morphologies, the proposed Monte Carlo method is generally able to identify lower free energy states that are often missed by traditional approaches. Importantly, the Monte Carlo method identifies such states from random initial configurations, thereby obviating the need for educated initial guesses that can be difficult to formulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of stochastically forced quasi-periodic attractors

A problem of the analysis of stochastically forced quasi-periodic auto-oscillations of nonlinear dynamic systems is considered. A stationary distribution of random trajectories in the neighborhood of the corresponding deterministic attractor (torus) is studied. A parametric description of quadratic approximation of the quasipotential based on the stochastic sensitivity functions (SSF) technique is given. Using this technique, we analyse a dispersion of stochastic flows near the torus. For the case of two-torus in three-dimensional space, the stochastic sensitivity function is constructed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗