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

Approach to the Steady State in Kinetic Models with Thermal Reservoirs at Different Temperatures

We continue the investigation of kinetic models of a system in contact via stochastic interactions with several spatially homogeneous thermal reservoirs at different temperatures. Considering models different from those investigated in Carlen et al. (Braz J Probab Stat 29:372–386, 2015), we explicitly compute the unique spatially uniform non-equilibrium steady state (NESS) and prove that it is approached exponentially fast from any uniform initial state. This leaves open the question of whether there exist NESS that are not spatially uniform. Making a further simplification of our models, we then prove non-existence of such NESS and exponential approach to the unique spatially uniform NESS (with a computably boundable rate). The method of proof relies on refined Doeblin estimates and other probabilistic techniques, and is quite different form the analysis in Carlen et al. (Braz J Probab Stat 29:372–386, 2015) that was based on contraction mapping methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical Properties for Flows with Unbounded Roof Function, Including the Lorenz Attractor

For geometric Lorenz attractors (including the classical Lorenz attractor) we obtain a greatly simplified proof of the central limit theorem which applies also to the more general class of codimension two singular hyperbolic attractors. We also obtain the functional central limit theorem and moment estimates, as well as iterated versions of these results. A consequence is deterministic homogenisation (convergence to a stochastic differential equation) for fast-slow dynamical systems whenever the fast dynamics is singularly hyperbolic of codimension two.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

On the Validity of Linear Response Theory in High-Dimensional Deterministic Dynamical Systems

This theoretical work considers the following conundrum: linear response theory is successfully used by scientists in numerous fields, but mathematicians have shown that typical low-dimensional dynamical systems violate the theory’s assumptions. Here we provide a proof of concept for the validity of linear response theory in high-dimensional deterministic systems for large-scale observables. We introduce an exemplary model in which observables of resolved degrees of freedom are weakly coupled to a large, inhomogeneous collection of unresolved chaotic degrees of freedom. By employing statistical limit laws we give conditions under which such systems obey linear response theory even if all the degrees of freedom individually violate linear response. We corroborate our result with numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Most Effective Model for Describing the Universal Behavior of a Noisy Kuramoto–Sivashinsky Equation as a Paradigmatic Model

We study a noisy Kuramoto–Sivashinsky (KS) equation which describes unstable surface growth and chemical turbulence. It has been conjectured that the universal long-wavelength behavior of the equation, which is characterized by scale-dependent parameters, is described by a Kardar–Parisi–Zhang (KPZ) equation. We consider this conjecture by analyzing a renormalization-group equation for a class of generalized KPZ equations. We then uniquely determine the parameter values of the KPZ equation that most effectively describes the universal long-wavelength behavior of the noisy KS equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Filtration of stationary Gaussian statistical experiments

The filtration of stationary Gaussian statistical experiments is determined by a solution of the equation of optimum filtration, which is characterized by the two-dimensional matrix of covariances. The parameters of a filtered signal are set by empiric covariances.

97 MATHEMATICS AND COMPUTING↗

Characteristic Functions and Compactness of Distributions of Sums of Independent Random Variables

The sequences of distributions of centered sums of independent random variables are considered within the framework of the series scheme, without assuming the classical conditions for uniform asymptotic smallness and uniform limit constancy. Necessary and sufficient conditions are obtained for relative and stochastic compactness of such sequences in terms of the characteristic functions of summable random variables and with using their τ-centers.

97 MATHEMATICS AND COMPUTING↗

Numerical and Analytical Methods of Study of Stochastic Systems with Delay

This paper is a brief review of methods of qualitative and quantitative study of various classes of stochastic systems with aftereffect. We describe schemes of analysis of stochastic ordinary and partial differential and integrodifferential equations, which are also applicable for their deterministic analogs. These numerical and analytical schemes are implemented in the software package Mathematica and in the language Intel Fortran.

97 MATHEMATICS AND COMPUTING↗

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↗

Spread of Values of a Cantor-Type Fractal Continuous Nonmonotone Function

By using the Q{sub 5}{sup ∗}-representation of numbers[0,1]∋x=β{sub α{sub 1(x)1}}+∑{sub k=2}{sup ∞}(β{sub α{sub k(x)k}}∏{sub j=1}{sup k−1}q{sub α{sub j(x)j}})=Δ{sub α{sub 1(x)α{sub 2(x)…α{sub k(x)…}{sup Q{sub 5}{sup ∗}}}}} determined by the quinary alphabet A{sub 5} ≡ {0, 1, 2, 3, 4} and an infinite stochastic matrix ‖q{sub ik}‖, i ∈ A{sub 5}, k ∈ N, with positive elements (q{sub 0k} + q{sub 1k} + q{sub 2k} + q{sub 3k} + q{sub 4k} = 1) such that ∏{sub k=1}{sup ∞}max{sub i}{q_i_k}=0 and β{sub 0k} = 0, β{sub i + 1, k} = β{sub ik} + q{sub ik}, i= (0,4)-bar , we define a continuous Cantor-type function by the equalityf(Δ{sub α{sub 1…α{sub k…}{sup Q{sub 5}{sup ∗}}}})=δ{sub α{sub 1(x)1}}+∑{sub k=2}{sup ∞}(δ{sub α{sub k(x)k}}∑{sub j=1}{sup k−1}g{sub α{sub j(x)j}})≡Δ{sub α{sub 1(x)…α{sub k(x)…}{sup G}}}, where δ{sub 0n} = 0, δ{sub 1n}=((2+ε{sub n})/4), δ{sub 2n}=(2/4)=δ{sub 3n}, and δ{sub 4n}=((2−ε{sub n})/4), i.e., δ{sub i + 1, n} = δ{sub in} + g{sub in}, n ∈ N, and (ε{sub n}) is a given sequence of real numbers such that 0 ≤ ε{sub n} ≤ 1. We prove that this function is well defined and continuous. Moreover, it does not have intervals of monotonicity, except the intervals where it is constant. A criterion of bounded variation of the function is also established. We are especially interested in the problem of level sets of the function and in the topological and metric properties of the images of Cantor-type sets.

97 MATHEMATICS AND COMPUTING↗

Stochastic Perturbations of Stable Dynamical Systems: Trajectory-Wise Approach

We study stochastic perturbations of a dynamical system with a locally stable fixed point. The perturbed system has the form of Ito stochastic differential equations. We assume that perturbations do not vanish at the equilibrium of the deterministic system. Using the approach based on consideration of trajectories to the analysis of stochastic differential equations, we find restrictions for perturbations under which the stability of the equilibrium is preserved with probability 1.

97 MATHEMATICS AND COMPUTING↗

A New Stochastic Modeling of 3-D Mud Drapes Inside Point Bar Sands in Meandering River Deposits

The environment of major sediments of eastern China oilfields is a meandering river where mud drapes inside point bar sand occur and are recognized as important factors for underground fluid flow and distribution of the remaining oil. The present detailed architectural analysis, and the related mud drapes' modeling inside a point bar, is practical work to enhance oil recovery. This paper illustrates a new stochastic modeling of mud drapes inside point bars. The method is a hierarchical strategy and composed of three nested steps. Firstly, the model of meandering channel bodies is established using the Fluvsim method. Each channel centerline obtained from the Fluvsim is preserved for the next simulation. Secondly, the curvature ratios of each meandering river at various positions are calculated to determine the occurrence of each point bar. The abandoned channel is used to characterize the geometry of each defined point bar. Finally, mud drapes inside each point bar are predicted through random sampling of various parameters, such as number, horizontal intervals, dip angle, and extended distance of mud drapes. A dataset, collected from a reservoir in the Shengli oilfield of China, was used to illustrate the mud drapes' building procedure proposed in this paper. The results show that the inner architectural elements of the meandering river are depicted fairly well in the model. More importantly, the high prediction precision from the cross validation of five drilled wells shows the practical value and significance of the proposed method.

42 ENGINEERING↗

Reservoir Modeling by Data Integration via Intermediate Spaces and Artificial Intelligence Tools in MPS Simulation Frameworks

Conditioning stochastic simulations are very important in many geostatistical applications that call for the introduction of nonlinear and multiple-point data in reservoir modeling. Here, a new methodology is proposed for the incorporation of different data types into multiple-point statistics (MPS) simulation frameworks. Unlike the previous techniques that call for an approximate forward model (filter) for integration of secondary data into geologically constructed models, the proposed approach develops an intermediate space where all the primary and secondary data are easily mapped onto. Definition of the intermediate space, as may be achieved via application of artificial intelligence tools like neural networks and fuzzy inference systems, eliminates the need for using filters as in previous techniques. The applicability of the proposed approach in conditioning MPS simulations to static and geologic data is verified by modeling a real example of discrete fracture networks using conventional well-log data. The training patterns are well reproduced in the realizations, while the model is also consistent with the map of secondary data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development Optimization and Uncertainty Analysis Methods for Oil and Gas Reservoirs

Uncertainty complicates the development optimization of oil and gas exploration and production projects, but methods have been devised to analyze uncertainty and its impact on optimal decision-making. This paper compares two methods for development optimization and uncertainty analysis: Monte Carlo (MC) simulation and stochastic programming. Two example problems for a gas field development and an oilfield development are solved and discussed to elaborate the advantages and disadvantages of each method. Development optimization involves decisions regarding the configuration of initial capital investment and subsequent operational decisions. Uncertainty analysis involves the quantification of the impact of uncertain parameters on the optimum design concept. The gas field development problem is designed to highlight the differences in the implementation of the two methods and to show that both methods yield the exact same optimum design. The results show that both MC optimization and stochastic programming provide unique benefits, and that the choice of method depends on the goal of the analysis. While the MC method generates more useful information, along with the optimum design configuration, the stochastic programming method is more computationally efficient in determining the optimal solution. Reservoirs comprise multiple compartments and layers with multiphase flow of oil, water, and gas. We present a workflow for development optimization under uncertainty for these reservoirs, and solve an example on the design optimization of a multicompartment, multilayer oilfield development.

02 PETROLEUM↗

Bayesian and LASSO Regressions for Comparative Permeability Modeling of Sandstone Reservoirs

In this research, Bayesian model averaging (BMA) and least absolute shrinkage and selection operator regression (LASSO) algorithms were adopted for the core permeability modeling as a function of well log and core measurements. More specifically, the core permeability (dependent factor) was modeled given the well log and core data (independent variables) and then predicted in non-cored intervals of a well in a sandstone formation. The BMA is a stochastic linear modeling and a Bayesian parameter selection. Among 50 linear models generated as a function of the independent variables and depicted in Occam’s window, the best model of the highest posterior probability is determined. The best model has the optimal subset parameters that most influences the response factor. The core permeability modeling was again conducted using the LASSO algorithm, which adopts a completely different way of modeling and subset selection by using the penalized least-squared equation. Both BMA and LASSO resulted in a very accurate prediction of the core permeability by achieving perfect measured and calculated permeability matching. Matches for the two algorithms were quantified by computing the adjusted R-squared and root-mean-squared prediction error. In addition, results of the BMA and LASSO were compared to the conventional multiple linear regression (LM). The LASSO algorithm led to an accurate matching similar to LM, but slightly better than BMA. Therefore, both BMA and LASSO represent integrated procedures for accurate permeability estimation similar to conventional regression analysis.

58 GEOSCIENCES↗

Pilot Point Optimization of Mining Boundaries for Lateritic Metal Deposits: Finding the Trade-off Between Dilution and Ore Loss

Geological contacts in lateritic metal deposits (footwall topographies) often delineate the orebody boundaries. Spatial variations seen in such contacts are frequently higher than those for the metal grades of the deposit. Therefore, borehole spacing chosen based on the grade variations cannot adequately capture the geological contact variability. Consequently, models created using such boreholes cause high volumetric uncertainties in the actual and targeted ore extraction volumes, which, in turn, lead to high unplanned dilution and ore losses. In this paper, a method to design optimum ore/mining boundaries for lateritic metal deposits is presented. The proposed approach minimizes the dilution/ore losses and comprises two main steps. First, the uncertainty on the orebody boundary is represented using a set of stochastic realizations generated with a multiple-point statistics algorithm. Then, the optimal orebody boundary is determined using an optimization technique inspired by a model calibration method called Pilot Points. The pilot points represent synthetic elevation values, and they are used to construct smooth mining boundaries using the multilevel B-spline technique. The performance of a generated surface is evaluated using the expected sum of losses in each of the stochastic realizations. The simulated annealing algorithm is used to iteratively determine the pilot point values which minimize the expected losses. The results show a significant reduction in the dilution volume as compared to those obtained from the actual mining operation.

58 GEOSCIENCES↗

Capturing Hidden Geochemical Anomalies in Scarce Data by Fractal Analysis and Stochastic Modeling

Fractal/multifractal modeling is a widely used geomathematical approach to capturing different populations in geochemical mapping. The rationale of this methodology is based on empirical frequency density functions attained from global or local distributions. This approach is quite popular because of its simplicity and versatility; it accounts for the frequency and spatial distribution of geochemical data considering self-similarity across a range of scales. Using this technique for detection of geochemical anomalies in scarce data, however, is problematic and can lead to systematic bias in the characterization of the underlying populations. In this paper, an innovative technique is presented that provides good results without a priori assumptions. A simulation approach is adopted for fractal analysis by generating different possible distribution scenarios for the variable under study to reveal the underlying populations that are frequently hidden due to lack of data. The proposed technique is called the global simulated size–number method, and it is validated in a case study with two synthetic datasets and another case study with real dataset from the Ushtagan gold deposit in northeast Kazakhstan.

58 GEOSCIENCES↗