Data Driven Surrogate Based Optimization in the Problem Solving Environment WBCSim

Large scale, multidisciplinary, engineering designs are always difficult due to the complexity and dimensionality of these problems. Direct coupling between the analysis codes and the optimization routines can be prohibitively time consuming due to the complexity of the underlying simulation codes. One way of tackling this problem is by constructing computationally cheap(er) approximations of the expensive simulations, that mimic the behavior of the simulation model as closely as possible. This paper presents a data driven, surrogate based optimization algorithm that uses a trust region based sequential approximate optimization (SAO) framework and a statistical sampling approach based on design of experiment (DOE) arrays. The algorithm is implemented using techniques from two packages—SURFPACK and SHEPPACK that provide a collection of approximation algorithms to build the surrogates and three different DOE techniques—full factorial (FF), Latin hypercube sampling (LHS), and central composite design (CCD)—are used to train the surrogates. The results are compared with the optimization results obtained by directly coupling an optimizer with the simulation code. The biggest concern in using the SAO framework based on statistical sampling is the generation of the required database. As the number of design variables grows, the computational cost of generating the required database grows rapidly. A data driven approach is proposed to tackle this situation, where the trick is to run the expensive simulation if and only if a nearby data point does not exist in the cumulatively growing database. Over time the database matures and is enriched as more and more optimizations are performed. Results show that the proposed methodology dramatically reduces the total number of calls to the expensive simulation runs during the optimization process

Fundamentals of Flakeboard Manufacture: Viscoelastic Behavior of the Wood Component

Theories of the viscoelastic behavior of amorphous polymers are reviewed and are used to describe the density gradient formation in flakeboard. This technique utilizes measured temperature and gas pressure at discrete locations inside a flake mat during hot pressing to predict the glass transition temperature of wood as a function of press time. The difference between the flake temperature and the predicted glass transition temperature is a relative indicator of the amount of flake deformation and stress relaxation at a location in the mat. A knowledge of the stress history imposed in the mat is then used to relate flake deformation and stress relaxation to the formation of a density gradient. This analysis allows for a significant portion of the density gradient to develop after the hot press has closed. Experimental data for various density gradients support the theories presented here

Approximations for many-body Green's functions: insights from the fundamental equations

Several widely used methods for the calculation of band structures and photo emission spectra, such as the GW approximation, rely on Many-Body Perturbation Theory. They can be obtained by iterating a set of functional differential equations relating the one-particle Green's function to its functional derivative with respect to an external perturbing potential. In the present work we apply a linear response expansion in order to obtain insights in various approximations for Green's functions calculations. The expansion leads to an effective screening, while keeping the effects of the interaction to all orders. In order to study various aspects of the resulting equations we discretize them, and retain only one point in space, spin, and time for all variables. Within this one-point model we obtain an explicit solution for the Green's function, which allows us to explore the structure of the general family of solutions, and to determine the specific solution that corresponds to the physical one. Moreover we analyze the performances of established approaches like $GW$ over the whole range of interaction strength, and we explore alternative approximations. Finally we link certain approximations for the exact solution to the corresponding manipulations for the differential equation which produce them. This link is crucial in view of a generalization of our findings to the real (multidimensional functional) case where only the differential equation is known.Comment: 17 pages, 7 figure

Matrix product approach for the asymmetric random average process

We consider the asymmetric random average process which is a one-dimensional stochastic lattice model with nearest neighbour interaction but continuous and unbounded state variables. First, the explicit functional representations, so-called beta densities, of all local interactions leading to steady states of product measure form are rigorously derived. This also completes an outstanding proof given in a previous publication. Then, we present an alternative solution for the processes with factorized stationary states by using a matrix product ansatz. Due to continuous state variables we obtain a matrix algebra in form of a functional equation which can be solved exactly.Comment: 17 pages, 1 figur

The Oscillating Universe: an Alternative to Inflation

The aim of this paper is to show, that the 'oscillating universe' is a viable alternative to inflation. We remind that this model provides a natural solution to the flatness or entropy and to the horizon problem of standard cosmology. We study the evolution of density perturbations and determine the power spectrum in a closed universe. The results lead to constraints of how a previous cycle might have looked like. We argue that most of the radiation entropy of the present universe may have originated from gravitational entropy produced in a previous cycle. We show that measurements of the power spectrum on very large scales could in principle decide whether our universe is closed, flat or open.Comment: revised version for publication in Classical and Quantum Gravity, 23 pages, uuencoded compressed tarred Latex file with 7 eps figures included, fig.8 upon reques

Regular spherical dust spacetimes

Physical (and weak) regularity conditions are used to determine and classify all the possible types of spherically symmetric dust spacetimes in general relativity. This work unifies and completes various earlier results. The junction conditions are described for general non-comoving (and non-null) surfaces, and the limits of kinematical quantities are given on all comoving surfaces where there is Darmois matching. We show that an inhomogeneous generalisation of the Kantowski-Sachs metric may be joined to the Lemaitre-Tolman-Bondi metric. All the possible spacetimes are explicitly divided into four groups according to topology, including a group in which the spatial sections have the topology of a 3-torus. The recollapse conjecture (for these spacetimes) follows naturally in this approach.Comment: Minor improvements, additional references. Accepted by GR

Green's function for a Schroedinger operator and some related summation formulas

Summation formulas are obtained for products of associated Lagurre polynomials by means of the Green's function K for the Hamiltonian H = -{d^2\over dx^2} + x^2 + Ax^{-2}, A > 0. K is constructed by an application of a Mercer type theorem that arises in connection with integral equations. The new approach introduced in this paper may be useful for the construction of wider classes of generating function.Comment: 14 page

Nonlinear modes for the Gross-Pitaevskii equation -- demonstrative computation approach

A method for the study of steady-state nonlinear modes for Gross-Pitaevskii equation (GPE) is described. It is based on exact statement about coding of the steady-state solutions of GPE which vanish as $x\to+\infty$ by reals. This allows to fulfill {\it demonstrative computation} of nonlinear modes of GPE i.e. the computation which allows to guarantee that {\it all} nonlinear modes within a given range of parameters have been found. The method has been applied to GPE with quadratic and double-well potential, for both, repulsive and attractive nonlinearities. The bifurcation diagrams of nonlinear modes in these cases are represented. The stability of these modes has been discussed.Comment: 21 pages, 6 figure

Statistics of Rare Events in Disordered Conductors

Asymptotic behavior of distribution functions of local quantities in disordered conductors is studied in the weak disorder limit by means of an optimal fluctuation method. It is argued that this method is more appropriate for the study of seldom occurring events than the approaches based on nonlinear $\sigma$-models because it is capable of correctly handling fluctuations of the random potential with large amplitude as well as the short-scale structure of the corresponding solutions of the Schr\"{o}dinger equation. For two- and three-dimensional conductors new asymptotics of the distribution functions are obtained which in some cases differ significantly from previously established results.Comment: 17 pages, REVTeX 3.0 and 1 Postscript figur

Particle trajectories in linearized irrotational shallow water flows

We investigate the particle trajectories in an irrotational shallow water flow over a flat bed as periodic waves propagate on the water's free surface. Within the linear water wave theory, we show that there are no closed orbits for the water particles beneath the irrotational shallow water waves. Depending on the strength of underlying uniform current, we obtain that some particle trajectories are undulating path to the right or to the left, some are looping curves with a drift to the right and others are parabolic curves or curves which have only one loop