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

The two-level atom laser: analytical results and the laser transition

The problem of the two-level atom laser is studied analytically. The steady-state solution is expressed as a continued fraction, and allows for accurate approximation by rational functions. Moreover, we show that the abrupt change observed in the pump dependence of the steady-state population is directly connected with the transition to the lasing regime. The condition for a sharp transition to Poissonian statistics is expressed as a scaling limit of vanishing cavity loss and light-matter coupling, $\kappa \to 0$, $g \to 0$, such that $g^2/\kappa$ stays finite and $g^2/\kappa > 2 \gamma$, where $\gamma$ is the rate of atomic losses. The same scaling procedure is also shown to describe a similar change to Poisson distribution in the Scully-Lamb laser model too, suggesting that the low-$\kappa$, low-$g$ asymptotics is of a more general significance for the laser transition.Comment: 23 pages, 3 figures. Extended discussion of the paper aim (in the Introduction) and of the results (Conclusions and Discussion). Results unchange

Effects of restricted basilar papillar lesions and hair cell regeneration on auditory forebrain frequency organization in adult European Starlings

The frequency organization of neurons in the forebrain Field L complex (FLC) of adult starlings was investigated to determine the effects of hair cell (HC) destruction in the basal portion of the basilar papilla (BP) and of subsequent HC regeneration. Conventional microelectrode mapping techniques were used in normal starlings and in lesioned starlings either 2 d or 6-10 weeks after aminoglycoside treatment. Histological examination of the BP and recordings of auditory brainstem evoked responses confirmed massive loss of HCs in the basal portion of the BP and hearing losses at frequencies >2 kHz in starlings tested 2 d after aminoglycoside treatment. In these birds, all neurons in the region of the FLC in which characteristic frequencies (CFs) normally increase from 2 to 6 kHz had CF in the range of 2-4 kHz. The significantly elevated thresholds of responses in this region of altered tonotopic organization indicated that they were the residue of prelesion responses and did not reflect CNS plasticity. In the long-term recovery birds, there was histological evidence of substantial HC regeneration. The tonotopic organization of the high-frequency region of the FLC did not differ from that in normal starlings, but the mean threshold at CF in this frequency range was intermediate between the values in the normal and lesioned short-recovery groups. The recovery of normal tonotopicity indicates considerable stability of the topography of neuronal connections in the avian auditory system, but the residual loss of sensitivity suggests deficiencies in high-frequency HC function

Exact Solutions of a (2+1)-Dimensional Nonlinear Klein-Gordon Equation

The purpose of this paper is to present a class of particular solutions of a C(2,1) conformally invariant nonlinear Klein-Gordon equation by symmetry reduction. Using the subgroups of similitude group reduced ordinary differential equations of second order and their solutions by a singularity analysis are classified. In particular, it has been shown that whenever they have the Painlev\'e property, they can be transformed to standard forms by Moebius transformations of dependent variable and arbitrary smooth transformations of independent variable whose solutions, depending on the values of parameters, are expressible in terms of either elementary functions or Jacobi elliptic functions.Comment: 16 pages, no figures, revised versio

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

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

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

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