A translation model for non-stationary, non-Gaussian random processes

A model for simulation of non-stationary, non-Gaussian processes based on non-linear translation of Gaussian random vectors is presented. This method is a generalization of traditional translation processes that includes the capability of simulating samples with spatially or temporally varying marginal probability density functions. A formal development of the properties of the resulting process includes joint probability density function, correlation distortion and lower and upper bounds that depend on the target marginal distributions. Examples indicate the possibility of exactly matching a wide range of marginal pdfs and second order moments through a simple interpolating algorithm. Furthermore, the application of the method in simulating statistically inhomogeneous random media is investigated, using the specific case of binary translation with stationary and non-stationary target correlations.     

A Monte Carlo simulation model for stationary non-Gaussian processes

A class of stationary non-Gaussian processes, referred to as the class of mixtures of translation processes, is defined by their finite dimensional distributions consisting of mixtures of finite dimensional distributions of translation processes. The class of mixtures of translation processes includes translation processes and is useful for both Monte Carlo simulation and analytical studies. As for translation processes, the mixture of translation processes can have a wide range of marginal distributions and correlation functions. Moreover, these processes can match a broader range of second order correlation functions than translation processes. The paper also develops an algorithm for generating samples of any non-Gaussian process in the class of mixtures of translation processes. The algorithm is based on the sampling representation theorem for stochastic processes and properties of the conditional distributions. Examples are presented to illustrate the proposed Monte Carlo algorithm and compare features of translation processes and mixture of translation processes.     

Translation models revisited

Translation models have been defined as memoryless mappings of Gaussian elements which match exactly/approximately target marginal distributions/correlations. We extend this class of translation models to include memoryless mappings of non-Gaussian elements. It is shown that quantities of interest inferred from equivalent translation models, i.e., models which share the same marginal distributions and have similar second moments, can differ significantly. It is suggested to construct families of equivalent translation models and select members of these families which are optimal for given quantities of interest.     

Model selection for a class of stochastic processes or random fields with bounded range

Methods are developed for finding an optimal model for a non-Gaussian stationary stochastic process or homogeneous random field under limited information. The available information consists of: (i) one or more finite length samples of the process or field; and (ii) knowledge that the process or field takes values in a bounded interval of the real line whose ends may or may not be known. The methods are developed and applied to the special case of non-Gaussian processes or fields belonging to the class of beta translation processes. Beta translation processes provide a flexible model for representing physical phenomena taking values in a bounded range, and are therefore useful for many applications. Numerical examples are presented to illustrate the utility of beta translation processes and the proposed methods for model selection.     

Extreme value distributions for nonlinear transformations of vector Gaussian processes

Approximations are developed for the marginal and joint probability distributions for the extreme values, associated with a vector of non-Gaussian random processes. The component non-Gaussian processes are obtained as nonlinear transformations of a vector of stationary, mutually correlated, Gaussian random processes and are thus, mutually dependent. The multivariate counting process, associated with the number of level crossings by the component non-Gaussian processes, is modelled as a multivariate Poisson point process. An analytical formulation is developed for determining the parameters of the multivariate Poisson process. This, in turn, leads to the joint probability distribution of the extreme values of the non-Gaussian processes, over a given time duration. For problems not amenable for analytical solutions, an algorithm is developed to determine these parameters numerically. The proposed extreme value distributions have applications in time-variant reliability analysis of randomly vibrating structural systems. The method is illustrated through three numerical examples and their accuracy is examined with respect to estimates from full scale Monte Carlo simulations of vector non-Gaussian processes.     

Non-Gaussian random processes in ocean engineering

This paper presents a state-of-the-art review on stochastic analysis and probabilistic prediction of non-Gaussian random processes in ocean engineering. The derivation of probability density functions which constitute the basis for stochastic analysis of non-Gaussian processes is discussed in detail, and then the probability distributions of peaks and troughs of non-Gaussian random process is discussed to provide information necessary for engineering design. As an example of application of these probability distribution functions, the procedure for predicting responses of an offshore structure which has substantial non-linear characteristics in random seas is presented.     

Modeling of atmospheric temperature fluctuations by translations of oscillatory random processes with application to spacecraft atmospheric re-entry

The presence of random fluctuations of air temperature within the Earth’s atmosphere is a well-documented phenomenon. During the past seventy years there have been numerous experimental efforts to accurately measure air temperature as a function of altitude and, through careful data analysis, provide statistics describing these fluctuations and the associated fluctuations in temperature gradients. In addition, several researchers suggest the presence of atmospheric layers or “sheets” where the statistics describing fluctuations in air temperature can vary significantly from layer to layer. Herein, we propose a model to represent fluctuations of air temperature within a layered atmosphere. The model is a special type of inhomogeneous non-Gaussian differentiable random process and can be calibrated to available data on the marginal statistics and spectral content of the fluctuating temperature field, as well as the associated first derivative of the process representing fluctuations in temperature gradients. Properties of the proposed model are presented, and statistical realizations of the fluctuating temperature field and its gradient are computed and presented for illustration. The random vibration response of a spacecraft falling to Earth through these fluctuating conditions is then considered to demonstrate the usefulness of the proposed model.     

Efficient stationary multivariate non-Gaussian simulation based on a Hermite PDF model

An efficient stationary multivariate non-Gaussian simulation method is developed using spectral representation and third order Hermite polynomial translation. An approximate closed form relationship is employed to identify the Hermite translation parameters based on target skewness and kurtosis. This preserves a high degree of accuracy over the entire admissible range of the Hermite translation, and eliminates the need for iterative solution of the translation parameters. The Hermite PDF model is suitable for a wide range of strongly non-Gaussian stochastic process. In addition, an explicit bidirectional relationship between the target non-Gaussian and Gaussian correlation is developed to eliminate the need for iteration or numerical integration to identify the underlying Gaussian correlation. Examples apply the simulation method to both theoretical targets and experimental wind pressure data.     

Simulation of strongly non-Gaussian processes using Karhunen–Loeve expansion

The non-Gaussian Karhunen–Loeve (K–L) expansion is very attractive because it can be extended readily to non-stationary and multi-dimensional fields in a unified way. However, for strongly non-Gaussian processes, the original procedure is unable to match the distribution tails well. This paper proposes an effective solution to this tail mismatch problem using a modified orthogonalization technique that reduces the degree of shuffling within columns containing empirical realizations of the K–L random variables. Numerical examples demonstrate that the present algorithm is capable of matching highly non-Gaussian marginal distributions and stationary/non-stationary covariance functions simultaneously to a very accurate degree. The ability to converge correctly to an abrupt lower bound in the target marginal distributions studied is noteworthy.     

An efficient simulation algorithm for non-Gaussian nonstationary processes

A novel algorithm is proposed for simulating univariate non-Gaussian nonstationary processes (NNP) with the specified evolutionary power spectral density (EPSD)/nonstationary auto-correlation function (NACF) and first four-order time-varying marginal moments (TVMMs). The sample realizations of the target NNP are generated as the outputs from a specific time-varying auto-regressive (TVAR) model via filtering the non-Gaussian and nonstationary white noise inputs. These white noise inputs are also non-Gaussian and nonstationary, and their first four-order TVMMs are predetermined using an approach developed herein according to the specified EPSD/NACF and first four-order TVMMs of the outputs. The conventional Johnson transformation is updated to accommodate the nonstationary cases for producing desired white noise inputs. This algorithm is developed from the linear filtering method (LFM), and inherits the simplicity and high efficiency from LFM. It fills the gaps in LFM-based algorithms for simulating NNP. Two numerical examples, i.e., a ground motion acceleration and a downburst velocity, are presented to fully demonstrate the capabilities of the proposed algorithm by comparing the simulation statistics with the targets.     

Reliability of dynamic systems under limited information

A method is developed for reliability analysis of dynamic systems under limited information. The available information includes one or more samples of the system output; any known information on features of the output can be used if available. The method is based on the theory of non-Gaussian translation processes and is shown to be particularly suitable for problems of practical interest. For illustration, we apply the proposed method to a series of relevant examples and compare with results given by traditional statistical estimators. It is demonstrated that the method delivers accurate results for the case of linear and nonlinear dynamic systems, and can be applied to analyze experimental data and/or mathematical model outputs.     

Fatigue damage assessment for a spectral model of non-Gaussian random loads

In this paper, a new model for random loads–the Laplace driven moving average–is presented. The model is second order, non-Gaussian, and strictly stationary. It shares with its Gaussian counterpart the ability to model any spectrum but has additional flexibility to model the skewness and kurtosis of the marginal distribution. Unlike most other non-Gaussian models proposed in the literature, such as the transformed Gaussian or Volterra series models, the new model is no longer derivable from Gaussian processes. In the paper, a summary of the properties of the new model is given and its upcrossing intensities are evaluated. Then it is used to estimate fatigue damage both from simulations and in terms of an upper bound that is of particular use for narrowband spectra.     

Simulation of stationary Gaussian/non-Gaussian stochastic processes based on stochastic harmonic functions

A new model is proposed to represent and simulate Gaussian/non-Gaussian stochastic processes. In the proposed model, stochastic harmonic function (SHF) is extended to represent multivariate Gaussian process firstly. Compared with the conventional spectral representation method (SRM), the SHF based model requires much fewer variables and Cholesky decompositions. Then, SHF based model is further extended to univariate/multivariate non-Gaussian stochastic process simulation. The target non-Gaussian process can be obtained from the corresponding underlying Gaussian processes by memoryless nonlinear transformation. For arbitrarily given marginal probability distribution function (PDF), the covariance function of the underlying multivariate Gaussian process can be determined easily by introducing the Mehler’s formula. And when the incompatibility between the target non-Gaussian power spectral density (PSD) or PSD matrix and marginal PDF exists, the calibration of the target non-Gaussian spectrum will be required. Hence, the proposed model can be regarded as SRM to efficiently generate Gaussian/non-Gaussian processes. Finally, several numerical examples are addressed to show the effectiveness of the proposed method.     

Linear models for non-Gaussian processes and applications to linear random vibration

Linear models are finite sums of specified deterministic, continuous functions of time with random coefficients. It is shown that linear models provide (i) accurate approximations for real-valued non-Gaussian processes with continuous samples defined on bounded time intervals, (ii) simple solutions for linear random vibration problems with non-Gaussian input, and (iii) efficient techniques for selecting optimal designs from collections of proposed alternatives. Theoretical arguments and numerical examples are presented to establish properties of linear models, illustrate the construction of linear models, solve linear random vibration with non-Gaussian input, and propose an approach for optimal design of linear dynamic systems. It is shown that the proposed linear model provides an efficient tool for analyzing linear systems in non-Gaussian environment.     

The geometry of random vibrations and solutions by FORM and SORM

The geometry of random vibration problems in the space of standard normal random variables obtained from discretization of the input process is investigated. For linear systems subjected to Gaussian excitation, the problems of interest are characterized by simple geometric forms, such as vectors, planes, half spaces, wedges and ellipsoids. For non-Gaussian responses, the problems of interest are generally characterized by non-linear geometric forms. Approximate solutions for such problems are obtained by use of the first- and second-order reliability methods (FORM and SORM). This article offers a new outlook to random vibration problems and an approximate method for their solution. Examples involving response to non-Gaussian excitation and out-crossing of a vector process from a non-linear domain are used to demonstrate the approach.     

Simulation of a non-Gaussian stochastic process based on a combined distribution of the UHPM and the GBD

To simulate non-Gaussian stochastic processes based on the first four moments, various simulation methods are presented, in which the determination of the transformation model and the calculation of the correlation coefficients between non-Gaussian stochastic processes and Gaussian stochastic processes are the critical procedures in these simulation methods. However, some existing simulation methods are limited to specific ranges. Furthermore, their practical applications are affected negatively due to the expensive cost of determining the transformation model and the correlation coefficients between non-Gaussian and Gaussian stochastic processes. Therefore, an accurate and efficient simulation method of a non-Gaussian stochastic process with a broader range is proposed in this article. Since the simulation of non-Gaussian processes and the Nataf transformation of non-Gaussian variables have some similar characteristics, a new combined distribution is proposed based on the unified Hermite polynomial model (UHPM) and the generalized beta distribution (GBD). Then, the combined distribution is employed in the simulation of non-Gaussian stochastic processes, in which the transformation model is deduced by the combined distribution. The correlation coefficient transformation function (CCTF) between the Gaussian and non-Gaussian stochastic processes can be evaluated through the interpolation method. Furthermore, numerical examples are presented to show the accuracy and effectiveness of the proposed simulation method for non-Gaussian stochastic processes.     

A simple and efficient methodology to approximate a general non-Gaussian stationary stochastic process by a translation process

Some widely used methodologies for simulation of non-Gaussian processes rely on translation process theory which imposes certain compatibility conditions between the non-Gaussian power spectral density function (PSDF) and the non-Gaussian probability density function (PDF) of the process. In many practical applications, the non-Gaussian PSDF and PDF are assigned arbitrarily; therefore, in general they can be incompatible. Several techniques to approximate such incompatible non-Gaussian PSDF/PDF pairs with a compatible pair have been proposed that involve either some iterative scheme on simulated sample functions or some general optimization approach. Although some of these techniques produce satisfactory results, they can be time consuming because of their nature. In this paper, a new iterative methodology is developed that estimates a non-Gaussian PSDF that: (a) is compatible with the prescribed non-Gaussian PDF, and (b) closely approximates the prescribed incompatible non-Gaussian PSDF. The corresponding underlying Gaussian PSDF is also determined. The basic idea is to iteratively upgrade the underlying Gaussian PSDF using the directly computed (through translation process theory) non-Gaussian PSDF at each iteration, rather than through expensive ensemble averaging of PSDFs computed from generated non-Gaussian sample functions. The proposed iterative scheme possesses two major advantages: it is conceptually very simple and it converges extremely fast with minimal computational effort. Once the underlying Gaussian PSDF is determined, generation of non-Gaussian sample functions is straightforward without any need for iterations. Numerical examples are provided demonstrating the capabilities of the methodology.     

Analysis of complex system reliability with correlated random vectors

The analysis of reliability of complex engineering systems remains a challenge in the field of reliability. It will be even more difficult if correlated random vectors are involved, which is generally the case as practical engineering systems invariably contain parameters that are mutually correlated. A new method for transforming correlated distributions, involving the Nataf transformation, is proposed that avoids the solution of integral equations; the method is based on the Taylor series expansion of the probability density function (PDF) of a bivariate normal distribution resulting in an explicit polynomial equation of the equivalent correlation coefficient. The required numerical results can be obtained efficiently and accurately.The proposed method for transformation of correlated random vectors is useful for developing a method for system reliability including complex systems with correlated random vectors. Based on the complete system failure process (originally defined as the development process of nonlinearity) and the fourth-moment method, the analysis of system reliability for elastic-plastic material avoids the identification of the potential failure modes of the system and their mutual correlations which are required in the traditional methods. Finally, four examples are presented – two examples to illustrate the potential of the new method for transformation of correlated random vectors, and two examples to illustrate the application of the proposed more effective method for system reliability.     

Simulation of non-Gaussian field applied to wind pressure fluctuations

A simulation algorithm to generate non-Gaussian wind pressure fields is proposed. This algorithm uses the correlation–distortion method based on translation vector processes. Conditions on the matrix of cross-covariance functions are given to assure the applicability of the model. The proposed method does not require iterative procedures and it is well suited when experimental data are available. In particular it requires cross-covariance functions and marginal distribution that can be directly estimated from data. To illustrate the procedure, the model is calibrated on experimental results obtained from wind tunnel tests on a tall building. The efficiency of the proposed methodology for reproducing the non-Gaussian nature of pressure fluctuations on separated flow regions is demonstrated.     

On the accuracy of the polynomial chaos approximation

Polynomial chaos representations for non-Gaussian random variables and stochastic processes are infinite series of Hermite polynomials of standard Gaussian random variables with deterministic coefficients. Finite truncations of these series are referred to as polynomial chaos (PC) approximations. This paper explores features and limitations of PC approximations. Metrics are developed to assess the accuracy of the PC approximation. A collection of simple, but relevant examples is examined in this paper. The number of terms in the PC approximations used in the examples exceeds the number of terms retained in most current applications. For the examples considered, it is demonstrated that (1) the accuracy of the PC approximation improves in some metrics as additional terms are retained, but does not exhibit this behavior in all metrics considered in the paper, (2) PC approximations for strictly stationary, non-Gaussian stochastic processes are initially nonstationary and gradually may approach weak stationarity as the number of terms retained increases, and (3) the development of PC approximations for certain processes may become computationally demanding, or even prohibitive, because of the large number of coefficients that need to be calculated. However, there have been many applications in which PC approximations have been successful.